------- Comment #6 from dominiq at lps dot ens dot fr 2009-05-08 15:36 ------- Adding the print *, ... in the subroutines and mixing C and Fortran, I get:
310.40363 0.0000000 Function 1: -------------------- result 1: 310.403625 result 2: 310.403625 23223.025 0.0000000 Function 2: -------------------- result 1: 23223.025391 result 2: 23223.025391 142007.84 1.56250000E-02 1.56250000E-02 Function 3: -------------------- result 1: 142007.843750 result 2: 142007.828125 > I think default-real-kind on win64 could be 8 It would as weird as it was with Crays. default-real-kind is usually the shortest float available, i.e. real(4) on intel procs. -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=40070