dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #20509
Re: [Question #139472]: a boundary condition for sphere geometry
Question #139472 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/139472
Status: Open => Answered
Garth Wells proposed the following answer:
There are some syntax errors in '1st code'. It should be:
from dolfin import *
sphere = UnitSphere(10)
Q = FunctionSpace(sphere, "CG", 1)
outer= compile_subdomains("x[0]*x[0] + x[1]*x[1] +x[2]*x[2] > 0.9*0.9 - 5*DOLFIN_EPS && on_boundary")
g = Constant(1.0)
BC = DirichletBC(Q, g, outer)
--
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.