Dear all
I have updated a code of mine to Petsc3.22, when trying to
uninterpolate a periodic mesh and get the following error
[0]PETSC ERROR: - Error Message
--
[0]PETSC ERROR: No support for this operation f
all is ok, but with
mpirun -np 2 ./saveDemo
it creates a 10x10 grid, but I get "output.vtu" with a total
of 120 cells. However the pointSF of the DMPlex seems correct.
Thanks
Matteo
Il 21/10/24 19:15, Matthew Knepley ha scri
you need (and do)
provide an appropriate mat.
Barry
On Dec 10, 2024, at 7:02 AM, Matteo Semplice via petsc-users
wrote:
Dear petsc-users,
I am trying with a student to modify the MG example 65 to use
mat-shells instead of assembled matrices. Our use case is for a
method that will
vec in the same way you need (and do)
provide an appropriate mat.
Barry
On Dec 10, 2024, at 7:02 AM, Matteo Semplice via petsc-users
wrote:
Dear petsc-users,
I am trying with a student to modify the MG example 65 to use
mat-shells instead of assembled matrices. Our use case is for a
Dear petsc-users,
I am trying with a student to modify the MG example 65 to use
mat-shells instead of assembled matrices. Our use case is for a method
that will use custom shell operators and shell interpolation/restrictions.
To start we have modified ex65 and tried to replace the standar
quot;output.vtu" with a total of
120 cells. However the pointSF of the DMPlex seems correct.
Thanks
Matteo
Il 21/10/24 19:15, Matthew Knepley ha scritto:
On Mon, Oct 21, 2024 at 12:22 PM Matteo Semplice via petsc-users
wrote:
Dear petsc-users,
pirun -np 2 ./saveDemo
it creates a 10x10 grid, but I get "output.vtu" with a total of 120
cells. However the pointSF of the DMPlex seems correct.
Thanks
Matteo
Il 21/10/24 19:15, Matthew Knepley ha scritto:
On Mon, Oct 21, 2024 at 12:22 PM Matteo Semplice via petsc-users
wr
Dear petsc-users,
I am having issues with output of parallel data attached to a
DMPlex (or maybe more fundamental ones about DMPlex...).
So I currently
1. create a DMPlex (DMPlexCreateGmshFromFile or DMPlexCreateBoxMesh)
2. partition it
3. and create a section for my data layout with
D
On 23/09/2024 17:16, Matthew Knepley wrote:
On Mon, Sep 23, 2024 at 5:33 AM Matteo Semplice via petsc-users
wrote:
Dear petsc,
I need to hand-code a jacobian and I can't figure out how to
translate the DMplex points/fields to matrix indices.
My DMPlex has a se
Dear petsc,
I need to hand-code a jacobian and I can't figure out how to
translate the DMplex points/fields to matrix indices.
My DMPlex has a section with m fields per cell (which makes for n>m dof
per cell since some fields are vector). Say I want to insert an nxn
block for the row cor
10 matches
Mail list logo
