Re: [petsc-users] DMView and DMLoad

2022-03-09 Thread Sagiyama, Koki
: [petsc-users] DMView and DMLoad Dear Koki, Many thanks for your help - that was very useful! We have been able to save/load non-periodic meshes. Unfortunately, the periodic ones are still not working for us, and I attach a small example code and output. As you suggested

Re: [petsc-users] DMView and DMLoad

2022-03-08 Thread Berend van Wachem
and see if it works. Thanks, Koki *From:* Berend van Wachem *Sent:* Thursday, February 17, 2022 9:06 AM *To:* Sagiyama, Koki ; Hapla Vaclav ; PETSc users list ; Lawrence Mitchell *Subject:* Re: [petsc-users] DMView and DMLoad Dear Koki, Many thanks for your help and sorry for the slow reply. I

Re: [petsc-users] DMView and DMLoad

2022-02-24 Thread Sagiyama, Koki
sion 2.0.0` option when saving and see if it works. Thanks, Koki From: Berend van Wachem Sent: Thursday, February 17, 2022 9:06 AM To: Sagiyama, Koki ; Hapla Vaclav ; PETSc users list ; Lawrence Mitchell Subject: Re: [petsc-users] DMView and DMLoad Dear Koki, M

Re: [petsc-users] DMView and DMLoad

2022-02-17 Thread Berend van Wachem
; Lawrence Mitchell *Subject:* Re: [petsc-users] DMView and DMLoad Dear Koki, Thanks for your email. In the example of your last email DMPlexCoordinatesLoad() takes sF0 (PetscSF) as a third argument. In our code this modification does not fix the error when loading a periodic dm. Are we doing something

Re: [petsc-users] DMView and DMLoad

2021-12-07 Thread Sagiyama, Koki
-users] DMView and DMLoad Dear Koki, Thanks for your email. In the example of your last email DMPlexCoordinatesLoad() takes sF0 (PetscSF) as a third argument. In our code this modification does not fix the error when loading a periodic dm. Are we doing something wrong? I've included an example code

Re: [petsc-users] DMView and DMLoad

2021-12-06 Thread Berend van Wachem
ing topology/coordinates. Thanks, Koki -------- *From:* Berend van Wachem *Sent:* Wednesday, November 17, 2021 3:16 PM *To:* Hapla Vaclav ; PETSc users list ; Lawrence Mitchell ; Sagiyama, Koki *Subject:* Re: [petsc-users] DMView and

Re: [petsc-users] DMView and DMLoad

2021-11-18 Thread Sagiyama, Koki
5_storage_version 2.0.0" option when saving topology/coordinates. Thanks, Koki From: Berend van Wachem Sent: Wednesday, November 17, 2021 3:16 PM To: Hapla Vaclav ; PETSc users list ; Lawrence Mitchell ; Sagiyama, Koki Subject: Re: [petsc-users] DMVi

Re: [petsc-users] DMView and DMLoad

2021-11-17 Thread Berend van Wachem
Dear Vaclav, Lawrence, Koki, Thanks for your help! Following your advice and following your example (https://petsc.org/main/docs/manual/dmplex/#saving-and-loading-data-with-hdf5) we are able to save and load the DM with a wrapped Vector in h5 format (PETSC_VIEWER_HDF5_PETSC) successfully.

Re: [petsc-users] DMView and DMLoad

2021-10-01 Thread Sagiyama, Koki
accordingly. Thanks, Koki From: petsc-users on behalf of Berend van Wachem Sent: Thursday, September 30, 2021 12:02 PM To: Hapla Vaclav Cc: PETSc users list Subject: Re: [petsc-users] DMView and DMLoad *** This email originates from outside

Re: [petsc-users] DMView and DMLoad

2021-09-30 Thread Berend van Wachem
Dear Vaclav, Lawrence, following your example, we have managed to save the DM with a wrapped Vector in h5 format (PETSC_VIEWER_HDF5_PETSC) with: DMPlexTopologyView(dm, viewer); DMClone(dm, ); ... DMPlexSectionView(dm, viewer, sdm); DMGetLocalVector(sdm, ); ...

Re: [petsc-users] DMView and DMLoad

2021-09-22 Thread Hapla Vaclav
To avoid confusions here, Berend seems to be specifically demanding XDMF (PETSC_VIEWER_HDF5_XDMF). The stuff we are now working on is parallel checkpointing in our own HDF5 format (PETSC_VIEWER_HDF5_PETSC), I will make a series of MRs on this topic in the following days. For XDMF, we are

Re: [petsc-users] DMView and DMLoad

2021-09-20 Thread Patrick Sanan
Thanks for reporting that! There are still some things in the "classic" docs build that we could make more robust but it's not clear if we should devote the energy to that or to continuing to replace those processes with new ones which are better integrated with the Sphinx build. MR that should

Re: [petsc-users] DMView and DMLoad

2021-09-20 Thread Pierre Jolivet
> On 20 Sep 2021, at 3:38 PM, Lawrence Mitchell wrote: > > Dear Sergio, > > (Added petsc-users back to cc), > >> On 20 Sep 2021, at 14:08, sergio.bengoec...@ovgu.de wrote: >> >> Dear Lawrence, >> >> thanks for the HDF5 saving and loading example. >> >> In the documentation you sent >>

Re: [petsc-users] DMView and DMLoad

2021-09-20 Thread Lawrence Mitchell
Dear Sergio, (Added petsc-users back to cc), > On 20 Sep 2021, at 14:08, sergio.bengoec...@ovgu.de wrote: > > Dear Lawrence, > > thanks for the HDF5 saving and loading example. > > In the documentation you sent > (https://petsc.org/main/docs/manual/dmplex/#saving-and-loading-data-with-hdf5)

Re: [petsc-users] DMView and DMLoad

2021-09-17 Thread Lawrence Mitchell
Hi Berend, > On 14 Sep 2021, at 12:23, Matthew Knepley wrote: > > On Tue, Sep 14, 2021 at 5:15 AM Berend van Wachem > wrote: > Dear PETSc-team, > > We are trying to save and load distributed DMPlex and its associated > physical fields (created with DMCreateGlobalVector) (Uvelocity, >

Re: [petsc-users] DMView and DMLoad

2021-09-14 Thread Matthew Knepley
On Tue, Sep 14, 2021 at 5:15 AM Berend van Wachem wrote: > Dear PETSc-team, > > We are trying to save and load distributed DMPlex and its associated > physical fields (created with DMCreateGlobalVector) (Uvelocity, > VVelocity, ...) in HDF5_XDMF format. To achieve this, we do the following: >

[petsc-users] DMView and DMLoad

2021-09-14 Thread Berend van Wachem
Dear PETSc-team, We are trying to save and load distributed DMPlex and its associated physical fields (created with DMCreateGlobalVector) (Uvelocity, VVelocity, ...) in HDF5_XDMF format. To achieve this, we do the following: 1) save in the same xdmf.h5 file: DMView( DM ,

Re: [petsc-users] DMView and DMLoad

2015-04-02 Thread Matthew Knepley
I have fixed this bug. The fix is merged to next, and will make it to master soon. Thanks, Matt On Tue, Mar 10, 2015 at 5:05 PM, Matthew Knepley knep...@gmail.com wrote: On Tue, Mar 10, 2015 at 5:02 PM, Ataollah Mesgarnejad ames...@tigers.lsu.edu wrote: Thanks Matt, I changed

Re: [petsc-users] DMView and DMLoad

2015-03-10 Thread Matthew Knepley
On Mon, Mar 9, 2015 at 11:28 AM, Ataollah Mesgarnejad ames...@tigers.lsu.edu wrote: Dear all, I'm trying to save and load distributed DMPlex in HDF5 format with DMView and DMLoad. The problem is even though there are no errors while loading the DM when I look at the DM loaded from file is

Re: [petsc-users] DMView and DMLoad

2015-03-10 Thread Ataollah Mesgarnejad
Thanks Matt, I changed your file to reproduce the error. If you run it on more than one processor you will get something similar to: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Invalid argument [0]PETSC

Re: [petsc-users] DMView and DMLoad

2015-03-10 Thread Matthew Knepley
On Tue, Mar 10, 2015 at 5:02 PM, Ataollah Mesgarnejad ames...@tigers.lsu.edu wrote: Thanks Matt, I changed your file to reproduce the error. If you run it on more than one processor you will get something similar to: [0]PETSC ERROR: - Error Message

[petsc-users] DMView and DMLoad

2015-03-09 Thread Ataollah Mesgarnejad
Dear all, I'm trying to save and load distributed DMPlex in HDF5 format with DMView and DMLoad. The problem is even though there are no errors while loading the DM when I look at the DM loaded from file is not equal to the DM that was saved. Also, when loading the DM if DMSetType is called before