-
assigned issue to
Screwed up geometry of ghosted plots
Plots of ghosted meshes seem to have correct topology but geometry of ghosted cells is rubish.
from dolfin import *
import sys
parameters['ghost_mode'] = sys.argv[1]
mesh = UnitSquareMesh(2, 2)
rank = MPI.rank(mesh.mpi_comm())
for opt in ['regular', 'ghost']:
for c in cells(mesh, opt):
print rank, opt, c.index(), c.global_index()
plot(mesh, title='mesh at rank %d' % rank)
interactive()
Run mpirun -n 2 python 1.py shared_facet
or mpirun -n 2 python 1.py shared_vertex
and push 'c' in plots.
Next step when this is resolved is checking whether Function
s are plotted correctly.
Comments (11)
-
reporter -
Is it really useful to use plot in parallel?
-
reporter Well, not for serious stuff. But debugging, checking topological properties is exactly thing where it is useful. Feel free to reassign to me. I guess that the issue has a good ratio priority/difficulty.
-
-
assigned issue to
I'm not into hacking around in VTK. I'd like to see more effort into getting some matplotlib plot routines working (that's what I used to debug ghost meshes).
-
assigned issue to
-
reporter - changed status to resolved
-
reporter -
OK - so that was easy.
-
reporter The sad part is that is another of several bugs introduced by omitting
MeshEntityIterator(mesh) -> MeshEntityIterator(mesh, opt)
in whole DOLFIN. -
I guess it was always a danger - where are the other problems?
-
reporter Other was
#501. -
- removed milestone
Removing milestone: 1.7 (automated comment)
- Log in to comment
Maybe easy for Chris?