Gmx rdf taking way too long reading frames – User discussions

GROMACS version: 2021
GROMACS modification: /No
Here post your question

I have a trajectory files of 200 ns simulation of a system with 488007 atoms, around 20000 frames and it takes about 10-20 minutes to read one frame on my HPC environment of 1 GPU, 10 CPUs per task.

I used a simple command something like :

gmx rdf -f trajectory.xtc -s topol.tpr -o rdf.xvg -ref Protein -sel SOL

If I read the same xtc file with e.g. gmx check, it takes about less than 30 minutes to read the whole frames.

Is something wrong with my input files or does it take usually so long for gmx rdf to read frames?

How can I improve the performance calculating rdf?

My guess is that this is not a reading issue, but an O(#atoms^2) operations issue for the RDF due to half a million atoms. Setting the -rmax option to something short sould give a significant speedup.



1 Like

thank you sir for your help. My mistak not considering the rmax. I set rmax to 3 nm and it speeded up already very much.

Read more here: Source link