Comment 2 for bug 1143977

Revision history for this message
Anders Logg (logg) wrote :

This is broken because someone broke the connectivity computation for 0-0 connectivity. I found the following in TopologyComputation.cpp. It just sets every vertex to be connected to itself. Why was this changed?

 if (d0 == 0 && d1 == 0)
  {
    std::vector<std::vector<std::size_t> >
      connectivity00(topology.size(0), std::vector<std::size_t>(1));
    for (MeshEntityIterator v(mesh, 0); !v.end(); ++v)
      connectivity00[v->index()][0] = v->index();
    topology(0, 0).set(connectivity00);
  }