Re: [petsc-users] [DMCreate from a Gmsh with mixed mesh crash for DMPlex]

On Thu, Jul 14, 2022 at 6:33 PM Mike Michell wrote: > VTK format itself supports pyramids for sure. But I do not think PETSc > output part through VTK format does not. Is there any possibility to debug > this issue? > VTK output now works with pyramids in https://gitlab.com/petsc/petsc/-/merg

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

Thank you for the corrections. It works. However, I still have a problem. I tested zero stratum for DMGetStratumIS() with the default depth label, since the graph of vertex is needed. Then, PetscSFGetGraph() from the Natural SF is crashed. I checked that pBcCompIS contains proper set of integers v

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

On Mon, Jul 18, 2022 at 12:14 PM Bora Jeong wrote: > Thank you for the corrections. It works. > > However, I still have a problem. > I tested zero stratum for DMGetStratumIS() with the default depth label, > since the graph of vertex is needed. Then, PetscSFGetGraph() from the > Natural SF is cra

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

On Mon, Jul 18, 2022 at 4:30 PM Bora Jeong wrote: > Actually, it seems that my problem is that I do not properly understand > what Natural SF means. > > What I want to achieve is to extract natural vertex index from Natural SF, > and make sure this natural vertex index is the same with what is wr

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

Thank you for the reply. I am actually putting 1dof on the vertex points. In the attached sample code, from line 63 to line 99, 1dof on all vertex points is defined for the section. But my problem is, if I print out Natural SF using the viewer, I only have 10 vertex points, instead of 15. Here 15

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

Alexis noticed a problem with the natural SF when constraints are defined. He has a MR coming.  @Alexis: Could it be that what Bora sees is exactly what you fixed? Blaise On Jul 18, 2022, at 8:35 PM, Bora Jeong wrote: Thank you for the reply.  I am act

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

On Mon, Jul 18, 2022 at 7:35 PM Bora Jeong wrote: > Thank you for the reply. > > I am actually putting 1dof on the vertex points. > In the attached sample code, from line 63 to line 99, 1dof on all vertex > points is defined for the section. But my problem is, if I print out > Natural SF using th

Re: [petsc-users] PETSc / AMRex

> On Jul 15, 2022, at 12:26 PM, Matthew Knepley wrote: > > On Fri, Jul 15, 2022 at 2:12 PM Randall Mackie > wrote: > > >> On Jul 15, 2022, at 11:58 AM, Matthew Knepley > > wrote: >> >> On Fri, Jul 15, 2022 at 1:46 PM Randall Mackie > <