MPI parallelisation is OFF _ \ | | ___| _ \ | | | | __ \ _` | __ \ __| _ \ __ `__ \ \___ \ | | | | ___/ | | | ( | | | | ( | | | | | ___/ ___ | _| _| |_|\__,_|_| _|\__|\___/ _| _| _|_____/ _| _| _| ___ , _ |) (| \ _, o _ |\ /|/ \ ,_ o _ _ |/\_| | | |/ | /|/| | |/ |/ |__// | | / |/ \/ \/|/ (\__/ \/|_/ | |_/|/|_/|_/ | |/|/\__/|_/ (| Version: 1.2.0 phantom_version not found; no git information to display Using NATIVE inverse sqrt Running on Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50GHz 24 cpus @ 2.50 Ghz, cache size 16384 KB Running in openMP on 16 threads --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1000000) 30.518 MB real(8) xyzh_soa(1000000,4) 30.518 MB real(8) vxyzu(4,1000000) 30.518 MB real(4) alphaind(2,1000000) 7.629 MB real(4) divcurlv(4,1000000) 15.259 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,1000000) 30.518 MB real(8) Bxyz(3,1000000) 22.888 MB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(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,1000000) 30.518 MB real(4) divBsumm(1000000) 3.815 MB real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,1000000) 30.518 MB 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: 585.573 MB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) checking wab0 = wab(0).................OK [max err = 0.000E+00, tol = 2.225-308] checking gradh0 = -3.*wab(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking dphidh0 = dphidh(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking radkern2 = radkern*radkern....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == wkern...........OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == grkern..........OK [max err = 0.000E+00, tol = 2.225-308] checking w gradient equal to gradw.....OK [max err = 1.499E-07, tol = 2.000E-05] checking potential gradient = force....OK [max err = 6.722E-08, tol = 2.000E-05] checking get_kernel_grav1 == wkern.....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel_grav1 == grkern....OK [max err = 0.000E+00, tol = 2.225-308] checking dphi/dh = phi - q*dphi/dq.....OK [max err = 8.105E-15, tol = 2.000E-07] --> TESTING LINKLIST / NEIGHBOUR FINDING ----------------------------------------------------------------------------------------- random seed = -43587 32768 particles set in uniform random distribution ----------------------------------------------------------------------------------------- Test 1 of 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 = 8945 min nneigh = 1 mean = 2206.93 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 = 8963 min nneigh = 1 mean = 1806.18 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 : 367.52, real neigh mean = 57.00 ratio try/act= 6.45 trial neigh max : 548, max real neigh = 57 n neighbour calls : 155520, mean per part = 1.156 n density calcs : 2000000, mean per part = 2.000 (wall: total = 5.62s force = 2.12s dens = 3.00s link = 0.50s) (cpu : total = 83.40s force = 33.58s dens = 46.90s link = 2.91s) (cpu/wall: total = 14.83 force = 15.80 dens = 15.63 link = 5.83) completed in 5.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking mean nneigh...................OK [max err = 0.000E+00, tol = 2.225-308] checking max nneigh....................OK [got 57 should be 57] checking n density calcs...............OK [is 2000000 should be 2000000] checking total nneigh..................OK [is 57000000 should be 57000000] checking ibin > 0......................OK [got T should be T] --> testing Hydro derivatives (on 1 active particles) trial neigh mean : 407.00, real neigh mean = 57.00 ratio try/act= 7.14 trial neigh max : 407, max real neigh = 57 n neighbour calls : 1, mean per part = 2.000 n density calcs : 2, mean per part = 2.000 (wall: total = 0.50s force = 0.00s dens = 0.12s link = 0.38s) (cpu : total = 4.41s force = 0.35s dens = 1.30s link = 2.75s) (cpu/wall: total = 8.82 force = --- dens = 10.41 link = 7.34) (moved 0.00010% of particles in 8.88889% of time, efficiency = 0.00112%) 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 : 394.00, real neigh mean = 57.00 ratio try/act= 6.91 trial neigh max : 502, max real neigh = 57 n neighbour calls : 6, mean per part = 1.600 n density calcs : 20, mean per part = 2.000 (wall: total = 0.62s force = 0.00s dens = 0.12s link = 0.50s) (cpu : total = 4.41s force = 0.36s dens = 1.31s link = 2.75s) (cpu/wall: total = 7.06 force = --- dens = 10.46 link = 5.49) (moved 0.00100% of particles in 11.11111% of time, efficiency = 0.00900%) 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 : 400.44, real neigh mean = 57.00 ratio try/act= 7.03 trial neigh max : 533, max real neigh = 57 n neighbour calls : 64, mean per part = 1.640 n density calcs : 200, mean per part = 2.000 (wall: total = 0.50s force = 0.00s dens = 0.00s link = 0.50s) (cpu : total = 4.43s force = 0.36s dens = 1.31s link = 2.76s) (cpu/wall: total = 8.86 force = --- dens = --- link = 5.51) (moved 0.01000% of particles in 8.88889% of time, efficiency = 0.11250%) 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 : 360.76, real neigh mean = 57.00 ratio try/act= 6.33 trial neigh max : 533, max real neigh = 57 n neighbour calls : 328, mean per part = 1.328 n density calcs : 2000, mean per part = 2.000 (wall: total = 0.62s force = 0.00s dens = 0.12s link = 0.50s) (cpu : total = 4.45s force = 0.38s dens = 1.31s link = 2.75s) (cpu/wall: total = 7.11 force = --- dens = 10.50 link = 5.50) (moved 0.10000% of particles in 11.11111% of time, efficiency = 0.90000%) 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 : 374.99, real neigh mean = 57.00 ratio try/act= 6.58 trial neigh max : 538, max real neigh = 57 n neighbour calls : 3312, mean per part = 1.331 n density calcs : 20000, mean per part = 2.000 (wall: total = 0.62s force = 0.00s dens = 0.12s link = 0.50s) (cpu : total = 4.95s force = 0.66s dens = 1.54s link = 2.75s) (cpu/wall: total = 7.92 force = --- dens = 12.29 link = 5.50) (moved 1.00000% of particles in 11.11111% of time, efficiency = 9.00000%) 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 : 381.20, real neigh mean = 57.00 ratio try/act= 6.69 trial neigh max : 548, max real neigh = 57 n neighbour calls : 17136, mean per part = 1.171 n density calcs : 200000, mean per part = 2.000 (wall: total = 1.00s force = 0.12s dens = 0.38s link = 0.50s) (cpu : total = 12.06s force = 3.64s dens = 5.67s link = 2.75s) (cpu/wall: total = 12.06 force = 29.11 dens = 15.12 link = 5.50) (moved 10.00000% of particles in 17.77778% of time, efficiency = 56.25000%) 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 (individual alpha) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 (wall: total = 4.00s force = 2.00s dens = 1.50s link = 0.50s) (cpu : total = 60.29s force = 33.06s dens = 24.47s link = 2.76s) (cpu/wall: total = 15.07 force = 16.53 dens = 16.31 link = 5.51) completed in 4.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.638E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking art. visc force(x)............OK [max err = 2.705E-02, tol = 3.200E-02] checking art. visc force(y)............OK [max err = 1.147E-02, tol = 2.400E-02] checking art. visc force(z)............OK [max err = 9.451E-03, tol = 2.400E-02] --> testing artificial viscosity terms (individual alpha) (on 10000 active particles) trial neigh mean : 372.13, real neigh mean = 57.00 ratio try/act= 6.53 trial neigh max : 534, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10000, mean per part = 1.000 (wall: total = 0.62s force = 0.12s dens = 0.00s link = 0.50s) (cpu : total = 4.74s force = 1.39s dens = 0.60s link = 2.75s) (cpu/wall: total = 7.59 force = 11.15 dens = --- link = 5.50) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.638E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking art. visc force(x)............OK [max err = 2.705E-02, tol = 3.200E-02] checking art. visc force(y)............OK [max err = 1.147E-02, tol = 2.400E-02] checking art. visc force(z)............OK [max err = 9.451E-03, tol = 2.400E-02] --> testing artificial viscosity terms (individual alpha) (on 100 active particles) trial neigh mean : 399.04, real neigh mean = 57.00 ratio try/act= 7.00 trial neigh max : 532, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 100, mean per part = 1.000 (wall: total = 0.62s force = 0.00s dens = 0.12s link = 0.50s) (cpu : total = 4.37s force = 0.36s dens = 1.26s link = 2.75s) (cpu/wall: total = 7.00 force = --- dens = 10.09 link = 5.50) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.638E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking art. visc force(x)............OK [max err = 2.705E-02, tol = 3.200E-02] checking art. visc force(y)............OK [max err = 1.147E-02, tol = 2.400E-02] checking art. visc force(z)............OK [max err = 9.451E-03, tol = 2.400E-02] --> testing artificial viscosity terms (individual alpha) (on 1 active particles) trial neigh mean : 407.00, real neigh mean = 57.00 ratio try/act= 7.14 trial neigh max : 407, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1, mean per part = 1.000 (wall: total = 0.62s force = 0.00s dens = 0.12s link = 0.50s) (cpu : total = 4.43s force = 0.36s dens = 1.31s link = 2.76s) (cpu/wall: total = 7.09 force = --- dens = 10.51 link = 5.52) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.638E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking art. visc force(x)............OK [max err = 2.705E-02, tol = 3.200E-02] checking art. visc force(y)............OK [max err = 1.147E-02, tol = 2.400E-02] checking art. visc force(z)............OK [max err = 9.451E-03, tol = 2.400E-02] --> testing ddivv/dt in Cullen & Dehnen switch trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 completed in 2.00 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.689E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking alphaloc......................OK [max err = 3.421E-04, tol = 3.500E-04] --> testing physical viscosity terms (two first derivatives) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 thread 0 WARNING: viscosity constraining timestep on 1000000 particles by factor 462754.43948004389 (wall: total = 4.12s force = 2.12s dens = 1.50s link = 0.50s) (cpu : total = 62.07s force = 34.93s dens = 24.39s link = 2.75s) (cpu/wall: total = 15.05 force = 16.44 dens = 16.26 link = 5.49) completed in 4.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 = 2.015E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvxdx.........................OK [max err = 8.214E-04, tol = 1.700E-03] checking dvxdy.........................OK [max err = 1.636E-16, tol = 2.500E-15] checking dvxdz.........................OK [max err = 1.636E-16, tol = 2.500E-15] checking dvydx.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvydy.........................OK [max err = 2.561E-16, tol = 2.500E-15] checking dvydz.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvzdx.........................OK [max err = 4.699E-17, tol = 2.500E-15] checking dvzdy.........................OK [max err = 6.568E-04, tol = 1.500E-03] checking dvzdz.........................OK [max err = 4.498E-17, tol = 2.500E-15] checking viscous force(x)..............OK [max err = 3.949E-02, tol = 4.000E-02] checking viscous force(y)..............OK [max err = 2.994E-02, tol = 3.000E-02] checking viscous force(z)..............OK [max err = 2.730E-02, tol = 3.100E-02] checking mean nneigh...................OK [max err = 0.000E+00, tol = 2.225-308] checking max nneigh....................OK [got 57 should be 57] checking n density calcs...............OK [is 1000000 should be 1000000] checking total nneigh..................OK [is 57000000 should be 57000000] checking \sum v.dv/dt + du/dt = 0......OK [max err = 4.499E-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 MHD derivatives (using B/rho directly) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 (wall: total = 4.12s force = 2.00s dens = 1.62s link = 0.50s) (cpu : total = 60.32s force = 33.12s dens = 24.45s link = 2.75s) (cpu/wall: total = 14.62 force = 16.56 dens = 15.05 link = 5.49) completed in 4.12 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB (symm)...................OK [max err = 8.214E-04, tol = 2.000E-03] checking dBx/dt........................OK [max err = 2.952E-04, tol = 2.000E-03] checking dBy/dt........................OK [max err = 1.161E-03, tol = 2.000E-03] checking dBz/dt........................OK [max err = 2.482E-04, tol = 2.000E-02] checking mhd force(x)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(y)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(z)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking div B (diff)..................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlB(x)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(y)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(z)......................OK [max err = 8.219E-04, tol = 1.000E-03] --> testing MHD derivatives (using B/rho directly) (on 10000 active particles) trial neigh mean : 372.13, real neigh mean = 57.00 ratio try/act= 6.53 trial neigh max : 534, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10000, mean per part = 1.000 (wall: total = 0.62s force = 0.00s dens = 0.12s link = 0.50s) (cpu : total = 4.73s force = 0.67s dens = 1.31s link = 2.75s) (cpu/wall: total = 7.57 force = --- dens = 10.52 link = 5.49) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB (symm)...................OK [max err = 8.214E-04, tol = 2.000E-03] checking dBx/dt........................OK [max err = 2.952E-04, tol = 2.000E-03] checking dBy/dt........................OK [max err = 1.161E-03, tol = 2.000E-03] checking dBz/dt........................OK [max err = 2.482E-04, tol = 2.000E-02] checking mhd force(x)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(y)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(z)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking div B (diff)..................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlB(x)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(y)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(z)......................OK [max err = 8.219E-04, tol = 1.000E-03] --> testing MHD derivatives (using B/rho directly) (on 100 active particles) trial neigh mean : 399.04, real neigh mean = 57.00 ratio try/act= 7.00 trial neigh max : 532, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 100, mean per part = 1.000 (wall: total = 0.62s force = 0.12s dens = 0.00s link = 0.50s) (cpu : total = 5.33s force = 1.27s dens = 1.32s link = 2.75s) (cpu/wall: total = 8.53 force = 10.14 dens = --- link = 5.49) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB (symm)...................OK [max err = 8.214E-04, tol = 2.000E-03] checking dBx/dt........................OK [max err = 2.952E-04, tol = 2.000E-03] checking dBy/dt........................OK [max err = 1.161E-03, tol = 2.000E-03] checking dBz/dt........................OK [max err = 2.482E-04, tol = 2.000E-02] checking mhd force(x)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(y)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(z)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking div B (diff)..................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlB(x)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(y)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(z)......................OK [max err = 8.219E-04, tol = 1.000E-03] --> testing MHD derivatives (using B/rho directly) (on 1 active particles) trial neigh mean : 407.00, real neigh mean = 57.00 ratio try/act= 7.14 trial neigh max : 407, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1, mean per part = 1.000 (wall: total = 0.62s force = 0.12s dens = 0.00s link = 0.50s) (cpu : total = 4.46s force = 0.38s dens = 1.32s link = 2.75s) (cpu/wall: total = 7.13 force = 3.07 dens = --- link = 5.51) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB (symm)...................OK [max err = 8.214E-04, tol = 2.000E-03] checking dBx/dt........................OK [max err = 2.952E-04, tol = 2.000E-03] checking dBy/dt........................OK [max err = 1.161E-03, tol = 2.000E-03] checking dBz/dt........................OK [max err = 2.482E-04, tol = 2.000E-02] checking mhd force(x)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(y)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(z)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking div B (diff)..................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlB(x)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(y)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(z)......................OK [max err = 8.219E-04, tol = 1.000E-03] --> testing artificial resistivity terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 (wall: total = 4.00s force = 2.00s dens = 1.62s link = 0.38s) (cpu : total = 59.74s force = 32.55s dens = 24.44s link = 2.75s) (cpu/wall: total = 14.94 force = 16.27 dens = 15.04 link = 7.34) completed in 4.00 s checking dBx/dt (resist)...............OK [max err = 0.000E+00, tol = 3.700E-02] checking dBy/dt (resist)...............OK [max err = 0.000E+00, tol = 3.400E-02] checking dBz/dt (resist)...............OK [max err = 0.000E+00, tol = 2.200E-01] checking \sum du/dt + B.dB/dt = 0......OK [max err = 0.000E+00, tol = 2.700E-03] --> testing artificial resistivity terms (on 10000 active particles) trial neigh mean : 372.13, real neigh mean = 57.00 ratio try/act= 6.53 trial neigh max : 534, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10000, mean per part = 1.000 (wall: total = 0.62s force = 0.12s dens = 0.00s link = 0.50s) (cpu : total = 4.74s force = 0.67s dens = 1.33s link = 2.75s) (cpu/wall: total = 7.59 force = 5.33 dens = --- link = 5.50) completed in 0.62 s checking dBx/dt (resist)...............OK [max err = 0.000E+00, tol = 3.700E-02] checking dBy/dt (resist)...............OK [max err = 0.000E+00, tol = 3.400E-02] checking dBz/dt (resist)...............OK [max err = 0.000E+00, tol = 2.200E-01] --> testing artificial resistivity terms (on 100 active particles) trial neigh mean : 399.04, real neigh mean = 57.00 ratio try/act= 7.00 trial neigh max : 532, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 100, mean per part = 1.000 (wall: total = 0.50s force = 0.00s dens = 0.12s link = 0.38s) (cpu : total = 4.65s force = 1.32s dens = 0.58s link = 2.75s) (cpu/wall: total = 9.30 force = --- dens = 4.66 link = 7.34) completed in 0.50 s checking dBx/dt (resist)...............OK [max err = 0.000E+00, tol = 3.700E-02] checking dBy/dt (resist)...............OK [max err = 0.000E+00, tol = 3.400E-02] checking dBz/dt (resist)...............OK [max err = 0.000E+00, tol = 2.200E-01] --> testing artificial resistivity terms (on 1 active particles) trial neigh mean : 407.00, real neigh mean = 57.00 ratio try/act= 7.14 trial neigh max : 407, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1, mean per part = 1.000 (wall: total = 0.62s force = 0.12s dens = 0.00s link = 0.50s) (cpu : total = 3.72s force = 0.58s dens = 0.39s link = 2.75s) (cpu/wall: total = 5.95 force = 4.64 dens = --- link = 5.50) completed in 0.62 s checking dBx/dt (resist)...............OK [max err = 0.000E+00, tol = 3.700E-02] checking dBy/dt (resist)...............OK [max err = 0.000E+00, tol = 3.400E-02] checking dBz/dt (resist)...............OK [max err = 0.000E+00, tol = 2.200E-01] --> testing div B cleaning terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 (wall: total = 4.00s force = 2.00s dens = 1.62s link = 0.38s) (cpu : total = 60.19s force = 32.95s dens = 24.48s link = 2.76s) (cpu/wall: total = 15.05 force = 16.47 dens = 15.07 link = 7.35) completed in 4.00 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradpsi_x.....................OK [max err = 4.261E-04, tol = 8.500E-04] checking gradpsi_y.....................OK [max err = 8.500E-04, tol = 8.500E-04] checking gradpsi_z.....................OK [max err = 3.729E-04, tol = 2.000E-03] --> testing div B cleaning terms (on 10000 active particles) trial neigh mean : 372.13, real neigh mean = 57.00 ratio try/act= 6.53 trial neigh max : 534, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10000, mean per part = 1.000 (wall: total = 0.62s force = 0.12s dens = 0.00s link = 0.50s) (cpu : total = 5.67s force = 1.34s dens = 1.58s link = 2.75s) (cpu/wall: total = 9.07 force = 10.71 dens = --- link = 5.50) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradpsi_x.....................OK [max err = 4.261E-04, tol = 8.500E-04] checking gradpsi_y.....................OK [max err = 8.500E-04, tol = 8.500E-04] checking gradpsi_z.....................OK [max err = 3.729E-04, tol = 2.000E-03] --> testing div B cleaning terms (on 100 active particles) trial neigh mean : 399.04, real neigh mean = 57.00 ratio try/act= 7.00 trial neigh max : 532, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 100, mean per part = 1.000 (wall: total = 0.62s force = 0.12s dens = 0.00s link = 0.50s) (cpu : total = 4.73s force = 0.39s dens = 1.55s link = 2.79s) (cpu/wall: total = 7.57 force = 3.12 dens = --- link = 5.57) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradpsi_x.....................OK [max err = 4.261E-04, tol = 8.500E-04] checking gradpsi_y.....................OK [max err = 8.500E-04, tol = 8.500E-04] checking gradpsi_z.....................OK [max err = 3.729E-04, tol = 2.000E-03] --> testing div B cleaning terms (on 1 active particles) trial neigh mean : 407.00, real neigh mean = 57.00 ratio try/act= 7.14 trial neigh max : 407, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1, mean per part = 1.000 (wall: total = 0.62s force = 0.00s dens = 0.12s link = 0.50s) (cpu : total = 4.75s force = 0.39s dens = 1.61s link = 2.75s) (cpu/wall: total = 7.60 force = --- dens = 12.90 link = 5.50) completed in 0.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradpsi_x.....................OK [max err = 4.261E-04, tol = 8.500E-04] checking gradpsi_y.....................OK [max err = 8.500E-04, tol = 8.500E-04] checking gradpsi_z.....................OK [max err = 3.729E-04, tol = 2.000E-03] --> testing Hydro derivs in setup with density contrast ---------------- particles set on 500 x 500 x 500 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-03 dy: 2.000E-03 dz: 2.000E-03 ----------------------------------------------------------------------------------------- ---------------- particles set on 500 x 500 x 500 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-03 dy: 2.000E-03 dz: 2.000E-03 ----------------------------------------------------------------------------------------- ---------------- particles set on 50 x 50 x 50 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-02 dy: 2.000E-02 dz: 2.000E-02 ----------------------------------------------------------------------------------------- thread 0 npart = 648432 in blob = 523984 to test = 73824 WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 430.91, real neigh mean = 57.47 ratio try/act= 7.50 trial neigh max : 88189, max real neigh = 988 n neighbour calls : 124015, mean per part = 1.191 n density calcs : 1452028, mean per part = 2.239 force controlling timestep on 20200 gas particles (wall: total = 5.25s force = 1.75s dens = 3.12s link = 0.38s) (cpu : total = 80.33s force = 28.30s dens = 50.13s link = 1.90s) (cpu/wall: total = 15.30 force = 16.17 dens = 16.04 link = 5.08) completed in 5.25 s checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] checking mean nneigh...................OK [max err = 0.000E+00, tol = 1.000E-17] checking max nneigh....................OK [got 988 should be 988] checking total nneigh..................OK [is 37263216 should be 37263216] --> 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 = 0.50s force = 0.00s dens = 0.12s link = 0.38s) (cpu : total = 3.96s force = 0.40s dens = 1.66s link = 1.90s) (cpu/wall: total = 7.92 force = --- dens = 13.24 link = 5.07) (moved 0.00154% of particles in 9.52381% of time, efficiency = 0.01619%) 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 = 0.50s force = 0.00s dens = 0.12s link = 0.38s) (cpu : total = 3.97s force = 0.40s dens = 1.67s link = 1.90s) (cpu/wall: total = 7.95 force = --- dens = 13.33 link = 5.08) (moved 0.01542% of particles in 9.52381% of time, efficiency = 0.16193%) 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 = 0.38s force = 0.00s dens = 0.12s link = 0.25s) (cpu : total = 4.00s force = 0.43s dens = 1.67s link = 1.90s) (cpu/wall: total = 10.67 force = --- dens = 13.34 link = 7.60) (moved 0.15422% of particles in 7.14286% of time, efficiency = 2.15905%) 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 = 0.50s force = 0.00s dens = 0.12s link = 0.38s) (cpu : total = 4.52s force = 0.77s dens = 1.85s link = 1.91s) (cpu/wall: total = 9.04 force = --- dens = 14.79 link = 5.09) (moved 1.54218% of particles in 9.52381% of time, efficiency = 16.19291%) 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 = 0.75s force = 0.12s dens = 0.25s link = 0.38s) (cpu : total = 8.19s force = 3.08s dens = 3.20s link = 1.90s) (cpu/wall: total = 10.92 force = 24.65 dens = 12.81 link = 5.07) (moved 11.38500% of particles in 14.28571% of time, efficiency = 79.69501%) 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.336E-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 : 611.48, real neigh mean = 50.23 ratio try/act= 12.17 trial neigh max : 1152, max real neigh = 123 n neighbour calls : 55963, mean per part = 1.448 n density calcs : 708021, mean per part = 5.664 force controlling timestep on 10147 gas particles (wall: total = 1.50s force = 0.38s dens = 1.12s link = 0.00s) (cpu : total = 24.24s force = 5.33s dens = 18.30s link = 0.60s) (cpu/wall: total = 16.16 force = 14.21 dens = 16.27 link = ---) evaluating derivs with 12500 particles active... trial neigh mean : 605.02, real neigh mean = 50.22 ratio try/act= 12.05 trial neigh max : 1119, 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 988 gas particles (wall: total = 0.12s force = 0.00s dens = 0.12s link = 0.00s) (cpu : total = 2.12s force = 0.88s dens = 0.64s link = 0.60s) (cpu/wall: total = 16.95 force = --- dens = 5.11 link = ---) checking force(x)......................OK [max err = 6.826E-05, L2 err = 3.180E-11 tol = 1.000E-04] checking force(y)......................OK [max err = 9.126E-05, L2 err = 2.893E-11 tol = 1.000E-04] checking force(z)......................OK [max err = 5.823E-05, L2 err = 2.049E-11 tol = 1.000E-04] checking du/dt.........................OK [max err = 5.419E-08, L2 err = 6.156E-10 tol = 1.000E-05] checking dBx/dt........................OK [max err = 1.912E-07, L2 err = 4.502E-11 tol = 1.000E-05] checking dBy/dt........................OK [max err = 1.058E-06, L2 err = 1.638E-10 tol = 1.000E-05] checking dBz/dt........................OK [max err = 1.568E-07, L2 err = 3.228E-11 tol = 1.000E-05] checking dpsi/dt.......................OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-05] checking div B (symm)..................OK [max err = 6.104E-05, L2 err = 6.196E-11 tol = 1.000E-03] <-- DERIVS TEST COMPLETE --> TESTING EQUATION OF STATE MODULE --> testing equation of state initialisation checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] Reading output_DE_z0.00x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading opacs.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading /opacs.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] Reading helm_data.tab in /home/dprice/phantom-nightly/phantom/data/eos/helmholtz/ (from PHANTOM_DIR setting) checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] --> testing barotropic equation of state Barotropic eq of state: cs_ld = 3.162E-01 code units = 3.643E+04 cm/s Barotropic eq of state: cs = 3.162E-01 code units = 3.643E+04 cm/s Barotropic eq of state: rhocritT == rhoT = 2.514E-04 code units = 5.000E-19 g/cm^3 Barotropic eq of state: rhocrit0 == rho0 = 5.027E-04 code units = 1.000E-18 g/cm^3 Barotropic eq of state: rhocrit1 == rho1 = 5.027E+00 code units = 1.000E-14 g/cm^3 Barotropic eq of state: rhocrit2 == rho2 = 5.027E+04 code units = 1.000E-10 g/cm^3 Barotropic eq of state: rhocrit3 == rho3 = 5.027E+11 code units = 1.000E-03 g/cm^3 Barotropic eq of state: Barotropic eq of state: P = cs_ld*rho for rho/(g/cm^3) < rhoT Barotropic eq of state: P = 10**(log10(cs_bg**2) + M*(log10(rhoT)-log10(rho))) for rhoT <= rho/(g/cm^3) < rho0 Barotropic eq of state: P = cs*rho for rho0 <= rho/(g/cm^3) < rho1 Barotropic eq of state: P = cs*rho1*(rho /rho1)^ 1.400 for rho1 <= rho/(g/cm^3) < rho2 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^ 1.400*(rho /rho2)^ 1.100 for rho2 <= rho/(g/cm^3) < rho3 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^ 1.400*(rho3/rho2)^ 1.100*(rho /rho3)^ 1.667 for rho3 <= rho/(g/cm^3) checking equation of state is continuous......... checking p/rho is continuous...........OK [checked 4999 of 5000 values] <-- EQUATION OF STATE TEST COMPLETE --> TESTING COOLING MODULE --> testing find_in_table routine checking table(i) < val < table(i+1)... <-- COOLING TEST COMPLETE --> TESTING KDTREE <-- KDTREE TEST COMPLETE --> SKIPPING SELF-GRAVITY TESTS (need -DGRAVITY) --> SKIPPING DUST TEST (REQUIRES -DDUST) --> SKIPPING NON-IDEAL MHD TEST (REQUIRES -DNONIDEALMHD) --> TESTING READ/WRITE from dump file --------> TIME = 20.00 : full dump written to file test.dump <-------- >>> reading setup from file: test.dump on unit 9 FT:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 06:17:37.2 n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.3000000000000000 External field found, Bext = 124.00000000000000 125.00000000000000 126.00000000000000 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. time = 20.000000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(2,1001) 7.820 KB real(4) divcurlv(4,1001) 15.641 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,1001) 31.281 KB real(8) Bxyz(3,1001) 23.461 KB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(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,1001) 31.281 KB real(4) divBsumm(1001) 3.910 KB real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,1001) 31.281 KB 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: 738.734 KB --------------------------------------------------------- number of array sizes = 4 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 (MHD) file --> checking header checking read without errors...........OK [got 0 should be 0] checking time..........................OK [max err = 0.000E+00, tol = 2.225-308] checking hfact.........................OK [max err = 0.000E+00, tol = 2.225-308] checking gamma.........................OK [max err = 0.000E+00, tol = 2.225-308] checking polyk.........................OK [max err = 0.000E+00, tol = 2.225-308] checking m(gas)........................OK [max err = 0.000E+00, tol = 2.225-308] checking m(dust).......................OK [max err = 0.000E+00, tol = 2.225-308] checking npartoftype(1)................OK [got 991 should be 991] checking npartoftype(2)................OK [got 10 should be 10] checking umass.........................OK [max err = 0.000E+00, tol = 2.225-308] checking udist.........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bexty.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextz.........................OK [max err = 0.000E+00, tol = 2.225-308] --> checking particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking vx............................OK [max err = 0.000E+00, tol = 2.225-308] checking vy............................OK [max err = 0.000E+00, tol = 2.225-308] checking vz............................OK [max err = 0.000E+00, tol = 2.225-308] checking u.............................OK [max err = 0.000E+00, tol = 2.225-308] checking alpha.........................OK [max err = 0.000E+00, tol = 1.175E-38] checking Bx............................OK [max err = 0.000E+00, tol = 1.175E-38] checking By............................OK [max err = 0.000E+00, tol = 1.175E-38] checking Bz............................OK [max err = 0.000E+00, tol = 1.175E-38] checking psi...........................OK [max err = 0.000E+00, tol = 1.175E-38] checking particle type 1...............OK [checked 991 values] checking particle type 2...............OK [checked 10 values] --> checking sink particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking m.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking hsoft.........................OK [max err = 0.000E+00, tol = 2.225-308] checking maccreted.....................OK [max err = 0.000E+00, tol = 2.225-308] checking spinx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spiny.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spinz.........................OK [max err = 0.000E+00, tol = 2.225-308] checking tlast.........................OK [max err = 0.000E+00, tol = 2.225-308] checking vx............................OK [max err = 0.000E+00, tol = 2.225-308] checking vy............................OK [max err = 0.000E+00, tol = 2.225-308] checking vz............................OK [max err = 0.000E+00, tol = 2.225-308] --> checking read of single array from file x ->y z vx vy vz u checking error flag....................OK [got 0 should be 0] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] --------> TIME = 20.00 : small dump written to file test.dump <-------- --> checking read_dump >>> reading setup from file: test.dump on unit 9 ERROR! read_dump: file is not a Phantom full dump checking read_dump returns is_small_dump error code........OK [got 1978 should be 1978] --> checking read_smalldump >>> reading small dump file: test.dump on unit 9 ST:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 06:17:37.2 WARNING! read_smalldump: *** VELOCITY WILL BE MISSING FROM SMALL DUMP FILES *** n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.2999999523162842 comment: hfact was 1.1230000257492065 now 1.1230000000000000 comment: tolh was 9.9999997473787516E-005 now 1.0000000000000000E-004 comment: C_cour was 0.30000001192092896 now 0.29999999999999999 External field found, Bext = 124.00000000000000 125.00000000000000 126.00000000000000 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. time = 20.000000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(2,1001) 7.820 KB real(4) divcurlv(4,1001) 15.641 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,1001) 31.281 KB real(8) Bxyz(3,1001) 23.461 KB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(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,1001) 31.281 KB real(4) divBsumm(1001) 3.910 KB real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,1001) 31.281 KB 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: 742.645 KB --------------------------------------------------------- number of array sizes = 4 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: alpha not found in file WARNING! sink particle velocities not found got 11 sink properties from 10 sink particles ----------------------------------------------- ID| Mass | Racc | Macc | hsoft | ----------------------------------------------- 1| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 2| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 3| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 4| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 5| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 6| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 7| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 8| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 9| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 10| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | ----------------------------------------------- WARNING! div B cleaning field (Psi) not found in Phantom dump file: assuming psi=0 <<< finished reading (MHD) file --> checking header checking read without errors...........OK [got 0 should be 0] checking time..........................OK [max err = 0.000E+00, tol = 1.192E-07] checking hfact.........................OK [max err = 2.293E-08, tol = 1.192E-07] checking gamma.........................OK [max err = 3.668E-08, tol = 1.192E-07] checking polyk.........................OK [max err = 1.445E-08, tol = 1.192E-07] checking m(gas)........................OK [max err = 0.000E+00, tol = 1.192E-07] checking m(dust).......................OK [max err = 0.000E+00, tol = 1.192E-07] checking npartoftype(1)................OK [got 991 should be 991] checking npartoftype(2)................OK [got 10 should be 10] checking umass.........................OK [max err = 0.000E+00, tol = 2.225-308] checking udist.........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bexty.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextz.........................OK [max err = 0.000E+00, tol = 2.225-308] --> checking particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking particle type 1...............OK [checked 991 values] checking particle type 2...............OK [checked 10 values] --> checking sink particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking m.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking hsoft.........................OK [max err = 0.000E+00, tol = 2.225-308] checking maccreted.....................OK [max err = 0.000E+00, tol = 2.225-308] checking spinx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spiny.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spinz.........................OK [max err = 0.000E+00, tol = 2.225-308] checking tlast.........................OK [max err = 0.000E+00, tol = 2.225-308] <-- READ/WRITE TEST COMPLETE --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1000000) 30.518 MB real(8) xyzh_soa(1000000,4) 30.518 MB real(8) vxyzu(4,1000000) 30.518 MB real(4) alphaind(2,1000000) 7.629 MB real(4) divcurlv(4,1000000) 15.259 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,1000000) 30.518 MB real(8) Bxyz(3,1000000) 22.888 MB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(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,1000000) 30.518 MB real(4) divBsumm(1000000) 3.815 MB real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,1000000) 30.518 MB 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: 558.870 MB --------------------------------------------------------- --> TESTING STEP MODULE / boundary crossing ---------------- particles set on 50 x 50 x 50 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-02 dy: 2.000E-02 dz: 2.000E-02 ----------------------------------------------------------------------------------------- t = 0.20000000000000001 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.40000000000000002 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.60000000000000009 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.80000000000000004 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.0000000000000000 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.2000000000000000 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.3999999999999999 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.5999999999999999 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.7999999999999998 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.9999999999999998 dt = 0.20000000000000001 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] <-- STEP TEST COMPLETE --> 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.2500000000000001E-002 checking potential energy..............OK [max err = 0.000E+00, tol = 2.220E-16] checking angular momentum..............OK [max err = 0.000E+00, tol = 2.220E-10] nsteps per orbit = 711 norbits = 100 checking angular momentum..............OK [max err = 6.594E-15, tol = 3.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] checking total energy..................OK [max err = 2.440E-08, tol = 3.000E-08] --> skipping circumbinary disc test (-DPERIODIC is set) --> skipping circumbinary disc test (-DPERIODIC is set) --> testing softening in sink particle binary checking potential energy..............OK [max err = 3.391E-16, tol = 4.441E-16] nsteps per orbit = 1099 norbits = 10 checking angular momentum..............OK [max err = 6.324E-15, tol = 2.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.225-308] checking total energy..................OK [max err = 1.021E-10, tol = 2.000E-09] --> testing accretion onto sink particles checking accretion flag................OK [got T should be T] checking isdead_or_accreted flag.......OK [got T should be T] checking x(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking y(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking z(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking vx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (uniform density) ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.500E-02 -2.500E-02 -2.500E-02 0.000E+00 ) by accreting 93 particles: M= 3.321E-01 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (sin) >>>>>> s t r e t c h m a p p i n g <<<<<< stretching to match density profile in r direction density at r = 0.0000000000000000 is 1.0000000000000000 total mass = 6.9604063530853081E-004 >>>>>> done ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.765E-18 6.370E-19 1.870E-18 0.000E+00 ) by accreting 280 particles: M= 1.000E+00 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] <-- PTMASS TEST COMPLETE --> TESTING GNEWTON MODULE --> testing relativistic orbit ---------- orbital parameters ---------- apocenter : 180.000 pericenter : 20.000 semi-major axis : 100.000 eccentricity : 0.800 initial velocity : 36.724 gravitational radius : 2.122 central mass : 1.000E+06 ---------------------------------------- integrating the orbit for one azimuthal period... ( C_force = 1.0000000000000000E-004 ) checking energy........................OK [max err = 9.701E-06, tol = 1.000E-03] checking angular momentum..............OK [max err = 5.863E-07, tol = 1.000E-03] checking precession angle..............OK [max err = 5.762E-07, tol = 1.000E-05] <-- GNEWTON TEST COMPLETE --> TESTING COROTATION MODULE --> testing centrifugal force with rotation about z axis checking fx=Omega^2 x..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=Omega^2 y..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] checking phi=1/2 Omega^2 R^2...........OK [max err = 2.321E-16, tol = 4.441E-16] --> testing coriolis force with rotation about z axis checking fx=2 Omega vy.................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=-2 Omega vx................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] <-- COROTATE TEST COMPLETE --> SKIPPING DISC SETUP TESTS --> TESTING GEOMETRY MODULE --> checking transforms for coords and vectors are reversible checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 8.882E-16, L2 err = 7.326E-17 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 8.882E-16, L2 err = 1.433E-16 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 4.441E-16, L2 err = 6.410E-17 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 2.538E-16, L2 err = 1.465E-16 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 4.441E-16, L2 err = 7.166E-17 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 8.882E-16, L2 err = 1.036E-16 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 8.882E-16, L2 err = 1.282E-16 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 8.882E-16, L2 err = 1.794E-16 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 8.882E-16, L2 err = 1.923E-16 tol = 1.000E-15] --> checking transforms for coords and vectors are reversible at r=0 checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 1.225E-16, L2 err = 1.500E-32 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 4.441E-16, L2 err = 9.065E-17 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 4.441E-16, L2 err = 7.166E-17 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 1.000E-15, L2 err = 1.000E-30 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] --> testing get_coord_limits matches coord_transform checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] <-- GEOMETRY TEST COMPLETE --> testing Sedov blast wave ---------------- particles set on 16 x 16 x 16 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 6.250E-02 dy: 6.250E-02 dz: 6.250E-02 ----------------------------------------------------------------------------------------- npart = 4096 particle mass = 2.4414062500000000E-004 > step 1 / 128 t = 0.7812500E-03 dt = 7.81E-04 moved 4096 in 0.26 cpu-s < | np = 4096 | > step 3 / 256 t = 0.1171875E-02 dt = 3.91E-04 moved 56 in 0.057 cpu-s < > step 4 / 256 t = 0.1562500E-02 dt = 3.91E-04 moved 4096 in 0.27 cpu-s < > step 5 / 256 t = 0.1953125E-02 dt = 3.91E-04 moved 56 in 0.057 cpu-s < > step 6 / 256 t = 0.2343750E-02 dt = 3.91E-04 moved 352 in 0.076 cpu-s < > step 7 / 256 t = 0.2734375E-02 dt = 3.91E-04 moved 56 in 0.057 cpu-s < > step 8 / 256 t = 0.3125000E-02 dt = 3.91E-04 moved 4096 in 0.28 cpu-s < > step 9 / 256 t = 0.3515625E-02 dt = 3.91E-04 moved 56 in 0.057 cpu-s < > step 10 / 256 t = 0.3906250E-02 dt = 3.91E-04 moved 400 in 0.083 cpu-s < > step 11 / 256 t = 0.4296875E-02 dt = 3.91E-04 moved 32 in 0.054 cpu-s < > step 12 / 256 t = 0.4687500E-02 dt = 3.91E-04 moved 1232 in 0.14 cpu-s < > step 13 / 256 t = 0.5078125E-02 dt = 3.91E-04 moved 32 in 0.054 cpu-s < > step 14 / 256 t = 0.5468750E-02 dt = 3.91E-04 moved 408 in 0.083 cpu-s < > step 15 / 256 t = 0.5859375E-02 dt = 3.91E-04 moved 8 in 0.052 cpu-s < > step 16 / 256 t = 0.6250000E-02 dt = 3.91E-04 moved 4096 in 0.28 cpu-s < > step 17 / 256 t = 0.6640625E-02 dt = 3.91E-04 moved 8 in 0.051 cpu-s < > step 18 / 256 t = 0.7031250E-02 dt = 3.91E-04 moved 408 in 0.082 cpu-s < > step 19 / 256 t = 0.7421875E-02 dt = 3.91E-04 moved 8 in 0.052 cpu-s < > step 20 / 256 t = 0.7812500E-02 dt = 3.91E-04 moved 1232 in 0.14 cpu-s < > step 21 / 256 t = 0.8203125E-02 dt = 3.91E-04 moved 8 in 0.052 cpu-s < > step 22 / 256 t = 0.8593750E-02 dt = 3.91E-04 moved 456 in 0.087 cpu-s < > step 12 / 128 t = 0.9375000E-02 dt = 7.81E-04 moved 2728 in 0.23 cpu-s < > step 13 / 128 t = 0.1015625E-01 dt = 7.81E-04 moved 480 in 0.090 cpu-s < > step 14 / 128 t = 0.1093750E-01 dt = 7.81E-04 moved 1520 in 0.16 cpu-s < > step 15 / 128 t = 0.1171875E-01 dt = 7.81E-04 moved 456 in 0.088 cpu-s < > step 16 / 128 t = 0.1250000E-01 dt = 7.81E-04 moved 4096 in 0.29 cpu-s < > step 17 / 128 t = 0.1328125E-01 dt = 7.81E-04 moved 456 in 0.089 cpu-s < > step 18 / 128 t = 0.1406250E-01 dt = 7.81E-04 moved 1520 in 0.16 cpu-s < > step 19 / 128 t = 0.1484375E-01 dt = 7.81E-04 moved 456 in 0.089 cpu-s < > step 20 / 128 t = 0.1562500E-01 dt = 7.81E-04 moved 2944 in 0.26 cpu-s < > step 21 / 128 t = 0.1640625E-01 dt = 7.81E-04 moved 408 in 0.086 cpu-s < > step 22 / 128 t = 0.1718750E-01 dt = 7.81E-04 moved 1568 in 0.17 cpu-s < > step 23 / 128 t = 0.1796875E-01 dt = 7.81E-04 moved 376 in 0.085 cpu-s < > step 24 / 128 t = 0.1875000E-01 dt = 7.81E-04 moved 3816 in 0.30 cpu-s < > step 25 / 128 t = 0.1953125E-01 dt = 7.81E-04 moved 328 in 0.084 cpu-s < > step 26 / 128 t = 0.2031250E-01 dt = 7.81E-04 moved 1616 in 0.17 cpu-s < > step 27 / 128 t = 0.2109375E-01 dt = 7.81E-04 moved 232 in 0.078 cpu-s < > step 28 / 128 t = 0.2187500E-01 dt = 7.81E-04 moved 3088 in 0.27 cpu-s < > step 29 / 128 t = 0.2265625E-01 dt = 7.81E-04 moved 256 in 0.081 cpu-s < > step 30 / 128 t = 0.2343750E-01 dt = 7.81E-04 moved 1784 in 0.19 cpu-s < > step 31 / 128 t = 0.2421875E-01 dt = 7.81E-04 moved 152 in 0.071 cpu-s < > step 32 / 128 t = 0.2500000E-01 dt = 7.81E-04 moved 4096 in 0.32 cpu-s < > step 33 / 128 t = 0.2578125E-01 dt = 7.81E-04 moved 80 in 0.065 cpu-s < > step 34 / 128 t = 0.2656250E-01 dt = 7.81E-04 moved 1760 in 0.19 cpu-s < > step 35 / 128 t = 0.2734375E-01 dt = 7.81E-04 moved 56 in 0.062 cpu-s < > step 36 / 128 t = 0.2812500E-01 dt = 7.81E-04 moved 3304 in 0.29 cpu-s < > step 37 / 128 t = 0.2890625E-01 dt = 7.81E-04 moved 32 in 0.058 cpu-s < > step 38 / 128 t = 0.2968750E-01 dt = 7.81E-04 moved 1328 in 0.16 cpu-s < > step 39 / 128 t = 0.3046875E-01 dt = 7.81E-04 moved 32 in 0.058 cpu-s < > step 40 / 128 t = 0.3125000E-01 dt = 7.81E-04 moved 3936 in 0.32 cpu-s < > step 41 / 128 t = 0.3203125E-01 dt = 7.81E-04 moved 8 in 0.053 cpu-s < > step 42 / 128 t = 0.3281250E-01 dt = 7.81E-04 moved 1424 in 0.17 cpu-s < > step 22 / 64 t = 0.3437500E-01 dt = 0.0016 moved 2944 in 0.27 cpu-s < > step 23 / 64 t = 0.3593750E-01 dt = 0.0016 moved 1472 in 0.17 cpu-s < > step 24 / 64 t = 0.3750000E-01 dt = 0.0016 moved 4088 in 0.33 cpu-s < > step 25 / 64 t = 0.3906250E-01 dt = 0.0016 moved 1568 in 0.18 cpu-s < > step 26 / 64 t = 0.4062500E-01 dt = 0.0016 moved 3016 in 0.28 cpu-s < > step 27 / 64 t = 0.4218750E-01 dt = 0.0016 moved 1616 in 0.18 cpu-s < > step 28 / 64 t = 0.4375000E-01 dt = 0.0016 moved 3816 in 0.33 cpu-s < > step 29 / 64 t = 0.4531250E-01 dt = 0.0016 moved 1640 in 0.18 cpu-s < > step 30 / 64 t = 0.4687500E-01 dt = 0.0016 moved 3232 in 0.29 cpu-s < > step 31 / 64 t = 0.4843750E-01 dt = 0.0016 moved 1688 in 0.19 cpu-s < > step 32 / 64 t = 0.5000000E-01 dt = 0.0016 moved 4096 in 0.35 cpu-s < > step 33 / 64 t = 0.5156250E-01 dt = 0.0016 moved 1808 in 0.20 cpu-s < > step 34 / 64 t = 0.5312500E-01 dt = 0.0016 moved 3352 in 0.31 cpu-s < > step 35 / 64 t = 0.5468750E-01 dt = 0.0016 moved 1640 in 0.19 cpu-s < > step 36 / 64 t = 0.5625000E-01 dt = 0.0016 moved 3936 in 0.34 cpu-s < > step 37 / 64 t = 0.5781250E-01 dt = 0.0016 moved 1776 in 0.20 cpu-s < > step 38 / 64 t = 0.5937500E-01 dt = 0.0016 moved 3400 in 0.32 cpu-s < > step 39 / 64 t = 0.6093750E-01 dt = 0.0016 moved 1728 in 0.20 cpu-s < > step 40 / 64 t = 0.6250000E-01 dt = 0.0016 moved 4088 in 0.36 cpu-s < > step 41 / 64 t = 0.6406250E-01 dt = 0.0016 moved 1680 in 0.20 cpu-s < > step 42 / 64 t = 0.6562500E-01 dt = 0.0016 moved 3600 in 0.34 cpu-s < > step 43 / 64 t = 0.6718750E-01 dt = 0.0016 moved 1832 in 0.21 cpu-s < > step 44 / 64 t = 0.6875000E-01 dt = 0.0016 moved 4016 in 0.36 cpu-s < > step 45 / 64 t = 0.7031250E-01 dt = 0.0016 moved 1832 in 0.21 cpu-s < > step 46 / 64 t = 0.7187500E-01 dt = 0.0016 moved 3600 in 0.34 cpu-s < > step 47 / 64 t = 0.7343750E-01 dt = 0.0016 moved 1784 in 0.21 cpu-s < > step 48 / 64 t = 0.7500000E-01 dt = 0.0016 moved 4096 in 0.38 cpu-s < > step 49 / 64 t = 0.7656250E-01 dt = 0.0016 moved 1544 in 0.20 cpu-s < > step 50 / 64 t = 0.7812500E-01 dt = 0.0016 moved 3744 in 0.36 cpu-s < > step 51 / 64 t = 0.7968750E-01 dt = 0.0016 moved 1232 in 0.17 cpu-s < > step 52 / 64 t = 0.8125000E-01 dt = 0.0016 moved 4064 in 0.38 cpu-s < > step 53 / 64 t = 0.8281250E-01 dt = 0.0016 moved 1096 in 0.16 cpu-s < > step 54 / 64 t = 0.8437500E-01 dt = 0.0016 moved 3744 in 0.36 cpu-s < > step 55 / 64 t = 0.8593750E-01 dt = 0.0016 moved 968 in 0.15 cpu-s < > step 56 / 64 t = 0.8750000E-01 dt = 0.0016 moved 4096 in 0.39 cpu-s < > step 57 / 64 t = 0.8906250E-01 dt = 0.0016 moved 992 in 0.15 cpu-s < > step 58 / 64 t = 0.9062500E-01 dt = 0.0016 moved 3696 in 0.37 cpu-s < > step 59 / 64 t = 0.9218750E-01 dt = 0.0016 moved 728 in 0.13 cpu-s < > step 60 / 64 t = 0.9375000E-01 dt = 0.0016 moved 4064 in 0.40 cpu-s < > step 61 / 64 t = 0.9531250E-01 dt = 0.0016 moved 704 in 0.12 cpu-s < > step 62 / 64 t = 0.9687500E-01 dt = 0.0016 moved 3600 in 0.37 cpu-s < > step 63 / 64 t = 0.9843750E-01 dt = 0.0016 moved 632 in 0.12 cpu-s < > step 64 / 64 t = 0.1000000 dt = 0.0016 moved 4096 in 0.40 cpu-s < --------> TIME = 0.1000 : full dump written to file test002 <-------- input file test.in written successfully. Since code start: 95 timesteps, wall: 1.1s cpu: 19s cpu/wall: 17 Since last dump : 95 timesteps, wall: 1.1s cpu: 19s cpu/wall: 16 wall cpu cpu/wall frac step : 1.12s 18.47s 16.42 90.00% step (force): 0.12s 6.54s 52.35 10.00% step (dens) : 0.88s 7.58s 8.67 70.00% step (link) : 0.12s 3.97s 31.79 10.00% write_dump : 0.12s 0.03s 0.21 10.00% ------------------------------------------------------------------------------ |** Number of steps since last summary: 95 **| |** Wall time since last summary: 1.125 seconds **| ------------------------------------------------------------------------------ |* particles woken *| | #steps | mean # part/step | max # part/step | | 24| 97.00| 216 | ------------------------------------------------------------------------------ ------------------------------------------------- | bin | dt | npart | frac | cpufrac | ------------------------------------------------- | 0 1.000E-01 0 0.00% 2.46% | | 1 5.000E-02 0 0.00% 2.18% | | 2 2.500E-02 0 0.00% 4.31% | | 3 1.250E-02 80 1.95% 8.55% | | 4 6.250E-03 488 11.91% 16.83% | | 5 3.125E-03 2920 71.29% 30.60% | | 6 1.563E-03 608 14.84% 35.06% | ------------------------------------------------- Etot= 9.668E-01, Ekin= 1.982E-01, Etherm= 7.686E-01, Epot= 0.000E+00 Emag= 0.000E+00 Linm= 1.266E-17, Angm= 1.349E-18 Centre of Mass = 1.626E-18, 3.334E-18, 0.000E+00 density (max)= 1.603E+00 (mean)= 1.159E+00 (max)= 9.456E+00 g/cm^3 alpha(max)= 1.000E+00 RMS Mach #= 6.243E-01 div B (max)= 0.000E+00, div B (mean)= 0.000E+00 h|div B|/B (max)= 0.000E+00, h|div B|/B (mean)= 0.000E+00 --> TESTING M_4 cubic kernel <-- KERNEL TEST COMPLETE checking total energy..................OK [max err = 2.831E-04, tol = 4.700E-04] checking linear momentum...............OK [max err = 1.266E-17, tol = 7.000E-15] <-- testing complete total wall time = 3 min, 39.00 s (= 2.1900E+02s) total cpu time = 19 min, 20.59 s (= 1.1606E+03s) SUMMARY OF ALL TESTS: PASSED: 151 of 151 100.0% FAILED: 0 of 151 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found