VELOCIraptor with cooling gets stuck in infinite loop
When I run SWIFT with VELOCIraptor it gets in an infinite loop of recomputing VELOCIraptor in the same snap shot over and over again. SWIFT was configured as follows:
./configure --disable-mpi --disable-parallel-hdf5 --with-velociraptor=/path/to/velociraptor/VELOCIraptor-STF/stf/lib --with-cooling=const-lambda
And the Small Cosmo Volume example was run using constant Lambda cooling and VELOCIraptor:
../swift -c -s -x -C -G -t 8 small_cosmo_volume.yml 2>&1 | tee output.log
In which the small_cosmo_volume.yml file was adjusted by adding two sections for the cooling and VELOCIraptor:
StructureFinding:
config_file_name: stf_input.cfg # Name of the STF config file.
basename: ./stf # Common part of the name of output files.
output_time_format: 1 # Specifies the frequency format of structure finding. 0 for simulation steps (delta_step) and 1 for simulation time intervals (delta_time).
scale_factor_first: 0.02 # Scale-factor of the first snaphot (cosmological run)
delta_time: 1.10 # Time difference between consecutive structure finding outputs (in internal units) in simulation time intervals.
# Constant lambda cooling function
LambdaCooling:
lambda_nH2_cgs: 1e-24 # Cooling rate divided by square Hydrogen number density (in cgs units [erg * s^-1 * cm^3])
cooling_tstep_mult: 1.0 # (Optional) Dimensionless pre-factor for the time-step condition.See the attached file for the complete yaml file.
If however the delta_time parameter is changed to a value equal to the one used for the snap shots (delta_time=1.02), the code runs fine and does not get stuck in an infinite loop.