Astropy checksum verification fails for fits produced by EleFits

[kabasset]

Mirror of https://euclid.roe.ac.uk/issues/26018

I'm trying to add the checksum to the fits file using EleFits, but I'm encountering troubles. Here is a minimal example:

auto obj = Euclid::Fits::MefFile ("myfile.fits", Euclid::Fits::FileMode::Create);
std::vector<long> naxes {3, 2};
auto shape = Euclid::Fits::Position<-1> (naxes.begin(), naxes.end());
std::vector<double> datavec {1, 2, 3, 4, 5, 6};
Euclid::Fits::PtrRaster<double, -1> ptrRaster(shape, datavec.data());
obj.primary().updateShape<double, -1>(shape);
obj.primary().writeRaster(ptrRaster);
obj.primary().update_checksums();

The file is written and it contains the fields CHECKSUM and DATASUM in the header. But when I open it, astropy.io.fits raises the following warning:

WARNING: Checksum verification failed for HDU ('PRIMARY', 1).

The EleFits version is 5.3.0.
Could you please help me with this?
Thanks
Emiliano