Matt,
If DMPlexBuildFromCellListParallel() is called with two edges that have
the same two vertices what will happen? It looks like it ends up with an
incorrect PetscSF if the two edges are passed on different ranks. Hence the
DMPLEX is not valid and produces garbage.
Neurons can be connected to themselves which seems to be breaking DMPLEX
and hence DMNETWORK.
Is there anything we can do to support having multiple edges between the
same two vertices? If not is there a way we can have
DMPlexBuildFromCellListParallel() generate an error automatically if there are
such extra edges in the input data.
Thanks
Barry
In this work, the neurons are represented by vertices in the network and each
synapse is a graph edge.
