[petsc-users] PLEX depth stratum bug ?

Dharmendar Reddy dharmareddy84 at gmail.com
Thu Jan 23 23:05:21 CST 2014


Hello,
         I am seeing inconsistent results using get stratum call.  For a 2d mesh

call DMPlexGetDepthStratum(dmMesh,2,globalidStart,globalidEnd,ierr)

prints: globalIdStart = 0 and globalIdEnd = 0


where as

call DMPlexGetHeightStratum(dmMesh,0,globalidStart,globalidEnd,ierr)

prints globalIdStart = 0 and globalIdEnd = 4

For a mesh with 4 cells. I have attached the test case. I am on petsc
version v3.4.3-2568-gb45093d
-------------- next part --------------
A non-text attachment was scrubbed...
Name: testDMDepthStrata.F90
Type: text/x-fortran
Size: 5857 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140123/8715a689/attachment-0001.bin>


More information about the petsc-users mailing list