[petsc-users] ISView() in PETSc 3.22

[Adrian Croucher]

hi, I am testing my (Fortran) code on PETSc 3.22 and have got it to 
build. However I am getting some unusual new behaviour when I write an 
IS to an HDF5 file using ISView().

The attached minimal example shows the issue. It creates a simple 
10-element IS and writes it to HDF5. With previous versions of PETSc 
this would give me a 10x1 dataset containing the values 0 - 9, as expected.

When I run it with PETSc 3.22 (in serial), I again get the expected 
values written on stdout, so it looks like the IS itself is correct. But 
in the HDF5 file I get a 1x3 dataset containing the values (10,1,0).

Has something changed here?

- Adrian

-- 
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email: a.croucher at auckland.ac.nz
tel: +64 (0)9 923 4611
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_is.F90
Type: text/x-fortran
Size: 817 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20241001/2431970e/attachment-0001.bin>