dolfin team mailing list archive

["B. Emek Abali" <question155482@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

B. Emek Abali proposed the following answer:
you are on a good way, please send the whole code if it is finished, it
looks good!

first, interior domain is cell_domain for the FEniCS jargon, so if you
mark your sub_channel and use dx(0) in your variational form, try to
assemble(a, cell_domains=subdomains, interior_...) so that compiler
knows where to look for the values which have 0. MeshFunction is a list
of cell ids and its values, value is just a property to mark, if you
use...mark(...,1) you change all the values to 1 and then in dx(1) you
call only the cell ids, listed in the MeshFunction X (which is linked:
assemble(...cell_domains=X...) ) ,

second, in your bound_mem (bounding membran? ) you call with ds (small
letter s) ONLY the boundary (exterior facets) lines, boundary to your
background mesh=Rectangle(...) you want to use dS (capital letter s) to
reach all the interior facets or dS(0) for the one with value (marked as
0 ) in your MeshFunction (assmble(...interior...=boundary_parts...)

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