Re: [petsc-users] dm_view of high-order geometry/solution

CGNS supports fourth order and it's coded in PETSc, but Paraview hasn't implemented reading it yet. I think it would not be much work for someone (maybe you) to add it to Paraview. I have lots of applications on cubics, but not much beyond that so it hasn't risen to top priority for me. There

Re: [petsc-users] dm_view of high-order geometry/solution

Dear Jed, Thank you for your help! Now I moved the line using "DMViewFromOptions" after the function "PetscDSSetObjective", and it works for "-dm_coord_petscspace_degree 3 -petscspace_degree 3". But when I tried degree 4: ./ex33 -dm_plex_simplex 0 -dm_plex_box_faces 1,1 -mesh_transform

Re: [petsc-users] dm_view of high-order geometry/solution

And here's an MR to do what you want without any code/arg changes. https://gitlab.com/petsc/petsc/-/merge_requests/6588 Jed Brown writes: > Duan Junming writes: > >> Dear Jed, >> >> >> Thank you for the suggestion. >> >> When I run tests/ex33.c with >> >> ./ex33 -dm_plex_simplex 0

Re: [petsc-users] dm_view of high-order geometry/solution

Duan Junming writes: > Dear Jed, > > > Thank you for the suggestion. > > When I run tests/ex33.c with > > ./ex33 -dm_plex_simplex 0 -dm_plex_box_faces 1,1 -mesh_transform annulus > -dm_coord_space 0 -dm_coord_petscspace_degree 3 -dm_refine 1 -dm_view > cgns:test.cgns > > and load it using

Re: [petsc-users] dm_view of high-order geometry/solution

Dear Jed, Thank you for the suggestion. When I run tests/ex33.c with ./ex33 -dm_plex_simplex 0 -dm_plex_box_faces 1,1 -mesh_transform annulus -dm_coord_space 0 -dm_coord_petscspace_degree 3 -dm_refine 1 -dm_view cgns:test.cgns and load it using Paraview, the mesh is still with straight

Re: [petsc-users] dm_view of high-order geometry/solution

Matthew Knepley writes: > On Mon, Jun 12, 2023 at 6:01 AM Duan Junming wrote: > >> Dear Matt, >> >> Thank you for the reply. I have a more specific question about the >> spectral element example. Do you have any suggestions that how to write >> all the nodes in each cell to .vtu? >> > It is the

[petsc-users] Bug Report TaoALMM class

Dear PETSc/Tao team, I think there might be a bug in the Tao ALMM class:  In the function TaoALMMComputeAugLagAndGradient_Private(), see, eg. https://petsc.org/release/src/tao/constrained/impls/almm/almm.c.html#TAOALMM line 648 the gradient seems to be wrong. The given function and gradient

Re: [petsc-users] dm_view of high-order geometry/solution

On Mon, Jun 12, 2023 at 6:01 AM Duan Junming wrote: > Dear Matt, > > Thank you for the reply. I have a more specific question about the > spectral element example. Do you have any suggestions that how to write > all the nodes in each cell to .vtu? > It is the same procedure. VTU is not a great

Re: [petsc-users] dm_view of high-order geometry/solution

Dear Matt, Thank you for the reply. I have a more specific question about the spectral element example. Do you have any suggestions that how to write all the nodes in each cell to .vtu? Thanks! Junming From: knep...@gmail.com Sent: Monday, June 12, 2023

Re: [petsc-users] dm_view of high-order geometry/solution

On Mon, Jun 12, 2023 at 4:13 AM Duan Junming via petsc-users < petsc-users@mcs.anl.gov> wrote: > Dear all, > > > I am playing with dm/impls/plex/tests/ex33.c and know how to set > high-order geometry. > > Is it possible to output the final mesh to vtu, e.g. annulus example? > The problem is that

[petsc-users] dm_view of high-order geometry/solution

Dear all, I am playing with dm/impls/plex/tests/ex33.c and know how to set high-order geometry. Is it possible to output the final mesh to vtu, e.g. annulus example? Thanks! Junming