<div dir="ltr"><div>Dear PETSc users,</div><div><br></div><div>I have a mesh where I have a label that contains some elements. If I understand correctly, after uniform refinement, the label also contains the refined faces and vertices. Since I'd like to have only cells in the label, I thought to remove the rest using something like:</div><div><br></div><div>DMPlexGetDepthStratum(dm, 3, &vStart, &vEnd);</div><div>DMLabelFilter(label, vStart, vEnd);</div><div><br></div><div>However, I don't fully understand the functioning of DMLabelFilter. When vStart>0, it works as I would expect, but when vStart=0, it simply removes all the points from the label. Why is this happening?</div><div><br></div><div>Since for the cells vStart usually is 0, I can't use this method. Am I missing something?</div><div><br></div><div>Attached is a minimal working example.<br></div><div><br></div><div>Best,</div><div>Emiel<br></div></div>