dolfin team mailing list archive

[Niclas Jansson <njansson@xxxxxx>]

logg <logg@xxxxxxxxx> writes:

> Public bug reported:
>
> When running the parallel assembly/solve with 9 processes, I get the
> following error message:
>
> dolfin/fem/Assembler.cpp:202: static void
> dolfin::Assembler::assemble_exterior_facets(dolfin::GenericTensor&,
> const dolfin::Form&, dolfin::UFC&, const dolfin::MeshFunction<unsigned
> int>*, std::vector<double, std::allocator<double> >*): Assertion
> `cell_map' failed.
>

It fails since the boundary meshes can be empty (some processes might
have partitions completely inside the domain).  

Another problem besides the failed assert is the CellIterator, it cant
handle empty meshes.


A quick, and somewhat ugly fix in assembly

Compute boundary
if (empty) return;
assert(....)

for(CellIterator
....


Niclas