I have a Dolfin mesh that defines two or more material/sub-domains like:
mesh = Mesh("test.xml") subd = MeshFunction("size_t", mesh, mesh.topology().dim(), mesh.domains() ) dx = Measure("dx", domain=mesh, subdomain_data=subd) E0 = FiniteElement("P", mesh.ufl_cell(), 1) V = FunctionSpace(mesh, E0) u = Function(V) v = TestFunction(V)
I can use this (dx measuse) to define PDE forms in different domains like:
c1*u*v*dx(0) + c2*u*v*dx(1)
However, I would like to interpolate/project different constants to
the different subdomains for the initial solution vector, something like
u0 = [ interpolate(Constant(0.0), V.dx(0)), \ interpolate(Constant(1.0), V.dx(1)) ]
Is this possible somehow (as interpolate(Constant(1.0),V) affects everything)?