[petsc-users] VecView to hdf5 broken for large (complex) vectors

[Matthew Knepley]

On Tue, Apr 16, 2019 at 10:58 AM Sajid Ali <sajidsyed2021 at u.northwestern.edu>
wrote:

> Hi Matt,
>
> I tried running the same example with a smaller grid on a workstation and
> I see that for a grid size of 8192x8192 (vector write dims 67108864, 2),
> the output file has a chunk size of (16777215, 2).
>
> I want to see what chunksize is when you get the failure.


> I can’t see HDF5_INT_MAX in the spack build-log (which includes
> configure). Is there a better way to look it up?
>
master *:~/TripReceipts/RiceTomDefense2019$ find
/PETSc3/petsc/petsc-dev/include/ -name "*.h" | xargs grep PETSC_HDF5

/PETSc3/petsc/petsc-dev/include//petscviewerhdf5.h:#define
PETSC_HDF5_INT_MAX  2147483647

/PETSc3/petsc/petsc-dev/include//petscviewerhdf5.h:#define
PETSC_HDF5_INT_MIN -2147483647

/PETSc3/petsc/petsc-dev/include//petscviewerhdf5.h:  if ((a) >
PETSC_HDF5_INT_MAX) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Array
too long for HDF5");

  Matt

> [sajid at xrmlite .spack]$ cat build.out | grep "HDF"
> #define PETSC_HAVE_HDF5 1
> #define PETSC_HAVE_LIBHDF5HL_FORTRAN 1
> #define PETSC_HAVE_LIBHDF5 1
> #define PETSC_HAVE_LIBHDF5_HL 1
> #define PETSC_HAVE_LIBHDF5_FORTRAN 1
> #define PETSC_HAVE_HDF5_RELEASE_VERSION 5
> #define PETSC_HAVE_HDF5_MINOR_VERSION 10
> #define PETSC_HAVE_HDF5_MAJOR_VERSION 1
>
> Thank You,
> Sajid Ali
> Applied Physics
> Northwestern University
>


-- 
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/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190416/aef8b8be/attachment.html>