dolfin team mailing list archive

[Anders Logg <question91944@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Anders Logg proposed the following answer:
On Fri, Nov 27, 2009 at 03:22:40PM -0000, Patrick Riesen wrote:
> New question #91944 on DOLFIN:
> https://answers.launchpad.net/dolfin/+question/91944
>
> Hello all,
>
> I'd like to extract the points of a 2d mesh boundary subdomain as a 1D mesh
> for separate computation of initial/boundary conditions.

You can do this simply by

  boundary_mesh = BoundaryMesh(mesh)

> Using the SubMesh class crashed in dolfin::TopologyComputation::compute_entities with a segfault.
>
> Is it supposed to work to get a 1D mesh from 2D with SubMesh?
> Or do i need to take another approach to extract a 1D mesh from 2D?
>
> thanks for your help

The problem with the above (BoundaryMesh) is that it does give you the
boundary mesh which is topologically 1D but it is embedded in R^2
which means you'll have trouble assembling over it since DOLFIN
assumes nD meshes to be embedded in R^n.

The SubMesh thing is intended for different use.

--
Anders

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