dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Thu, Sep 10, 2009 at 12:45:02PM +0200, Niclas Jansson wrote:
> 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
> ....

I think that fix looks good. And it also works. :-)

I'll make the fix and close the bug.

--
Anders

Attachment: signature.asc
Description: Digital signature