Re: [petsc-users] How to find the map between the high order coordinates of DMPlex and vertex numbering?

2022-06-16 Thread Zongze Yang
I tried the projection operation. However, it seems that the projection gives the wrong solution. After projection, the bounding box is changed! See logs below. First, I patch the petsc4py by adding `DMProjectCoordinates`: ``` diff --git a/src/binding/petsc4py/src/PETSc/DM.pyx b/src/binding/petsc4

Re: [petsc-users] Customizing NASM subsnes

2022-06-16 Thread Matthew Knepley
On Thu, Jun 16, 2022 at 5:57 PM tt73 wrote: > > Hi, > > I am using NASM as the outer solver for a nonlinear problem. For one of > the subdomains, I want to run the local solve with a different set of > options form the others. Is there any way to set options for each > subdomain? > I can see tw

[petsc-users] Customizing NASM subsnes

2022-06-16 Thread tt73
Hi, I am using  NASM as the outer solver for a nonlinear problem. For one of the subdomains, I want to run the local solve with a different set of options form the others. Is there any way to set options for each subdomain? 

Re: [petsc-users] How to find the map between the high order coordinates of DMPlex and vertex numbering?

2022-06-16 Thread Matthew Knepley
On Thu, Jun 16, 2022 at 12:06 PM Zongze Yang wrote: > > > 在 2022年6月16日,23:22,Matthew Knepley 写道: > >  > On Thu, Jun 16, 2022 at 11:11 AM Zongze Yang wrote: > >> Hi, if I load a `gmsh` file with second-order elements, the coordinates >> will be stored in a DG-P2 space. After obtaining the coord

Re: [petsc-users] How to find the map between the high order coordinates of DMPlex and vertex numbering?

2022-06-16 Thread Zongze Yang
> 在 2022年6月16日,23:22,Matthew Knepley 写道: > >  >> On Thu, Jun 16, 2022 at 11:11 AM Zongze Yang wrote: > >> Hi, if I load a `gmsh` file with second-order elements, the coordinates will >> be stored in a DG-P2 space. After obtaining the coordinates of a cell, how >> can I map the coordinates

Re: [petsc-users] How to find the map between the high order coordinates of DMPlex and vertex numbering?

2022-06-16 Thread Matthew Knepley
On Thu, Jun 16, 2022 at 11:11 AM Zongze Yang wrote: > Hi, if I load a `gmsh` file with second-order elements, the coordinates > will be stored in a DG-P2 space. After obtaining the coordinates of a cell, > how can I map the coordinates to vertex and edge? > By default, they are stored as P2, not

Re: [petsc-users] SLEPc EPSGD: too much time in single iteration

2022-06-16 Thread Matthew Knepley
On Wed, Jun 15, 2022 at 9:32 PM Runfeng Jin wrote: > Hi! Thank you for your reply. > > I am a little confused about the problem of machine. These two matrices > solved in the same cluster, if there are some problems about the machine, > why the low performance just happen to the matrix B? > The