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

2023-05-14 Thread Zongze Yang
Could you try to project the coordinates into the continuity space by enabling the option `-dm_plex_gmsh_project_petscdualspace_lagrange_continuity true`? Best wishes, Zongze On Mon, 15 May 2023 at 04:24, Matthew Knepley wrote: > On Sun, May 14, 2023 at 12:27 PM Zongze Yang wrote: > >> >> >>

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

2023-05-14 Thread Matthew Knepley
On Sun, May 14, 2023 at 12:27 PM Zongze Yang wrote: > > > > On Sun, 14 May 2023 at 23:54, Matthew Knepley wrote: > >> On Sun, May 14, 2023 at 9:21 AM Zongze Yang wrote: >> >>> Hi, Matt, >>> >>> The issue has been resolved while testing on the latest version of >>> PETSc. It seems that the probl

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

2023-05-14 Thread Zongze Yang
On Sun, 14 May 2023 at 23:54, Matthew Knepley wrote: > On Sun, May 14, 2023 at 9:21 AM Zongze Yang wrote: > >> Hi, Matt, >> >> The issue has been resolved while testing on the latest version of PETSc. >> It seems that the problem has been fixed in the following merge request: >> https://gitlab.c

Re: [petsc-users] SNESDMDASNESSetFunctionLocal in Fortran

2023-05-14 Thread Matthew Knepley
On Sun, May 14, 2023 at 12:06 PM Khaled Nabil Shar Abdelaziz < kabde...@purdue.edu> wrote: > Hey there, > > > > I'm having a problem with the DMDASNESSetFunctionLocal() function in C and > its Fortran counterpart. The thing is, in C, you can pass a bunch of > variables using the ctx parameter, but

[petsc-users] SNESDMDASNESSetFunctionLocal in Fortran

2023-05-14 Thread Khaled Nabil Shar Abdelaziz
Hey there, I'm having a problem with the DMDASNESSetFunctionLocal() function in C and its Fortran counterpart. The thing is, in C, you can pass a bunch of variables using the ctx parameter, but in Fortran, it only seems to accept one variable. What's weird is that the SNESSetFunction() function

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

2023-05-14 Thread Matthew Knepley
On Sun, May 14, 2023 at 9:21 AM Zongze Yang wrote: > Hi, Matt, > > The issue has been resolved while testing on the latest version of PETSc. > It seems that the problem has been fixed in the following merge request: > https://gitlab.com/petsc/petsc/-/merge_requests/5970 > No problem. Glad it is

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

2023-05-14 Thread Zongze Yang
Yes, you are correct. I have conducted tests using high-order 3D meshes of a unit cube, and regrettably, the tests have failed. I have attached the files for your reference. Kindly review the output provided below: ( The volume should be 1) ``` $ mpiexec -n 3 ./ex33 -coord_space 0 -dm_plex_filena

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

2023-05-14 Thread Jed Brown
Good to hear this works for you. I believe there is still a problem with high order tetrahedral elements (we've been coping with it for months and someone asked last week) and plan to look at it as soon as possible now that my semester finished. Zongze Yang writes: > Hi, Matt, > > The issue h

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

2023-05-14 Thread Zongze Yang
Hi, Matt, The issue has been resolved while testing on the latest version of PETSc. It seems that the problem has been fixed in the following merge request: https://gitlab.com/petsc/petsc/-/merge_requests/5970 I sincerely apologize for any inconvenience caused by my previous message. However, I w

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

2023-05-14 Thread Matthew Knepley
On Sat, May 13, 2023 at 6:08 AM Zongze Yang wrote: > Hi, Matt, > > There seem to be ongoing issues with projecting high-order coordinates > from a gmsh file to other spaces. I would like to inquire whether there are > any plans to resolve this problem. > > Thank you for your attention to this mat