[petsc-users] hdf5 and FILE_MODE_APPEND

[Mohamad M. Nasr-Azadani]

I know, it is getting very confusing. But this sounds very simple but got
this complicated so far.
I think you have seen my emails yesterday.
I want to write a vector (3D DMDA, structured) at each time to an hdf5
file.
But, I want to also add the coodinates to the same *.h5 file as well.
The grid is structured, so I only need 3 1-D arrays representing the
coordinates (and not 1D*1D*1D coordinates cell coordinates).
How would you go about it?

I thought about creating a "parallel" vector which has the size of 0 on all
processors and N on processor zero and then dump the coordinates into that
vector and into the same hdf5 file. No luck, since hdf5 writer does not
like zero-sized local vectors (if I could have done this, then I could use
the same hdf5 viewer and not close the file). But now, I am stuck.
I think about closing the hdf5 viewer after I dumped the field data, create
a vector using PETSC_COMM_SELF including the coordinates on all processor
and then add that vector to the end of same hdf5 file (that's why I need to
open it and use append mode). Of course, I would do the writing only on
processor zero.

I hope I did not confuse you.
Thanks,
Mohamad





On Fri, Jan 20, 2012 at 1:37 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> On Fri, Jan 20, 2012 at 15:33, Mohamad M. Nasr-Azadani <mmnasr at gmail.com>wrote:
>
>> At this point, I can't since the first time the viewer is created via
>>  PETSC_COMM_WORLD and the second time, PETSC_COMM_SELF.
>> (I am not even 100% sure if that could cause any troubles with the hdf5
>> file?)
>>
>
> Well this sounds pretty confusing. What are you actually trying to do?
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120120/c6a0441a/attachment.htm>