SZ compression of 3D single precision data

Created by: tleger26

I’m trying to use SZ through ADIOS2 to compress some single precision data and have run into an odd issue. When passing the data in 3D matrix form, the decompressed data for some portions of the domain are significantly (orders of magnitude) different from the original data being compressed, regardless of the error bound mode or level specified. If I convert the data into a vector (1D array) first and then pass it into ADIOS2, then the SZ compression works as expected. I’ve also tried using ZFP 0.5.5 and it works as expected for data in both 3D matrix and vector form. I’ve tried version 2.4.0 and 2.5.0 of ADIOS2, along with SZ 2.1.5.3, 2.1.6.0, and 2.1.7.0, all with the same results. Has anyone else tried and/or experienced this behaviour?