# Dirichlet boundary function not defined over entire face

**URL:** <https://community.freefem.org/t/dirichlet-boundary-function-not-defined-over-entire-face/2226>\
**Category:** General Discussion\
**Created:** [January 16, 2023, 10:55pm UTC](https://community.freefem.org/t/dirichlet-boundary-function-not-defined-over-entire-face/2226 "2023-01-16T22:55:14Z")\
**Posts on this page:** 1\
**Showing post:** 6

<div class="post-metadata">

**Author:** ![jsmacd](https://avatars.discourse-cdn.com/v4/letter/j/df788c/32.png) [@jsmacd](https://community.freefem.org/u/jsmacd)\
**Post date:** [January 24, 2023, 7:33pm UTC](https://community.freefem.org/t/dirichlet-boundary-function-not-defined-over-entire-face/2226/6 "2023-01-24T19:33:08Z")

</div>

Hi @cmd, I’ve been working on a solution to my problem and I think I’ve isolated where my issues are. Would you mind taking a look at this post?

> [@Piecewise function to eliminate singularity but singularity still causing problem](https://community.freefem.org/t/piecewise-function-to-eliminate-singularity-but-singularity-still-causing-problem/2235):
>
> Dear community, I am trying to eliminate a singularity within a cube by replacing the function with the singularity with a piecewise function that has no singularity. However, I’m getting some funky errors including an attempt to divide by 0. But I’m trying very hard not to divide by zero. Can someone see where I’ve gone wrong? Here is the code that’s giving me the error. load "msh3" // to use cube load "MUMPS" // real myRadius = 0.3; func uBC = (1/(((x - 0.5)^2 + (y - 0.5)^2)^0.5))\*(((x - 0…

---

_[View the full topic](https://community.freefem.org/t/dirichlet-boundary-function-not-defined-over-entire-face/2226)._
