As commented under #952 (closed), the accuracy problems seem to come from I/O, not from core numerics.