This message finally appeared in the spam filters for the libmesh-devel
mailing list, but I believe it has already been discussed.
--
John
On Mon, Jan 9, 2017 at 2:57 PM, Nestola Maria Giuseppina Chiara <
maria.giuseppina.chiara.nest...@usi.ch> wrote:
> Dear all,
> we are trying to build parallel meshes. What i need is that each
> processor has just its own mesh.
>
> We are coupling a fem code with FD code thus we need to guarantee also the
> same partition.
>
> We have already implemented a routine to do this but we have some
> questions.
>
> 1)
>
> If our setup of elements already includes elements across process
> boundaries, do we need to still call gather_neighboring_elements()?
> I added the following line to our code:
>
> MeshCommunication().gather_neighboring_elements(cast_ref<DistributedMesh
> &>(mesh));
>
> but it seems to take forever to execute that in dbg mode (Even with very
> few grid points i.e. ~100 [8 procs]).
>
> 2)
>
> The mesh looks ok if I output it even without the above line added but I
> trip different asserts in dbg mode depending on the type of element I use
> (I have both implemented).
>
> If I use Tri3 I trip:
>
> Assertion `oldvar == static_cast<Told>(static_cast<Tnew>(oldvar))' failed.
> oldvar = 18446744073709551615
> static_cast<Told>(static_cast<Tnew>(oldvar)) = 4294967295
>
> If I use Quad4 I trip:
>
> Assertion `i < _points.size()' failed.
> i = 0
> _points.size() = 0
>
> The messages are identical on each process.
>
>
> Best
> Barna & Maria
>
> ------------------------------------------------------------
> ------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
> _______________________________________________
> Libmesh-devel mailing list
> Libmesh-devel@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/libmesh-devel
>
>
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, SlashDot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
