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.

  Thanks,

    Matt 

Thanks!

Junming


From: knepley@gmail.com <knepley@gmail.com>
Sent: Monday, June 12, 2023 1:34:24 PM
To: Duan Junming
Cc: petsc-users@mcs.anl.gov
Subject: 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 VTK has no nice, standard way to talk about higher order polynomials. Thus, the
best way to do this is to call DMRefine() once or twice and output that mesh.

  Thanks,

    Matt
 

Thanks!

Junming



--
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener



--
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener