MPI parallelisation is OFF _ \ | | ___| _ \ | | | | __ \ _` | __ \ __| _ \ __ `__ \ \___ \ | | | | ___/ | | | ( | | | | ( | | | | | ___/ ___ | _| _| |_|\__,_|_| _|\__|\___/ _| _| _|_____/ _| _| _| ___ , _ |) (| \ _, o _ |\ /|/ \ ,_ o _ _ |/\_| | | |/ | /|/| | |/ |/ |__// | | / |/ \/ \/|/ (\__/ \/|_/ | |_/|/|_/|_/ | |/|/\__/|_/ (| Version: 1.2.0 phantom_version not found; no git information to display Using NATIVE inverse sqrt ERROR! compile settings: -DDUST currently not compatible with magnetic fields (-DMHD) 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(4,1000000) 15.259 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,1000000) 30.518 MB real(8) Bxyz(3,1000000) 22.888 MB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B 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,0) 0.000 B 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,1000000) 30.518 MB real(4) divBsumm(1000000) 3.815 MB 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,1000000) 30.518 MB real(8) dustproppred(4,0) 0.000 B 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.370 GB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) checking wab0 = wab(0).................OK [max err = 0.000E+00, tol = 2.225-308] checking gradh0 = -3.*wab(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking dphidh0 = dphidh(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking radkern2 = radkern*radkern....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == wkern...........OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == grkern..........OK [max err = 0.000E+00, tol = 2.225-308] checking w gradient equal to gradw.....OK [max err = 1.499E-07, tol = 2.000E-05] checking potential gradient = force....OK [max err = 6.722E-08, tol = 2.000E-05] checking get_kernel_grav1 == wkern.....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel_grav1 == grkern....OK [max err = 0.000E+00, tol = 2.225-308] checking dphi/dh = phi - q*dphi/dq.....OK [max err = 8.105E-15, tol = 2.000E-07] --> TESTING LINKLIST / NEIGHBOUR FINDING ----------------------------------------------------------------------------------------- random seed = -43587 32768 particles set in uniform random distribution ----------------------------------------------------------------------------------------- Test 1 of 2: building linked list... checking ncells > 0....................OK [got T should be T] checking neighbour number.............. checking nneigh (cached)...............OK [checked 32768 values] checking nneigh (no cache).............OK [checked 32768 of 32768 values] max nneigh = 8945 min nneigh = 1 mean = 2206.93 Test 2 of 2: building linked list... checking ncells > 0....................OK [got T should be T] Test 3: building linked list... particles in a line in x direction checking ncells > 0....................OK [got T should be T] particles in a line in y direction checking ncells > 0....................OK [got T should be T] particles in a line in z direction checking ncells > 0....................OK [got T should be T] <-- LINKLIST TEST COMPLETE --> TESTING DERIVS MODULE ---------------- particles set on 100 x 100 x 100 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 1.000E-02 dy: 1.000E-02 dz: 1.000E-02 ----------------------------------------------------------------------------------------- thread 0 npart = 1000000 hfact = 1.2 --> testing Hydro derivatives trial neigh mean : 367.52, real neigh mean = 57.00 ratio try/act= 6.45 trial neigh max : 548, max real neigh = 57 n neighbour calls : 155520, mean per part = 1.156 n density calcs : 2000000, mean per part = 2.000 dtmin = 1.0483190326299730E-003 dtmax = 1.9552730199381144E-002 dtmax/dtmin = 18.651507404503516 dtcour/dtf = 0.16444910687077402 (wall: total = 2.38s force = 1.12s dens = 1.12s link = 0.12s) (cpu : total = 35.38s force = 17.35s dens = 17.04s link = 0.98s) (cpu/wall: total = 14.90 force = 15.42 dens = 15.15 link = 7.87) completed in 2.38 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking mean nneigh...................OK [max err = 0.000E+00, tol = 2.225-308] checking max nneigh....................OK [got 57 should be 57] checking n density calcs...............OK [is 2000000 should be 2000000] checking total nneigh..................OK [is 57000000 should be 57000000] --> testing artificial viscosity terms (individual alpha) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 5.1214812385350550E-004 dtmax = 5.1214812385350853E-004 dtmax/dtmin = 0.99999999999987588 dtcour/dtf = 1.6818595648850399E-002 (wall: total = 1.75s force = 1.12s dens = 0.50s link = 0.12s) (cpu : total = 26.81s force = 17.08s dens = 8.82s link = 0.92s) (cpu/wall: total = 15.32 force = 15.18 dens = 17.63 link = 7.32) completed in 1.75 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.669E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking art. visc force(x)............OK [max err = 2.705E-02, tol = 3.200E-02] checking art. visc force(y)............OK [max err = 1.147E-02, tol = 2.400E-02] checking art. visc force(z)............OK [max err = 9.451E-03, tol = 2.400E-02] --> testing ddivv/dt in Cullen & Dehnen switch trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 completed in 0.75 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.689E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking alphaloc......................OK [max err = 3.421E-04, tol = 3.500E-04] --> testing physical viscosity terms (two first derivatives) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 0.29999999999999999 dtmax = 0.29999999999999999 dtmax/dtmin = 0.99999999999999978 dtcour/dtf = 447.73485835607624 WARNING! force: viscosity constraining Courant timestep (wall: total = 1.75s force = 1.12s dens = 0.50s link = 0.12s) (cpu : total = 27.35s force = 17.45s dens = 8.99s link = 0.91s) (cpu/wall: total = 15.63 force = 15.52 dens = 17.97 link = 7.28) completed in 1.75 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 2.015E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvxdx.........................OK [max err = 8.214E-04, tol = 1.700E-03] checking dvxdy.........................OK [max err = 1.636E-16, tol = 2.500E-15] checking dvxdz.........................OK [max err = 1.636E-16, tol = 2.500E-15] checking dvydx.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvydy.........................OK [max err = 2.561E-16, tol = 2.500E-15] checking dvydz.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvzdx.........................OK [max err = 4.699E-17, tol = 2.500E-15] checking dvzdy.........................OK [max err = 6.568E-04, tol = 1.500E-03] checking dvzdz.........................OK [max err = 4.498E-17, tol = 2.500E-15] checking viscous force(x)..............OK [max err = 3.949E-02, tol = 4.000E-02] checking viscous force(y)..............OK [max err = 2.994E-02, tol = 3.000E-02] checking viscous force(z)..............OK [max err = 2.730E-02, tol = 3.100E-02] checking mean nneigh...................OK [max err = 0.000E+00, tol = 2.225-308] checking max nneigh....................OK [got 57 should be 57] checking n density calcs...............OK [is 1000000 should be 1000000] checking total nneigh..................OK [is 57000000 should be 57000000] checking \sum v.dv/dt + du/dt = 0......OK [max err = 4.500E-07, tol = 1.520E-06] checking du/dt >= 0 for all particles..OK [got T should be T] --> testing dust evolution terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 1.0483190357062076E-003 dtmax = 1.9552730256757554E-002 dtmax/dtmin = 18.651507404503516 dtcour/dtf = 0.15783470237374561 (wall: total = 3.00s force = 2.38s dens = 0.50s link = 0.12s) (cpu : total = 46.24s force = 36.42s dens = 8.91s link = 0.91s) (cpu/wall: total = 15.41 force = 15.34 dens = 17.81 link = 7.29) completed in 3.00 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking deps/dt.......................OK [max err = 1.826E-05, tol = 4.000E-05] checking du/dt.........................OK [max err = 4.482E-04, tol = 5.000E-04] checking deltavx.......................OK [max err = 2.282E-05, tol = 2.300E-05] checking energy conservation (dE=0)....OK [max err = 6.073E-15, tol = 6.500E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] checking dust mass conservation........OK [max err = 9.108E-16, tol = 1.000E-15] --> testing MHD derivatives (using B/rho directly) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 4.0724710578236395E-004 dtmax = 3.6278231850391555E-003 dtmax/dtmin = 8.9081619820701672 dtcour/dtf = 0.26984546408093529 (wall: total = 2.38s force = 1.75s dens = 0.50s link = 0.12s) (cpu : total = 37.36s force = 27.60s dens = 8.84s link = 0.91s) (cpu/wall: total = 15.73 force = 15.77 dens = 17.68 link = 7.30) completed in 2.38 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB (symm)...................OK [max err = 8.214E-04, tol = 2.000E-03] checking dBx/dt........................OK [max err = 2.952E-04, tol = 2.000E-03] checking dBy/dt........................OK [max err = 1.161E-03, tol = 2.000E-03] checking dBz/dt........................OK [max err = 2.482E-04, tol = 2.000E-02] checking mhd force(x)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(y)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(z)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking div B (diff)..................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlB(x)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(y)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(z)......................OK [max err = 8.219E-04, tol = 1.000E-03] --> testing artificial resistivity terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 4.5348182413328494E-004 dtmax = 3.1672007673647666E-003 dtmax/dtmin = 6.9841845886934921 dtcour/dtf = 0.27809867308719438 (wall: total = 2.38s force = 1.62s dens = 0.62s link = 0.12s) (cpu : total = 37.14s force = 27.38s dens = 8.84s link = 0.92s) (cpu/wall: total = 15.64 force = 16.85 dens = 14.14 link = 7.33) completed in 2.38 s checking dBx/dt (resist)...............OK [max err = 0.000E+00, tol = 3.700E-02] checking dBy/dt (resist)...............OK [max err = 0.000E+00, tol = 3.400E-02] checking dBz/dt (resist)...............OK [max err = 0.000E+00, tol = 2.200E-01] checking \sum du/dt + B.dB/dt = 0......OK [max err = 0.000E+00, tol = 2.700E-03] --> testing div B cleaning terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 1.9494471122519190E-004 dtmax = 3.8121849255141557E-004 dtmax/dtmin = 1.9555210816205104 dtcour/dtf = 0.11955025015767445 (wall: total = 2.38s force = 1.75s dens = 0.50s link = 0.12s) (cpu : total = 37.45s force = 27.70s dens = 8.83s link = 0.91s) (cpu/wall: total = 15.77 force = 15.83 dens = 17.67 link = 7.32) completed in 2.38 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradpsi_x.....................OK [max err = 4.261E-04, tol = 8.500E-04] checking gradpsi_y.....................OK [max err = 8.500E-04, tol = 8.500E-04] checking gradpsi_z.....................OK [max err = 3.729E-04, tol = 2.000E-03] --> testing Hydro derivs in setup with density contrast ---------------- particles set on 500 x 500 x 500 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-03 dy: 2.000E-03 dz: 2.000E-03 ----------------------------------------------------------------------------------------- ---------------- particles set on 500 x 500 x 500 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-03 dy: 2.000E-03 dz: 2.000E-03 ----------------------------------------------------------------------------------------- ---------------- particles set on 50 x 50 x 50 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-02 dy: 2.000E-02 dz: 2.000E-02 ----------------------------------------------------------------------------------------- thread 0 npart = 648432 in blob = 523984 to test = 73824 WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 430.91, real neigh mean = 57.47 ratio try/act= 7.50 trial neigh max : 88189, max real neigh = 988 n neighbour calls : 124015, mean per part = 1.191 n density calcs : 1452028, mean per part = 2.239 dtmin = 2.8817634064656154E-004 dtmax = 2.1540775519298411E-002 dtmax/dtmin = 74.748591334610850 dtcour/dtf = 1.3274325079742730 (wall: total = 2.62s force = 1.38s dens = 1.12s link = 0.12s) (cpu : total = 40.64s force = 21.60s dens = 18.33s link = 0.70s) (cpu/wall: total = 15.48 force = 15.71 dens = 16.29 link = 5.63) completed in 2.62 s checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] checking mean nneigh...................OK [max err = 0.000E+00, tol = 1.000E-17] checking max nneigh....................OK [got 988 should be 988] checking total nneigh..................OK [is 37263216 should be 37263216] <-- DERIVS TEST COMPLETE --> TESTING EQUATION OF STATE MODULE --> testing equation of state initialisation checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] Reading output_DE_z0.00x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading opacs.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading /opacs.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] Reading helm_data.tab in /home/dprice/phantom-nightly/phantom/data/eos/helmholtz/ (from PHANTOM_DIR setting) checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] --> testing barotropic equation of state Barotropic eq of state: cs_ld = 3.162E-01 code units = 3.643E+04 cm/s Barotropic eq of state: cs = 3.162E-01 code units = 3.643E+04 cm/s Barotropic eq of state: rhocritT == rhoT = 2.514E-04 code units = 5.000E-19 g/cm^3 Barotropic eq of state: rhocrit0 == rho0 = 5.027E-04 code units = 1.000E-18 g/cm^3 Barotropic eq of state: rhocrit1 == rho1 = 5.027E+00 code units = 1.000E-14 g/cm^3 Barotropic eq of state: rhocrit2 == rho2 = 5.027E+04 code units = 1.000E-10 g/cm^3 Barotropic eq of state: rhocrit3 == rho3 = 5.027E+11 code units = 1.000E-03 g/cm^3 Barotropic eq of state: Barotropic eq of state: P = cs_ld*rho for rho/(g/cm^3) < rhoT Barotropic eq of state: P = 10**(log10(cs_bg**2) + M*(log10(rhoT)-log10(rho))) for rhoT <= rho/(g/cm^3) < rho0 Barotropic eq of state: P = cs*rho for rho0 <= rho/(g/cm^3) < rho1 Barotropic eq of state: P = cs*rho1*(rho /rho1)^ 1.400 for rho1 <= rho/(g/cm^3) < rho2 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^ 1.400*(rho /rho2)^ 1.100 for rho2 <= rho/(g/cm^3) < rho3 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^ 1.400*(rho3/rho2)^ 1.100*(rho /rho3)^ 1.667 for rho3 <= rho/(g/cm^3) checking equation of state is continuous......... checking p/rho is continuous...........OK [checked 4999 of 5000 values] <-- EQUATION OF STATE TEST COMPLETE --> TESTING COOLING MODULE --> testing find_in_table routine checking table(i) < val < table(i+1)... <-- COOLING TEST COMPLETE --> TESTING KDTREE <-- KDTREE TEST COMPLETE --> SKIPPING SELF-GRAVITY TESTS (need -DGRAVITY) --> 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 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 : 5438.64, real neigh mean = 50.29 ratio try/act= 108.15 trial neigh max : 12675, max real neigh = 123 n neighbour calls : 46204, mean per part = 1.258 n density calcs : 785537, mean per part = 4.387 dtmin = 3.5904894719452420E-004 dtmax = 3.2977318517523487E-003 dtmax/dtmin = 9.1846303338834296 dtcour/dtf = 63904.867021328188 (wall: total = 5.38s force = 2.62s dens = 2.75s link = 0.00s) (cpu : total = 85.65s force = 43.12s dens = 42.21s link = 0.32s) (cpu/wall: total = 15.94 force = 16.43 dens = 15.35 link = ---) checking acceleration from drag conserves momentum(x)......OK [max err = 4.290E-11, tol = 7.000E-07] checking acceleration from drag conserves momentum(y)......OK [max err = 2.012E-10, tol = 7.000E-07] checking acceleration from drag conserves momentum(z)......OK [max err = 3.152E-10, tol = 7.000E-07] checking acceleration from drag conserves energy...........OK [max err = 1.080E-07, tol = 7.000E-07] --> testing DUSTYBOX checking dust velocities match exact solution....OK [max err = 9.090E-05, tol = 1.000E-04] checking dust accel matches exact solution.......OK [max err = 2.848E-03, tol = 3.000E-03] checking gas velocities match exact solution.....OK [max err = 9.090E-05, tol = 1.000E-04] checking gas accel matches exact solution........OK [max err = 2.848E-03, tol = 3.000E-03] checking kinetic energy decay matches exact......OK [max err = 1.302E-05, tol = 1.000E-04] --> 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.430E-03, tol = 2.600E-03] <-- DUST TEST COMPLETE --> SKIPPING DUSTGROWTH TEST (REQUIRES -DDUST -DDUSTGROWTH) --> SKIPPING NON-IDEAL MHD TEST (REQUIRES -DNONIDEALMHD) --> TESTING READ/WRITE from dump file --------> TIME = 20.00 : full dump written to file test.dump <-------- >>> reading setup from file: test.dump on unit 9 FT:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 05:41:38.4 n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.3000000000000000 External field found, Bext = 124.00000000000000 125.00000000000000 126.00000000000000 ERROR reading values to verify conservation laws. Resetting initial values. ERROR reading grain size/density from file header time = 20.000000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(2,1001) 7.820 KB real(4) divcurlv(4,1001) 15.641 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,1001) 31.281 KB real(8) Bxyz(3,1001) 23.461 KB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1001) 172.047 KB real(8) dustevol(11,1001) 86.023 KB real(8) ddustevol(11,1001) 86.023 KB real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,11,1001) 258.070 KB 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,1001) 31.281 KB real(8) dBevol(4,1001) 31.281 KB real(4) divBsumm(1001) 3.910 KB real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(11,1001) 86.023 KB real(8) Bpred(4,1001) 31.281 KB real(8) dustproppred(4,0) 0.000 B integer(1) iphase(1001) 1001.000 B integer(1) iphase_soa(1001) 1001.000 B real(4) gradh(1,1001) 3.910 KB real(8) tstop(22,1001) 172.047 KB integer(4) ll(1001) 3.910 KB integer(4) ibelong(1001) 3.910 KB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,2003) 15.648 KB integer(4) inodeparts(1001) 3.910 KB real(8) xyzh_swap(1001,4) 31.281 KB integer(4) inodeparts_swap(1001) 3.910 KB integer(1) iphase_swap(1001) 1001.000 B integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) cellatid(2003) 7.824 KB integer(4) ifirstincell(2003) 7.824 KB kdnode nodeglobal(2003) 101.715 KB kdnode node(2003) 101.715 KB integer(4) nodemap(2003) 7.824 KB --------------------------------------------------------- Total memory allocated to arrays: 1.540 MB --------------------------------------------------------- reading particles 1: 1001, from block 1 lims= 1- 1001 got 11 sink properties from 10 sink particles ----------------------------------------------- ID| Mass | Racc | Macc | hsoft | ----------------------------------------------- 1| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 2| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 3| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 4| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 5| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 6| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 7| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 8| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 9| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 10| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | ----------------------------------------------- <<< finished reading (MHD) file --> checking header checking read without errors...........OK [got 0 should be 0] checking time..........................OK [max err = 0.000E+00, tol = 2.225-308] checking hfact.........................OK [max err = 0.000E+00, tol = 2.225-308] checking gamma.........................OK [max err = 0.000E+00, tol = 2.225-308] checking polyk.........................OK [max err = 0.000E+00, tol = 2.225-308] checking m(gas)........................OK [max err = 0.000E+00, tol = 2.225-308] checking m(dust).......................OK [max err = 0.000E+00, tol = 2.225-308] checking npartoftype(1)................OK [got 991 should be 991] checking npartoftype(2)................OK [got 10 should be 10] checking umass.........................OK [max err = 0.000E+00, tol = 2.225-308] checking udist.........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bexty.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextz.........................OK [max err = 0.000E+00, tol = 2.225-308] --> checking particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking vx............................OK [max err = 0.000E+00, tol = 2.225-308] checking vy............................OK [max err = 0.000E+00, tol = 2.225-308] checking vz............................OK [max err = 0.000E+00, tol = 2.225-308] checking u.............................OK [max err = 0.000E+00, tol = 2.225-308] checking alpha.........................OK [max err = 0.000E+00, tol = 1.175E-38] checking Bx............................OK [max err = 0.000E+00, tol = 1.175E-38] checking By............................OK [max err = 0.000E+00, tol = 1.175E-38] checking Bz............................OK [max err = 0.000E+00, tol = 1.175E-38] checking psi...........................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking particle type 1...............OK [checked 991 values] checking particle type 2...............OK [checked 10 values] --> checking sink particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking m.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking hsoft.........................OK [max err = 0.000E+00, tol = 2.225-308] checking maccreted.....................OK [max err = 0.000E+00, tol = 2.225-308] checking spinx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spiny.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spinz.........................OK [max err = 0.000E+00, tol = 2.225-308] checking tlast.........................OK [max err = 0.000E+00, tol = 2.225-308] checking vx............................OK [max err = 0.000E+00, tol = 2.225-308] checking vy............................OK [max err = 0.000E+00, tol = 2.225-308] checking vz............................OK [max err = 0.000E+00, tol = 2.225-308] --> checking read of single array from file x ->y z vx vy vz u dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop checking error flag....................OK [got 0 should be 0] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] --------> TIME = 20.00 : small dump written to file test.dump <-------- --> checking read_dump >>> reading setup from file: test.dump on unit 9 ERROR! read_dump: file is not a Phantom full dump checking read_dump returns is_small_dump error code........OK [got 1978 should be 1978] --> checking read_smalldump >>> reading small dump file: test.dump on unit 9 ST:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 05:41:38.6 WARNING! read_smalldump: *** VELOCITY WILL BE MISSING FROM SMALL DUMP FILES *** n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.2999999523162842 comment: hfact was 1.1230000257492065 now 1.1230000000000000 comment: tolh was 9.9999997473787516E-005 now 1.0000000000000000E-004 comment: C_cour was 0.30000001192092896 now 0.29999999999999999 External field found, Bext = 124.00000000000000 125.00000000000000 126.00000000000000 ERROR reading values to verify conservation laws. Resetting initial values. ERROR reading grain size/density from file header time = 20.000000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(2,1001) 7.820 KB real(4) divcurlv(4,1001) 15.641 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,1001) 31.281 KB real(8) Bxyz(3,1001) 23.461 KB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1001) 172.047 KB real(8) dustevol(11,1001) 86.023 KB real(8) ddustevol(11,1001) 86.023 KB real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,11,1001) 258.070 KB 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,1001) 31.281 KB real(8) dBevol(4,1001) 31.281 KB real(4) divBsumm(1001) 3.910 KB real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(11,1001) 86.023 KB real(8) Bpred(4,1001) 31.281 KB real(8) dustproppred(4,0) 0.000 B integer(1) iphase(1001) 1001.000 B integer(1) iphase_soa(1001) 1001.000 B real(4) gradh(1,1001) 3.910 KB real(8) tstop(22,1001) 172.047 KB integer(4) ll(1001) 3.910 KB integer(4) ibelong(1001) 3.910 KB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,2003) 15.648 KB integer(4) inodeparts(1001) 3.910 KB real(8) xyzh_swap(1001,4) 31.281 KB integer(4) inodeparts_swap(1001) 3.910 KB integer(1) iphase_swap(1001) 1001.000 B integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) cellatid(2003) 7.824 KB integer(4) ifirstincell(2003) 7.824 KB kdnode nodeglobal(2003) 101.715 KB kdnode node(2003) 101.715 KB integer(4) nodemap(2003) 7.824 KB --------------------------------------------------------- Total memory allocated to arrays: 1.540 MB --------------------------------------------------------- reading particles 1: 1001, from block 1 lims= 1- 1001 ERROR: missing velocity information from file WARNING: u not in file but setting u = (K*rho**(gamma-1))/(gamma-1) WARNING: alpha not found in file WARNING! sink particle velocities not found got 11 sink properties from 10 sink particles ----------------------------------------------- ID| Mass | Racc | Macc | hsoft | ----------------------------------------------- 1| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 2| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 3| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 4| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 5| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 6| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 7| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 8| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 9| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 10| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | ----------------------------------------------- WARNING! div B cleaning field (Psi) not found in Phantom dump file: assuming psi=0 <<< finished reading (MHD) file --> checking header checking read without errors...........OK [got 0 should be 0] checking time..........................OK [max err = 0.000E+00, tol = 1.192E-07] checking hfact.........................OK [max err = 2.293E-08, tol = 1.192E-07] checking gamma.........................OK [max err = 3.668E-08, tol = 1.192E-07] checking polyk.........................OK [max err = 1.445E-08, tol = 1.192E-07] checking m(gas)........................OK [max err = 0.000E+00, tol = 1.192E-07] checking m(dust).......................OK [max err = 0.000E+00, tol = 1.192E-07] checking npartoftype(1)................OK [got 991 should be 991] checking npartoftype(2)................OK [got 10 should be 10] checking umass.........................OK [max err = 0.000E+00, tol = 2.225-308] checking udist.........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bexty.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextz.........................OK [max err = 0.000E+00, tol = 2.225-308] --> checking particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking particle type 1...............OK [checked 991 values] checking particle type 2...............OK [checked 10 values] --> checking sink particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking m.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking hsoft.........................OK [max err = 0.000E+00, tol = 2.225-308] checking maccreted.....................OK [max err = 0.000E+00, tol = 2.225-308] checking spinx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spiny.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spinz.........................OK [max err = 0.000E+00, tol = 2.225-308] checking tlast.........................OK [max err = 0.000E+00, tol = 2.225-308] <-- READ/WRITE TEST COMPLETE --> 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(4,1000000) 15.259 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,1000000) 30.518 MB real(8) Bxyz(3,1000000) 22.888 MB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B 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,0) 0.000 B 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,1000000) 30.518 MB real(4) divBsumm(1000000) 3.815 MB 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,1000000) 30.518 MB real(8) dustproppred(4,0) 0.000 B 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.370 GB --------------------------------------------------------- --> TESTING STEP MODULE / boundary crossing ---------------- particles set on 50 x 50 x 50 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-02 dy: 2.000E-02 dz: 2.000E-02 ----------------------------------------------------------------------------------------- t = 0.20000000000000001 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.40000000000000002 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.60000000000000009 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.80000000000000004 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.0000000000000000 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.2000000000000000 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.3999999999999999 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.5999999999999999 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.7999999999999998 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.9999999999999998 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] <-- STEP TEST COMPLETE --> SKIPPING TEST OF IND TIMESTEP UTILS (need -DIND_TIMESTEPS) --> TESTING EXTERNAL FORCES MODULE ----------------------------------------------------------------------------------------- random seed = -43587 1000 particles set in uniform random distribution ----------------------------------------------------------------------------------------- checking star external force initialisation......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking corotate external force initialisation..OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 3.537E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 3.552E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking binary external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 5.620E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 4.658E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 5.188E-11, tol = 1.500E-03] checking prdrag external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking torus external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 2.704E-10, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 1.051E-09, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking toystar external force initialisation...OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking external B field external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] No spiral arm spheroids set. No potential to be read in. checking spiral external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.338E-05, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 8.375E-05, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 6.864E-05, tol = 1.500E-03] checking Lense-Thirring external force initialisation......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] ERROR! extern_neutronstar: Error opening ns-rdensity.tab checking neutronstar external force initialisation.........OK [got 1 should be 1] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking Einstein-prec external force initialisation.......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 1.619E-10, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 1.526E-10, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 9.963E-11, tol = 1.500E-03] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking static sinusoid external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.196E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] ERROR! externalforces: Require number of particles per star : iexternalforce = 14 checking grav. wave inspiral external force initialisation...........OK [got 2 should be 2] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking disc gravity external force initialisation........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 8.134E-10, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking corotating binary external force initialisation...OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 4.960E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 4.744E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 1.824E-06, tol = 1.500E-03] --> testing accrete_particles routine checking accreted=was_accreted for all externf...OK [checked 16 values] --> testing velocity-dependent external force solvers checking corotate external force initialisation..OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking prdrag external force initialisation....OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking Lense-Thirring external force initialisation......OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.417E-16, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.184E-16, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 3.533E-16, tol = 1.000E-14] checking Einstein-prec external force initialisation.......OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.417E-16, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.184E-16, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 3.533E-16, tol = 1.000E-14] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking corotating binary external force initialisation...OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] <-- EXTERNAL FORCE TESTS COMPLETE --> TESTING PTMASS MODULE --> testing integration of binary orbit dt for sinks = 1.2500000000000001E-002 checking potential energy..............OK [max err = 0.000E+00, tol = 2.220E-16] checking angular momentum..............OK [max err = 0.000E+00, tol = 2.220E-10] nsteps per orbit = 711 norbits = 100 checking angular momentum..............OK [max err = 6.594E-15, tol = 3.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] checking total energy..................OK [max err = 2.440E-08, tol = 3.000E-08] --> skipping circumbinary disc test (-DPERIODIC is set) --> skipping circumbinary disc test (-DPERIODIC is set) --> testing softening in sink particle binary checking potential energy..............OK [max err = 3.391E-16, tol = 4.441E-16] nsteps per orbit = 1099 norbits = 10 checking angular momentum..............OK [max err = 6.324E-15, tol = 2.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.225-308] checking total energy..................OK [max err = 1.021E-10, tol = 2.000E-09] --> testing accretion onto sink particles checking accretion flag................OK [got T should be T] checking isdead_or_accreted flag.......OK [got T should be T] checking x(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking y(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking z(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking vx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (uniform density) ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.500E-02 -2.500E-02 -2.500E-02 0.000E+00 ) by accreting 93 particles: M= 3.321E-01 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (sin) >>>>>> s t r e t c h m a p p i n g <<<<<< stretching to match density profile in r direction density at r = 0.0000000000000000 is 1.0000000000000000 total mass = 6.9604063530853081E-004 >>>>>> done ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.765E-18 6.370E-19 1.870E-18 0.000E+00 ) by accreting 280 particles: M= 1.000E+00 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] <-- PTMASS TEST COMPLETE --> TESTING GNEWTON MODULE --> testing relativistic orbit ---------- orbital parameters ---------- apocenter : 180.000 pericenter : 20.000 semi-major axis : 100.000 eccentricity : 0.800 initial velocity : 36.724 gravitational radius : 2.122 central mass : 1.000E+06 ---------------------------------------- integrating the orbit for one azimuthal period... ( C_force = 1.0000000000000000E-004 ) checking energy........................OK [max err = 9.701E-06, tol = 1.000E-03] checking angular momentum..............OK [max err = 5.863E-07, tol = 1.000E-03] checking precession angle..............OK [max err = 5.762E-07, tol = 1.000E-05] <-- GNEWTON TEST COMPLETE --> TESTING COROTATION MODULE --> testing centrifugal force with rotation about z axis checking fx=Omega^2 x..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=Omega^2 y..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] checking phi=1/2 Omega^2 R^2...........OK [max err = 2.321E-16, tol = 4.441E-16] --> testing coriolis force with rotation about z axis checking fx=2 Omega vy.................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=-2 Omega vx................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] <-- COROTATE TEST COMPLETE --> SKIPPING DISC SETUP TESTS --> TESTING GEOMETRY MODULE --> checking transforms for coords and vectors are reversible checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 8.882E-16, L2 err = 7.326E-17 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 8.882E-16, L2 err = 1.433E-16 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 4.441E-16, L2 err = 6.410E-17 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 2.538E-16, L2 err = 1.465E-16 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 4.441E-16, L2 err = 7.166E-17 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 8.882E-16, L2 err = 1.036E-16 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 8.882E-16, L2 err = 1.282E-16 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 8.882E-16, L2 err = 1.794E-16 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 8.882E-16, L2 err = 1.923E-16 tol = 1.000E-15] --> checking transforms for coords and vectors are reversible at r=0 checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 1.225E-16, L2 err = 1.500E-32 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 4.441E-16, L2 err = 9.065E-17 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 4.441E-16, L2 err = 7.166E-17 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 1.000E-15, L2 err = 1.000E-30 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] --> testing get_coord_limits matches coord_transform checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] <-- GEOMETRY TEST COMPLETE --> testing Sedov blast wave ---------------- particles set on 16 x 16 x 16 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 6.250E-02 dy: 6.250E-02 dz: 6.250E-02 ----------------------------------------------------------------------------------------- npart = 4096 particle mass = 2.4414062500000000E-004 t = 0.38059E-03 dt = 2.822E-04 (tolv), np = 4096 t = 0.66275E-03 dt = 3.979E-04 (tolv) t = 0.10606E-02 dt = 4.818E-04 (tolv) t = 0.15424E-02 dt = 5.881E-04 (tolv) t = 0.21305E-02 dt = 6.514E-04 (courant) t = 0.27819E-02 dt = 6.608E-04 (courant) t = 0.34427E-02 dt = 6.721E-04 (courant) t = 0.41147E-02 dt = 6.851E-04 (courant) t = 0.47999E-02 dt = 7.000E-04 (courant) t = 0.54998E-02 dt = 7.166E-04 (courant) t = 0.62164E-02 dt = 7.349E-04 (courant) t = 0.69513E-02 dt = 7.550E-04 (courant) t = 0.77063E-02 dt = 7.768E-04 (courant) t = 0.84831E-02 dt = 8.002E-04 (courant) t = 0.92833E-02 dt = 8.254E-04 (courant) t = 0.10109E-01 dt = 8.521E-04 (courant) t = 0.10961E-01 dt = 8.805E-04 (courant) t = 0.11841E-01 dt = 9.104E-04 (courant) t = 0.12752E-01 dt = 9.418E-04 (courant) t = 0.13693E-01 dt = 9.746E-04 (courant) t = 0.14668E-01 dt = 1.009E-03 (courant) t = 0.15677E-01 dt = 1.044E-03 (courant) t = 0.16721E-01 dt = 1.080E-03 (courant) t = 0.17801E-01 dt = 1.118E-03 (courant) t = 0.18919E-01 dt = 1.156E-03 (courant) t = 0.20074E-01 dt = 1.194E-03 (courant) t = 0.21268E-01 dt = 1.233E-03 (courant) t = 0.22501E-01 dt = 1.271E-03 (courant) t = 0.23773E-01 dt = 1.311E-03 (courant) t = 0.25084E-01 dt = 1.352E-03 (courant) t = 0.26436E-01 dt = 1.392E-03 (courant) t = 0.27828E-01 dt = 1.432E-03 (courant) t = 0.29260E-01 dt = 1.473E-03 (courant) t = 0.30732E-01 dt = 1.513E-03 (courant) t = 0.32245E-01 dt = 1.553E-03 (courant) t = 0.33799E-01 dt = 1.593E-03 (courant) t = 0.35392E-01 dt = 1.632E-03 (courant) t = 0.37023E-01 dt = 1.670E-03 (courant) t = 0.38693E-01 dt = 1.707E-03 (courant) t = 0.40400E-01 dt = 1.742E-03 (courant) t = 0.42142E-01 dt = 1.776E-03 (courant) t = 0.43918E-01 dt = 1.809E-03 (courant) t = 0.45726E-01 dt = 1.842E-03 (courant) t = 0.47569E-01 dt = 1.876E-03 (courant) t = 0.49445E-01 dt = 1.909E-03 (courant) t = 0.51353E-01 dt = 1.942E-03 (courant) t = 0.53295E-01 dt = 1.974E-03 (courant) t = 0.55269E-01 dt = 2.006E-03 (courant) t = 0.57275E-01 dt = 2.038E-03 (courant) t = 0.59313E-01 dt = 2.069E-03 (courant) t = 0.61382E-01 dt = 2.100E-03 (courant) t = 0.63482E-01 dt = 2.131E-03 (courant) t = 0.65613E-01 dt = 2.162E-03 (courant) t = 0.67774E-01 dt = 2.192E-03 (courant) t = 0.69966E-01 dt = 2.223E-03 (courant) t = 0.72189E-01 dt = 2.253E-03 (courant) t = 0.74442E-01 dt = 2.284E-03 (courant) t = 0.76726E-01 dt = 2.291E-03 (courant) t = 0.79016E-01 dt = 2.268E-03 (courant) t = 0.81284E-01 dt = 2.251E-03 (courant) t = 0.83535E-01 dt = 2.241E-03 (courant) t = 0.85776E-01 dt = 2.054E-03 (courant) t = 0.87830E-01 dt = 2.032E-03 (courant) t = 0.89863E-01 dt = 2.017E-03 (courant) t = 0.91880E-01 dt = 2.009E-03 (courant) t = 0.93890E-01 dt = 2.009E-03 (courant) t = 0.95898E-01 dt = 2.015E-03 (courant) t = 0.97913E-01 dt = 2.028E-03 (courant) t = 0.99942E-01 dt = 5.849E-05 (dtprint) t = 0.10000 dt = 1.823E-03 (courant) --------> TIME = 0.1000 : full dump written to file test002 <-------- input file test.in written successfully. Since code start: 70 timesteps, wall: 0.62s cpu: 9.4s cpu/wall: 15 Since last dump : 70 timesteps, wall: 0.62s cpu: 9.4s cpu/wall: 15 wall cpu cpu/wall frac step : 0.62s 9.34s 14.94 100.00% step (force): 0.50s 4.64s 9.27 80.00% step (dens) : 0.12s 3.53s 28.20 20.00% Etot= 9.669E-01, Ekin= 1.982E-01, Etherm= 7.686E-01, Epot= 0.000E+00 Emag= 0.000E+00 Linm= 1.410E-17, Angm= 3.001E-19 Centre of Mass = -2.033E-18, -4.716E-18, -2.602E-18 density (max)= 1.602E+00 (mean)= 1.159E+00 (max)= 9.455E+00 g/cm^3 Mgas = 1.000E+00 Mdust1 = 0.000E+00 Mdust2 = 0.000E+00 Mdust3 = 0.000E+00 Mdust4 = 0.000E+00 Mdust5 = 0.000E+00 Mdust6 = 0.000E+00 Mdust7 = 0.000E+00 Mdust8 = 0.000E+00 Mdust9 = 0.000E+00 Mdust10 = 0.000E+00 Mdust11 = 0.000E+00 alpha(max)= 1.000E+00 RMS Mach #= 6.243E-01 div B (max)= 0.000E+00, div B (mean)= 0.000E+00 h|div B|/B (max)= 0.000E+00, h|div B|/B (mean)= 0.000E+00 --> TESTING M_4 cubic kernel <-- KERNEL TEST COMPLETE checking total energy..................OK [max err = 2.332E-04, tol = 4.700E-04] checking linear momentum...............OK [max err = 1.410E-17, tol = 7.000E-15] <-- testing complete total wall time = 1 min, 36.50 s (= 9.6500E+01s) total cpu time = 16 min, 54.63 s (= 1.0146E+03s) SUMMARY OF ALL TESTS: PASSED: 137 of 137 100.0% FAILED: 0 of 137 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found