dolfin team mailing list archive

[Melanie Jahny <question157341@xxxxxxxxxxxxxxxxxxxxx>]

Question #157341 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/157341

Melanie Jahny gave more information on the question:
I tried several tests with boundary integrals.
Is it possible that sometimes the compiler cannot find the boundary?
Because for some (outer) boundary it works, for some it doesn't.
I use the ffc version 0.9.4.
Here are the ways I tried to define the boundary:

1.
Bound_Test = compile_subdomains('x[1]==10.0 && on_boundary')
Bound_Test.mark(boundary_parts,1)

2.

class Out_test(SubDomain):  
   def inside(self,x,on_boundary):
       return True if x[1] > 10.0-DOLFIN_EPS else False

test = Out_test()
test.mark(boundary_parts,1)

3.

test_pt = Point(0.0,10.0)
near_test = Mesh.closest_point(mesh,test_pt)
y_test = near_test[1]

class Out_test(SubDomain):  
   def inside(self,x,on_boundary):
       return True if x[0] > y_test-DOLFIN_EPS else False

test = Out_test()
test.mark(boundary_parts,2)

Is there any other way I can try? Do I make any mistake?
Thanks for any help
Melanie

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.