MPI parallelisation is OFF _ \ | | ___| _ \ | | | | __ \ _` | __ \ __| _ \ __ `__ \ \___ \ | | | | ___/ | | | ( | | | | ( | | | | | ___/ ___ | _| _| |_|\__,_|_| _|\__|\___/ _| _| _|_____/ _| _| _| ___ , _ |) (| \ _, o _ |\ /|/ \ ,_ o _ _ |/\_| | | |/ | /|/| | |/ |/ |__// | | / |/ \/ \/|/ (\__/ \/|_/ | |_/|/|_/|_/ | |/|/\__/|_/ (| Version: 1.2.0 phantom_version not found; no git information to display Using NATIVE inverse sqrt Running on Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50GHz 24 cpus @ 2.50 Ghz, cache size 16384 KB Running in openMP on 16 threads --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1000000) 30.518 MB real(8) xyzh_soa(1000000,4) 30.518 MB real(8) vxyzu(4,1000000) 30.518 MB real(4) alphaind(1,0) 0.000 B real(4) divcurlv(1,1000000) 3.815 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,0) 0.000 B real(8) Bxyz(3,0) 0.000 B real(8) dustprop(4,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(2,0) 0.000 B real(8) dustevol(1,0) 0.000 B real(8) ddustevol(1,0) 0.000 B real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,1,0) 0.000 B real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(0) 0.000 B real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1000000) 30.518 MB real(8) dBevol(4,0) 0.000 B real(4) divBsumm(0) 0.000 B real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,0) 0.000 B 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(2,1000000) 15.259 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: 437.753 MB --------------------------------------------------------- --> 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 = 184 min nneigh = 1 mean = 32.33 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 : 353.08, real neigh mean = 56.33 ratio try/act= 6.27 trial neigh max : 548, max real neigh = 57 n neighbour calls : 191944, mean per part = 1.192 n density calcs : 2230472, mean per part = 2.230 dtmin = 1.0483190326299730E-003 dtmax = 1.9552730199381144E-002 dtmax/dtmin = 18.651507404503516 dtcour/dtf = 1.3131832163368353 (wall: total = 2.00s force = 0.62s dens = 1.25s link = 0.12s) (cpu : total = 30.25s force = 9.35s dens = 19.97s link = 0.93s) (cpu/wall: total = 15.12 force = 14.96 dens = 15.97 link = 7.45) completed in 2.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 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] --> testing artificial viscosity terms (disc viscosity) trial neigh mean : 349.69, real neigh mean = 56.33 ratio try/act= 6.21 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 = 7.2744556021618536E-004 dtmax/dtmin = 1.4203811872679410 dtcour/dtf = 1.4848776275147300 (wall: total = 1.00s force = 0.62s dens = 0.25s link = 0.12s) (cpu : total = 14.01s force = 9.20s dens = 3.91s link = 0.90s) (cpu/wall: total = 14.01 force = 14.71 dens = 15.66 link = 7.22) completed in 1.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 art. visc force(x)............OK [max err = 3.079E-03, tol = 3.200E-02] checking art. visc force(y)............OK [max err = 3.104E-03, tol = 2.400E-02] checking art. visc force(z)............OK [max err = 6.242E-04, tol = 2.400E-02] --> SKIPPING Cullen-Dehnen terms (need nalpha=2) --> testing physical viscosity terms (two first derivatives) trial neigh mean : 349.69, real neigh mean = 56.33 ratio try/act= 6.21 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 = 2129.2332745283952 WARNING! force: viscosity constraining Courant timestep (wall: total = 1.12s force = 0.75s dens = 0.25s link = 0.12s) (cpu : total = 16.09s force = 10.04s dens = 5.15s link = 0.90s) (cpu/wall: total = 14.31 force = 13.39 dens = 20.61 link = 7.21) completed in 1.12 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 dvxdx.........................OK [max err = 8.214E-04, tol = 1.700E-03] checking dvxdy.........................OK [max err = 1.460E-16, tol = 2.500E-15] checking dvxdz.........................OK [max err = 1.460E-16, tol = 2.500E-15] checking dvydx.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvydy.........................OK [max err = 2.376E-16, tol = 2.500E-15] checking dvydz.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvzdx.........................OK [max err = 4.323E-17, tol = 2.500E-15] checking dvzdy.........................OK [max err = 6.568E-04, tol = 1.500E-03] checking dvzdz.........................OK [max err = 4.346E-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 \sum v.dv/dt + du/dt = 0......OK [max err = 5.014E-07, tol = 1.520E-06] checking du/dt >= 0 for all particles..OK [got T should be T] --> SKIPPING dust evolution terms (need -DDUST) --> 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 : 426.73, real neigh mean = 57.20 ratio try/act= 7.46 trial neigh max : 88189, max real neigh = 988 n neighbour calls : 133141, mean per part = 1.205 n density calcs : 1504111, mean per part = 2.320 dtmin = 2.8817634064656154E-004 dtmax = 2.1540775519298411E-002 dtmax/dtmin = 74.748591334610850 dtcour/dtf = 1.3274325079742730 (wall: total = 1.50s force = 0.50s dens = 0.88s link = 0.12s) (cpu : total = 23.40s force = 8.26s dens = 14.44s link = 0.69s) (cpu/wall: total = 15.60 force = 16.52 dens = 16.51 link = 5.55) completed in 1.50 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 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] <-- 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) --> SKIPPING DUST TEST (REQUIRES -DDUST) --> 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: (hydro+dust): 21/02/2019 05:47:58.0 n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.3000000000000000 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. 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(1,0) 0.000 B real(4) divcurlv(1,1001) 3.910 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,0) 0.000 B real(8) Bxyz(3,0) 0.000 B 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(2,0) 0.000 B real(8) dustevol(1,0) 0.000 B real(8) ddustevol(1,0) 0.000 B real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,1,0) 0.000 B 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,0) 0.000 B real(4) divBsumm(0) 0.000 B real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,0) 0.000 B 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(2,1001) 15.641 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: 591.126 KB --------------------------------------------------------- number of array sizes = 2 number of blocks = 1 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 (hydro) 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 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 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 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: (hydro+dust): 21/02/2019 05:47:58.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.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 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. 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(1,0) 0.000 B real(4) divcurlv(1,1001) 3.910 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,0) 0.000 B real(8) Bxyz(3,0) 0.000 B 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(2,0) 0.000 B real(8) dustevol(1,0) 0.000 B real(8) ddustevol(1,0) 0.000 B real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,1,0) 0.000 B 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,0) 0.000 B real(4) divBsumm(0) 0.000 B real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,0) 0.000 B 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(2,1001) 15.641 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: 591.126 KB --------------------------------------------------------- number of array sizes = 2 number of blocks = 1 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! 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 | ----------------------------------------------- <<< finished reading (hydro) 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 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(1,0) 0.000 B real(4) divcurlv(1,1000000) 3.815 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,0) 0.000 B real(8) Bxyz(3,0) 0.000 B real(8) dustprop(4,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(2,0) 0.000 B real(8) dustevol(1,0) 0.000 B real(8) ddustevol(1,0) 0.000 B real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,1,0) 0.000 B real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(0) 0.000 B real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1000000) 30.518 MB real(8) dBevol(4,0) 0.000 B real(4) divBsumm(0) 0.000 B real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,0) 0.000 B 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(2,1000000) 15.259 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: 426.309 MB --------------------------------------------------------- --> SKIPPING TEST OF STEP MODULE (need -DPERIODIC) --> 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 = 8.195E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.159E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 5.017E-11, tol = 1.500E-03] checking corotate external force initialisation..OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 3.300E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 3.074E-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 = 3.866E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 5.082E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 7.494E-11, tol = 1.500E-03] checking prdrag external force initialisation....OK [got 0 should be 0] updating prdrag at t = 0.0000000000000000 Mdot = 0.0000000000000000 LAcc = 0.0000000000000000 checking fextx = -grad phi.............OK [max err = 8.195E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.159E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 5.017E-11, tol = 1.500E-03] checking torus external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 1.598E-09, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 5.828E-09, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 5.017E-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.552E-05, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.543E-05, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 7.052E-05, tol = 1.500E-03] checking Lense-Thirring external force initialisation......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 8.195E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.159E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 5.017E-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.154E-10, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 1.263E-10, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 1.000E-10, tol = 1.500E-03] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 8.195E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.159E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 5.017E-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.220E-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 = 9.219E-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 = 3.848E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 3.938E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 1.571E-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] updating prdrag at t = 0.0000000000000000 Mdot = 0.0000000000000000 LAcc = 0.0000000000000000 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] --> testing integration of circumbinary disc Centre of mass is at (x,y,z) = ( 0.000E+00, 1.224E-17, -1.027E-18) Particle setup OK checking no errors during disc setup...OK [got 0 should be 0] dt for sinks = 1.2500000000000001E-002 nsteps per orbit = 711 norbits = 10 checking angular momentum..............OK [max err = 1.053E-07, tol = 3.000E-07] checking linear momentum...............OK [max err = 2.447E-14, tol = 6.000E-14] checking total energy..................OK [max err = 9.685E-04, tol = 2.000E-03] --> testing integration of disc around eccentric binary Centre of mass is at (x,y,z) = ( 0.000E+00, 1.090E-17, 2.840E-18) Particle setup OK checking no errors during disc setup...OK [got 0 should be 0] dt for sinks = 2.2963966338592295E-002 nsteps per orbit = 387 norbits = 10 checking angular momentum..............OK [max err = 7.100E-07, tol = 1.100E-06] checking linear momentum...............OK [max err = 2.237E-14, tol = 4.000E-14] checking total energy..................OK [max err = 2.603E-03, tol = 1.200E-02] --> 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 --> TESTING DISC SETUP --> testing set_disc Phantom: general disc setup (see .discparams file for details) Setting up disc containing 1000000 gas particles 100000 200000 300000 400000 500000 600000 700000 800000 900000 1000000 (/H) per particle...: 0.8035 (/H) mean...........: 0.8642 (/H) minimum........: 0.7837 (/H) maximum........: 1.1942 e.g. for alpha_SS = 0.1, use alpha_AV = 1.15716 using alpha_AV = 1.0 gives alpha_SS ~ 0.0784-> 0.1194 # gas disc parameters - this file is NOT read by setup R_in = 0.500 ! inner disc boundary R_ref = 0.500 ! reference radius R_out = 10. ! outer disc boundary H/R_in = 0.020 ! disc aspect ratio H/R at R=R_in H/R_ref = 0.020 ! disc aspect ratio H/R at R=R_ref H/R_out = 0.00945741609 ! disc aspect ratio H/R at R=R_out sig_in = 0.000 ! surface density (g/cm^2) at R=R_in sig_ref = 0.000 ! surface density (g/cm^2) at R=R_ref sig_out = 1.399E-07 ! surface density (g/cm^2) at R=R_out sig_max = 1.700E-06 ! maximum surface density (g/cm^2) Qmin = 215.122813 ! minimum Toomre Q parameter n = 1000000 ! number of particles in the disc p_index = 1.500 ! power law index of surface density profile q_index = 0.750 ! power law index of sound speed profile M_star = 1.000 ! mass of central star M_disc = 1.000E-04 ! disc mass M_disc/M_star = 1.000E-04 ! relative disc mass cs0 = 0.0168179283051 ! sound speed at R=1 T_in = 1.540E-18 ! temperature (K) at R=R_in T_ref = 1.540E-18 ! temperature (K) at R=R_ref T_out = 1.721E-20 ! temperature (K) at R=R_out inc.deg = 0.000 ! disc inclination in degrees = 0.86418656 ! approx. mean smoothing length over disc scale height alphaSS_min = 0.0783737666392 ! minimum Shakura-Sunyaev alpha viscosity in disc alphaSS_max = 0.119415242 ! maximum Shakura-Sunyaev alpha viscosity in disc udist = 1.000 ! distance units (cgs) umass = 1.000 ! mass units (cgs) utime = 3871.42319 ! time units (cgs) G = 1.000 ! in code units c = 1.161E+14 ! in code units # Temperature profile = 0.0K (R/ 0.5)^(-1.5) # Surface density = 1.61E-05 g/cm^2 (R/ 0.5)^(-1.5) (1 - sqrt( 0.5/R)) # Disc total angular momentum = 5.57E-08 g*cm^2/sec Centre of mass is at (x,y,z) = (-1.647E-14, 7.583E-14, 3.018E-15) Particle setup OK checking setup errors..................OK [got 0 should be 0] checking setup warnings................OK [got 0 should be 0] completed in 2.88 s checking vphi**2/r = force_r...........OK [max err = 2.407E-01, tol = 3.100E-01] <-- DISC SETUP TESTS COMPLETE --> 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 --> SKIPPING Sedov blast wave (needs -DPERIODIC) <-- testing complete total wall time = 48.62 s total cpu time = 10 min, 28.69 s (= 6.2869E+02s) --> TESTING M_4 cubic kernel <-- KERNEL TEST COMPLETE SUMMARY OF ALL TESTS: PASSED: 112 of 112 100.0% FAILED: 0 of 112 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found