On Wed, Apr 20, 2011 at 10:38:36AM -0000, B. Emek Abali wrote: > Question #153567 on DOLFIN changed: > https://answers.launchpad.net/dolfin/+question/153567 > > Status: Open => Answered > > B. Emek Abali proposed the following answer: > instead of > > x[0] > (20.0 - DOLFIN_EPS) > > try to use > > fabs(x[0] - 20.0) < DOLFIN_EPS > > or an easier one > right = compile_subdomains('x[0] == 20.0') > ... > bc2 = DirichletBC(system.sub(1), zero, right)
You can also use near(x[0], 20.0) -- Anders _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
