Floating point exception when running 'triangle' set up
If I run the following in examples/self_gravity_testing/triangle...
cosma-b > python makeIC.py 3.0
UnitMass_in_cgs: 1.98855e+33
UnitLength_in_cgs: 1.495978707e+13
UnitVelocity_in_cgs: 474372.0
UnitTime_in_cgs: 31535982.4568
G= 39.4120520353
[[ 6.8 6.5 6.5 ]
[ 6.35 6.75980762 6.5 ]
[ 6.35 6.24019238 6.5 ]
[ 9.5 6.5 6.5 ]]
cosma-b > ../../swift -G -e -t 1 triangle_theta_0p7.yml
Welcome to the cosmological hydrodynamical code
______ _________________
/ ___/ | / / _/ ___/_ __/
\__ \| | /| / // // /_ / /
___/ /| |/ |/ // // __/ / /
/____/ |__/|__/___/_/ /_/
SPH With Inter-dependent Fine-grained Tasking
Version : 0.6.0
Revision: v0.6.0-214-g05313121-dirty, Branch: gravity_testing, Date: 2017-08-24 16:33:20 +0100
Webpage : www.swiftsim.com
Config. options: '--disable-vec --enable-no-gravity-below-id=4 --with-multipole-order=2'
Compiler: ICC, Version: 16.0.20160415
CFLAGS : '-idirafter /usr/include/linux -O3 -ansi_alias -xAVX -no-vec -no-simd -w2 -Wunused-variable -Werror'
HDF5 library version: 1.8.14
FFTW library version: 3.x (details not available)
[00000.0] main: CPU frequency used for tick conversion: 2399966719 Hz
[00000.0] main: Running on: cosma-b
[00000.0] main: WARNING: Floating point exceptions will be reported.
[00000.0] main: sizeof(part) is 128 bytes.
[00000.0] main: sizeof(xpart) is 64 bytes.
[00000.0] main: sizeof(spart) is 64 bytes.
[00000.0] main: sizeof(gpart) is 96 bytes.
[00000.0] main: sizeof(multipole) is 52 bytes.
[00000.0] main: sizeof(grav_tensor) is 80 bytes.
[00000.0] main: sizeof(task) is 64 bytes.
[00000.0] main: sizeof(cell) is 640 bytes.
[00000.0] main: Reading runtime parameters from file 'triangle_theta_0p7.yml'
[00000.0] main: Reading ICs from file 'triangle.hdf5'
[00000.0] read_ic_single: IC and internal units match. No conversion needed.
[00000.0] main: Reading initial conditions took 3.144 ms.
[00000.0] main: Read 0 gas particles, 0 star particles and 4 gparts from the ICs.
[00000.0] main: space_init took 1.585 ms.
[00000.0] main: space dimensions are [ 12.000 12.000 12.000 ].
[00000.0] main: space is periodic.
[00000.0] main: highest-level cell dimensions are [ 12 12 12 ].
[00000.0] main: 0 parts in 1728 cells.
[00000.0] main: 4 gparts in 1728 cells.
[00000.0] main: 0 sparts in 1728 cells.
[00000.0] main: maximum depth is 0.
[00000.0] engine_init: no processor affinity used
[00000.0] engine_policy: engine policies are [ steal keep numa_affinity self_gravity ]
[00000.0] gravity_props_print: Self-gravity scheme: FMM-MM with m-poles of order 2
[00000.0] gravity_props_print: Self-gravity time integration: eta=0.0100
[00000.0] gravity_props_print: Self-gravity opening angle: theta=0.7000
[00000.0] gravity_props_print: Self-gravity softening: epsilon=0.0300 (Plummer equivalent: 0.0100)
[00000.0] gravity_props_print: Self-gravity mesh smoothing-scale: a_smooth=1.250000
[00000.0] gravity_props_print: Self-gravity tree cut-off: r_cut_max=4.500000
[00000.0] gravity_props_print: Self-gravity truncation cut-off: r_cut_min=0.100000
[00000.0] engine_init: Absolute minimal timestep size: 6.938894e-16
[00000.0] engine_init: Minimal timestep size (on time-line): 6.938894e-16
[00000.0] engine_init: Maximal timestep size (on time-line): 9.765625e-02
[00000.0] main: engine_init took 2.682 ms.
[00000.0] main: Running on 0 gas particles, 0 star particles and 4 DM particles (4 gravity particles)
[00000.0] main: from t=0.000e+00 until t=1.000e+02 with 1 threads and 1 queues (dt_min=1.000e-15, dt_max=1.000e-01)...
[00000.0] engine_init_particles: Computing initial gas densities.
[00000.0] engine_init_particles: Converting internal energy variable.
[00000.0] engine_init_particles: Running initial fake time-step.
Floating exception (core dumped)
On the other hand, if we run the 'particle line' with two particles (examples/self_gravity_testing/particle_line, remove the -d) there is no issue
cosma-b > python makeIC.py 2 2.0
UnitMass_in_cgs: 1.98855e+33
UnitLength_in_cgs: 1.495978707e+13
UnitVelocity_in_cgs: 474372.0
UnitTime_in_cgs: 31535982.4568
G= 39.4120520353
Coordinates of particles are:
[[ 6.25 6. 6. ]
[ 6.75 6. 6. ]
[ 8.5 6. 6. ]]
Masses of particles are:
[ 5.00000000e-01 5.00000000e-01 1.00000000e-06]
Velocities of particles are:
[[ 0. 0. 0. ]
[ 0. 0. 0. ]
[ 0. 4.43914699 0. ]]
cosma-b > ../../swift -G -e -d -t 1 particle_line_theta_0p7.yml
Welcome to the cosmological hydrodynamical code
______ _________________
/ ___/ | / / _/ ___/_ __/
\__ \| | /| / // // /_ / /
___/ /| |/ |/ // // __/ / /
/____/ |__/|__/___/_/ /_/
SPH With Inter-dependent Fine-grained Tasking
Version : 0.6.0
Revision: v0.6.0-215-gc44f5af7-dirty, Branch: gravity_testing, Date: 2017-08-24 19:26:46 +0100
Webpage : www.swiftsim.com
Config. options: '--disable-vec --enable-no-gravity-below-id=3 --with-multipole-order=2'
Compiler: ICC, Version: 16.0.20160415
CFLAGS : '-idirafter /usr/include/linux -O3 -ansi_alias -xAVX -no-vec -no-simd -w2 -Wunused-variable -Werror'
HDF5 library version: 1.8.14
FFTW library version: 3.x (details not available)
[00000.0] main: Executing a dry run. No i/o or time integration will be performed.
[00000.0] main: CPU frequency used for tick conversion: 2399966187 Hz
[00000.0] main: Running on: cosma-b
[00000.0] main: WARNING: Floating point exceptions will be reported.
[00000.0] main: sizeof(part) is 128 bytes.
[00000.0] main: sizeof(xpart) is 64 bytes.
[00000.0] main: sizeof(spart) is 64 bytes.
[00000.0] main: sizeof(gpart) is 96 bytes.
[00000.0] main: sizeof(multipole) is 52 bytes.
[00000.0] main: sizeof(grav_tensor) is 80 bytes.
[00000.0] main: sizeof(task) is 64 bytes.
[00000.0] main: sizeof(cell) is 640 bytes.
[00000.0] main: Reading runtime parameters from file 'particle_line_theta_0p7.yml'
[00000.0] main: Reading ICs from file 'particle_line.hdf5'
[00000.0] read_ic_single: IC and internal units match. No conversion needed.
[00000.0] main: Reading initial conditions took 2.615 ms.
[00000.0] main: Read 0 gas particles, 0 star particles and 3 gparts from the ICs.
[00000.0] main: space_init took 0.017 ms.
[00000.0] main: space dimensions are [ 12.000 12.000 12.000 ].
[00000.0] main: space isn't periodic.
[00000.0] main: highest-level cell dimensions are [ 0 0 0 ].
[00000.0] main: 0 parts in 0 cells.
[00000.0] main: 3 gparts in 0 cells.
[00000.0] main: 0 sparts in 0 cells.
[00000.0] main: maximum depth is 0.
[00000.0] engine_init: no processor affinity used
[00000.0] engine_policy: engine policies are [ steal keep numa_affinity self_gravity ]
[00000.0] gravity_props_print: Self-gravity scheme: FMM-MM with m-poles of order 2
[00000.0] gravity_props_print: Self-gravity time integration: eta=0.0100
[00000.0] gravity_props_print: Self-gravity opening angle: theta=1.0000
[00000.0] gravity_props_print: Self-gravity softening: epsilon=0.0300 (Plummer equivalent: 0.0100)
[00000.0] gravity_props_print: Self-gravity mesh smoothing-scale: a_smooth=1.250000
[00000.0] gravity_props_print: Self-gravity tree cut-off: r_cut_max=4.500000
[00000.0] gravity_props_print: Self-gravity truncation cut-off: r_cut_min=0.100000
[00000.0] engine_init: Absolute minimal timestep size: 6.938894e-16
[00000.0] engine_init: Minimal timestep size (on time-line): 8.881784e-14
[00000.0] engine_init: Maximal timestep size (on time-line): 7.812500e-01
[00000.0] main: engine_init took 1.205 ms.
[00000.0] main: Running on 0 gas particles, 0 star particles and 3 DM particles (3 gravity particles)
[00000.0] main: from t=0.000e+00 until t=1.000e+02 with 1 threads and 1 queues (dt_min=1.000e-13, dt_max=1.000e+00)...
[00000.0] main: Time integration ready to start. End of dry-run.