with variable block timesteps
rungravidy in= outdir= ...
gravidy is not installed by default in NEMO. See INSTALL below. The rungravidy NEMO frontend is available for those who want to work with NEMO’s snapshot(5NEMO) files.
Note the current version of rungravidy does not support the -p (screen output) and -r (resume an old run) options.
GraviDy: -h [ --help ] Display this message Required options: -i [ --input ] <filename> Input data filename -t [ --time ] <value> (=1) Integration time (In N-body units) -r [ --resume ] <filename>.info Resume a simulation with an .info file Elective options: -o [ --output ] <filename> Output data filename -s [ --softening ] <value> Softening parameter (default 1e-4) -e [ --eta ] <value> ETA of time-step calculation (default 0.01) -p [ --screen ] Print summary in the screen instead of a file -z [ --interval ] <value> Output time interval (default 0.125) Extra options: -l [ --lagrange ] Print information of the Lagrange Radii in every integration time -a [ --all ] Print all the information of N-particles in every integration time -g [ --gpu ] <value> (=0) GPUs to use, by default is the maximum available devices (use even numbers)
mkplummer p128 128 snapprint p128 i,m,x,y,z,vx,vy,vz format=%.15f > p128.tab gravidy --input p128.tab --output p128 --interval 0.1 --time 10 --eta 0.01 --softening 0.05 cat p128.out.snapshot_0* | sed ’s/# Time://’ | tabtos - - time skip,m,pos,vel nbody=128 | csf - p128.out SnapShotalthough note that the last snapshot (time=10.0 in this case) might be replicated in versions of gravidy before June 2023. The rungravidy frontend would simplify this to
mkplummer p128 128 mkdir p128.d rungravidy p128 p128.d deltat=0.1 tcrit=10 eta=0.01 eps=0.05 snapplot p128.d/OUT3.snapwith the caveat that rungravidy enforces the use of a run directory.
cd $NEMO/usr/gravidy make bench1and some example timings
k2-bat 21.9 i7-1260P lma 35.0? jansky 12.0? i9-12xxxK
cd $NEMO/usr/gravidy make clean install make cmp34this will install the vanilla "cpu" version. The "gpu" and "mpi" version will need more guidance.
The "cmp34" target will compare the results of a direct run of gravidy with the same one of rungravidy. There should be no differences.
NEMO/usr/gravidy - root directory rungravidy.c - example NEMO frontend
Maureira-Fredes, C and Amaro-Seoane, P. "GraviDy, a GPU modular, parallel direct-summation N-body integrator: Dynamics with softening", MNRAS, vol. 473, pp. 3113-3127, Jan . 2018. https://ui.adsabs.harvard.edu/abs/2018MNRAS.473.3113M
http://github.com/cmaureir/gravidy.git
12-feb-2017 example added PJT 6-apr-2023 add install note PJT 16-jun-2023 (rungravidy v0.4 now using gravidy default parameters PJT 21-jun-2023 fixed gravidy last duplicated snapshot PJT