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 <jed@jedbrown.org>
Sent: Monday, June 12, 2023 2:52:48 PM
To: Matthew Knepley; Duan Junming
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] dm_view of high-order geometry/solution
 
Matthew Knepley <knepley@gmail.com> writes:

> On Mon, Jun 12, 2023 at 6:01 AM Duan Junming <junming.duan@epfl.ch> 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.