Re: [petsc-users] Including constrained dofs in a matrix/vector created with DMCreateXXXX()

2019-02-20 Thread Jordan Wagner via petsc-users
On 2/19/19 7:40 AM, Matthew Knepley wrote: On Tue, Feb 19, 2019 at 1:13 AM Jordan Wagner via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Hi, Over the past few months, I have implemented dmplex/section in my preexisting finite element code. We already had o

[petsc-users] Including constrained dofs in a matrix/vector created with DMCreateXXXX()

2019-02-18 Thread Jordan Wagner via petsc-users
Hi, Over the past few months, I have implemented dmplex/section in my preexisting finite element code. We already had our own implementations for things like the finite element space, boundary conditions, projections, etc, so I am mostly using dmplex/section at its lowest level; no DT/DS

[petsc-users] Can PetscSections of a DMPlex be written out and read in with HDF5?

2018-07-13 Thread Jordan Wagner
Hi, In my code, I have a separate preprocessing script that I usually run in serial to set up various things for my solver. One of the things I want the preprocessor to do is setup everything related to the DMPlex, and then I write it out to an HDF5 file to be read in the main solver. In the

[petsc-users] Is there a partition section in a distributed DM object?

2018-01-23 Thread Jordan Wagner
Hi, At some point, I would like to export my distributed DM along with an associated scalar field corresponding to the process rank for which each point is owned. I know that what I need is the partition section, which is more or less identical to Fig. 6 in Knepley and Karpeev (2009). I was

Re: [petsc-users] Function to convert a dense matrix holding the cell connectivity to a MPIADJ for use with MatMeshToCellGraph

2017-12-12 Thread Jordan Wagner
03:54 PM, Matthew Knepley wrote: Barry wrote this, so he probably knows how to fix it. Another option is to use DMPlex for your mesh. It will give you the dual if you want.   Thanks,      Matt On Tue, Dec 12, 2017 at 3:44 PM, Jordan Wagner <j.wag...@rice.edu <mailto:j.wag...@rice.edu&g

[petsc-users] Function to convert a dense matrix holding the cell connectivity to a MPIADJ for use with MatMeshToCellGraph

2017-12-12 Thread Jordan Wagner
Hi, I am trying to use the function MatMeshToCellGraph. I currently have a matrix that holds the cell connectivity of simplex elements. So it is a numCells x 3 matrix where the row corresponds to the cell number and the column is a vertex of that cell. To use this function, it appears I need