Re: [petsc-users] List of points with dof>0 in a PetscSection

2022-06-28 Thread Matthew Knepley
On Fri, Jun 10, 2022 at 4:06 PM Blaise Bourdin wrote: > Hi, > > Given a PetscSection, is there an easy way to get a list of point at which > the number of dof is >0? > For instance, when projecting over a FE space, I’d rather do a loop over > such points than do a loop over all points in a DM, ge

Re: [petsc-users] Load mesh as DMPlex along with Solution Fields obtained from External Codes

2022-06-28 Thread Matthew Knepley
On Wed, Jun 8, 2022 at 12:15 AM Mike Michell wrote: > Dear PETSc developer team, > > I am a user of PETSc DMPlex for a finite-volume solver. So far, I have > loaded a mesh file made by Gmsh as a DMPlex object without pre-computed > solution field. > But what if I need to load the mesh as well as