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,0) 0.000 B real(8) Bxyz(3,0) 0.000 B real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(2,0) 0.000 B real(8) dustevol(1,0) 0.000 B real(8) ddustevol(1,0) 0.000 B real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,1,0) 0.000 B real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(1000000) 3.815 MB real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1000000) 30.518 MB real(8) dBevol(4,0) 0.000 B real(4) divBsumm(0) 0.000 B real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(1,0) 0.000 B real(8) Bpred(4,0) 0.000 B real(8) dustproppred(4,0) 0.000 B integer(1) iphase(1000000) 976.562 KB integer(1) iphase_soa(1000000) 976.562 KB real(4) gradh(2,1000000) 7.629 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) 102.997 MB kdnode node(1000001) 102.997 MB integer(4) nodemap(1000001) 3.815 MB --------------------------------------------------------- Total memory allocated to arrays: 575.083 MB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) --> TESTING SELF-GRAVITY --> testing taylor series expansion about current node checking fx taylor series about f0.....OK [max err = 2.939E-04, tol = 3.000E-04] checking fy taylor series about f0.....OK [max err = 1.069E-04, tol = 1.100E-04] checking fz taylor series about f0.....OK [max err = 8.815E-05, tol = 9.000E-05] checking phi taylor series about f0....OK [max err = 7.685E-04, tol = 8.000E-04] --> testing taylor series expansion about distant node checking fx taylor series about f0.....OK [max err = 8.625E-05, tol = 8.700E-05] checking fy taylor series about f0.....OK [max err = 1.418E-06, tol = 1.500E-06] checking fz taylor series about f0.....OK [max err = 1.512E-05, tol = 1.600E-05] checking phi taylor series about f0....OK [max err = 5.827E-06, tol = 5.900E-06] --> checking results of compute_fnode routine checking fx=-dphi/dx...................OK [max err = 2.279E-10, tol = 1.000E-06] checking fy=-dphi/dy...................OK [max err = 2.797E-10, tol = 1.000E-06] checking fz=-dphi/dz...................OK [max err = 1.714E-10, tol = 1.000E-06] checking dfx/dx........................OK [max err = 2.003E-09, tol = 1.000E-06] checking dfx/dy........................OK [max err = 9.794E-11, tol = 1.000E-06] checking dfx/dz........................OK [max err = 1.771E-10, tol = 1.000E-06] checking dfy/dy........................OK [max err = 6.952E-09, tol = 1.000E-06] checking dfx/dz........................OK [max err = 1.768E-10, tol = 1.000E-06] checking dfz/dz........................OK [max err = 1.450E-08, tol = 1.000E-06] checking d^2fx/dx^2....................OK [max err = 1.669E-04, tol = 1.000E-03] checking d^2fx/dy^2....................OK [max err = 1.035E-03, tol = 1.100E-03] checking d^2fx/dz^2....................OK [max err = 4.838E-04, tol = 1.000E-03] checking d^2fy/dx^2....................OK [max err = 4.748E-05, tol = 1.000E-03] checking d^2fy/dy^2....................OK [max err = 6.492E-04, tol = 1.000E-03] checking d^2fy/dz^2....................OK [max err = 7.659E-04, tol = 1.000E-03] checking d^2fz/dx^2....................OK [max err = 2.254E-05, tol = 1.000E-03] checking d^2fz/dy^2....................OK [max err = 1.185E-03, tol = 1.200E-03] checking d^2fz/dz^2....................OK [max err = 3.533E-04, tol = 1.000E-03] --> testing taylor series expansion about both current and distant nodes checking fx taylor series about f0.....OK [max err = 4.216E-05, tol = 4.300E-05] checking fy taylor series about f0.....OK [max err = 1.358E-04, tol = 1.400E-04] checking fz taylor series about f0.....OK [max err = 3.163E-04, tol = 3.200E-04] checking phi taylor series about f0....OK [max err = 9.619E-04, tol = 9.700E-04] --> testing gravity force in densityforce for gas particles trial neigh mean : 259.35, real neigh mean = 49.90 ratio try/act= 5.20 trial neigh max : 406, max real neigh = 57 n neighbour calls : 331, mean per part = 1.448 n density calcs : 2989, mean per part = 4.045 dtmin = 0.29999999999999999 dtmax = 0.29999999999999999 dtmax/dtmin = 0.99999999999999978 dtcour/dtf = 9.1416825968195159 (wall: total = 0.00s force = 0.00s dens = 0.00s link = 0.00s) (cpu : total = 0.05s force = 0.01s dens = 0.02s link = 0.02s) (cpu/wall: total = --- force = --- dens = --- link = ---) completed in 0.00 s checking fgrav(x)......................OK [max err = 3.612E-03, L2 err = 1.066E-03 tol = 5.000E-03] checking fgrav(y)......................OK [max err = 4.294E-03, L2 err = 1.200E-03 tol = 6.000E-03] checking fgrav(z)......................OK [max err = 4.419E-03, L2 err = 1.149E-03 tol = 9.400E-03] checking potential.....................OK [max err = 4.391E-04, tol = 4.800E-04] --> testing gravity force in densityforce for empty particles trial neigh mean : 259.35, real neigh mean = 49.90 ratio try/act= 5.20 trial neigh max : 406, max real neigh = 57 n neighbour calls : 331, mean per part = 1.448 n density calcs : 2989, mean per part = 4.045 dtmin = 0.29999999999999999 dtmax = 0.29999999999999999 dtmax/dtmin = 0.99999999999999978 dtcour/dtf = 9.1416825968195159 (wall: total = 0.00s force = 0.00s dens = 0.00s link = 0.00s) (cpu : total = 0.03s force = 0.01s dens = 0.02s link = 0.00s) (cpu/wall: total = --- force = --- dens = --- link = ---) completed in 0.00 s checking fgrav(x)......................OK [max err = 3.612E-03, L2 err = 1.066E-03 tol = 5.000E-03] checking fgrav(y)......................OK [max err = 4.294E-03, L2 err = 1.200E-03 tol = 6.000E-03] checking fgrav(z)......................OK [max err = 4.419E-03, L2 err = 1.149E-03 tol = 9.400E-03] checking potential.....................OK [max err = 4.391E-04, tol = 4.800E-04] --> testing gravity force in densityforce for star particles trial neigh mean : 259.35, real neigh mean = 49.90 ratio try/act= 5.20 trial neigh max : 406, max real neigh = 57 n neighbour calls : 331, mean per part = 1.448 n density calcs : 2989, mean per part = 4.045 dtmin = 0.29999999999999999 dtmax = 0.29999999999999999 dtmax/dtmin = 0.99999999999999978 dtcour/dtf = 9.1416825968195159 (wall: total = 0.00s force = 0.00s dens = 0.00s link = 0.00s) (cpu : total = 0.04s force = 0.01s dens = 0.02s link = 0.01s) (cpu/wall: total = --- force = --- dens = --- link = ---) completed in 0.00 s checking fgrav(x)......................OK [max err = 3.612E-03, L2 err = 1.066E-03 tol = 5.000E-03] checking fgrav(y)......................OK [max err = 4.294E-03, L2 err = 1.200E-03 tol = 6.000E-03] checking fgrav(z)......................OK [max err = 4.419E-03, L2 err = 1.149E-03 tol = 9.400E-03] checking potential.....................OK [max err = 4.391E-04, tol = 4.800E-04] --> testing gravity force in densityforce for darkm particles trial neigh mean : 259.35, real neigh mean = 49.90 ratio try/act= 5.20 trial neigh max : 406, max real neigh = 57 n neighbour calls : 331, mean per part = 1.448 n density calcs : 2989, mean per part = 4.045 dtmin = 0.29999999999999999 dtmax = 0.29999999999999999 dtmax/dtmin = 0.99999999999999978 dtcour/dtf = 9.1416825968195159 (wall: total = 0.00s force = 0.00s dens = 0.00s link = 0.00s) (cpu : total = 0.03s force = 0.01s dens = 0.02s link = 0.01s) (cpu/wall: total = --- force = --- dens = --- link = ---) completed in 0.00 s checking fgrav(x)......................OK [max err = 3.612E-03, L2 err = 1.066E-03 tol = 5.000E-03] checking fgrav(y)......................OK [max err = 4.294E-03, L2 err = 1.200E-03 tol = 6.000E-03] checking fgrav(z)......................OK [max err = 4.419E-03, L2 err = 1.149E-03 tol = 9.400E-03] checking potential.....................OK [max err = 4.391E-04, tol = 4.800E-04] --> testing gravity force in densityforce for bulge particles trial neigh mean : 259.35, real neigh mean = 49.90 ratio try/act= 5.20 trial neigh max : 406, max real neigh = 57 n neighbour calls : 331, mean per part = 1.448 n density calcs : 2989, mean per part = 4.045 dtmin = 0.29999999999999999 dtmax = 0.29999999999999999 dtmax/dtmin = 0.99999999999999978 dtcour/dtf = 9.1416825968195159 (wall: total = 0.00s force = 0.00s dens = 0.00s link = 0.00s) (cpu : total = 0.03s force = 0.01s dens = 0.02s link = 0.00s) (cpu/wall: total = --- force = --- dens = --- link = ---) completed in 0.00 s checking fgrav(x)......................OK [max err = 3.612E-03, L2 err = 1.066E-03 tol = 5.000E-03] checking fgrav(y)......................OK [max err = 4.294E-03, L2 err = 1.200E-03 tol = 6.000E-03] checking fgrav(z)......................OK [max err = 4.419E-03, L2 err = 1.149E-03 tol = 9.400E-03] checking potential.....................OK [max err = 4.391E-04, tol = 4.800E-04] <-- SELF-GRAVITY TESTS COMPLETE --> TESTING PTMASS MODULE --> testing integration of binary orbit dt for sinks = 1.2500000000000001E-002 checking potential energy..............OK [max err = 0.000E+00, tol = 2.220E-16] checking angular momentum..............OK [max err = 0.000E+00, tol = 2.220E-10] nsteps per orbit = 711 norbits = 100 checking angular momentum..............OK [max err = 6.594E-15, tol = 3.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] checking total energy..................OK [max err = 2.440E-08, tol = 3.000E-08] --> testing integration of circumbinary disc Centre of mass is at (x,y,z) = ( 0.000E+00, -2.411E-17, 2.068E-18) Particle setup OK checking no errors during disc setup...OK [got 0 should be 0] dt for sinks = 1.2500000000000001E-002 nsteps per orbit = 711 norbits = 10 checking angular momentum..............OK [max err = 7.814E-09, tol = 3.000E-07] checking linear momentum...............OK [max err = 1.896E-14, tol = 6.000E-14] checking total energy..................OK [max err = 2.617E-03, tol = 3.100E-03] --> testing integration of disc around eccentric binary Centre of mass is at (x,y,z) = ( 0.000E+00, 1.640E-17, -3.223E-18) Particle setup OK checking no errors during disc setup...OK [got 0 should be 0] dt for sinks = 2.2963966338592295E-002 nsteps per orbit = 387 norbits = 10 checking angular momentum..............OK [max err = 2.952E-07, tol = 1.100E-06] checking linear momentum...............OK [max err = 1.112E-14, tol = 4.000E-14] checking total energy..................OK [max err = 2.446E-03, tol = 1.200E-02] --> testing softening in sink particle binary checking potential energy..............OK [max err = 3.391E-16, tol = 4.441E-16] nsteps per orbit = 1099 norbits = 10 checking angular momentum..............OK [max err = 6.324E-15, tol = 2.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.225-308] checking total energy..................OK [max err = 1.021E-10, tol = 2.000E-09] --> testing accretion onto sink particles checking accretion flag................OK [got T should be T] checking isdead_or_accreted flag.......OK [got T should be T] checking x(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking y(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking z(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking vx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (uniform density) got rhomax = 28.595009469164072 on particle 179 ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.500E-02 -2.500E-02 -2.500E-02 0.000E+00 ) by accreting 93 particles: M= 3.321E-01 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (sin) >>>>>> s t r e t c h m a p p i n g <<<<<< stretching to match density profile in r direction density at r = 0.0000000000000000 is 1.0000000000000000 total mass = 6.9604063530853081E-004 >>>>>> done got rhomax = 1176.9049829475812 on particle 163 checking rhomax........................OK [max err = 0.000E+00, tol = 2.220E-16] 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 complete total wall time = 49.12 s total cpu time = 11 min, 23.37 s (= 6.8337E+02s) SUMMARY OF ALL TESTS: PASSED: 35 of 35 100.0% FAILED: 0 of 35 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found