Dirichlet boundary function not defined over entire face

[jsmacd]

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

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…