❯ MPI=yes ./run-benchmarks.sh polar
Checking polar benchmark...
Building polar benchmark...
Running polar benchmark...
*** results differ from reference ***
output written to opt-status-gfortran.html
❯ ./diffdumps polar_00001 polar_00001.ref
diffdumps: we welcome you
>>> reading setup from file: polar_00001 on unit 9
FT:Phantom:2021.0.0:7034209 (hydro): 27/08/2021 15:41:53.7
npart(total) = 1000000
setting isothermal sound speed^2 (polyk) = 9.4060406122873990E-004 gamma = 1.0000000000000000
time = 2.5000000000000000
--> ALLOCATING ALL ARRAYS
------------------------------------------------------------
Total memory allocated to arrays: 588.631 MB n = 1000000
------------------------------------------------------------
reading particles 1: 249493, from block 1 lims= 1- 249493
got 17 sink properties from 2 sink particles
-----------------------------------------------
ID| Mass | Racc | Macc | hsoft |
-----------------------------------------------
1| 0.70 | 0.50 | 0.0 | 0.0 |
2| 0.58 | 0.50 | 0.0 | 0.0 |
-----------------------------------------------
reading particles 1: 250808, from block 2 lims= 249494- 500301
reading particles 1: 250711, from block 3 lims= 500302- 751012
reading particles 1: 248988, from block 4 lims= 751013- 1000000
<<< finished reading (hydro) file
> sorting particles by ID...
copying arrays...
> sort completed in 1.4217190000000000 s
>>> reading setup from file: polar_00001.ref on unit 10
FT:Phantom:1.2.0:aec70bb (hydro): 18/02/2019 18:30:24.3
npart = 1000000
npart(total) = 1000000
setting isothermal sound speed^2 (polyk) = 9.4060406122873990E-004 gamma = 1.0000000000000000
time = 2.5000000000000000
reading particles 1: 1000000, from block 1 lims= 1- 1000000
WARNING! sink particle lum not found
WARNING! sink particle Teff not found
WARNING! sink particle Reff not found
WARNING! sink particle mdotloss not found
WARNING! sink particle mdotav not found
WARNING! sink particle mprev not found
got 17 sink properties from 2 sink particles
-----------------------------------------------
ID| Mass | Racc | Macc | hsoft |
-----------------------------------------------
1| 0.70 | 0.50 | 0.0 | 0.0 |
2| 0.58 | 0.50 | 0.0 | 0.0 |
-----------------------------------------------
WARNING: Particle IDs not in dump; resetting IDs
<<< finished reading (hydro) file
> sorting particles by ID...
copying arrays...
> sort completed in 0.33386000000000005 s
checking id............................OK [checked 1000000 values]
checking x.............................
x = 4.408E+00 should be 4.408E+00 ratio = 1.000E+00 err = 1.036E-05 ( 1)
x = -7.282E+00 should be -7.282E+00 ratio = 1.000E+00 err = 1.063E-07 ( 2)
x = -6.551E+00 should be -6.551E+00 ratio = 1.000E+00 err = 2.824E-07 ( 3)
x = 8.881E+00 should be 8.881E+00 ratio = 1.000E+00 err = 7.466E-07 ( 4)
x = 1.487E+00 should be 1.482E+00 ratio = 1.003E+00 err = 3.395E-03 ( 5)
x = 7.888E+00 should be 7.888E+00 ratio = 1.000E+00 err = 1.372E-06 ( 6)
x = -5.488E+00 should be -5.488E+00 ratio = 1.000E+00 err = 2.163E-06 ( 7)
x = 1.200E+01 should be 1.200E+01 ratio = 1.000E+00 err = 9.958E-08 ( 8)
x = -5.295E+00 should be -5.295E+00 ratio = 1.000E+00 err = 5.727E-07 ( 9)
x = -6.719E-02 should be -6.640E-02 ratio = 1.012E+00 err = 1.188E-02 ( 136)
x = 6.410E-02 should be 6.152E-02 ratio = 1.042E+00 err = 4.193E-02 ( 1407)
x = 6.819E-03 should be 5.893E-03 ratio = 1.157E+00 err = 1.571E-01 ( 2178)
x = 3.644E-03 should be 5.459E-03 ratio = 6.674E-01 err = 3.326E-01 ( 10057)
x = 1.851E-04 should be -8.559E-04 ratio =-2.163E-01 err = 1.216E+00 ( 48273)
x = -1.951E-03 should be -4.143E-05 ratio = 4.710E+01 err = 4.610E+01 ( 503279)
FAILED [on 1000000 of 1000000 values, max err = 4.610E+01, L2 err = 4.448E-05 tol = 2.225-308]
checking y.............................
y = -2.844E+00 should be -2.844E+00 ratio = 1.000E+00 err = 9.836E-06 ( 1)
y = 1.092E+01 should be 1.092E+01 ratio = 1.000E+00 err = 1.554E-07 ( 2)
y = -9.386E+00 should be -9.386E+00 ratio = 1.000E+00 err = 2.339E-07 ( 3)
y = 5.582E-01 should be 5.582E-01 ratio = 1.000E+00 err = 2.031E-06 ( 4)
y = -1.116E+00 should be -1.116E+00 ratio = 1.000E+00 err = 4.976E-04 ( 5)
y = -2.317E+00 should be -2.317E+00 ratio = 1.000E+00 err = 7.427E-07 ( 6)
y = -4.900E+00 should be -4.900E+00 ratio = 1.000E+00 err = 2.127E-06 ( 7)
y = 3.729E+00 should be 3.729E+00 ratio = 1.000E+00 err = 1.133E-07 ( 8)
y = 8.482E+00 should be 8.482E+00 ratio = 1.000E+00 err = 6.554E-07 ( 9)
y = 1.733E+00 should be 1.727E+00 ratio = 1.003E+00 err = 3.175E-03 ( 29)
y = 2.590E-03 should be 4.097E-04 ratio = 6.323E+00 err = 5.323E+00 ( 76)
y = 5.372E-04 should be 1.136E-05 ratio = 4.730E+01 err = 4.630E+01 ( 37822)
FAILED [on 1000000 of 1000000 values, max err = 4.630E+01, L2 err = 4.609E-05 tol = 2.225-308]
checking z.............................
z = 2.162E+00 should be 2.162E+00 ratio = 1.000E+00 err = 1.503E-05 ( 1)
z = -2.757E+00 should be -2.757E+00 ratio = 1.000E+00 err = 3.285E-07 ( 2)
z = -4.439E+00 should be -4.439E+00 ratio = 1.000E+00 err = 8.068E-08 ( 3)
z = 4.390E+00 should be 4.390E+00 ratio = 1.000E+00 err = 3.756E-07 ( 4)
z = 4.642E-01 should be 4.620E-01 ratio = 1.005E+00 err = 4.694E-03 ( 5)
z = 3.111E+00 should be 3.111E+00 ratio = 1.000E+00 err = 8.023E-07 ( 6)
z = -3.012E+00 should be -3.012E+00 ratio = 1.000E+00 err = 9.185E-07 ( 7)
z = 7.950E+00 should be 7.950E+00 ratio = 1.000E+00 err = 7.632E-08 ( 8)
z = -1.196E+00 should be -1.196E+00 ratio = 1.000E+00 err = 1.517E-06 ( 9)
z = 5.096E-02 should be 4.944E-02 ratio = 1.031E+00 err = 3.083E-02 ( 93)
z = -1.413E-02 should be -1.312E-02 ratio = 1.076E+00 err = 7.643E-02 ( 743)
z = 5.615E-03 should be 3.715E-03 ratio = 1.512E+00 err = 5.117E-01 ( 2453)
z = -1.073E-03 should be -1.283E-05 ratio = 8.363E+01 err = 8.263E+01 ( 9178)
FAILED [on 1000000 of 1000000 values, max err = 8.357E+01, L2 err = 2.808E-05 tol = 2.225-308]
checking h.............................
h = 1.210E-01 should be 1.210E-01 ratio = 1.000E+00 err = 4.927E-07 ( 1)
h = 1.583E-01 should be 1.583E-01 ratio = 1.000E+00 err = 9.411E-08 ( 2)
h = 1.144E-01 should be 1.144E-01 ratio = 1.000E+00 err = 6.513E-08 ( 4)
h = 4.206E-02 should be 4.208E-02 ratio = 9.994E-01 err = 5.547E-04 ( 5)
h = 1.093E-01 should be 1.093E-01 ratio = 1.000E+00 err = 1.363E-07 ( 6)
h = 1.024E-01 should be 1.024E-01 ratio = 1.000E+00 err = 4.364E-07 ( 7)
h = 1.383E-01 should be 1.383E-01 ratio = 1.000E+00 err = 3.232E-07 ( 9)
h = 5.619E-02 should be 5.620E-02 ratio = 9.999E-01 err = 1.031E-04 ( 11)
h = 1.869E-01 should be 1.869E-01 ratio = 1.000E+00 err = 7.975E-08 ( 12)
FAILED [on 811384 of 1000000 values, max err = 3.338E-03, L2 err = 7.098E-06 tol = 2.225-308]
checking vx............................
vx = -1.943E-01 should be -1.944E-01 ratio = 9.998E-01 err = 2.487E-04 ( 1)
vx = 2.178E-01 should be 2.178E-01 ratio = 1.000E+00 err = 2.600E-06 ( 2)
vx = -2.326E-01 should be -2.326E-01 ratio = 1.000E+00 err = 8.443E-06 ( 3)
vx = 3.162E-02 should be 3.161E-02 ratio = 1.000E+00 err = 2.323E-04 ( 4)
vx = -2.241E-01 should be -2.296E-01 ratio = 9.761E-01 err = 2.394E-02 ( 5)
vx = -7.021E-02 should be -7.022E-02 ratio = 9.998E-01 err = 1.708E-04 ( 6)
vx = -2.375E-01 should be -2.375E-01 ratio = 1.000E+00 err = 5.086E-05 ( 7)
vx = 8.716E-02 should be 8.716E-02 ratio = 1.000E+00 err = 1.563E-05 ( 8)
vx = 2.615E-01 should be 2.615E-01 ratio = 1.000E+00 err = 9.407E-06 ( 9)
vx = -8.935E-03 should be -1.248E-02 ratio = 7.158E-01 err = 2.842E-01 ( 39)
vx = 2.598E-03 should be 1.622E-04 ratio = 1.602E+01 err = 1.502E+01 ( 11355)
vx = -3.238E-04 should be -4.611E-06 ratio = 7.023E+01 err = 6.923E+01 ( 114103)
vx = -1.266E-03 should be -1.006E-06 ratio = 1.259E+03 err = 1.258E+03 ( 853245)
FAILED [on 1000000 of 1000000 values, max err = 1.258E+03, L2 err = 9.764E-04 tol = 2.225-308]
checking vy............................
vy = -4.104E-01 should be -4.103E-01 ratio = 1.000E+00 err = 8.467E-05 ( 1)
vy = 1.733E-01 should be 1.733E-01 ratio = 1.000E+00 err = 1.080E-05 ( 2)
vy = 2.026E-01 should be 2.026E-01 ratio = 1.000E+00 err = 1.060E-05 ( 3)
vy = -3.543E-01 should be -3.543E-01 ratio = 1.000E+00 err = 4.344E-06 ( 4)
vy = -6.705E-01 should be -6.685E-01 ratio = 1.003E+00 err = 2.852E-03 ( 5)
vy = -3.662E-01 should be -3.662E-01 ratio = 1.000E+00 err = 4.979E-06 ( 6)
vy = 3.196E-01 should be 3.196E-01 ratio = 1.000E+00 err = 3.329E-05 ( 7)
vy = -2.777E-01 should be -2.777E-01 ratio = 1.000E+00 err = 1.902E-06 ( 8)
vy = 1.845E-01 should be 1.845E-01 ratio = 1.000E+00 err = 3.385E-05 ( 9)
vy = -2.027E-01 should be -2.085E-01 ratio = 9.725E-01 err = 2.747E-02 ( 29)
vy = 1.670E-04 should be 1.322E-04 ratio = 1.264E+00 err = 2.638E-01 ( 197)
vy = 3.260E-05 should be -4.006E-05 ratio =-8.139E-01 err = 1.814E+00 ( 1138)
vy = -6.021E-04 should be -1.536E-05 ratio = 3.919E+01 err = 3.819E+01 ( 3077)
vy = 2.586E-03 should be -6.401E-06 ratio =-4.041E+02 err = 4.051E+02 ( 220900)
FAILED [on 1000000 of 1000000 values, max err = 4.051E+02, L2 err = 1.029E-03 tol = 2.225-308]
checking vz............................
vz = -1.505E-01 should be -1.505E-01 ratio = 9.998E-01 err = 2.418E-04 ( 1)
vz = 1.264E-01 should be 1.264E-01 ratio = 1.000E+00 err = 8.866E-06 ( 2)
vz = -8.027E-02 should be -8.028E-02 ratio = 1.000E+00 err = 1.358E-05 ( 3)
vz = -3.426E-02 should be -3.426E-02 ratio = 1.000E+00 err = 3.592E-05 ( 4)
vz = -2.210E-01 should be -2.233E-01 ratio = 9.899E-01 err = 1.012E-02 ( 5)
vz = -8.032E-02 should be -8.033E-02 ratio = 1.000E+00 err = 3.317E-05 ( 6)
vz = -7.645E-02 should be -7.646E-02 ratio = 9.999E-01 err = 8.577E-05 ( 7)
vz = 6.677E-03 should be 6.677E-03 ratio = 1.000E+00 err = 7.322E-05 ( 8)
vz = 1.470E-01 should be 1.470E-01 ratio = 1.000E+00 err = 1.757E-05 ( 9)
vz = 2.507E-02 should be 2.412E-02 ratio = 1.040E+00 err = 3.965E-02 ( 305)
vz = -4.960E-03 should be -5.784E-03 ratio = 8.575E-01 err = 1.425E-01 ( 489)
vz = 7.926E-04 should be -1.636E-04 ratio =-4.845E+00 err = 5.845E+00 ( 3030)
vz = 7.751E-04 should be -6.237E-05 ratio =-1.243E+01 err = 1.343E+01 ( 267966)
vz = 4.225E-04 should be -9.144E-06 ratio =-4.621E+01 err = 4.721E+01 ( 560427)
FAILED [on 1000000 of 1000000 values, max err = 5.821E+01, L2 err = 8.758E-04 tol = 2.225-308]
diffdumps: we wish you a pleasant journey
particle IDs differ 0 times
positions differ 1000000 times
smoothing lengths differ 811384 times
velocities differ 1000000 times
MAX RMS ERROR: 1.0287E-03
FILES DIFFER
To specify a tolerance, use:
diffdumps polar_00001 polar_00001.ref 1.e-15
