MPI parallelisation is OFF _ \ | | ___| _ \ | | | | __ \ _` | __ \ __| _ \ __ `__ \ \___ \ | | | | ___/ | | | ( | | | | ( | | | | | ___/ ___ | _| _| |_|\__,_|_| _|\__|\___/ _| _| _|_____/ _| _| _| ___ , _ |) (| \ _, o _ |\ /|/ \ ,_ o _ _ |/\_| | | |/ | /|/| | |/ |/ |__// | | / |/ \/ \/|/ (\__/ \/|_/ | |_/|/|_/|_/ | |/|/\__/|_/ (| Version: 1.3.0 Compilation details: Date compiled: Fri Aug 16 10:03:12 AEST 2019 Pre-processor flags used: -DPERIODIC -DMHD -DDUST -DCONST_ARTRES -DCURLV Version details: Date this version was committed to the repository: Fri Aug 16 09:13:36 2019 +1000 Commit tag: ok-20190815-6-g45b00f3 Full commit ID: 45b00f3519e293fa0858a3d29d67f224ea413284 Locally modified files: 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,1200000) 36.621 MB real(8) xyzh_soa(1200000,4) 36.621 MB real(8) vxyzu(4,1200000) 36.621 MB real(4) alphaind(2,1200000) 9.155 MB real(4) divcurlv(4,1200000) 18.311 MB real(4) divcurlB(4,1200000) 18.311 MB real(8) Bevol(4,1200000) 36.621 MB real(8) Bxyz(3,1200000) 27.466 MB real(8) dustprop(2,0) 0.000 B real(8) dustgasprop(4,0) 0.000 B real(4) VrelVf(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1200000) 201.416 MB real(8) dustevol(11,1200000) 100.708 MB real(8) ddustevol(11,1200000) 100.708 MB real(8) ddustprop(2,0) 0.000 B real(8) deltav(3,11,1200000) 302.124 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,1200000) 36.621 MB real(8) dBevol(4,1200000) 36.621 MB real(4) divBsumm(1200000) 4.578 MB real(8) fext(3,1200000) 27.466 MB real(8) vpred(4,1200000) 36.621 MB real(8) dustpred(11,1200000) 100.708 MB real(8) Bpred(4,1200000) 36.621 MB real(8) dustproppred(2,0) 0.000 B integer(1) iphase(1200000) 1.144 MB integer(1) iphase_soa(1200000) 1.144 MB real(4) gradh(1,1200000) 4.578 MB real(8) tstop(22,1200000) 201.416 MB integer(4) ll(1200000) 4.578 MB integer(4) ibelong(1200000) 4.578 MB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,2400001) 18.311 MB integer(4) inodeparts(1200000) 4.578 MB real(8) xyzh_swap(1200000,4) 36.621 MB integer(4) inodeparts_swap(1200000) 4.578 MB integer(1) iphase_swap(1200000) 1.144 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) cellatid(2400001) 9.155 MB integer(4) ifirstincell(2400001) 9.155 MB kdnode nodeglobal(2400001) 119.019 MB kdnode node(2400001) 119.019 MB integer(4) nodemap(2400001) 9.155 MB --------------------------------------------------------- Total memory allocated to arrays: 1.783 GB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) --> TESTING M_4 cubic kernel 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.944E-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 = 5.218E-15, tol = 2.000E-07] <-- KERNEL TEST COMPLETE --> 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.048319032629973E-003 dtmax = 1.955273019938107E-002 dtmax/dtmin = 18.6515074045035 dtcour/dtf = 0.164449106870774 (wall: total = 2.50s force = 1.00s dens = 1.25s link = 0.25s) (cpu : total = 39.95s force = 16.53s dens = 20.61s link = 2.80s) (cpu/wall: total = 15.98 force = 16.53 dens = 16.49 link = 11.19) completed in 2.50 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.746E-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.121481238535053E-004 dtmax = 5.121481238535085E-004 dtmax/dtmin = 0.999999999999876 dtcour/dtf = 1.681859564884744E-002 (wall: total = 1.75s force = 1.00s dens = 0.62s link = 0.12s) (cpu : total = 26.50s force = 14.87s dens = 9.72s link = 1.91s) (cpu/wall: total = 15.14 force = 14.87 dens = 15.55 link = 15.30) 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.617E-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 = 4.360E-03, tol = 5.700E-03] checking art. visc force(y)............OK [max err = 1.334E-02, tol = 1.400E-02] checking art. visc force(z)............OK [max err = 1.209E-02, tol = 1.300E-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.62 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.300000000000000 dtmax = 0.300000000000000 dtmax/dtmin = 1.00000000000000 dtcour/dtf = 447.734858356076 WARNING! force: viscosity constraining Courant timestep (wall: total = 1.75s force = 1.00s dens = 0.62s link = 0.12s) (cpu : total = 27.33s force = 15.66s dens = 9.76s link = 1.91s) (cpu/wall: total = 15.62 force = 15.66 dens = 15.61 link = 15.27) 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.529E-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.154E-16, tol = 2.500E-15] checking dvxdz.........................OK [max err = 1.155E-16, tol = 2.500E-15] checking dvydx.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvydy.........................OK [max err = 2.244E-16, tol = 2.500E-15] checking dvydz.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvzdx.........................OK [max err = 4.134E-17, tol = 2.500E-15] checking dvzdy.........................OK [max err = 6.568E-04, tol = 1.500E-03] checking dvzdz.........................OK [max err = 4.179E-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 n density calcs...............OK [is 1000000 should be 1000000] checking mean nneigh...................OK [max err = 0.000E+00, tol = 2.225-308] checking max nneigh....................OK [got 57 should be 57] checking total nneigh..................OK [is 57000000 should be 57000000] checking \sum v.dv/dt + du/dt = 0......OK [max err = 4.501E-07, tol = 1.700E-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.048319035706207E-003 dtmax = 1.955273025675749E-002 dtmax/dtmin = 18.6515074045035 dtcour/dtf = 0.157834702373746 (wall: total = 2.88s force = 2.12s dens = 0.62s link = 0.12s) (cpu : total = 45.43s force = 33.24s dens = 9.61s link = 2.59s) (cpu/wall: total = 15.80 force = 15.64 dens = 15.37 link = 20.72) completed in 2.88 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 = 1.000E-03] checking deltavx.......................OK [max err = 2.282E-05, tol = 1.010E-03] checking energy conservation (dE=0)....OK [max err = 1.050E-15, tol = 6.500E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-15] checking dust mass conservation........OK [max err = 1.202E-16, tol = 1.100E-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.072471057823643E-004 dtmax = 3.627823185039157E-003 dtmax/dtmin = 8.90816198207016 dtcour/dtf = 0.269845464080934 (wall: total = 2.25s force = 1.50s dens = 0.62s link = 0.12s) (cpu : total = 35.77s force = 23.92s dens = 9.59s link = 2.26s) (cpu/wall: total = 15.90 force = 15.95 dens = 15.35 link = 18.09) completed in 2.25 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.534818241332850E-004 dtmax = 3.167200767364767E-003 dtmax/dtmin = 6.98418458869349 dtcour/dtf = 0.278098673087194 (wall: total = 2.25s force = 1.50s dens = 0.62s link = 0.12s) (cpu : total = 35.05s force = 23.53s dens = 9.59s link = 1.93s) (cpu/wall: total = 15.58 force = 15.69 dens = 15.35 link = 15.43) completed in 2.25 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.949447112251919E-004 dtmax = 3.812184925514154E-004 dtmax/dtmin = 1.95552108162051 dtcour/dtf = 0.119550250157675 (wall: total = 2.25s force = 1.50s dens = 0.62s link = 0.12s) (cpu : total = 35.44s force = 23.91s dens = 9.60s link = 1.93s) (cpu/wall: total = 15.75 force = 15.94 dens = 15.36 link = 15.44) completed in 2.25 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 = 9.300E-04, tol = 9.300E-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 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.881763406465615E-004 dtmax = 2.154077551929839E-002 dtmax/dtmin = 74.7485913346108 dtcour/dtf = 1.32743250797427 (wall: total = 2.62s force = 1.25s dens = 1.25s link = 0.12s) (cpu : total = 41.44s force = 19.13s dens = 20.79s link = 1.52s) (cpu/wall: total = 15.79 force = 15.30 dens = 16.63 link = 12.16) 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 = 9.520E-17, 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] Reading eos_binary_table.dat in /home/dprice/phantom-nightly/phantom/data/eos/shen/ (from PHANTOM_DIR setting) /home/dprice/phantom-nightly/phantom/data/eos/shen/eos_binary_table.dat 0 Reading eos_binary_table.dat in /home/dprice/phantom-nightly/phantom/data/eos/shen/ (from PHANTOM_DIR setting) 1 0.000000000000000E+000 2 0.158489320000000 3 0.167880400000000 4 0.177827940000000 5 0.188364910000000 6 0.199526230000000 7 0.211348900000000 8 0.223872110000000 9 0.237137370000000 10 0.251188640000000 11 0.266072510000000 12 0.281838290000000 13 0.298538260000000 14 0.316227770000000 15 0.334965440000000 16 0.354813390000000 17 0.375837410000000 18 0.398107170000000 19 0.421696510000000 20 0.446683590000000 21 0.473151260000000 22 0.501187240000000 23 0.530884450000000 24 0.562341330000000 25 0.595662150000000 26 0.630957350000000 27 0.668343920000000 28 0.707945790000000 29 0.749894210000000 30 0.794328240000000 31 0.841395150000000 32 0.891250940000000 33 0.944060880000000 34 1.00000000000000 35 1.05925370000000 36 1.12201850000000 37 1.18850220000000 38 1.25892540000000 39 1.33352140000000 40 1.41253760000000 41 1.49623570000000 42 1.58489320000000 43 1.67880400000000 44 1.77827940000000 45 1.88364910000000 46 1.99526230000000 47 2.11348900000000 48 2.23872110000000 49 2.37137370000000 50 2.51188640000000 51 2.66072510000000 52 2.81838290000000 53 2.98538260000000 54 3.16227770000000 55 3.34965440000000 56 3.54813390000000 57 3.75837410000000 58 3.98107170000000 59 4.21696510000000 60 4.46683590000000 61 4.73151260000000 62 5.01187240000000 63 5.30884450000000 64 5.62341330000000 65 5.95662150000000 66 6.30957350000000 67 6.68343920000000 68 7.07945790000000 69 7.49894210000000 70 7.94328240000000 71 8.41395150000000 72 8.91250940000000 73 9.44060880000000 74 10.0000000000000 75 10.5925370000000 76 11.2201850000000 77 11.8850220000000 78 12.5892540000000 79 13.3352140000000 80 14.1253760000000 81 14.9623570000000 82 15.8489320000000 83 16.7880400000000 84 17.7827940000000 85 18.8364910000000 86 19.9526230000000 87 21.1348900000000 88 22.3872110000000 89 23.7137370000000 90 25.1188640000000 91 26.6072510000000 92 28.1838290000000 93 29.8538260000000 94 31.6227770000000 95 33.4965440000000 96 35.4813390000000 97 37.5837410000000 98 39.8107170000000 99 42.1696510000000 100 44.6683590000000 101 47.3151260000000 102 50.1187240000000 103 53.0884450000000 104 56.2341330000000 105 59.5662150000000 106 63.0957350000000 107 66.8343920000000 108 70.7945790000000 109 74.9894210000000 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 < rhoT Barotropic eq of state: P = 10**(log10(cs_ld**2) + M*(log10(rhoT)-log10(rho))) for rhoT <= rho < rho0 Barotropic eq of state: P = cs*rho for rho0 <= rho < rho1 Barotropic eq of state: P = cs*rho1*(rho /rho1)^1.400 for rho1 <= rho < rho2 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^1.400*(rho /rho2)^1.100 for rho2 <= rho < rho3 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^1.400*(rho3/rho2)^1.100*(rho /rho3)^1.667 for rho3 <= rho 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.590489471945241E-004 dtmax = 3.297731851752349E-003 dtmax/dtmin = 9.18463033388343 dtcour/dtf = 63904.8670213282 (wall: total = 5.38s force = 2.12s dens = 3.25s link = 0.00s) (cpu : total = 87.03s force = 35.21s dens = 51.45s link = 0.37s) (cpu/wall: total = 16.19 force = 16.57 dens = 15.83 link = ---) checking acceleration from drag conserves momentum(x)......OK [max err = 1.940E-10, tol = 7.000E-07] checking acceleration from drag conserves momentum(y)......OK [max err = 1.283E-10, tol = 7.000E-07] checking acceleration from drag conserves momentum(z)......OK [max err = 3.923E-10, tol = 7.000E-07] checking acceleration from drag conserves energy...........OK [max err = 1.062E-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.3.0:45b00f3 (mhd+clean+dust+1dust) : 16/08/2019 10:04:24.9 n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.30000000000000 External field found, Bext = 124.000000000000 125.000000000000 126.000000000000 time = 20.0000000000000 --> 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(2,0) 0.000 B real(8) dustgasprop(4,0) 0.000 B real(4) VrelVf(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(2,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(2,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.536 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 = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] checking dustfrac......................OK [max err = 2.235E-08, tol = 1.192E-07] 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 dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz deltavx deltavy deltavz 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.3.0:45b00f3 (mhd+clean+dust+1dust) : 16/08/2019 10:04:25.0 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.29999995231628 External field found, Bext = 124.000000000000 125.000000000000 126.000000000000 time = 20.0000000000000 --> 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(2,0) 0.000 B real(8) dustgasprop(4,0) 0.000 B real(4) VrelVf(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(2,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(2,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,1200000) 36.621 MB real(8) xyzh_soa(1200000,4) 36.621 MB real(8) vxyzu(4,1200000) 36.621 MB real(4) alphaind(2,1200000) 9.155 MB real(4) divcurlv(4,1200000) 18.311 MB real(4) divcurlB(4,1200000) 18.311 MB real(8) Bevol(4,1200000) 36.621 MB real(8) Bxyz(3,1200000) 27.466 MB real(8) dustprop(2,0) 0.000 B real(8) dustgasprop(4,0) 0.000 B real(4) VrelVf(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1200000) 201.416 MB real(8) dustevol(11,1200000) 100.708 MB real(8) ddustevol(11,1200000) 100.708 MB real(8) ddustprop(2,0) 0.000 B real(8) deltav(3,11,1200000) 302.124 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,1200000) 36.621 MB real(8) dBevol(4,1200000) 36.621 MB real(4) divBsumm(1200000) 4.578 MB real(8) fext(3,1200000) 27.466 MB real(8) vpred(4,1200000) 36.621 MB real(8) dustpred(11,1200000) 100.708 MB real(8) Bpred(4,1200000) 36.621 MB real(8) dustproppred(2,0) 0.000 B integer(1) iphase(1200000) 1.144 MB integer(1) iphase_soa(1200000) 1.144 MB real(4) gradh(1,1200000) 4.578 MB real(8) tstop(22,1200000) 201.416 MB integer(4) ll(1200000) 4.578 MB integer(4) ibelong(1200000) 4.578 MB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,2400001) 18.311 MB integer(4) inodeparts(1200000) 4.578 MB real(8) xyzh_swap(1200000,4) 36.621 MB integer(4) inodeparts_swap(1200000) 4.578 MB integer(1) iphase_swap(1200000) 1.144 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) list(1200000) 4.578 MB integer(4) cellatid(2400001) 9.155 MB integer(4) ifirstincell(2400001) 9.155 MB kdnode nodeglobal(2400001) 119.019 MB kdnode node(2400001) 119.019 MB integer(4) nodemap(2400001) 9.155 MB --------------------------------------------------------- Total memory allocated to arrays: 1.783 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.200000000000000 dt = 0.200000000000000 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.400000000000000 dt = 0.200000000000000 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.600000000000000 dt = 0.200000000000000 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.800000000000000 dt = 0.200000000000000 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.00000000000000 dt = 0.200000000000000 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.20000000000000 dt = 0.200000000000000 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.40000000000000 dt = 0.200000000000000 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.60000000000000 dt = 0.200000000000000 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.80000000000000 dt = 0.200000000000000 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 = 2.00000000000000 dt = 0.200000000000000 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 = 5.474E-48, 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.250000000000000E-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 = 2.387E-14, 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 = 1.130E-16, tol = 4.441E-16] nsteps per orbit = 1099 norbits = 10 checking angular momentum..............OK [max err = 8.717E-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.000000000000000E+000 is 1.00000000000000 total mass = 6.960406353085308E-004 >>>>>> done ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.548E-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.000000000000000E-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 = 8.882E-16, L2 err = 7.326E-17 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 = 8.882E-16, L2 err = 1.638E-16 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 8.882E-16, L2 err = 1.433E-16 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 = 0.000E+00, L2 err = 0.000E+00 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 = 4.441E-16, L2 err = 9.065E-17 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 = 8.882E-16, L2 err = 1.433E-16 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.441406250000000E-004 t = 0.38059E-03 dt = 2.822E-04 (tolv), np = 4096 t = 0.66275E-03 dt = 3.979E-04 (tolv) WARNING! evolve: N gas particles with energy = 0 : N = 1664 t = 0.10606E-02 dt = 4.818E-04 (tolv) t = 0.15424E-02 dt = 5.881E-04 (tolv) WARNING! evolve: N gas particles with energy = 0 : N = 160 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.50s cpu: 9.0s cpu/wall: 18 Since last dump : 70 timesteps, wall: 0.50s cpu: 9.0s cpu/wall: 18 wall cpu cpu/wall frac step : 0.50s 8.95s 17.90 80.00% step (force): 0.25s 3.95s 15.79 40.00% step (dens) : 0.25s 3.81s 15.23 40.00% write_dump : 0.12s 1.59s 12.73 20.00% Etot= 9.669E-01, Ekin= 1.982E-01, Etherm= 7.686E-01, Epot= 0.000E+00 Emag= 0.000E+00 Linm= 1.976E-18, Angm= 7.930E-19 Centre of Mass = -1.816E-18, -4.825E-18, -2.776E-17 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 Mdust12 = 0.000E+00 Mdust13 = 0.000E+00 Mdust14 = 0.000E+00 Mdust15 = 0.000E+00 Mdust16 = 0.000E+00 Mdust17 = 0.000E+00 Mdust18 = 0.000E+00 Mdust19 = 0.000E+00 Mdust20 = 0.000E+00 Mdust21 = 0.000E+00 Mdust22 = 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 checking total energy..................OK [max err = 2.332E-04, tol = 4.700E-04] checking linear momentum...............OK [max err = 1.976E-18, tol = 7.000E-15] <-- testing complete total wall time = 1 min, 24.00 s (= 8.4000E+01s) total cpu time = 14 min, 56.56 s (= 8.9656E+02s) SUMMARY OF ALL TESTS: PASSED: 136 of 136 100.0% FAILED: 0 of 136 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found