← Back to team overview

dolfin team mailing list archive

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.