Commits

Matt Knepley committed e19c037

DMPlex: Fixed small bug
- Was not initializing cone size

Comments (0)

Files changed (1)

src/dm/impls/plex/plexrefine.c

         supportRef[0] = fStartNew + (f - fStart)*4 + r;
         supportRef[1] = fStartNew + (f - fStart)*4 + (r+1)%4;
         for (s = 0; s < supportSize; ++s) {
+          ierr = DMPlexGetConeSize(dm, f, &coneSize);CHKERRQ(ierr);
           ierr = DMPlexGetCone(dm, f, &coneCell);CHKERRQ(ierr);
           for (c = 0; c < coneSize; ++c) if (coneCell[c] == f) break;
           supportRef[2+s] = fStartNew + (f - fStart)*4 + newFaces[c*4+r];