The reached precision is for constructing the writer in the IO routine. I guess this should be calculated outside of the solver.