Re: [petsc-users] Adding more dofs to the coordinate DM crashes with either PetscSectionSetDof or PetscSectionFieldDof

2022-07-25 Thread Duan Junming via petsc-users
Thank you very much! From: knep...@gmail.com Sent: Sunday, July 24, 2022 5:32:34 PM To: Duan Junming Cc: PETSc Subject: Re: [petsc-users] Adding more dofs to the coordinate DM crashes with either PetscSectionSetDof or PetscSectionFieldDof On Sun, Jul 24, 2022

Re: [petsc-users] Adding more dofs to the coordinate DM crashes with either PetscSectionSetDof or PetscSectionFieldDof

2022-07-24 Thread Matthew Knepley
On Sun, Jul 24, 2022 at 8:46 AM Duan Junming via petsc-users < petsc-users@mcs.anl.gov> wrote: > Dear all, > > > I want to add more dofs to the coordinate DM to represent a curved mesh. > > I first create a simple box mesh with one cell: > > > PetscCall(DMCreate(comm, )); >

[petsc-users] Adding more dofs to the coordinate DM crashes with either PetscSectionSetDof or PetscSectionFieldDof

2022-07-24 Thread Duan Junming via petsc-users
Dear all, I want to add more dofs to the coordinate DM to represent a curved mesh. I first create a simple box mesh with one cell: PetscCall(DMCreate(comm, )); PetscCall(DMSetType(dm, DMPLEX)); dim = 2; PetscInt n_faces[2] = {1, 1}; DMBoundaryType periodicity[2] =