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: 589.388 MB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) --> TESTING M_4 cubic kernel checking wab0 = wab(0).................OK [max err = 0.000E+00, tol = 2.225-308] checking gradh0 = -3.*wab(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking dphidh0 = dphidh(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking radkern2 = radkern*radkern....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == wkern...........OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == grkern..........OK [max err = 0.000E+00, tol = 2.225-308] checking w gradient equal to gradw.....OK [max err = 1.499E-07, tol = 2.000E-05] checking potential gradient = force....OK [max err = 6.722E-08, tol = 2.000E-05] checking get_kernel_grav1 == wkern.....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel_grav1 == grkern....OK [max err = 0.000E+00, tol = 2.225-308] checking dphi/dh = phi - q*dphi/dq.....OK [max err = 8.105E-15, tol = 2.000E-07] <-- KERNEL TEST COMPLETE --> TESTING LINKLIST / NEIGHBOUR FINDING ----------------------------------------------------------------------------------------- random seed = -43587 32768 particles set in uniform random distribution ----------------------------------------------------------------------------------------- Test 1 of 3: building linked list... checking ncells > 0....................OK [got T should be T] checking active/inactive cells......... checking active cells have at least one active particle....OK [checked 4326 of 32768 values] checking neighbour number.............. checking nneigh (cached)...............OK [checked 32768 values] checking nneigh (no cache).............OK [checked 32768 of 32768 values] max nneigh = 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 = 18.50s force = 6.50s dens = 9.25s link = 2.75s) (cpu : total = 288.37s force = 104.04s dens = 149.79s link = 34.54s) (cpu/wall: total = 15.59 force = 16.01 dens = 16.19 link = 12.56) completed in 18.50 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.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 = 2.75s force = 0.00s dens = 0.00s link = 2.75s) (cpu : total = 35.63s force = 0.57s dens = 0.58s link = 34.49s) (cpu/wall: total = 12.96 force = --- dens = --- link = 12.54) (moved 0.00010% of particles in 14.86486% of time, efficiency = 0.00067%) 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 = 2.75s force = 0.00s dens = 0.00s link = 2.75s) (cpu : total = 35.65s force = 0.56s dens = 0.58s link = 34.51s) (cpu/wall: total = 12.96 force = --- dens = --- link = 12.55) (moved 0.00100% of particles in 14.86486% of time, efficiency = 0.00673%) 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 = 2.88s force = 0.00s dens = 0.12s link = 2.75s) (cpu : total = 35.64s force = 0.56s dens = 0.58s link = 34.50s) (cpu/wall: total = 12.40 force = --- dens = 4.60 link = 12.55) (moved 0.01000% of particles in 15.54054% of time, efficiency = 0.06435%) 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 = 2.88s force = 0.00s dens = 0.12s link = 2.75s) (cpu : total = 35.81s force = 0.64s dens = 0.70s link = 34.47s) (cpu/wall: total = 12.46 force = --- dens = 5.61 link = 12.53) (moved 0.10000% of particles in 15.54054% of time, efficiency = 0.64348%) 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 = 3.00s force = 0.12s dens = 0.12s link = 2.75s) (cpu : total = 38.25s force = 1.59s dens = 2.16s link = 34.49s) (cpu/wall: total = 12.75 force = 12.75 dens = 17.30 link = 12.54) (moved 1.00000% of particles in 16.21622% of time, efficiency = 6.16667%) 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 = 4.50s force = 0.75s dens = 1.00s link = 2.75s) (cpu : total = 61.05s force = 10.92s dens = 15.64s link = 34.48s) (cpu/wall: total = 13.57 force = 14.56 dens = 15.64 link = 12.54) (moved 10.00000% of particles in 24.32432% of time, efficiency = 41.11111%) 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 = 14.00s force = 6.38s dens = 4.88s link = 2.75s) (cpu : total = 215.18s force = 103.04s dens = 77.41s link = 34.73s) (cpu/wall: total = 15.37 force = 16.16 dens = 15.88 link = 12.63) completed in 14.38 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.670E-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 = 2.88s force = 0.12s dens = 0.00s link = 2.75s) (cpu : total = 37.47s force = 1.60s dens = 1.34s link = 34.54s) (cpu/wall: total = 13.03 force = 12.79 dens = --- link = 12.56) completed in 2.88 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.670E-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 = 2.88s force = 0.12s dens = 0.00s link = 2.75s) (cpu : total = 35.61s force = 0.56s dens = 0.57s link = 34.47s) (cpu/wall: total = 12.39 force = 4.50 dens = --- link = 12.54) completed in 2.88 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.670E-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 = 2.75s force = 0.00s dens = 0.00s link = 2.75s) (cpu : total = 35.63s force = 0.57s dens = 0.58s link = 34.49s) (cpu/wall: total = 12.96 force = --- dens = --- link = 12.54) completed in 2.75 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.670E-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 7.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.730E-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.439480044 (wall: total = 14.75s force = 7.12s dens = 4.88s link = 2.75s) (cpu : total = 226.13s force = 113.88s dens = 77.55s link = 34.69s) (cpu/wall: total = 15.33 force = 15.98 dens = 15.91 link = 12.62) completed in 14.75 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.621E-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.426E-16, tol = 2.500E-15] checking dvxdz.........................OK [max err = 1.426E-16, tol = 2.500E-15] checking dvydx.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvydy.........................OK [max err = 2.342E-16, tol = 2.500E-15] checking dvydz.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvzdx.........................OK [max err = 4.243E-17, tol = 2.500E-15] checking dvzdy.........................OK [max err = 6.568E-04, tol = 1.500E-03] checking dvzdz.........................OK [max err = 4.328E-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.501E-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 = 14.00s force = 6.38s dens = 4.88s link = 2.75s) (cpu : total = 214.99s force = 103.03s dens = 77.44s link = 34.52s) (cpu/wall: total = 15.36 force = 16.16 dens = 15.89 link = 12.55) completed in 14.00 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 = 2.88s force = 0.00s dens = 0.12s link = 2.75s) (cpu : total = 37.43s force = 1.59s dens = 1.34s link = 34.50s) (cpu/wall: total = 13.02 force = --- dens = 10.68 link = 12.55) completed in 2.88 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 = 2.75s force = 0.00s dens = 0.12s link = 2.62s) (cpu : total = 35.66s force = 0.56s dens = 0.58s link = 34.51s) (cpu/wall: total = 12.97 force = --- dens = 4.66 link = 13.15) completed in 2.75 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 = 2.75s force = 0.00s dens = 0.00s link = 2.75s) (cpu : total = 35.90s force = 0.57s dens = 0.58s link = 34.75s) (cpu/wall: total = 13.05 force = --- dens = --- link = 12.63) completed in 2.75 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 = 14.00s force = 6.50s dens = 4.75s link = 2.75s) (cpu : total = 214.54s force = 102.48s dens = 77.56s link = 34.50s) (cpu/wall: total = 15.32 force = 15.77 dens = 16.33 link = 12.55) completed in 14.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 = 2.88s force = 0.00s dens = 0.12s link = 2.75s) (cpu : total = 37.43s force = 1.59s dens = 1.34s link = 34.50s) (cpu/wall: total = 13.02 force = --- dens = 10.71 link = 12.55) completed in 2.88 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 = 2.88s force = 0.00s dens = 0.12s link = 2.75s) (cpu : total = 35.64s force = 0.57s dens = 0.58s link = 34.49s) (cpu/wall: total = 12.40 force = --- dens = 4.67 link = 12.54) completed in 2.88 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 = 2.88s force = 0.12s dens = 0.00s link = 2.75s) (cpu : total = 35.81s force = 0.56s dens = 0.59s link = 34.66s) (cpu/wall: total = 12.45 force = 4.51 dens = --- link = 12.60) completed in 2.88 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 = 14.00s force = 6.50s dens = 4.75s link = 2.75s) (cpu : total = 214.81s force = 102.91s dens = 77.41s link = 34.49s) (cpu/wall: total = 15.34 force = 15.83 dens = 16.30 link = 12.54) completed in 14.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 = 2.88s force = 0.12s dens = 0.00s link = 2.75s) (cpu : total = 37.42s force = 1.59s dens = 1.34s link = 34.49s) (cpu/wall: total = 13.01 force = 12.70 dens = --- link = 12.54) completed in 2.88 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 = 2.88s force = 0.12s dens = 0.00s link = 2.75s) (cpu : total = 35.60s force = 0.56s dens = 0.58s link = 34.47s) (cpu/wall: total = 12.38 force = 4.48 dens = --- link = 12.53) completed in 2.88 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 = 2.88s force = 0.00s dens = 0.00s link = 2.88s) (cpu : total = 35.81s force = 0.57s dens = 0.57s link = 34.67s) (cpu/wall: total = 12.45 force = --- dens = --- link = 12.06) completed in 2.88 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 = 15.75s force = 5.00s dens = 8.75s link = 2.00s) (cpu : total = 245.51s force = 79.58s dens = 140.48s link = 25.45s) (cpu/wall: total = 15.59 force = 15.92 dens = 16.05 link = 12.73) completed in 15.75 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 = 2.12s force = 0.00s dens = 0.12s link = 2.00s) (cpu : total = 27.70s force = 0.58s dens = 0.59s link = 26.53s) (cpu/wall: total = 13.03 force = --- dens = 4.74 link = 13.26) (moved 0.00154% of particles in 13.49206% of time, efficiency = 0.01143%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.202E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] --> testing Hydro derivs in setup with density contrast (nactive= 100) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 351.83, real neigh mean = 57.00 ratio try/act= 6.17 trial neigh max : 489, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 100, mean per part = 1.000 (wall: total = 2.00s force = 0.00s dens = 0.00s link = 2.00s) (cpu : total = 27.68s force = 0.57s dens = 0.59s link = 26.52s) (cpu/wall: total = 13.84 force = --- dens = --- link = 13.26) (moved 0.01542% of particles in 12.69841% of time, efficiency = 0.12145%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.205E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] --> testing Hydro derivs in setup with density contrast (nactive= 1000) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 354.22, real neigh mean = 57.00 ratio try/act= 6.21 trial neigh max : 500, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000, mean per part = 1.000 (wall: total = 2.12s force = 0.00s dens = 0.12s link = 2.00s) (cpu : total = 27.84s force = 0.66s dens = 0.65s link = 26.53s) (cpu/wall: total = 13.10 force = --- dens = 5.20 link = 13.27) (moved 0.15422% of particles in 13.49206% of time, efficiency = 1.14303%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.206E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] --> testing Hydro derivs in setup with density contrast (nactive= 10000) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 353.60, real neigh mean = 57.00 ratio try/act= 6.20 trial neigh max : 500, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 10000, mean per part = 1.000 (wall: total = 2.25s force = 0.12s dens = 0.12s link = 2.00s) (cpu : total = 29.56s force = 1.65s dens = 1.36s link = 26.55s) (cpu/wall: total = 13.14 force = 13.20 dens = 10.88 link = 13.28) (moved 1.54218% of particles in 14.28571% of time, efficiency = 10.79527%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.206E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] --> testing Hydro derivs in setup with density contrast (nactive= 73824) WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 344.82, real neigh mean = 57.00 ratio try/act= 6.05 trial neigh max : 500, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 73824, mean per part = 1.000 (wall: total = 2.88s force = 0.50s dens = 0.38s link = 2.00s) (cpu : total = 41.16s force = 8.40s dens = 6.22s link = 26.54s) (cpu/wall: total = 14.32 force = 16.79 dens = 16.60 link = 13.27) (moved 11.38500% of particles in 18.25397% of time, efficiency = 62.37001%) checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.260E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.206E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 7.283E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.205E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.289E-05, tol = 1.000E-03] --> testing force evaluation with ind_timesteps ----------------------------------------------------------------------------------------- random seed = -43587 125000 particles set in uniform random distribution ----------------------------------------------------------------------------------------- evaluating derivs with all particles active... trial neigh mean : 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 = 4.62s force = 0.88s dens = 3.38s link = 0.38s) (cpu : total = 74.48s force = 15.37s dens = 54.32s link = 4.78s) (cpu/wall: total = 16.10 force = 17.57 dens = 16.10 link = 12.75) 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.62s force = 0.12s dens = 0.12s link = 0.38s) (cpu : total = 10.12s force = 2.74s dens = 1.99s link = 5.39s) (cpu/wall: total = 16.19 force = 21.94 dens = 15.89 link = 14.37) 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 <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004947BA dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000966D52 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000497F17 dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000966D52 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000499191 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000966D52 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000049A40F dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000966D52 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000049CEF0 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 0000000000966D52 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004A1FFB dump_utils_mp_wri 1681 utils_dumpfiles.f90 phantom 000000000072EC9A readwrite_dumps_m 539 readwrite_dumps.F90 phantom 0000000000966D52 testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown >>> reading setup from file: test.dump on unit 9 FT:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 00:46:33.8 forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000047C3CF dump_utils_mp_rea 1225 utils_dumpfiles.f90 phantom 00000000007382AB readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000968B18 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000481BD3 dump_utils_mp_rea 1246 utils_dumpfiles.f90 phantom 00000000007382AB readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000968B18 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000048393E dump_utils_mp_rea 1253 utils_dumpfiles.f90 phantom 00000000007382AB readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000968B18 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004856AD dump_utils_mp_rea 1260 utils_dumpfiles.f90 phantom 00000000007382AB readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000968B18 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000489CC7 dump_utils_mp_rea 1281 utils_dumpfiles.f90 phantom 00000000007382AB readwrite_dumps_m 782 readwrite_dumps.F90 phantom 0000000000968B18 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.30000000000000 External field found, Bext = 124.000000000000 125.000000000000 126.000000000000 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. time = 20.0000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(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 forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004A7C78 dump_utils_mp_rea 1884 utils_dumpfiles.f90 phantom 0000000000740B27 readwrite_dumps_m 1288 readwrite_dumps.F90 phantom 00000000007398EB readwrite_dumps_m 898 readwrite_dumps.F90 phantom 0000000000968B18 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004ABA26 dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000741176 readwrite_dumps_m 1305 readwrite_dumps.F90 phantom 00000000007398EB readwrite_dumps_m 898 readwrite_dumps.F90 phantom 0000000000968B18 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown 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 forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004B1630 dump_utils_mp_rea 2103 utils_dumpfiles.f90 phantom 000000000096DC8E testrwdump_mp_tes 309 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown ->y z vx vy vz u checking error flag....................OK [got 0 should be 0] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] --------> TIME = 20.00 : small dump written to file test.dump <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004947BA dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 0000000000735253 readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000966CF1 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000497F17 dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 0000000000735253 readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000966CF1 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000499191 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 0000000000735253 readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000966CF1 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000049A40F dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 0000000000735253 readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000966CF1 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000049CEF0 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 0000000000735253 readwrite_dumps_m 631 readwrite_dumps.F90 phantom 0000000000966CF1 testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown --> checking read_dump >>> reading setup from file: test.dump on unit 9 ERROR! read_dump: file is not a Phantom full dump checking read_dump returns is_small_dump error code........OK [got 1978 should be 1978] --> checking read_smalldump >>> reading small dump file: test.dump on unit 9 ST:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 00:46:35.2 WARNING! read_smalldump: *** VELOCITY WILL BE MISSING FROM SMALL DUMP FILES *** forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000047C3CF dump_utils_mp_rea 1225 utils_dumpfiles.f90 phantom 000000000073BE48 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000968A2E testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000481BD3 dump_utils_mp_rea 1246 utils_dumpfiles.f90 phantom 000000000073BE48 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000968A2E testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000048393E dump_utils_mp_rea 1253 utils_dumpfiles.f90 phantom 000000000073BE48 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000968A2E testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004856AD dump_utils_mp_rea 1260 utils_dumpfiles.f90 phantom 000000000073BE48 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000968A2E testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000489CC7 dump_utils_mp_rea 1281 utils_dumpfiles.f90 phantom 000000000073BE48 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 0000000000968A2E testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C3C62 test_mp_testsuite 251 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.29999995231628 comment: hfact was 1.12300002574921 now 1.12300000000000 comment: tolh was 9.999999747378752E-005 now 1.000000000000000E-004 comment: C_cour was 0.300000011920929 now 0.300000000000000 External field found, Bext = 124.000000000000 125.000000000000 126.000000000000 ERROR: could not find mdust_in in header ERROR reading values to verify conservation laws. Resetting initial values. time = 20.0000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(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 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: 585.573 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.200000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.400000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.600000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.800000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.00000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.20000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.40000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.60000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.80000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 2.00000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] <-- STEP TEST COMPLETE --> TESTING IND TIMESTEP UTILS checking bin for dt = 1.500E+00.......OK [got 0 should be 0] checking bin for dt = 7.500E-01.......OK [got 1 should be 1] checking bin for dt = 3.750E-01.......OK [got 2 should be 2] checking bin for dt = 1.875E-01.......OK [got 3 should be 3] checking bin for dt = 9.375E-02.......OK [got 4 should be 4] checking bin for dt = 4.688E-02.......OK [got 5 should be 5] checking bin for dt = 2.344E-02.......OK [got 6 should be 6] checking bin for dt = 1.172E-02.......OK [got 7 should be 7] checking bin for dt = 5.859E-03.......OK [got 8 should be 8] checking bin for dt = 2.930E-03.......OK [got 9 should be 9] checking bin for dt = 1.465E-03.......OK [got 10 should be 10] checking bin for dt = 7.324E-04.......OK [got 11 should be 11] checking bin for dt = 3.662E-04.......OK [got 12 should be 12] checking bin for dt = 1.831E-04.......OK [got 13 should be 13] checking bin for dt = 9.155E-05.......OK [got 14 should be 14] checking bin for dt = 4.578E-05.......OK [got 15 should be 15] checking bin for dt = 2.289E-05.......OK [got 16 should be 16] checking bin for dt = 1.144E-05.......OK [got 17 should be 17] checking bin for dt = 5.722E-06.......OK [got 18 should be 18] checking bin for dt = 2.861E-06.......OK [got 19 should be 19] checking bin for dt = 1.431E-06.......OK [got 20 should be 20] checking bin for dt = 7.153E-07.......OK [got 21 should be 21] checking bin for dt = 3.576E-07.......OK [got 22 should be 22] checking bin for dt = 1.788E-07.......OK [got 23 should be 23] checking bin for dt = 8.941E-08.......OK [got 24 should be 24] checking bin for dt = 4.470E-08.......OK [got 25 should be 25] checking bin for dt = 2.235E-08.......OK [got 26 should be 26] checking bin for dt = 1.118E-08.......OK [got 27 should be 27] checking bin for dt = 5.588E-09.......OK [got 28 should be 28] checking bin for dt = 2.794E-09.......OK [got 29 should be 29] checking bin for dt = huge.............OK [got 0 should be 0] checking bin for dt = 0................OK [got 30 should be 30] ---> checking change_nbinmax routine checking nbinmax = 1 checking nbinmax = 2 checking nbinmax = 3 checking nbinmax = 4 checking nbinmax = 5 checking nbinmax = 6 checking nbinmax = 7 checking nbinmax = 8 checking nbinmax = 9 checking nbinmax = 10 checking nbinmax = 11 checking nbinmax = 12 checking nbinmax = 13 checking nbinmax = 14 checking nbinmax = 15 checking nbinmax = 16 checking nbinmax = 17 checking nbinmax = 18 checking nbinmax = 19 checking nbinmax = 20 checking nbinmax = 21 checking nbinmax = 22 checking nbinmax = 23 checking nbinmax = 24 checking increasing nbins keeps step fraction....OK [max err = 0.000E+00, tol = 2.225-308] checking decreasing nbins keeps step fraction....OK [max err = 0.000E+00, tol = 2.225-308] <-- IND TIMESTEP UTILS TEST COMPLETE --> TESTING EXTERNAL FORCES MODULE ----------------------------------------------------------------------------------------- random seed = -43587 1000 particles set in uniform random distribution ----------------------------------------------------------------------------------------- checking star external force initialisation......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-11, tol = 1.500E-03] checking corotate external force initialisation..OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 2.237E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 2.812E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking binary external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 1.012E-10, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 5.524E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 6.598E-11, tol = 1.500E-03] checking prdrag external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-11, tol = 1.500E-03] checking torus external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 3.904E-07, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 4.945E-07, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.459E-10, tol = 1.500E-03] checking toystar external force initialisation...OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking external B field external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] No spiral arm spheroids set. No potential to be read in. checking spiral external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.271E-05, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 7.191E-05, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 7.013E-05, tol = 1.500E-03] checking Lense-Thirring external force initialisation......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-11, tol = 1.500E-03] ERROR! extern_neutronstar: Error opening ns-rdensity.tab checking neutronstar external force initialisation.........OK [got 1 should be 1] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking Einstein-prec external force initialisation.......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 9.697E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 8.048E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 1.058E-10, tol = 1.500E-03] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.838E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 6.038E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 3.769E-11, tol = 1.500E-03] checking static sinusoid external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.214E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] ERROR! externalforces: Require number of particles per star : iexternalforce = 14 checking grav. wave inspiral external force initialisation...........OK [got 2 should be 2] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking disc gravity external force initialisation........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.109E-10, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking corotating binary external force initialisation...OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 3.458E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 3.145E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 2.126E-06, tol = 1.500E-03] --> testing accrete_particles routine checking accreted=was_accreted for all externf...OK [checked 16 values] --> testing velocity-dependent external force solvers checking corotate external force initialisation..OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking prdrag external force initialisation....OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking Lense-Thirring external force initialisation......OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.417E-16, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.184E-16, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.766E-16, tol = 1.000E-14] checking Einstein-prec external force initialisation.......OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.417E-16, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.184E-16, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.766E-16, tol = 1.000E-14] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking corotating binary external force initialisation...OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] <-- EXTERNAL FORCE TESTS COMPLETE --> TESTING PTMASS MODULE --> testing integration of binary orbit dt for sinks = 1.250000000000000E-002 checking potential energy..............OK [max err = 0.000E+00, tol = 2.220E-16] checking angular momentum..............OK [max err = 0.000E+00, tol = 2.220E-10] nsteps per orbit = 711 norbits = 100 checking angular momentum..............OK [max err = 6.594E-15, tol = 3.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] checking total energy..................OK [max err = 2.440E-08, tol = 3.000E-08] --> 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) forrtl: warning (406): fort: (1): In call to GETNEIGH_POS, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000780AF0 ptmass_mp_ptmass_ 964 ptmass.F90 phantom 00000000009A4A37 testptmass_mp_tes 594 test_ptmass.F90 phantom 00000000009C3D2E test_mp_testsuite 283 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.500E-02 -2.500E-02 -2.500E-02 0.000E+00 ) by accreting 93 particles: M= 3.321E-01 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (sin) >>>>>> s t r e t c h m a p p i n g <<<<<< stretching to match density profile in r direction density at r = 0.000000000000000E+000 is 1.00000000000000 total mass = 6.960406353085308E-004 >>>>>> done forrtl: warning (406): fort: (1): In call to GETNEIGH_POS, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000780AF0 ptmass_mp_ptmass_ 964 ptmass.F90 phantom 00000000009A4A37 testptmass_mp_tes 594 test_ptmass.F90 phantom 00000000009C3D2E test_mp_testsuite 283 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.765E-18 6.370E-19 1.870E-18 0.000E+00 ) by accreting 280 particles: M= 1.000E+00 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] <-- PTMASS TEST COMPLETE --> TESTING GNEWTON MODULE --> testing relativistic orbit ---------- orbital parameters ---------- apocenter : 180.000 pericenter : 20.000 semi-major axis : 100.000 eccentricity : 0.800 initial velocity : 36.724 gravitational radius : 2.122 central mass : 1.000E+06 ---------------------------------------- integrating the orbit for one azimuthal period... ( C_force = 1.000000000000000E-004 ) checking energy........................OK [max err = 9.701E-06, tol = 1.000E-03] checking angular momentum..............OK [max err = 5.863E-07, tol = 1.000E-03] checking precession angle..............OK [max err = 5.762E-07, tol = 1.000E-05] <-- GNEWTON TEST COMPLETE --> TESTING COROTATION MODULE --> testing centrifugal force with rotation about z axis checking fx=Omega^2 x..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=Omega^2 y..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] checking phi=1/2 Omega^2 R^2...........OK [max err = 2.321E-16, tol = 4.441E-16] --> testing coriolis force with rotation about z axis checking fx=2 Omega vy.................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=-2 Omega vx................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] <-- COROTATE TEST COMPLETE --> 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 = 4.441E-16, L2 err = 1.110E-16 tol = 1.000E-15] --> checking transforms for coords and vectors are reversible at r=0 checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 1.225E-16, L2 err = 1.500E-32 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 4.441E-16, L2 err = 9.065E-17 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 4.441E-16, L2 err = 7.166E-17 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 1.000E-15, L2 err = 1.000E-30 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] --> testing get_coord_limits matches coord_transform checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] <-- GEOMETRY TEST COMPLETE --> testing Sedov blast wave ---------------- particles set on 16 x 16 x 16 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 6.250E-02 dy: 6.250E-02 dz: 6.250E-02 ----------------------------------------------------------------------------------------- npart = 4096 particle mass = 2.441406250000000E-004 > step 1 / 128 t = 0.7812500E-03 dt = 7.81E-04 moved 4096 in 1.1 cpu-s < | np = 4096 | > step 3 / 256 t = 0.1171875E-02 dt = 3.91E-04 moved 56 in 0.26 cpu-s < > step 4 / 256 t = 0.1562500E-02 dt = 3.91E-04 moved 4096 in 1.2 cpu-s < > step 5 / 256 t = 0.1953125E-02 dt = 3.91E-04 moved 56 in 0.26 cpu-s < > step 6 / 256 t = 0.2343750E-02 dt = 3.91E-04 moved 352 in 0.34 cpu-s < > step 7 / 256 t = 0.2734375E-02 dt = 3.91E-04 moved 56 in 0.26 cpu-s < > step 8 / 256 t = 0.3125000E-02 dt = 3.91E-04 moved 4096 in 1.2 cpu-s < > step 9 / 256 t = 0.3515625E-02 dt = 3.91E-04 moved 56 in 0.26 cpu-s < > step 10 / 256 t = 0.3906250E-02 dt = 3.91E-04 moved 400 in 0.35 cpu-s < > step 11 / 256 t = 0.4296875E-02 dt = 3.91E-04 moved 32 in 0.25 cpu-s < > step 12 / 256 t = 0.4687500E-02 dt = 3.91E-04 moved 1232 in 0.61 cpu-s < > step 13 / 256 t = 0.5078125E-02 dt = 3.91E-04 moved 32 in 0.25 cpu-s < > step 14 / 256 t = 0.5468750E-02 dt = 3.91E-04 moved 408 in 0.36 cpu-s < > step 15 / 256 t = 0.5859375E-02 dt = 3.91E-04 moved 8 in 0.24 cpu-s < > step 16 / 256 t = 0.6250000E-02 dt = 3.91E-04 moved 4096 in 1.2 cpu-s < > step 17 / 256 t = 0.6640625E-02 dt = 3.91E-04 moved 8 in 0.24 cpu-s < > step 18 / 256 t = 0.7031250E-02 dt = 3.91E-04 moved 408 in 0.36 cpu-s < > step 19 / 256 t = 0.7421875E-02 dt = 3.91E-04 moved 8 in 0.24 cpu-s < > step 20 / 256 t = 0.7812500E-02 dt = 3.91E-04 moved 1232 in 0.61 cpu-s < > step 21 / 256 t = 0.8203125E-02 dt = 3.91E-04 moved 8 in 0.24 cpu-s < > step 22 / 256 t = 0.8593750E-02 dt = 3.91E-04 moved 456 in 0.38 cpu-s < > step 12 / 128 t = 0.9375000E-02 dt = 7.81E-04 moved 2728 in 0.98 cpu-s < > step 13 / 128 t = 0.1015625E-01 dt = 7.81E-04 moved 480 in 0.38 cpu-s < > step 14 / 128 t = 0.1093750E-01 dt = 7.81E-04 moved 1520 in 0.70 cpu-s < > step 15 / 128 t = 0.1171875E-01 dt = 7.81E-04 moved 456 in 0.38 cpu-s < > step 16 / 128 t = 0.1250000E-01 dt = 7.81E-04 moved 4096 in 1.3 cpu-s < > step 17 / 128 t = 0.1328125E-01 dt = 7.81E-04 moved 456 in 0.38 cpu-s < > step 18 / 128 t = 0.1406250E-01 dt = 7.81E-04 moved 1520 in 0.70 cpu-s < > step 19 / 128 t = 0.1484375E-01 dt = 7.81E-04 moved 456 in 0.38 cpu-s < > step 20 / 128 t = 0.1562500E-01 dt = 7.81E-04 moved 2944 in 1.1 cpu-s < > step 21 / 128 t = 0.1640625E-01 dt = 7.81E-04 moved 408 in 0.37 cpu-s < > step 22 / 128 t = 0.1718750E-01 dt = 7.81E-04 moved 1568 in 0.73 cpu-s < > step 23 / 128 t = 0.1796875E-01 dt = 7.81E-04 moved 376 in 0.37 cpu-s < > step 24 / 128 t = 0.1875000E-01 dt = 7.81E-04 moved 3816 in 1.3 cpu-s < > step 25 / 128 t = 0.1953125E-01 dt = 7.81E-04 moved 328 in 0.35 cpu-s < > step 26 / 128 t = 0.2031250E-01 dt = 7.81E-04 moved 1616 in 0.74 cpu-s < > step 27 / 128 t = 0.2109375E-01 dt = 7.81E-04 moved 232 in 0.33 cpu-s < > step 28 / 128 t = 0.2187500E-01 dt = 7.81E-04 moved 3088 in 1.1 cpu-s < > step 29 / 128 t = 0.2265625E-01 dt = 7.81E-04 moved 256 in 0.33 cpu-s < > step 30 / 128 t = 0.2343750E-01 dt = 7.81E-04 moved 1784 in 0.79 cpu-s < > step 31 / 128 t = 0.2421875E-01 dt = 7.81E-04 moved 152 in 0.30 cpu-s < > step 32 / 128 t = 0.2500000E-01 dt = 7.81E-04 moved 4096 in 1.4 cpu-s < > step 33 / 128 t = 0.2578125E-01 dt = 7.81E-04 moved 80 in 0.28 cpu-s < > step 34 / 128 t = 0.2656250E-01 dt = 7.81E-04 moved 1760 in 0.79 cpu-s < > step 35 / 128 t = 0.2734375E-01 dt = 7.81E-04 moved 56 in 0.27 cpu-s < > step 36 / 128 t = 0.2812500E-01 dt = 7.81E-04 moved 3304 in 1.2 cpu-s < > step 37 / 128 t = 0.2890625E-01 dt = 7.81E-04 moved 32 in 0.26 cpu-s < > step 38 / 128 t = 0.2968750E-01 dt = 7.81E-04 moved 1328 in 0.65 cpu-s < > step 39 / 128 t = 0.3046875E-01 dt = 7.81E-04 moved 32 in 0.26 cpu-s < > step 40 / 128 t = 0.3125000E-01 dt = 7.81E-04 moved 3936 in 1.4 cpu-s < > step 41 / 128 t = 0.3203125E-01 dt = 7.81E-04 moved 8 in 0.25 cpu-s < > step 42 / 128 t = 0.3281250E-01 dt = 7.81E-04 moved 1424 in 0.69 cpu-s < > step 22 / 64 t = 0.3437500E-01 dt = 0.0016 moved 2944 in 1.2 cpu-s < > step 23 / 64 t = 0.3593750E-01 dt = 0.0016 moved 1472 in 0.71 cpu-s < > step 24 / 64 t = 0.3750000E-01 dt = 0.0016 moved 4088 in 1.4 cpu-s < > step 25 / 64 t = 0.3906250E-01 dt = 0.0016 moved 1568 in 0.73 cpu-s < > step 26 / 64 t = 0.4062500E-01 dt = 0.0016 moved 3016 in 1.2 cpu-s < > step 27 / 64 t = 0.4218750E-01 dt = 0.0016 moved 1616 in 0.75 cpu-s < > step 28 / 64 t = 0.4375000E-01 dt = 0.0016 moved 3816 in 1.4 cpu-s < > step 29 / 64 t = 0.4531250E-01 dt = 0.0016 moved 1640 in 0.76 cpu-s < > step 30 / 64 t = 0.4687500E-01 dt = 0.0016 moved 3232 in 1.2 cpu-s < > step 31 / 64 t = 0.4843750E-01 dt = 0.0016 moved 1688 in 0.77 cpu-s < > step 32 / 64 t = 0.5000000E-01 dt = 0.0016 moved 4096 in 1.5 cpu-s < > step 33 / 64 t = 0.5156250E-01 dt = 0.0016 moved 1808 in 0.80 cpu-s < > step 34 / 64 t = 0.5312500E-01 dt = 0.0016 moved 3352 in 1.3 cpu-s < > step 35 / 64 t = 0.5468750E-01 dt = 0.0016 moved 1640 in 0.76 cpu-s < > step 36 / 64 t = 0.5625000E-01 dt = 0.0016 moved 3936 in 1.4 cpu-s < > step 37 / 64 t = 0.5781250E-01 dt = 0.0016 moved 1776 in 0.80 cpu-s < > step 38 / 64 t = 0.5937500E-01 dt = 0.0016 moved 3400 in 1.3 cpu-s < > step 39 / 64 t = 0.6093750E-01 dt = 0.0016 moved 1728 in 0.78 cpu-s < > step 40 / 64 t = 0.6250000E-01 dt = 0.0016 moved 4088 in 1.5 cpu-s < > step 41 / 64 t = 0.6406250E-01 dt = 0.0016 moved 1680 in 0.77 cpu-s < > step 42 / 64 t = 0.6562500E-01 dt = 0.0016 moved 3600 in 1.4 cpu-s < > step 43 / 64 t = 0.6718750E-01 dt = 0.0016 moved 1832 in 0.83 cpu-s < > step 44 / 64 t = 0.6875000E-01 dt = 0.0016 moved 4016 in 1.5 cpu-s < > step 45 / 64 t = 0.7031250E-01 dt = 0.0016 moved 1832 in 0.83 cpu-s < > step 46 / 64 t = 0.7187500E-01 dt = 0.0016 moved 3600 in 1.4 cpu-s < > step 47 / 64 t = 0.7343750E-01 dt = 0.0016 moved 1784 in 0.82 cpu-s < > step 48 / 64 t = 0.7500000E-01 dt = 0.0016 moved 4096 in 1.5 cpu-s < > step 49 / 64 t = 0.7656250E-01 dt = 0.0016 moved 1544 in 0.75 cpu-s < > step 50 / 64 t = 0.7812500E-01 dt = 0.0016 moved 3744 in 1.4 cpu-s < > step 51 / 64 t = 0.7968750E-01 dt = 0.0016 moved 1232 in 0.65 cpu-s < > step 52 / 64 t = 0.8125000E-01 dt = 0.0016 moved 4064 in 1.5 cpu-s < > step 53 / 64 t = 0.8281250E-01 dt = 0.0016 moved 1096 in 0.61 cpu-s < > step 54 / 64 t = 0.8437500E-01 dt = 0.0016 moved 3744 in 1.4 cpu-s < > step 55 / 64 t = 0.8593750E-01 dt = 0.0016 moved 968 in 0.56 cpu-s < > step 56 / 64 t = 0.8750000E-01 dt = 0.0016 moved 4096 in 1.5 cpu-s < > step 57 / 64 t = 0.8906250E-01 dt = 0.0016 moved 992 in 0.58 cpu-s < > step 58 / 64 t = 0.9062500E-01 dt = 0.0016 moved 3696 in 1.4 cpu-s < > step 59 / 64 t = 0.9218750E-01 dt = 0.0016 moved 728 in 0.49 cpu-s < > step 60 / 64 t = 0.9375000E-01 dt = 0.0016 moved 4064 in 1.5 cpu-s < > step 61 / 64 t = 0.9531250E-01 dt = 0.0016 moved 704 in 0.48 cpu-s < > step 62 / 64 t = 0.9687500E-01 dt = 0.0016 moved 3600 in 1.4 cpu-s < > step 63 / 64 t = 0.9843750E-01 dt = 0.0016 moved 632 in 0.47 cpu-s < > step 64 / 64 t = 0.1000000 dt = 0.0016 moved 4096 in 1.5 cpu-s < --------> TIME = 0.1000 : full dump written to file test002 <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004947BA dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000087DBC2 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C1A32 testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C3E2D test_mp_testsuite 323 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000497F17 dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000087DBC2 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C1A32 testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C3E2D test_mp_testsuite 323 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 0000000000499191 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000087DBC2 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C1A32 testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C3E2D test_mp_testsuite 323 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000049A40F dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000087DBC2 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C1A32 testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C3E2D test_mp_testsuite 323 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 000000000049CEF0 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 000000000072B4DC readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000087DBC2 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C1A32 testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C3E2D test_mp_testsuite 323 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002B0601CB35A3 for_emit_diagnost Unknown Unknown phantom 00000000004A1FFB dump_utils_mp_wri 1681 utils_dumpfiles.f90 phantom 000000000072EC9A readwrite_dumps_m 539 readwrite_dumps.F90 phantom 000000000087DBC2 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C1A32 testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C3E2D test_mp_testsuite 323 testsuite.F90 phantom 00000000009E18CF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002B0604C31C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown input file test.in written successfully. Since code start: 95 timesteps, wall: 4.9s cpu: 77s cpu/wall: 16 Since last dump : 95 timesteps, wall: 4.9s cpu: 77s cpu/wall: 16 wall cpu cpu/wall frac step : 4.88s 76.54s 15.70 100.00% step (force): 1.62s 23.50s 14.46 33.33% step (dens) : 2.00s 31.12s 15.56 41.03% step (link) : 1.12s 20.90s 18.58 23.08% ------------------------------------------------------------------------------ |** Number of steps since last summary: 95 **| |** Wall time since last summary: 4.875 seconds **| ------------------------------------------------------------------------------ |* particles woken *| | #steps | mean # part/step | max # part/step | | 24| 97.00| 216 | ------------------------------------------------------------------------------ IND TIMESTEPS efficiency = 111.35% wall time per particle (last full step) : 3.05E-05s wall time per particle (ave. all steps) : 2.74E-05s ------------------------------------------------- | bin | dt | npart | frac | cpufrac | ------------------------------------------------- | 0 1.000E-01 0 0.00% 2.33% | | 1 5.000E-02 0 0.00% 2.23% | | 2 2.500E-02 0 0.00% 4.40% | | 3 1.250E-02 80 1.95% 8.72% | | 4 6.250E-03 488 11.91% 17.08% | | 5 3.125E-03 2920 71.29% 30.62% | | 6 1.563E-03 608 14.84% 34.63% | ------------------------------------------------- Etot= 9.668E-01, Ekin= 1.982E-01, Etherm= 7.686E-01, Epot= 0.000E+00 Emag= 0.000E+00 Linm= 1.647E-17, Angm= 8.260E-19 Centre of Mass = 1.084E-18, 2.467E-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 checking total energy..................OK [max err = 2.831E-04, tol = 4.700E-04] checking linear momentum...............OK [max err = 1.647E-17, tol = 7.000E-15] <-- testing complete total wall time = 11 min, 27.25 s (= 6.8725E+02s) total cpu time = 1 hr, 5 min, 24.88 s (= 3.9249E+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