diff --git a/Template-Config.sh b/Template-Config.sh index 2269516793bce4e873a7f2f715d9f2b854d049ca..d2f1be5102e9cead3b0f79348bed9c4d45e7f1c8 100644 --- a/Template-Config.sh +++ b/Template-Config.sh @@ -1,202 +1,221 @@ -################################################## -# Enable/Disable compile-time options as needed # -################################################## +######################################### +# Pick compile-time options as needed # +# from this template. # +######################################### - -#----------------------------------------Parallelization options - -IMPOSE_PINNING -#IMPOSE_PINNING_OVERRIDE_MODE -#EXPLICIT_VECTORIZATION # This uses AVX at selected places through the vectorclass C++ library -#PRESERVE_SHMEM_BINARY_INVARIANCE -#SIMPLE_DOMAIN_AGGREGATION - #--------------------------------------- Basic operation mode of code -PERIODIC -#TWODIMS -#ONEDIMS -#LONG_X_BITS=2 -#LONG_Y_BITS=2 -#LONG_Z_BITS=1 -#NTYPES=6 # Number of particle types. Note that this may only be changed from the default value of 6 if - # HDF5 snapshot files are used. +#PERIODIC # enables periodic boundary condistions +#NTYPES=6 # number of particle types +#RANDOMIZE_DOMAINCENTER # shifts the particle distribution randomly each step to reduce correlations of force errors in time +#LEAN # selects a special 'lean' mode of code operation, which is optimized for aggressive memory saving +#LONG_X_BITS=2 # can be used to reduce periodic box-dimension in x-direction relative to nominal box size +#LONG_Y_BITS=2 # can be used to reduce periodic box-dimension in y-direction relative to nominal box size +#LONG_Z_BITS=1 # can be used to reduce periodic box-dimension in z-direction relative to nominal box size +#GRAVITY_TALLBOX=2 # this can be used to set-up gravity with two periodicity only in two of three dimensions +#TWODIMS # restricts SPH formulation to two dimensions (for tests) +#ONEDIMS # restricts SPH formulation to one dimension (for tests) +#GADGET2_HEADER # allows reading of snapshots with GADGET-2/3 header format +#SECOND_ORDER_LPT_ICS # treats second order LPT ICs generated with Adrian Jenkin's code -#GADGET2_HEADER # allows reading in of Gadget2/3 snapshots by using this code's header format for snaphot file formats 1 and 2 -#SECOND_ORDER_LPT_ICS # treats second order LPT ICs generated with Adrian Jenkin's code -#LEAN # selects a special 'lean' mode of code operation, which is optimized for extreme memory saving #--------------------------------------- Gravity calculation -SELFGRAVITY # switch on for self-gravity -HIERARCHICAL_GRAVITY -#FMM # enables FMM tree algorithm -ALLOW_DIRECT_SUMMATION # allows calculation of small number of particles with direct summation -#EXTERNALGRAVITY # switch on for external potential -#EVALPOTENTIAL # computes gravitational potential -#EXTRAPOTTERM # if this is set, the extra multipole term needed in principle for the potential is computed even though it does not enter the force -#EXTRA_HIGH_EWALD_ACCURACY # if this is activate, a third order Taylor expansion is used for the Ewald corrections + +SELFGRAVITY # switch to enable self-gravity of particles (typically always on) +#FMM # enables Fast Multipole Method instead of one-sided tree algorithm +#MULTIPOLE_ORDER=2 # sets the multipole order of Tree or FMM computations +#EVALPOTENTIAL # computes gravitational potential besides force +#EXTRAPOTTERM # this computes an extra multipole term for the potential which is not needed for the forces +#EXTRA_HIGH_EWALD_ACCURACY # this uses third-order instead of second-order Taylor expansion to interpolate Ewald corrections from table +#ALLOW_DIRECT_SUMMATION # allows calculation of direct summation gravity force if only a tiny number of particles as active +#EXTERNALGRAVITY # switches on inclusion of external gravitational potential -#MULTIPOLE_ORDER=2 # enables in tree and/or FMM for a value of 2 monopoles+dipoles (this is default), 3 gives quadrupoles, 4 octupoles, and 5 hexadecupoles -#RANDOMIZE_DOMAINCENTER #--------------------------------------- TreePM Options -PMGRID=512 -#ASMTH=1.25 -#RCUT=6.0 -#NTAB=128 # size of short-range look-up table -#TREEPM_NOTIMESPLIT # if this is activated, long-range and short-range gravity are time-integrated on a single timestep -#PLACEHIGHRESREGION=2 -#HRPMGRID=512 # High-res PM grid size (optional, default is HRPMGRID=PMGRID) - -#FFT_COLUMN_BASED -#PM_ZOOM_OPTIMIZED -#GRAVITY_TALLBOX=2 # this can be used to set-up gravity with two-dimensional boxes -#TREE_NUM_BEFORE_NODESPLIT=10 # Optional number that tells how many particles are allowed in a tree node before it is split + +#PMGRID=512 # basic mesh size for TreePM calculations +#ASMTH=1.25 # sets the smoothing of the PM force and thus the split-scale between Tree/FMM force and PM force +#RCUT=6.0 # cut-off radius beyond which Tree or FMM evaluations are stopped in TreePM / FMM-PM +#NTAB=128 # size of short-range look-up table +#PLACEHIGHRESREGION=2 # bitmask selecting the particle types that are used to define the location of a secondary PM mesh +#HRPMGRID=512 # dimension of high-res PM grid (optional, default is HRPMGRID=PMGRID) +#FFT_COLUMN_BASED # uses a column-based FFT algorithm instead of the default slab-based one +#PM_ZOOM_OPTIMIZED # selects a communication strategy in the PM code that is better balanced for zoom simulations +#TREE_NUM_BEFORE_NODESPLIT=4 # number of particles are are at most allowed in a tree node before it is split (can be 1) + + +#--------------------------------------- Time integration options + +#TREEPM_NOTIMESPLIT # if this is activated, long-range and short-range gravity are time-integrated on a common timestep +#HIERARCHICAL_GRAVITY # enables hierarchical time integration of the gravity +#FORCE_EQUAL_TIMESTEPS # this chooses a global timestep for all particles + #--------------------------------------- Treatment of gravitational softening -#INDIVIDUAL_GRAVITY_SOFTENING=4+8+16+32 -#NSOFTCLASSES=4 # Number of different softening values. Normally equal to number of particle types, but can be chosen differently if DECOUPLE_TYPE_AND_SOFTTYPE is set. -#ADAPTIVE_HYDRO_SOFTENING +#INDIVIDUAL_GRAVITY_SOFTENING=4+8+16+32 # bitmasks which selects the particle type(s) which pick their softening class based on particle mass +#NSOFTCLASSES=4 # number of different softening classes +#ADAPTIVE_HYDRO_SOFTENING # makes SPH gas particles pick an adaptive gravitational softening proportional to their SPH smoothing lengths #--------------------------------------- SPH treatmeant and formulation -#PRESSURE_ENTROPY_SPH # enables the Hopkins 2013 pressure-entropy formulation -#OUTPUT_PRESSURE_SPH_DENSITY # output also density computed in the pressure-entropy forumulation -#GAMMA=1.4 -#ISOTHERM_EQS -#REUSE_HYDRO_ACCELERATIONS_FROM_PREVIOUS_STEP -IMPROVED_VELOCITY_GRADIENTS # use higher-order gradients of the velocity according to Hu et. al (2014) -VISCOSITY_LIMITER_FOR_LARGE_TIMESTEPS #limits the acceleration due to the viscosity + +#REUSE_HYDRO_ACCELERATIONS_FROM_PREVIOUS_STEP # does not recompute the pressure forces after application of source functions +#VISCOSITY_LIMITER_FOR_LARGE_TIMESTEPS # limits the acceleration due to the viscosity +#PRESSURE_ENTROPY_SPH # enables the Hopkins (2013) pressure-entropy formulation, other density-entropy is used +#GAMMA=1.4 # sets the adiabatic index +#ISOTHERM_EQS # selects an isothermal equation of state +#IMPROVED_VELOCITY_GRADIENTS # use higher-order gradients of the velocities according to Hu et. al (2014) + #--------------------------------------- SPH kernels -CUBIC_SPLINE_KERNEL # uses the cubic spline kernel -#WENDLAND_C2_KERNEL # the Wendland C2 kernel from Dehnen & Aly 2012 -#WENDLAND_C4_KERNEL # the Wendland C4 kernel from Dehnen & Aly 2012 -#WENDLAND_C6_KERNEL # the Wendland C6 kernel from Dehnen & Aly 2012 -WENDLAND_BIAS_CORRECTION # reduces self-contribution for Wendland kernels according to Dehnen & Aly 2012 +#CUBIC_SPLINE_KERNEL # uses the cubic spline kernel (default) +#WENDLAND_C2_KERNEL # use the Wendland C2 kernel from Dehnen & Aly 2012 +#WENDLAND_C4_KERNEL # use the Wendland C4 kernel from Dehnen & Aly 2012 +#WENDLAND_C6_KERNEL # use the Wendland C6 kernel from Dehnen & Aly 2012 +#WENDLAND_BIAS_CORRECTION # reduces self-contribution for Wendland kernels according to Dehnen & Aly 2012 #--------------------------------------- SPH viscosity options -#TIMEDEP_ART_VISC # Enables time dependend viscosity -#HIGH_ART_VISC_START # Start with high rather than low viscosity -#NO_SHEAR_VISCOSITY_LIMITER # Turns of the shear viscosity supression -OUTPUT_VISCOSITY_PARAMETER + +#TIMEDEP_ART_VISC # enables time dependend viscosity +#HIGH_ART_VISC_START # start with high rather than low viscosity +#NO_SHEAR_VISCOSITY_LIMITER # turns of the shear viscosity supression #--------------------------------------- Extra physics -#COOLING -#STARFORMATION - -#--------------------------------------- Time integration options -#FORCE_EQUAL_TIMESTEPS # this chooses a variable but global timestep +#COOLING # Enables radiative atomic cooling by hydrogen and helium +#STARFORMATION # Enables star formation with the Springel & Hernquist (2003) model #---------------------------------------- Single/double precision and data types -POSITIONS_IN_32BIT # if set, use 32-integers for positions -POSITIONS_IN_64BIT # if set, use 64-integers for positions -POSITIONS_IN_128BIT # if set, use 128-integers for positions -DOUBLEPRECISION=1 -DOUBLEPRECISION_FFTW -#OUTPUT_IN_DOUBLEPRECISION # output files will be written in double precision -#ENLARGE_DYNAMIC_RANGE_IN_TIME # This extends the dynamic range of the integer timeline from 32 to 64 bit -#IDS_32BIT # Selects 32-bit IDs for internal storage (default) -#IDS_48BIT # Selects 48-bit IDs for internal storage -#IDS_64BIT # Selects 64-bit IDs for internal storage -#USE_SINGLEPRECISION_INTERNALLY # reduces default double precision for most internal computations to single precision + +#POSITIONS_IN_32BIT # if set, use 32-integers for positions (default for single precision) +#POSITIONS_IN_64BIT # if set, use 64-integers for positions (default for double precision) +#POSITIONS_IN_128BIT # if set, use 128-integers for positions (may only work on some systems) +DOUBLEPRECISION=1 # if activated and set to 1, use double precision internally, for 2 use mixed precision, otherwise single precision +#DOUBLEPRECISION_FFTW # if set, carries out FFTs in double precision, otherwise single precision +#USE_SINGLEPRECISION_INTERNALLY # reduces default double precision for most internal computations to single precision +#ENLARGE_DYNAMIC_RANGE_IN_TIME # eextend the dynamic range of the integer timeline from 32 to 64 bits +#IDS_32BIT # selects 32-bit IDs for internal storage (default) +#IDS_48BIT # selects 48-bit IDs for internal storage +#IDS_64BIT # selects 64-bit IDs for internal storage +#OUTPUT_IN_DOUBLEPRECISION # output files will be written in double precision + #---------------------------------------- Output/Input options -#REDUCE_FLUSH -#OUTPUT_VELOCITY_GRADIENT -#OUTPUT_PRESSURE -#OUTPUT_ENTROPY -#OUTPUT_CHANGEOFENTROPY -#OUTPUT_POTENTIAL -#OUTPUT_ACCELERATION -#OUTPUT_TIMESTEP -#OUTPUT_DIVVEL # output velocity divergence -#OUTPUT_CURLVEL # output velocity curl -#OUTPUT_COOLHEAT # output actual energy loss/gain in cooling/heating routine -#POWERSPEC_ON_OUTPUT -#OUTPUT_NON_SYNCHRONIZED_ALLOWED -#OUTPUT_VELOCITIES_IN_HALF_PRECISION -#OUTPUT_ACCELERATIONS_IN_HALF_PRECISION -#OUTPUT_COORDINATES_AS_INTEGERS -#ALLOW_HDF5_COMPRESSION + +#OUTPUT_VELOCITY_GRADIENT # output velocity gradients +#OUTPUT_PRESSURE # output gas pressure +#OUTPUT_ENTROPY # output gas entropy +#OUTPUT_CHANGEOFENTROPY # output rate of change of entropy +#OUTPUT_POTENTIAL # output gravitational potential +#OUTPUT_ACCELERATION # output total acceleration +#OUTPUT_TIMESTEP # output timesteps +#OUTPUT_DIVVEL # output velocity divergence +#OUTPUT_CURLVEL # output velocity curl +#OUTPUT_COOLHEAT # output actual energy loss/gain in cooling/heating routine +#OUTPUT_PRESSURE_SPH_DENSITY # output also density computed in the pressure-entropy forumulation +#OUTPUT_VISCOSITY_PARAMETER # output current viscosity parameter in SPH particle output (only for time-dependent viscosity) +#OUTPUT_NON_SYNCHRONIZED_ALLOWED # allow snapshot creation also at steps that are not fully synchronized +#OUTPUT_VELOCITIES_IN_HALF_PRECISION # special option to store velocities in reduced precision +#OUTPUT_ACCELERATIONS_IN_HALF_PRECISION # special option to store acclerations in reduced precision +#OUTPUT_COORDINATES_AS_INTEGERS # special option to store coordinates as integers that are used internally +#POWERSPEC_ON_OUTPUT # computes a matter power spectrum when the code writes a snapshot output +#ALLOW_HDF5_COMPRESSION # applies HDF5 loss-less compression to selected output fields +#REDUCE_FLUSH # do not flush the I/O streams of the log-files every system step + #---------------------------------------- On the fly FOF groupfinder -#FOF # enable FoF output -#FOF_PRIMARY_LINK_TYPES=2 # 2^type for the primary dark matter type -#FOF_SECONDARY_LINK_TYPES=1+16+32 # 2^type for the types linked to nearest primaries -##FOF_GROUP_MIN_LEN=32 # default is 32 -##FOF_LINKLENGTH=0.16 # Linkinglength for FoF (default=0.2) -#FOF_ALLOW_HUGE_GROUPLENGTH # if this is set, groups and subhalos may have more than 2 billion particles in length + +#FOF # enable FoF output +#FOF_PRIMARY_LINK_TYPES=2 # bitmask, 2^type for the primary dark matter type +#FOF_SECONDARY_LINK_TYPES=1+16+32 # bitmask, 2^type for the type(s) linked to nearest primary particle +#FOF_GROUP_MIN_LEN=32 # minimum group length (default is 32) +#FOF_LINKLENGTH=0.2 # linking length for FoF (default=0.2) +#FOF_ALLOW_HUGE_GROUPLENGTH # set if there are groups and subhalos with more than 2 billion particles in length + #---------------------------------------- Subfind -SUBFIND # enables substructure finder -#SUBFIND_STORE_LOCAL_DENSITY # will calculate local densities and velocity dispersions for *all* particles - # selected by FOF primary and seconday link types (not only for particles in gropus), - # and store them in snapshots -#SUBFIND_ORPHAN_TREATMENT # creates special snapshots with formerly most bound particles -#SUBFIND_HBT # use previous subhalo catalogue instead of density excursions to define subhalo candidates + +#SUBFIND # enables substructure finder +#SUBFIND_HBT # use previous subhalo catalogue instead of density excursions to define subhalo candidates +#SUBFIND_STORE_LOCAL_DENSITY # calculates local densities and velocity dispersions for all particles and stores them in snapshots +#SUBFIND_ORPHAN_TREATMENT # creates special snapshots with formerly most bound particles + #---------------------------------------- Merger tree code -MERGERTREE # enables on-the-fly calculation of descendants, and merger tree construction from group catalogues + +#MERGERTREE # enables on-the-fly calculation of descendants, and merger tree construction from group catalogues + #---------------------------------------- On-the-fly lightcone creation -#LIGHTCONE -#LIGHTCONE_PARTICLES -#LIGHTCONE_MASSMAPS -#LIGHTCONE_PARTICLES_GROUPS -#LIGHTCONE_OUTPUT_ACCELERATIONS -#LIGHTCONE_IMAGE_COMP_HSML_VELDISP - -#REARRANGE_OPTION - -#--------------------------------------- IC generation -#NGENIC=256 # generate cosmological ICs, set NGENIC to the FFT grid size used for IC generation -#CREATE_GRID # start with a regular cartesian DM particle grid, instead of reading a glass file (for NGENIC) -#GENERATE_GAS_IN_ICS # add SPH particles to dark matter only ICs -#SPLIT_PARTICLE_TYPE=4+8 # particle types to be split if GENERATE_GAS_IN_ICS is activated -#NGENIC_FIX_MODE_AMPLITUDES # when activated, this leaves the mode amplitudes at sqrt(P(k)), instead of sampling from a Rayleigh distribution -#NGENIC_MIRROR_PHASES # if this is activated, all phases are turned by 180 degrees -#NGENIC_2LPT # applies 2LPT instead of just Zeldovich -#NGENIC_TEST # can be used to create ICs, measure the power spectrum, and then stop + +#LIGHTCONE # master option for lightcone output creation +#LIGHTCONE_PARTICLES # produces particle lightcones +#LIGHTCONE_MASSMAPS # produces mass shells on the lightcone +#LIGHTCONE_PARTICLES_GROUPS # computes groups for particles buffered on the lightcone +#LIGHTCONE_OUTPUT_ACCELERATIONS # stores accelerations for particles on lightcone +#LIGHTCONE_IMAGE_COMP_HSML_VELDISP # option for computing densities and smoothing length for lightcones in postprocessing +#REARRANGE_OPTION # special option to reorder lightcone data in mergertree order + + +#--------------------------------------- IC creation + +#NGENIC=256 # generate cosmological ICs, set NGENIC to the FFT grid size used for IC generation +#NGENIC_2LPT # applies 2LPT instead of just Zeldovich approximation +#CREATE_GRID # start with a regular Cartesian DM particle grid, instead of reading a glass file (for NGENIC) +#GENERATE_GAS_IN_ICS # add SPH particles to created or read dark matter only ICs +#SPLIT_PARTICLE_TYPE=4+8 # specifies particle types to be split if GENERATE_GAS_IN_ICS is activated +#NGENIC_FIX_MODE_AMPLITUDES # when activated, this leaves the mode amplitudes at sqrt(P(k)), instead of sampling from a Rayleigh distribution +#NGENIC_MIRROR_PHASES # if this is activated, all phases are turned by 180 degrees +#NGENIC_TEST # can be used to create ICs, measure the power spectrum, and then stop + + +#----------------------------------------Parallelization options + +#IMPOSE_PINNING # enables pinning of MPI processes to CPU cores +#IMPOSE_PINNING_OVERRIDE_MODE # tries to do the pinning even if a prior pinning is detected +#PRESERVE_SHMEM_BINARY_INVARIANCE # preserve binary invariance of results despite machine weather, at the price of more tree walk overhead +#EXPLICIT_VECTORIZATION # use AVX at selected places in SPH kernels through the vectorclass C++ library +#SIMPLE_DOMAIN_AGGREGATION # this is an experimental modification of the domain decomposition algorithm (can either help or harm performance) + #---------------------------------------- MPI related settings -#ISEND_IRECV_IN_DOMAIN -#USE_MPIALLTOALLV_IN_DOMAINDECOMP -#MPI_HYPERCUBE_ALLGATHERV # some MPI-libraries may use quite a bit of internal storage for MPI_Allgatherv. This uses hypercubes instead as a work-around -#MPI_MESSAGE_SIZELIMIT_IN_MB=200 -#NUMPART_PER_TASK_LARGE # Set this if the number of particle per task is very large (so that more than 2GB data is comprised just by the particle data) -#MAX_NUMBER_OF_RANKS_WITH_SHARED_MEMORY=64 -#NUMBER_OF_MPI_LISTENERS_PER_NODE=1 +#NUMBER_OF_MPI_LISTENERS_PER_NODE=1 # set such that the number of MPI-ranks per node and listener is maller than MAX_NUMBER_OF_RANKS_WITH_SHARED_MEMORY +#MAX_NUMBER_OF_RANKS_WITH_SHARED_MEMORY=64 # default is 64, but can also be set to 32 +#NUMPART_PER_TASK_LARGE # set this if the number of particles per task is so large that more than 2 GB are comprised just by particle data +#USE_MPIALLTOALLV_IN_DOMAINDECOMP # replaces hypercube communication in domain particle exchance with a single MPI_Allgatherv (can be less stable) +#MPI_HYPERCUBE_ALLGATHERV # if your MPI-library uses too much internal storage for MPI_Allgatherv, this uses a hypercube as a work-around +#MPI_MESSAGE_SIZELIMIT_IN_MB=200 # limit the message size of very large MPI transfers +#ISEND_IRECV_IN_DOMAIN # uses asynchronous communication instead of synchronous communication in hypercube pattern (can be less stable) + #---------------------------------------- Testing and Debugging options -DEBUG # enables core-dumps -#DEBUG_ENABLE_FPU_EXCEPTIONS # tries to enable FPU exceptions -#DEBUG_SYMTENSORS # carries out a few unit tests related to the tensor algebra routines of the code -HOST_MEMORY_REPORTING # reports after start-up the available system memory by analyzing /proc/meminfo -#FORCETEST=0.001 # calculates for given fraction of particles direct summation forces to check accuracy of tree force -#FORCETEST_FIXEDPARTICLESET # if this set, always the same particle receive a force accuracy check during a run -#FORCETEST_TESTFORCELAW=1 # Special option for measuring the effective force law, can be set to 1 or 2 for TreePM+HighRes -#VTUNE_INSTRUMENT -#DEBUG_MD5 # make MD5 check sum routines available for debugging -#SQUASH_TEST -#DOMAIN_SPECIAL_CHECK -#RECREATE_UNIQUE_IDS # Before carrying out the ID uniqueness test, this sets new IDs. Can be used if IC files contain broken IDs. -#NO_STOP_BELOW_MINTIMESTEP # do not stop when the code wants to go below minimum timestep -#TILING=2 # duplicated an IC in each dimension -#DO_NOT_PRODUCE_BIG_OUTPUT # for scaling tests, one may disable the writing of restart files and snapshot dumps -#MEASURE_TOTAL_MOMENTUM -#ENABLE_HEALTHTEST -#EWALD_TEST # a development test for the Ewald tables -#STOP_AFTER_STEP=10 # ends a simulation after the specified timestep (to simplify scalability tests) - -#TREE_NO_SAFETY_BOX # when set, this disables the geometric 'near node' protection + +#HOST_MEMORY_REPORTING # reports after start-up the available system memory by analyzing /proc/meminfo (active by default on Linux) +#ENABLE_HEALTHTEST # can be used to enable some auto-checking of the machine you're running on +#DEBUG # enables core-dumps in case MPI_Init() has disabled them +#DEBUG_ENABLE_FPU_EXCEPTIONS # tries to enable FPU exceptions so that a bad floating point operation triggers a core dump +#DEBUG_SYMTENSORS # carries out a few unit tests related to the tensor algebra routines of the code +#FORCETEST=0.001 # calculates for given fraction of particles direct summation forces to check accuracy of tree force +#FORCETEST_FIXEDPARTICLESET # if this set, always the same particle receive a force accuracy check during a run +#FORCETEST_TESTFORCELAW=1 # special option for measuring the effective force law, can be set to 1 or 2 for TreePM+HighRes +#VTUNE_INSTRUMENT # outputs auxiliary info to simplify Vtune performance analysis +#DEBUG_MD5 # make MD5 check sum routines available for debugging and code development +#SQUASH_TEST # special check-option for code development +#DOMAIN_SPECIAL_CHECK # special check-option for code development +#RECREATE_UNIQUE_IDS # before carrying out the ID uniqueness test, this sets new IDs (e.g. to fix broken ICs) +#NO_STOP_BELOW_MINTIMESTEP # do not stop when the code wants to go below a prescribed minimum timestep +#TILING=2 # duplicate an IC in each dimension for scaling tests +#DO_NOT_PRODUCE_BIG_OUTPUT # for scaling tests, disable the writing of restart files and snapshot dumps +#MEASURE_TOTAL_MOMENTUM # special test option for code development +#EWALD_TEST # a development test for the Ewald tables +#STOP_AFTER_STEP=10 # ends a simulation after the specified timestep (to simplify scalability tests) +#TREE_NO_SAFETY_BOX # when set, this disables the geometric 'near node' protection