{"type":"doc","content":[{"type":"paragraph","content":[{"text":"GROMACS is a versatile package to perform molecular dynamics, i.e. simulate the Newtonian equations of motion for systems with hundreds to millions of particles. It is primarily designed for biochemical molecules like proteins, lipids and nucleic acids that have a lot of complicated bonded interactions, but since GROMACS is extremely fast at calculating the nonbonded interactions many groups are also using it for research on non-biological systems.","type":"text"}]},{"type":"paragraph","content":[{"text":"Some test cases to start with can be found at ","type":"text"},{"text":"ftp://ftp.gromacs.org/pub/benchmarks/gmxbench-3.0.tar.gz","type":"text","marks":[{"type":"link","attrs":{"href":"ftp://ftp.gromacs.org/pub/benchmarks/gmxbench-3.0.tar.gz"}}]},{"text":".","type":"text"},{"type":"hardBreak"}]},{"type":"extension","attrs":{"layout":"default","extensionType":"com.atlassian.confluence.macro.core","extensionKey":"toc","parameters":{"macroParams":{},"macroMetadata":{"macroId":{"value":"065ecc75-aaf3-40f3-9e0e-a44c612a6b99"},"schemaVersion":{"value":"1"},"title":"Table of Contents"}}}},{"type":"heading","attrs":{"level":1},"content":[{"text":"Refereneces ","type":"text"}]},{"type":"bulletList","content":[{"type":"listItem","content":[{"type":"paragraph","content":[{"text":"http://manual.gromacs.org/documentation/current/install-guide/index.html","type":"text","marks":[{"type":"link","attrs":{"href":"http://manual.gromacs.org/documentation/current/install-guide/index.html"}}]}]}]}]},{"type":"paragraph"},{"type":"heading","attrs":{"level":1},"content":[{"text":"Build and Install Gromacs","type":"text"}]},{"type":"heading","attrs":{"level":2},"content":[{"text":"Download the source code","type":"text"}]},{"type":"table","attrs":{"layout":"default"},"content":[{"type":"tableRow","content":[{"type":"tableCell","attrs":{"colspan":1,"background":"initial","rowspan":1},"content":[{"type":"codeBlock","content":[{"text":"wget ftp://ftp.gromacs.org/pub/gromacs/gromacs-2020.2.tar.gz","type":"text"}]}]}]}]},{"type":"heading","attrs":{"level":2},"content":[{"text":"Build Gromacs for CPU only","type":"text"}]},{"type":"table","attrs":{"layout":"default"},"content":[{"type":"tableRow","content":[{"type":"tableCell","attrs":{"colspan":1,"background":"initial","rowspan":1},"content":[{"type":"codeBlock","content":[{"text":"tar xfz gromacs-2020.2.tar.gz\ncd gromacs-2020.2\nmodule load intel/2019.5.281\nmodule load mkl/2019.5.281\nmodule load gcc/8.4.0\nmodule load cmake/3.13.4\nmodule load hpcx/2.6.0\nmkdir build install run\ncd build\ncmake .. -DGMX_FFT_LIBRARY=mkl -DMKL_LIBRARIES=-mkl \\\n -DMKL_INCLUDE_DIR=$MKLROOT/include \\\n -DGMX_SIMD=AVX2_256 \\\n -DGMX_MPI=ON \\\n -DGMX_BUILD_MDRUN_ONLY=on \\\n -DBUILD_SHARED_LIBS=on \\\n -DGMX_HWLOC=off \\\n -DCMAKE_INSTALL_PREFIX=../install \\\n -DCMAKE_C_COMPILER=mpicc -DCMAKE_CXX_COMPILER=mpicxx\nmake -j 16 install","type":"text"}]}]}]}]},{"type":"paragraph","content":[{"text":"Check the install directory for the file ","type":"text"},{"text":"mdrun_mpi","type":"text","marks":[{"type":"strong"}]}]},{"type":"codeBlock","content":[{"text":"$ ls ../install/\nbin include lib64 share\n$ ls ../install/bin \ngmx-completion-mdrun_mpi.bash mdrun_mpi","type":"text"}]},{"type":"paragraph"},{"type":"heading","attrs":{"level":2},"content":[{"text":"Build Gromacs with GPU support","type":"text"}]},{"type":"table","attrs":{"layout":"default"},"content":[{"type":"tableRow","content":[{"type":"tableCell","attrs":{"colspan":1,"background":"initial","rowspan":1},"content":[{"type":"codeBlock","content":[{"text":"tar xfz gromacs-2020.2.tar.gz\ncd gromacs-2020.2\nmodule load intel/2019.5.281\nmodule load mkl/2019.5.281\nmodule load gcc/8.4.0\nmodule load cmake/3.13.4\nmodule load hpcx/2.6.0\nmodule load cuda/10.1\nmkdir build install run\ncd build\ncmake .. -DGMX_FFT_LIBRARY=mkl -DMKL_LIBRARIES=-mkl \\\n -DMKL_INCLUDE_DIR=$MKLROOT/include \\\n -DGMX_SIMD=AVX2_256 \\\n -DGMX_MPI=ON \\\n -DGMX_GPU=ON \\\n -DGMX_BUILD_MDRUN_ONLY=on \\\n -DBUILD_SHARED_LIBS=on \\\n -DGMX_HWLOC=off \\\n -DCMAKE_INSTALL_PREFIX=../install \\ \n -DCMAKE_C_COMPILER=mpicc -DCMAKE_CXX_COMPILER=mpicxx\nmake -j 16 install","type":"text"}]}]}]}]},{"type":"paragraph"},{"type":"heading","attrs":{"level":2},"content":[{"text":"Run Gromacs ","type":"text"}]},{"type":"paragraph"},{"type":"paragraph","content":[{"text":"To run Gromacs with CPU only using stmv case","type":"text","marks":[{"type":"strong"}]}]},{"type":"codeBlock","content":[{"text":"% mpirun -np 40 -x UCX_NET_DEVICES=mlx5_0:1 -bind-to core -report-bindings \\\nmdrun_mpi -v -s stmv.tpr -nsteps 10000 -noconfout -nb cpu -pin on\n\nCommand line:\n mdrun_mpi -v -s stmv.tpr -nsteps 10000 -noconfout -nb cpu -pin on\n\nReading file stmv.tpr, VERSION 2018.1 (single precision)\nNote: file tpx version 112, software tpx version 119\nOverriding nsteps with value passed on the command line: 10000 steps, 20 ps\nChanging nstlist from 10 to 80, rlist from 1.2 to 1.316\n\n\nUsing 40 MPI processes\nUsing 1 OpenMP thread per MPI process\n...\nstep 9900, remaining wall clock time: 6 s\nvol 0.96 imb F 1% pme/F 0.81 step 10000, remaining wall clock time: 0 s\n\n\nDynamic load balancing report:\n DLB was turned on during the run due to measured imbalance.\n Average load imbalance: 0.9%.\n The balanceable part of the MD step is 94%, load imbalance is computed from this.\n Part of the total run time spent waiting due to load imbalance: 0.8%.\n Steps where the load balancing was limited by -rdd, -rcon and/or -dds: X 0 % Y 0 % Z 0 %\n Average PME mesh/force load: 0.790\n Part of the total run time spent waiting due to PP/PME imbalance: 2.0 %\n\n\n Core t (s) Wall t (s) (%)\n Time: 27556.777 688.921 4000.0\n (ns/day) (hour/ns)\nPerformance: 2.509 9.567","type":"text"}]},{"type":"paragraph"},{"type":"paragraph","content":[{"text":"To run Gromacs with GPU","type":"text","marks":[{"type":"strong"}]}]},{"type":"codeBlock","content":[{"text":"% export OMP_NUM_THREADS=2\n% export KMP_AFFINITY=verbose,compact\n% mpirun -np 4 -x UCX_NET_DEVICES=mlx5_0:1 -map-by node:PE=$OMP_NUM_THREADS \\\nmdrun_mpi -v -s stmv.tpr -nsteps 10000 -noconfout -nb gpu -pin on \\\n-ntomp $OMP_NUM_THREADS \n\nCommand line:\n mdrun_mpi -v -s stmv.tpr -nsteps 10000 -noconfout -nb gpu -pin on\n\nReading file stmv.tpr, VERSION 2018.1 (single precision)\nNote: file tpx version 112, software tpx version 119\nOverriding nsteps with value passed on the command line: 10000 steps, 20 ps\nChanging nstlist from 10 to 100, rlist from 1.2 to 1.339\n\n\nOn host ops003.hpcadvisorycouncil.com 2 GPUs selected for this run.\nMapping of GPU IDs to the 2 GPU tasks in the 2 ranks on this node:\n PP:0,PP:1\nPP tasks will do (non-perturbed) short-ranged and most bonded interactions on the GPU\nPP task will update and constrain coordinates on the CPU\nUsing 4 MPI processes\nUsing 2 OpenMP threads per MPI process\n...\nimb F 0% step 9900, remaining wall clock time: 3 s\nimb F 0% step 10000, remaining wall clock time: 0 s\n\n\nDynamic load balancing report:\n DLB was off during the run due to low measured imbalance.\n Average load imbalance: 0.2%.\n The balanceable part of the MD step is 59%, load imbalance is computed from this.\n Part of the total run time spent waiting due to load imbalance: 0.1%.\n\n\n Core t (s) Wall t (s) (%)\n Time: 2581.428 322.681 800.0\n (ns/day) (hour/ns)\nPerformance: 5.356 4.481","type":"text"}]},{"type":"paragraph"},{"type":"heading","attrs":{"level":1},"content":[{"text":"Other command options ","type":"text"}]},{"type":"codeBlock","content":[{"text":"OPTIONS\n\nOptions to specify input files:\n\n -s [<.tpr>] (topol.tpr)\n Portable xdr run input file\n -cpi [<.cpt>] (state.cpt) (Opt.)\n Checkpoint file\n -table [<.xvg>] (table.xvg) (Opt.)\n xvgr/xmgr file\n -tablep [<.xvg>] (tablep.xvg) (Opt.)\n xvgr/xmgr file\n -tableb [<.xvg> [...]] (table.xvg) (Opt.)\n xvgr/xmgr file\n -rerun [<.xtc/.trr/...>] (rerun.xtc) (Opt.)\n Trajectory: xtc trr cpt gro g96 pdb tng\n -ei [<.edi>] (sam.edi) (Opt.)\n ED sampling input\n -multidir [ [...]] (rundir) (Opt.)\n Run directory\n -awh [<.xvg>] (awhinit.xvg) (Opt.)\n xvgr/xmgr file\n -membed [<.dat>] (membed.dat) (Opt.)\n Generic data file\n -mp [<.top>] (membed.top) (Opt.)\n Topology file\n -mn [<.ndx>] (membed.ndx) (Opt.)\n Index file\n\nOptions to specify output files:\n\n -o [<.trr/.cpt/...>] (traj.trr)\n Full precision trajectory: trr cpt tng\n -x [<.xtc/.tng>] (traj_comp.xtc) (Opt.)\n Compressed trajectory (tng format or portable xdr format)\n -cpo [<.cpt>] (state.cpt) (Opt.)\n Checkpoint file\n -c [<.gro/.g96/...>] (confout.gro)\n Structure file: gro g96 pdb brk ent esp\n -e [<.edr>] (ener.edr)\n Energy file\n -g [<.log>] (md.log)\n Log file\n -dhdl [<.xvg>] (dhdl.xvg) (Opt.)\n xvgr/xmgr file\n -field [<.xvg>] (field.xvg) (Opt.)\n xvgr/xmgr file\n -tpi [<.xvg>] (tpi.xvg) (Opt.)\n xvgr/xmgr file\n -tpid [<.xvg>] (tpidist.xvg) (Opt.)\n xvgr/xmgr file\n -eo [<.xvg>] (edsam.xvg) (Opt.)\n xvgr/xmgr file\n -px [<.xvg>] (pullx.xvg) (Opt.)\n xvgr/xmgr file\n -pf [<.xvg>] (pullf.xvg) (Opt.)\n xvgr/xmgr file\n -ro [<.xvg>] (rotation.xvg) (Opt.)\n xvgr/xmgr file\n -ra [<.log>] (rotangles.log) (Opt.)\n Log file\n -rs [<.log>] (rotslabs.log) (Opt.)\n Log file\n -rt [<.log>] (rottorque.log) (Opt.)\n Log file\n -mtx [<.mtx>] (nm.mtx) (Opt.)\n Hessian matrix\n -if [<.xvg>] (imdforces.xvg) (Opt.)\n xvgr/xmgr file\n -swap [<.xvg>] (swapions.xvg) (Opt.)\n xvgr/xmgr file\n\nOther options:\n\n -deffnm \n Set the default filename for all file options\n -xvg (xmgrace)\n xvg plot formatting: xmgrace, xmgr, none\n -dd (0 0 0)\n Domain decomposition grid, 0 is optimize\n -ddorder (interleave)\n DD rank order: interleave, pp_pme, cartesian\n -npme (-1)\n Number of separate ranks to be used for PME, -1 is guess\n -nt (0)\n Total number of threads to start (0 is guess)\n -ntmpi (0)\n Number of thread-MPI ranks to start (0 is guess)\n -ntomp (0)\n Number of OpenMP threads per MPI rank to start (0 is guess)\n -ntomp_pme (0)\n Number of OpenMP threads per MPI rank to start (0 is -ntomp)\n -pin (auto)\n Whether mdrun should try to set thread affinities: auto, on, off\n -pinoffset (0)\n The lowest logical core number to which mdrun should pin the first\n thread\n -pinstride (0)\n Pinning distance in logical cores for threads, use 0 to minimize\n the number of threads per physical core\n -gpu_id \n List of unique GPU device IDs available to use\n -gputasks \n List of GPU device IDs, mapping each PP task on each node to a\n device\n -[no]ddcheck (yes)\n Check for all bonded interactions with DD\n -rdd (0)\n The maximum distance for bonded interactions with DD (nm), 0 is\n determine from initial coordinates\n -rcon (0)\n Maximum distance for P-LINCS (nm), 0 is estimate\n -dlb (auto)\n Dynamic load balancing (with DD): auto, no, yes\n -dds (0.8)\n Fraction in (0,1) by whose reciprocal the initial DD cell size will\n be increased in order to provide a margin in which dynamic load\n balancing can act while preserving the minimum cell size.\n -nb (auto)\n Calculate non-bonded interactions on: auto, cpu, gpu\n -nstlist (0)\n Set nstlist when using a Verlet buffer tolerance (0 is guess)\n -[no]tunepme (yes)\n Optimize PME load between PP/PME ranks or GPU/CPU\n -pme (auto)\n Perform PME calculations on: auto, cpu, gpu\n -pmefft (auto)\n Perform PME FFT calculations on: auto, cpu, gpu\n -bonded (auto)\n Perform bonded calculations on: auto, cpu, gpu\n -update (auto)\n Perform update and constraints on: auto, cpu, gpu\n -[no]v (no)\n Be loud and noisy\n -pforce (-1)\n Print all forces larger than this (kJ/mol nm)\n -[no]reprod (no)\n Try to avoid optimizations that affect binary reproducibility\n -cpt (15)\n Checkpoint interval (minutes)\n -[no]cpnum (no)\n Keep and number checkpoint files\n -[no]append (yes)\n Append to previous output files when continuing from checkpoint\n instead of adding the simulation part number to all file names\n -nsteps (-2)\n Run this number of steps (-1 means infinite, -2 means use mdp\n option, smaller is invalid)\n -maxh (-1)\n Terminate after 0.99 times this time (hours)\n -replex (0)\n Attempt replica exchange periodically with this period (steps)\n -nex (0)\n Number of random exchanges to carry out each exchange interval (N^3\n is one suggestion). -nex zero or not specified gives neighbor\n replica exchange.\n -reseed (-1)\n Seed for replica exchange, -1 is generate a seed","type":"text"}]}],"version":1}

Browser not supported