The code in puma/src/legini.f90 is changed as workaround.
The optimization flag "-ffast-math" generated wrong code for the old version