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(4,1000000) 15.259 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) ibin(1000000) 976.562 KB integer(1) ibin_old(1000000) 976.562 KB integer(1) ibin_wake(1000000) 976.562 KB real(4) dt_in(1000000) 3.815 MB real(4) twas(1000000) 3.815 MB integer(1) iphase(1000000) 976.562 KB integer(1) iphase_soa(1000000) 976.562 KB real(4) gradh(1,1000000) 3.815 MB real(8) tstop(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: 463.503 MB --------------------------------------------------------- --> 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.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] <-- KERNEL TEST COMPLETE --> TESTING LINKLIST / NEIGHBOUR FINDING ----------------------------------------------------------------------------------------- random seed = -43587 32768 particles set in uniform random distribution ----------------------------------------------------------------------------------------- Test 1 of 3: building linked list... checking ncells > 0....................OK [got T should be T] checking active/inactive cells......... checking active cells have at least one active particle....OK [checked 4326 of 32768 values] 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 3: building linked list... checking ncells > 0....................OK [got T should be T] Test 3 of 3: building linked list... checking ncells > 0....................OK [got T should be T] checking active/inactive cells......... checking inactive cells have no active particles...........OK [checked 5842 of 32768 values] checking active cells have at least one active particle....OK [checked 3499 of 32768 values] checking neighbour number.............. checking nneigh (cached)...............OK [checked 32768 values] checking nneigh (no cache).............OK [checked 32768 of 32768 values] max nneigh = 165 min nneigh = 0 mean = 26.45 Test 4: 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 force controlling timestep on 58808 gas particles (wall: total = 15.12s force = 5.75s dens = 6.62s link = 2.75s) (cpu : total = 232.92s force = 90.92s dens = 107.37s link = 34.62s) (cpu/wall: total = 15.40 force = 15.81 dens = 16.21 link = 12.59) completed in 15.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 curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking ibin > 0......................OK [got T should be T] --> testing Hydro derivatives (on 1 active particles) trial neigh mean : 173.00, real neigh mean = 38.00 ratio try/act= 4.55 trial neigh max : 173, max real neigh = 38 n neighbour calls : 4, mean per part = 5.000 n density calcs : 6, mean per part = 6.000 force controlling timestep on 1 gas particles (wall: total = 2.75s force = 0.12s dens = 0.00s link = 2.62s) (cpu : total = 35.67s force = 0.54s dens = 0.55s link = 34.57s) (cpu/wall: total = 12.97 force = 4.34 dens = --- link = 13.17) (moved 0.00010% of particles in 18.18182% of time, efficiency = 0.00055%) checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] --> testing Hydro derivatives (on 10 active particles) trial neigh mean : 217.30, real neigh mean = 36.50 ratio try/act= 5.95 trial neigh max : 254, max real neigh = 38 n neighbour calls : 24, mean per part = 3.400 n density calcs : 52, mean per part = 5.200 force controlling timestep on 10 gas particles (wall: total = 2.75s force = 0.00s dens = 0.00s link = 2.75s) (cpu : total = 35.53s force = 0.55s dens = 0.55s link = 34.42s) (cpu/wall: total = 12.92 force = --- dens = --- link = 12.52) (moved 0.00100% of particles in 18.18182% of time, efficiency = 0.00550%) checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] --> testing Hydro derivatives (on 100 active particles) trial neigh mean : 244.69, real neigh mean = 36.90 ratio try/act= 6.63 trial neigh max : 328, max real neigh = 38 n neighbour calls : 256, mean per part = 3.560 n density calcs : 504, mean per part = 5.040 force controlling timestep on 100 gas particles (wall: total = 2.75s force = 0.00s dens = 0.12s link = 2.62s) (cpu : total = 35.68s force = 0.55s dens = 0.55s link = 34.58s) (cpu/wall: total = 12.98 force = --- dens = 4.44 link = 13.17) (moved 0.01000% of particles in 18.18182% of time, efficiency = 0.05500%) checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] --> testing Hydro derivatives (on 1000 active particles) trial neigh mean : 252.06, real neigh mean = 48.63 ratio try/act= 5.18 trial neigh max : 377, max real neigh = 51 n neighbour calls : 1058, mean per part = 2.058 n density calcs : 4206, mean per part = 4.206 force controlling timestep on 1000 gas particles (wall: total = 2.62s force = 0.00s dens = 0.00s link = 2.62s) (cpu : total = 35.74s force = 0.60s dens = 0.70s link = 34.43s) (cpu/wall: total = 13.61 force = --- dens = --- link = 13.12) (moved 0.10000% of particles in 17.35537% of time, efficiency = 0.57619%) checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] --> testing Hydro derivatives (on 10000 active particles) trial neigh mean : 277.39, real neigh mean = 50.21 ratio try/act= 5.53 trial neigh max : 398, max real neigh = 51 n neighbour calls : 10188, mean per part = 2.019 n density calcs : 40784, mean per part = 4.078 force controlling timestep on 10000 gas particles (wall: total = 2.88s force = 0.00s dens = 0.25s link = 2.62s) (cpu : total = 38.14s force = 1.31s dens = 2.40s link = 34.44s) (cpu/wall: total = 13.27 force = --- dens = 9.58 link = 13.12) (moved 1.00000% of particles in 19.00826% of time, efficiency = 5.26087%) checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] --> testing Hydro derivatives (on 100000 active particles) trial neigh mean : 348.16, real neigh mean = 55.42 ratio try/act= 6.28 trial neigh max : 548, max real neigh = 57 n neighbour calls : 25990, mean per part = 1.260 n density calcs : 260340, mean per part = 2.603 force controlling timestep on 13564 gas particles (wall: total = 4.00s force = 0.62s dens = 0.75s link = 2.62s) (cpu : total = 56.35s force = 9.43s dens = 12.47s link = 34.44s) (cpu/wall: total = 14.09 force = 15.09 dens = 16.63 link = 13.12) (moved 10.00000% of particles in 26.44628% of time, efficiency = 37.81250%) checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] --> 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 force controlling timestep on 58808 gas particles (wall: total = 11.38s force = 5.62s dens = 3.12s link = 2.62s) (cpu : total = 174.82s force = 89.83s dens = 50.57s link = 34.42s) (cpu/wall: total = 15.37 force = 15.97 dens = 16.18 link = 13.11) completed in 11.38 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.658E-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 = 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] --> testing artificial viscosity terms (disc viscosity) (on 10000 active particles) trial neigh mean : 276.49, real neigh mean = 50.21 ratio try/act= 5.51 trial neigh max : 398, max real neigh = 51 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10000, mean per part = 1.000 force controlling timestep on 10000 gas particles (wall: total = 2.88s force = 0.12s dens = 0.12s link = 2.62s) (cpu : total = 36.68s force = 1.30s dens = 0.96s link = 34.42s) (cpu/wall: total = 12.76 force = 10.39 dens = 7.66 link = 13.11) 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 curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.658E-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 = 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] --> testing artificial viscosity terms (disc viscosity) (on 100 active particles) trial neigh mean : 232.55, real neigh mean = 36.90 ratio try/act= 6.30 trial neigh max : 321, max real neigh = 38 n neighbour calls : 0, mean per part = 1.000 n density calcs : 100, mean per part = 1.000 force controlling timestep on 100 gas particles (wall: total = 2.75s force = 0.00s dens = 0.12s link = 2.62s) (cpu : total = 35.67s force = 0.55s dens = 0.56s link = 34.57s) (cpu/wall: total = 12.97 force = --- dens = 4.44 link = 13.17) completed in 2.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.658E-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 = 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] --> testing artificial viscosity terms (disc viscosity) (on 1 active particles) trial neigh mean : 173.00, real neigh mean = 38.00 ratio try/act= 4.55 trial neigh max : 173, max real neigh = 38 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1, mean per part = 1.000 force controlling timestep on 1 gas particles (wall: total = 2.75s force = 0.00s dens = 0.12s link = 2.62s) (cpu : total = 35.54s force = 0.55s dens = 0.55s link = 34.44s) (cpu/wall: total = 12.92 force = --- dens = 4.43 link = 13.12) completed in 2.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.658E-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 = 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 thread 0 WARNING: viscosity constraining timestep on 1000000 particles by factor 458160.336753278 (wall: total = 12.00s force = 6.25s dens = 3.12s link = 2.62s) (cpu : total = 184.87s force = 99.69s dens = 50.74s link = 34.43s) (cpu/wall: total = 15.41 force = 15.95 dens = 16.24 link = 13.12) completed in 12.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 curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.629E-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.324E-16, tol = 2.500E-15] checking dvxdz.........................OK [max err = 1.358E-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.325E-17, tol = 2.500E-15] checking dvzdy.........................OK [max err = 6.568E-04, tol = 1.500E-03] checking dvzdz.........................OK [max err = 4.283E-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.019E-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 force controlling timestep on 34608 gas particles (wall: total = 12.88s force = 4.50s dens = 6.38s link = 2.00s) (cpu : total = 198.88s force = 71.16s dens = 102.31s link = 25.40s) (cpu/wall: total = 15.45 force = 15.81 dens = 16.05 link = 12.70) completed in 12.88 s checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] --> testing Hydro derivs in setup with density contrast (nactive= 10) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 358.20, real neigh mean = 57.00 ratio try/act= 6.28 trial neigh max : 389, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10, mean per part = 1.000 (wall: total = 2.00s force = 0.00s dens = 0.12s link = 1.88s) (cpu : total = 27.93s force = 0.57s dens = 0.56s link = 26.80s) (cpu/wall: total = 13.97 force = --- dens = 4.48 link = 14.29) (moved 0.00154% of particles in 15.53398% of time, efficiency = 0.00993%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.202E-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] --> testing Hydro derivs in setup with density contrast (nactive= 100) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 351.83, real neigh mean = 57.00 ratio try/act= 6.17 trial neigh max : 489, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 100, mean per part = 1.000 (wall: total = 2.00s force = 0.00s dens = 0.00s link = 2.00s) (cpu : total = 27.52s force = 0.57s dens = 0.56s link = 26.39s) (cpu/wall: total = 13.76 force = --- dens = --- link = 13.19) (moved 0.01542% of particles in 15.53398% of time, efficiency = 0.09928%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.205E-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] --> testing Hydro derivs in setup with density contrast (nactive= 1000) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 354.22, real neigh mean = 57.00 ratio try/act= 6.21 trial neigh max : 500, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000, mean per part = 1.000 (wall: total = 2.00s force = 0.00s dens = 0.00s link = 2.00s) (cpu : total = 27.66s force = 0.64s dens = 0.61s link = 26.41s) (cpu/wall: total = 13.83 force = --- dens = --- link = 13.20) (moved 0.15422% of particles in 15.53398% of time, efficiency = 0.99278%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.206E-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] --> testing Hydro derivs in setup with density contrast (nactive= 10000) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 353.60, real neigh mean = 57.00 ratio try/act= 6.20 trial neigh max : 500, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10000, mean per part = 1.000 (wall: total = 2.12s force = 0.12s dens = 0.00s link = 2.00s) (cpu : total = 29.01s force = 1.52s dens = 1.09s link = 26.40s) (cpu/wall: total = 13.65 force = 12.18 dens = --- link = 13.20) (moved 1.54218% of particles in 16.50485% of time, efficiency = 9.34381%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.206E-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] --> testing Hydro derivs in setup with density contrast (nactive= 73824) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 344.82, real neigh mean = 57.00 ratio try/act= 6.05 trial neigh max : 500, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 73824, mean per part = 1.000 (wall: total = 2.75s force = 0.50s dens = 0.25s link = 2.00s) (cpu : total = 38.24s force = 7.50s dens = 4.37s link = 26.38s) (cpu/wall: total = 13.91 force = 14.99 dens = 17.47 link = 13.19) (moved 11.38500% of particles in 21.35922% of time, efficiency = 53.30251%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.260E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.206E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 7.283E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.205E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.289E-05, tol = 1.000E-03] --> testing force evaluation with ind_timesteps ----------------------------------------------------------------------------------------- random seed = -43587 125000 particles set in uniform random distribution ----------------------------------------------------------------------------------------- evaluating derivs with all particles active... trial neigh mean : 557.05, real neigh mean = 49.04 ratio try/act= 11.36 trial neigh max : 1131, max real neigh = 123 n neighbour calls : 59433, mean per part = 1.475 n density calcs : 712691, mean per part = 5.702 force controlling timestep on 51687 gas particles (wall: total = 3.50s force = 0.88s dens = 2.25s link = 0.38s) (cpu : total = 55.33s force = 13.15s dens = 37.45s link = 4.72s) (cpu/wall: total = 15.81 force = 15.03 dens = 16.64 link = 12.59) evaluating derivs with 12500 particles active... trial neigh mean : 551.76, real neigh mean = 49.06 ratio try/act= 11.25 trial neigh max : 1114, max real neigh = 103 n neighbour calls : 0, mean per part = 1.000 n density calcs : 12500, mean per part = 1.000 force controlling timestep on 5236 gas particles (wall: total = 0.62s force = 0.12s dens = 0.12s link = 0.38s) (cpu : total = 9.11s force = 2.42s dens = 1.55s link = 5.14s) (cpu/wall: total = 14.58 force = 19.33 dens = 12.44 link = 13.71) checking force(x)......................OK [max err = 2.323E-07, L2 err = 4.641E-11 tol = 1.000E-04] checking force(y)......................OK [max err = 7.352E-07, L2 err = 8.588E-11 tol = 1.000E-04] checking force(z)......................OK [max err = 5.010E-07, L2 err = 4.360E-11 tol = 1.000E-04] checking du/dt.........................OK [max err = 3.603E-08, L2 err = 4.223E-10 tol = 1.000E-05] <-- 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 <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000049478A dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 000000000072AE61 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000964EB1 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000497EE7 dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 000000000072AE61 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000964EB1 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000499161 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 000000000072AE61 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000964EB1 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000049A3DF dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 000000000072AE61 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000964EB1 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000049CEC0 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 000000000072AE61 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000964EB1 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown >>> reading setup from file: test.dump on unit 9 FT:Phantom:1.2.0: (hydro+dust): 21/02/2019 00:58:46.7 forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000047C39F dump_utils_mp_rea 1225 utils_dumpfiles.f90 phantom 0000000000737BD6 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000966C76 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000481BA3 dump_utils_mp_rea 1246 utils_dumpfiles.f90 phantom 0000000000737BD6 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000966C76 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000048390E dump_utils_mp_rea 1253 utils_dumpfiles.f90 phantom 0000000000737BD6 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000966C76 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000048567D dump_utils_mp_rea 1260 utils_dumpfiles.f90 phantom 0000000000737BD6 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000966C76 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000489C97 dump_utils_mp_rea 1281 utils_dumpfiles.f90 phantom 0000000000737BD6 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000966C76 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.30000000000000 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. 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(1,0) 0.000 B real(4) divcurlv(4,1001) 15.641 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) ibin(1001) 1001.000 B integer(1) ibin_old(1001) 1001.000 B integer(1) ibin_wake(1001) 1001.000 B real(4) dt_in(1001) 3.910 KB real(4) twas(1001) 3.910 KB 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: 609.699 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 forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 00000000004B1600 dump_utils_mp_rea 2103 utils_dumpfiles.f90 phantom 000000000096BDE9 testrwdump_mp_tes 309 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown ->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 <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000049478A dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 0000000000734BAD readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000964E50 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000497EE7 dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 0000000000734BAD readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000964E50 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000499161 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 0000000000734BAD readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000964E50 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000049A3DF dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 0000000000734BAD readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000964E50 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000049CEC0 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 0000000000734BAD readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000964E50 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown --> 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 00:58:47.8 WARNING! read_smalldump: *** VELOCITY WILL BE MISSING FROM SMALL DUMP FILES *** forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000047C39F dump_utils_mp_rea 1225 utils_dumpfiles.f90 phantom 000000000073B766 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000966B8C testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000481BA3 dump_utils_mp_rea 1246 utils_dumpfiles.f90 phantom 000000000073B766 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000966B8C testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000048390E dump_utils_mp_rea 1253 utils_dumpfiles.f90 phantom 000000000073B766 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000966B8C testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 000000000048567D dump_utils_mp_rea 1260 utils_dumpfiles.f90 phantom 000000000073B766 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000966B8C testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000489C97 dump_utils_mp_rea 1281 utils_dumpfiles.f90 phantom 000000000073B766 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000966B8C testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009B8312 test_mp_testsuite 251 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.29999995231628 comment: hfact was 1.12300002574921 now 1.12300000000000 comment: tolh was 9.999999747378752E-005 now 1.000000000000000E-004 comment: C_cour was 0.300000011920929 now 0.300000000000000 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. 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(1,0) 0.000 B real(4) divcurlv(4,1001) 15.641 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) ibin(1001) 1001.000 B integer(1) ibin_old(1001) 1001.000 B integer(1) ibin_wake(1001) 1001.000 B real(4) dt_in(1001) 3.910 KB real(4) twas(1001) 3.910 KB 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: 609.699 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(4,1000000) 15.259 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) ibin(1000000) 976.562 KB integer(1) ibin_old(1000000) 976.562 KB integer(1) ibin_wake(1000000) 976.562 KB real(4) dt_in(1000000) 3.815 MB real(4) twas(1000000) 3.815 MB integer(1) iphase(1000000) 976.562 KB integer(1) iphase_soa(1000000) 976.562 KB real(4) gradh(1,1000000) 3.815 MB real(8) tstop(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: 463.503 MB --------------------------------------------------------- --> SKIPPING TEST OF STEP MODULE (need -DPERIODIC) --> TESTING IND TIMESTEP UTILS checking bin for dt = 1.500E+00.......OK [got 0 should be 0] checking bin for dt = 7.500E-01.......OK [got 1 should be 1] checking bin for dt = 3.750E-01.......OK [got 2 should be 2] checking bin for dt = 1.875E-01.......OK [got 3 should be 3] checking bin for dt = 9.375E-02.......OK [got 4 should be 4] checking bin for dt = 4.688E-02.......OK [got 5 should be 5] checking bin for dt = 2.344E-02.......OK [got 6 should be 6] checking bin for dt = 1.172E-02.......OK [got 7 should be 7] checking bin for dt = 5.859E-03.......OK [got 8 should be 8] checking bin for dt = 2.930E-03.......OK [got 9 should be 9] checking bin for dt = 1.465E-03.......OK [got 10 should be 10] checking bin for dt = 7.324E-04.......OK [got 11 should be 11] checking bin for dt = 3.662E-04.......OK [got 12 should be 12] checking bin for dt = 1.831E-04.......OK [got 13 should be 13] checking bin for dt = 9.155E-05.......OK [got 14 should be 14] checking bin for dt = 4.578E-05.......OK [got 15 should be 15] checking bin for dt = 2.289E-05.......OK [got 16 should be 16] checking bin for dt = 1.144E-05.......OK [got 17 should be 17] checking bin for dt = 5.722E-06.......OK [got 18 should be 18] checking bin for dt = 2.861E-06.......OK [got 19 should be 19] checking bin for dt = 1.431E-06.......OK [got 20 should be 20] checking bin for dt = 7.153E-07.......OK [got 21 should be 21] checking bin for dt = 3.576E-07.......OK [got 22 should be 22] checking bin for dt = 1.788E-07.......OK [got 23 should be 23] checking bin for dt = 8.941E-08.......OK [got 24 should be 24] checking bin for dt = 4.470E-08.......OK [got 25 should be 25] checking bin for dt = 2.235E-08.......OK [got 26 should be 26] checking bin for dt = 1.118E-08.......OK [got 27 should be 27] checking bin for dt = 5.588E-09.......OK [got 28 should be 28] checking bin for dt = 2.794E-09.......OK [got 29 should be 29] checking bin for dt = huge.............OK [got 0 should be 0] checking bin for dt = 0................OK [got 30 should be 30] ---> checking change_nbinmax routine checking nbinmax = 1 checking nbinmax = 2 checking nbinmax = 3 checking nbinmax = 4 checking nbinmax = 5 checking nbinmax = 6 checking nbinmax = 7 checking nbinmax = 8 checking nbinmax = 9 checking nbinmax = 10 checking nbinmax = 11 checking nbinmax = 12 checking nbinmax = 13 checking nbinmax = 14 checking nbinmax = 15 checking nbinmax = 16 checking nbinmax = 17 checking nbinmax = 18 checking nbinmax = 19 checking nbinmax = 20 checking nbinmax = 21 checking nbinmax = 22 checking nbinmax = 23 checking nbinmax = 24 checking increasing nbins keeps step fraction....OK [max err = 0.000E+00, tol = 2.225-308] checking decreasing nbins keeps step fraction....OK [max err = 0.000E+00, tol = 2.225-308] <-- IND TIMESTEP UTILS TEST COMPLETE --> 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.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-11, tol = 1.500E-03] checking corotate external force initialisation..OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 2.237E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 2.812E-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 = 1.012E-10, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 5.524E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 6.598E-11, tol = 1.500E-03] checking prdrag external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-11, tol = 1.500E-03] checking torus external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 3.904E-07, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 4.945E-07, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.459E-10, 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.271E-05, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 7.191E-05, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 7.013E-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.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-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 = 9.697E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 8.048E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 1.058E-10, tol = 1.500E-03] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-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.214E-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.109E-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.458E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 3.145E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 2.126E-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 = 1.766E-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 = 1.766E-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.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 = 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) = (-5.523E-17, -1.017E-17, -5.192E-19) Particle setup OK checking no errors during disc setup...OK [got 0 should be 0] dt for sinks = 1.250000000000000E-002 nsteps per orbit = 711 norbits = 10 checking angular momentum..............OK [max err = 8.867E-08, tol = 3.000E-07] checking linear momentum...............OK [max err = 1.020E-14, tol = 6.000E-14] checking total energy..................OK [max err = 1.063E-03, tol = 2.000E-03] --> testing integration of disc around eccentric binary Centre of mass is at (x,y,z) = ( 0.000E+00, 1.505E-17, 1.551E-19) Particle setup OK checking no errors during disc setup...OK [got 0 should be 0] dt for sinks = 2.296396633859230E-002 nsteps per orbit = 387 norbits = 10 checking angular momentum..............OK [max err = 1.436E-07, tol = 2.100E-06] checking linear momentum...............OK [max err = 7.067E-15, tol = 5.000E-06] checking total energy..................OK [max err = 2.473E-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) forrtl: warning (406): fort: (1): In call to GETNEIGH_POS, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000780400 ptmass_mp_ptmass_ 964 ptmass.F90 phantom 000000000099D732 testptmass_mp_tes 594 test_ptmass.F90 phantom 00000000009B83DE test_mp_testsuite 283 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown 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 forrtl: warning (406): fort: (1): In call to GETNEIGH_POS, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B008F79E5A3 for_emit_diagnost Unknown Unknown phantom 0000000000780400 ptmass_mp_ptmass_ 964 ptmass.F90 phantom 000000000099D732 testptmass_mp_tes 594 test_ptmass.F90 phantom 00000000009B83DE test_mp_testsuite 283 testsuite.F90 phantom 00000000009D5F8F MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B009271CC05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown 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.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 --> 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.8033 (/H) mean...........: 0.8641 (/H) minimum........: 0.7810 (/H) maximum........: 1.2024 e.g. for alpha_SS = 0.1, use alpha_AV = 1.15733 using alpha_AV = 1.0 gives alpha_SS ~ 0.0781-> 0.1202 # 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.864060119 ! approx. mean smoothing length over disc scale height alphaSS_min = 0.0781017121566 ! minimum Shakura-Sunyaev alpha viscosity in disc alphaSS_max = 0.120235057 ! 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.58E-08 g*cm^2/sec Centre of mass is at (x,y,z) = (-2.281E-14, 6.682E-15, 7.408E-16) Particle setup OK checking setup errors..................OK [got 0 should be 0] checking setup warnings................OK [got 0 should be 0] completed in 26.38 s checking vphi**2/r = force_r...........OK [max err = 2.867E-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 = 4.441E-16, L2 err = 1.110E-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 = 8 min, 32.38 s (= 5.1238E+02s) total cpu time = 1 hr, 18 min, 58.63 s (= 4.7386E+03s) SUMMARY OF ALL TESTS: PASSED: 137 of 137 100.0% FAILED: 0 of 137 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found