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 lines. Should I modify the code to make cgns work? Or any other examples for me to start? Thanks! Junming ________________________________ From: Jed Brown <[email protected]> Sent: Monday, June 12, 2023 2:52:48 PM To: Matthew Knepley; Duan Junming Cc: [email protected] Subject: Re: [petsc-users] dm_view of high-order geometry/solution Matthew Knepley <[email protected]> writes:
On Mon, Jun 12, 2023 at 6:01 AM Duan Junming <[email protected]> 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 format for this. It wants everything at first order.
I would recommend configuring with --download-cgns and running with -dm_view cgns:output.cgns. This format has efficient parallel IO and curved elements work for moderate order in Paraview.