MPI parallelisation is OFF _ \ | | ___| _ \ | | | | __ \ _` | __ \ __| _ \ __ `__ \ \___ \ | | | | ___/ | | | ( | | | | ( | | | | | ___/ ___ | _| _| |_|\__,_|_| _|\__|\___/ _| _| _|_____/ _| _| _| ___ , _ |) (| \ _, o _ |\ /|/ \ ,_ o _ _ |/\_| | | |/ | /|/| | |/ |/ |__// | | / |/ \/ \/|/ (\__/ \/|_/ | |_/|/|_/|_/ | |/|/\__/|_/ (| Version: 1.2.0 phantom_version not found; no git information to display Using NATIVE inverse sqrt Running on Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50GHz 24 cpus @ 2.50 Ghz, cache size 16384 KB Running in openMP on 16 threads --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1000000) 30.518 MB real(8) xyzh_soa(1000000,4) 30.518 MB real(8) vxyzu(4,1000000) 30.518 MB real(4) alphaind(2,1000000) 7.629 MB real(4) divcurlv(1,1000000) 3.815 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,0) 0.000 B real(8) Bxyz(3,0) 0.000 B real(8) dustprop(4,1000000) 30.518 MB real(4) St(1000000) 3.815 MB real(4) csound(1000000) 3.815 MB real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1000000) 167.847 MB real(8) dustevol(11,1000000) 83.923 MB real(8) ddustevol(11,1000000) 83.923 MB real(8) ddustprop(4,1000000) 30.518 MB real(8) deltav(3,11,1000000) 251.770 MB real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(0) 0.000 B real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1000000) 30.518 MB real(8) dBevol(4,0) 0.000 B real(4) divBsumm(0) 0.000 B real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(11,1000000) 83.923 MB real(8) Bpred(4,0) 0.000 B real(8) dustproppred(4,1000000) 30.518 MB integer(1) iphase(1000000) 976.562 KB integer(1) iphase_soa(1000000) 976.562 KB real(4) gradh(1,1000000) 3.815 MB real(8) tstop(22,1000000) 167.847 MB integer(4) ll(1000000) 3.815 MB integer(4) ibelong(1000000) 3.815 MB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,1000001) 7.629 MB integer(4) inodeparts(1000000) 3.815 MB real(8) xyzh_swap(1000000,4) 30.518 MB integer(4) inodeparts_swap(1000000) 3.815 MB integer(1) iphase_swap(1000000) 976.562 KB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) cellatid(1000001) 3.815 MB integer(4) ifirstincell(1000001) 3.815 MB kdnode nodeglobal(1000001) 49.591 MB kdnode node(1000001) 49.591 MB integer(4) nodemap(1000001) 3.815 MB --------------------------------------------------------- Total memory allocated to arrays: 1.340 GB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) --> TESTING DUST MODULE --> testing drag initialisation checking drag initialisation...........OK [got 0 should be 0] checking drag initialisation...........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking deltav=0 gives Epstein drag...OK [got 1 should be 1] --> testing Epstein/Stokes drag transition checking ts is continuous into Stokes regime: rho=1.00E-19...........OK [max err = 3.172E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-18...........OK [max err = 3.172E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-17...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-16...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-15...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-14...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-13...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-12...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-11...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-10...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-09...........OK [max err = 6.244E-02, tol = 6.300E-02] checking Epstein drag formula matches non-linear solution..OK [max err = 3.865E-03, tol = 3.900E-03] --> testing DUST DRAG trial neigh mean : 8918.67, real neigh mean = 98.05 ratio try/act= 90.96 trial neigh max : 19855, max real neigh = 243 n neighbour calls : 45776, mean per part = 1.256 n density calcs : 776338, mean per part = 4.336 dtmin = 2.999301335176877E-004 dtmax = 2.505922581319789E-003 dtmax/dtmin = 8.35502105750055 dtcour/dtf = 47137.6310106911 (wall: total = 65.25s force = 27.25s dens = 37.50s link = 0.50s) (cpu : total = 1043.76s force = 436.49s dens = 599.74s link = 7.33s) (cpu/wall: total = 16.00 force = 16.02 dens = 15.99 link = 14.66) checking acceleration from drag conserves momentum(x)......OK [max err = 1.078E-10, tol = 7.000E-07] checking acceleration from drag conserves momentum(y)......OK [max err = 4.152E-10, tol = 7.000E-07] checking acceleration from drag conserves momentum(z)......OK [max err = 2.346E-09, tol = 7.000E-07] checking acceleration from drag conserves energy...........OK [max err = 4.470E-08, tol = 7.000E-07] --> testing DUSTYBOX checking dust velocities match exact solution....OK [max err = 7.268E-06, tol = 1.000E-05] checking dust accel matches exact solution.......OK [max err = 1.087E-04, tol = 3.300E-04] checking gas velocities match exact solution.....OK [max err = 7.268E-06, tol = 2.500E-05] checking gas accel matches exact solution........OK [max err = 1.087E-04, tol = 3.300E-04] checking kinetic energy decay matches exact......OK [max err = 1.037E-06, tol = 1.000E-05] --> testing DUSTYDIFFUSE checking sum(epsilon_k) = epsilon......OK [max err = 0.000E+00, tol = 1.000E-14] checking dust diffusion matches exact solution...OK [max err = 2.553E-03, tol = 2.600E-03] <-- DUST TEST COMPLETE --> TESTING DUSTGROWTH MODULE --> testing growth initialisation checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] checking growth initialisation.........OK [got 0 should be 0] --> testing GROWINGBOX ------------------ pure growth (ifrag = 0, St = const) ------------------ checking size match exact solution.....OK [max err = 1.330E-05, tol = 5.000E-04] ------------------ ifrag = 0, St=f(size) inspired from Laibe et al. (2008) ------------------ checking size match exact solution.....OK [max err = 1.039E-05, tol = 5.000E-04] ------------------ pure fragmentation (ifrag = 1, St = const) ------------------ checking size match exact solution.....OK [max err = 1.783E-05, tol = 5.000E-04] ------------------ pure fragmentation (ifrag = 2, St = const) ------------------ checking size match exact solution.....OK [max err = 1.763E-05, tol = 5.000E-04] ------------------ position based snow line ------------------ checking size match exact solution (in)..........OK [max err = 1.330E-05, tol = 5.000E-03] checking size match exact solution (out).........OK [max err = 1.706E-03, tol = 5.000E-03] ------------------ temperature based snow line ------------------ checking size match exact solution (in)..........OK [max err = 5.334E-05, tol = 5.000E-03] checking size match exact solution (out).........OK [max err = 3.115E-03, tol = 5.000E-03] --> testing INTERPOLATIONS -------- St ~ 1.6E-03 , cs = 3.16E+01 m/s -------- checking Stokes number interpolation match exact solution..OK [max err = 1.437E-04, tol = 5.000E-03] checking sound speed interpolation match exact solution....OK [max err = 3.906E-05, tol = 5.000E-03] checking size integration match exact solution...OK [max err = 1.241E-05, tol = 5.000E-03] -------- St ~ 1.6E-02 , cs = 1.00E+02 m/s -------- checking Stokes number interpolation match exact solution..OK [max err = 1.386E-04, tol = 5.000E-03] checking sound speed interpolation match exact solution....OK [max err = 2.432E-04, tol = 5.000E-03] checking size integration match exact solution...OK [max err = 1.216E-07, tol = 5.000E-03] -------- St ~ 1.6E-01 , cs = 3.16E+02 m/s -------- checking Stokes number interpolation match exact solution..OK [max err = 1.478E-04, tol = 5.000E-03] checking sound speed interpolation match exact solution....OK [max err = 2.432E-04, tol = 5.000E-03] checking size integration match exact solution...OK [max err = 4.086E-07, tol = 5.000E-03] -------- St ~ 1.6E+00 , cs = 1.00E+03 m/s -------- checking Stokes number interpolation match exact solution..OK [max err = 1.476E-04, tol = 5.000E-03] checking sound speed interpolation match exact solution....OK [max err = 2.432E-04, tol = 5.000E-03] checking size integration match exact solution...OK [max err = 5.593E-07, tol = 5.000E-03] -------- St ~ 1.6E+01 , cs = 3.16E+03 m/s -------- checking Stokes number interpolation match exact solution..OK [max err = 1.479E-04, tol = 5.000E-03] checking sound speed interpolation match exact solution....OK [max err = 2.432E-04, tol = 5.000E-03] checking size integration match exact solution...OK [max err = 5.810E-07, tol = 5.000E-03] -------- St ~ 1.6E+02 , cs = 1.00E+04 m/s -------- checking Stokes number interpolation match exact solution..OK [max err = 1.478E-04, tol = 5.000E-03] checking sound speed interpolation match exact solution....OK [max err = 2.432E-04, tol = 5.000E-03] checking size integration match exact solution...OK [max err = 5.863E-07, tol = 5.000E-03] <-- DUSTGROWTH TEST COMPLETE <-- testing complete total wall time = 11 min, 4.50 s (= 6.6450E+02s) total cpu time = 2 hr, 57 min, 0.32 s (= 1.0620E+04s) SUMMARY OF ALL TESTS: PASSED: 21 of 21 100.0% FAILED: 0 of 21 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found