Dear Gromacs users, I'm having problems with g_covar because sometimes it produces eigenvectors with nan values. This is related with a previous post (http://lists.gromacs.org/pipermail/gmx-users/2012-July/073492.html). I checked the integrity of my xtc and tpr/pdb input files, and everything seems fine, as the following works OK: echo System System | /usr/local/gromacs455/bin/g_covar455 -s mainchain_init.tpr -f mainchain_apo.xtc -o eigenvalues.xvg -v eigenvectors.trr -xpma covara.xpm -xpm covar.xpm -l covar.log
but running g_covar with parts of the same trajectory gives eigenvectors with nan values: echo System System | /usr/local/gromacs455/bin/g_covar455 -s mainchain_init.tpr -f mainchain_apo.xtc -b 2000 -e 22000 -o eigenvalues_1.xvg -v eigenvectors_1.trr -xpma covara_1.xpm -xpm covar_1.xpm -l covar_1.log -av av_1.pdb -ascii covar_1.dat echo System System | /usr/local/gromacs455/bin/g_covar455 -s mainchain_init.tpr -f mainchain_apo.xtc -b 22000 -e 42000 -o eigenvalues_2.xvg -v eigenvectors_2.trr -xpma covara_2.xpm -xpm covar_2.xpm -l covar_2.log -av av_2.pdb -ascii covar_2.dat eigenvectors_1.trr has one frame with all coordinates set to nan; eigenvectors_2.trr has 17. I don't see anything obviously wrong with the other files produced (xvg, xpm, pdb). What could be the problem? Thanks for your time, Rui Rodrigues-- gmx-users mailing list [email protected] http://lists.gromacs.org/mailman/listinfo/gmx-users * Only plain text messages are allowed! * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to [email protected]. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
