[petsc-dev] VecView HDF5 read / write partial array

[Matthew Knepley]

On Fri, Oct 20, 2017 at 4:32 PM, Blaise A Bourdin <bourdin at lsu.edu> wrote:

> Hi,
>
> Recent versions of exodus are implemented with netcdf4, which is based on
> hdf5.
> Instead of going through the process of rewriting exodus viewers, I would
> like to use PETSc’s HDF5 I/O operations.
> Problem is that as far as I can understand, exodus / netcdf does not use
> hdf5 time steps. Instead, it writes a single vector containing all values
> of a field at all time steps.
> (i.e. if vij denotes the i^th component of v at time step j, v is written
> as a single array v11, v21, v31, . . ., v12, v22, v32, . . ., v13, v23,  .
> . .) and so on).
> Is there a way to use VecView hdf5 to read/write only entries kn:(k+1)n-1
> in an hdf5 file?


This is exactly what HDF5 does. You write a 'hyperslab' of the array
corresponding to that time step.

   Matt


>
> Blaise
>
> --
> Department of Mathematics and Center for Computation & Technology
> Louisiana State University, Baton Rouge, LA 70803, USA
> Tel. +1 (225) 578 1612, Fax  +1 (225) 578 4276 http://www.math.lsu.edu/~
> bourdin
>
>
>
>
>
>
>
>


-- 
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.caam.rice.edu/~mk51/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171020/efc2c644/attachment-0001.html>