./ex1 -dm_view ascii:ex1.vtu cat ex1.vtu DM Object: 1 MPI processes type: da Processor [0] M 10 N 8 m 1 n 1 w 2 s 1 X range of indices: 0 10, Y range of indices: 0 8 The output as about, Is that need add some c code or use command option? Thanks.

On Apr 21, 2021, at 3:28 AM, Matthew Knepley <[email protected]> wrote:

On Tue, Apr 20, 2021 at 3:25 PM hu abel <[email protected] <mailto:[email protected]>> wrote: Hi everyone!

Could use DMDA create dm write to gmsh format file? Or vtk? I tried dm/tutorials/ex1.c ./ex1 -dm_view ascii:ex1.msh cat ex1.msh DM Object: 1 MPI processes type: da Processor [0] M 10 N 8 m 1 n 1 w 2 s 1 X range of indices: 0 10, Y range of indices: 0 8

You can write to .vtu, which should be what you want.

Thanks,

Matt

Thanks Abel

-- 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

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>

I tried ./ex1 -dm_view vtk:ex1.vtu But get nothing, I using Petsc Release Version 3.15.0, Mar 30, 2021 Configured as ./configure --prefix=$HOME/opt/petsc3 --with-zlib=1 \ --with-hypre=1 --download-hypre=yes \ --with-viennacl=1 --download-viennacl=yes \ --with-hdf5=1 --download-hdf5=yes \ --with-metis=1 --download-metis=yes \ --with-parmetis=1 --download-parmetis=yes \ --with-superlu=1 --download-superlu=yes \ --with-slepc=1 --download-slepc=yes \ --with-suitesparse=1 --download-suitesparse=yes \ --with-hpddm=1 --download-hpddm=yes \ --with-p4est=1 --download-p4est=yes \ --with-tetgen=1 --download-tetgen=yes \ --with-triangle=1 --download-triangle=yes \ --with-mumps=1 --download-mumps=yes \ --with-opengl=1 \ --with-scalapack=1 --download-scalapack=yes \ --with-debugging=no \ && make -j8 && make install Need add enable vtk in configure ? Thanks.

On Apr 21, 2021, at 3:37 AM, Matthew Knepley <[email protected]> wrote:

./ex1 -dm_view vtk:ex1.vtu

OK, Thanks abel

On Apr 21, 2021, at 4:08 AM, Matthew Knepley <[email protected]> wrote:

On Tue, Apr 20, 2021 at 3:47 PM hu abel <[email protected] <mailto:[email protected]>> wrote: I tried ./ex1 -dm_view vtk:ex1.vtu But get nothing, I using Petsc Release Version 3.15.0, Mar 30, 2021 Configured as

./configure --prefix=$HOME/opt/petsc3 --with-zlib=1 \ --with-hypre=1 --download-hypre=yes \ --with-viennacl=1 --download-viennacl=yes \ --with-hdf5=1 --download-hdf5=yes \ --with-metis=1 --download-metis=yes \ --with-parmetis=1 --download-parmetis=yes \ --with-superlu=1 --download-superlu=yes \ --with-slepc=1 --download-slepc=yes \ --with-suitesparse=1 --download-suitesparse=yes \ --with-hpddm=1 --download-hpddm=yes \ --with-p4est=1 --download-p4est=yes \ --with-tetgen=1 --download-tetgen=yes \ --with-triangle=1 --download-triangle=yes \ --with-mumps=1 --download-mumps=yes \ --with-opengl=1 \ --with-scalapack=1 --download-scalapack=yes \ --with-debugging=no \ && make -j8 && make install

Need add enable vtk in configure ?

Thanks.

Well, I know very little about VTK. Here is a working demo:

$ cd $PETSC_DIR/src/snes/tutorials $ make ex5 $ ./ex5 -dm_view vtk:ex5.vts -sol_view vtk:ex5.vts -mms 1 N: 16 error L2 4.24792e-13 inf 8.49584e-13

Thanks,

Matt

...

On Apr 21, 2021, at 3:37 AM, Matthew Knepley <[email protected] <mailto:[email protected]>> wrote:

./ex1 -dm_view vtk:ex1.vtu

-- 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

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>