Thanks for your help, Justin. I ran gmxcheck -c topol.tpr ("12284 atoms in file") and gmxcheck -f traj.xtc ("# Atoms 29"). The difference is the SOL atoms. I made another .tpr file with xtcgroups = TDR SOL, and now gmxcheck shows that the number of atoms in topol.tpr is 12284, but I still can't get g_covar to work. Maybe these are unrelated problems, I'm not sure.<br>
~Dayle<br>
<br><div class="gmail_quote">On Wed, Apr 1, 2009 at 2:49 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div><br>
<br>
Dayle Smith wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Tsjerk,<br>
Thank you for your quick and helpful response. I defined "xtc_grps = TDR" in my .mdp file (then I use grompp -f mdpme.mdp -c confout.gro -p topol.top -o topol.tpr -np 8 -n index.ndx, index contains TDR and SOL) and run it with -np 8. I'm afraid that I don't understand what "shuffling" or "matching series" refers to (a clue that I'm doing something wrong). This g_covar error message is probably another clue:<br>
<br>
WARNING: number of atoms in tpx (29) and trajectory (29) do not match<br>
-------------------------------------------------------<br>
Program g_covar, VERSION 3.3.3<br>
Source code file: nrjac.c, line: 129<br>
Fatal error:<br>
Error: Too many iterations in routine JACOBI<br>
<br>
So the # of atoms is the same, but some other key ingredient doesn't match. Can you please enlighten me?!<br>
</blockquote>
<br></div>
I think the output error message is bizarre, but you still have a number of atoms in the .tpr that does not match the .xtc. If your simulation is of TDR and SOL, then those groups will be in your topol.tpr. If your xtc-grps specify only TDR, then there will be a coordinate mismatch. Run the following:<br>
<br>
gmxcheck -c topol.tpr<br>
gmxcheck -f traj.xtc<br>
<br>
and see if gmxcheck reports the same number of atoms in both files. This is the quickest way to know for sure.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Thanks a lot,<br>
Dayle<div><br>
<br>
<br>
<br>
<br>
On Wed, Apr 1, 2009 at 11:12 AM, Tsjerk Wassenaar <<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a> <mailto:<a href="mailto:tsjerkw@gmail.com" target="_blank">tsjerkw@gmail.com</a>>> wrote:<br>
<br>
Hi Dayle,<br>
<br>
Errm, really, the only cases I know of this error to occur is when I<br>
had a mismatch between the reference and trajectory. Did you specify<br>
xtc-groups? Did you shuffle the system? How did you assert that you<br>
have matching series? Have you tried using the reference and the<br>
trajectory to convert (part of) the trajectory to .pdb and visualize?<br>
If all else fails, can you send (a link to) an archive containing a<br>
single frame from the trajectory and the reference?<br>
<br>
Cheers,<br>
<br>
Tsjerk<br>
<br>
2009/4/1 Dayle Smith <<a href="mailto:daylemariesmith@gmail.com" target="_blank">daylemariesmith@gmail.com</a><br></div>
<mailto:<a href="mailto:daylemariesmith@gmail.com" target="_blank">daylemariesmith@gmail.com</a>>>:<div><br>
> Greetings---<br>
> I'm working with a DNA system, and all of the routines I've<br>
worked with that<br>
> require Jacobi diagonalization (g_covar, g_rms, etc) fail with<br>
the "Too many<br>
> iterations in routine JACOBI" error. I'm using gromacs-3.3.3 with<br>
ffamber99<br>
> on the NCSA Mercury cluster. I've searched the archives, and I've<br>
found<br>
> several entries in which users are advised to check that the<br>
coordinates in<br>
> the trajectory and structure files match (mine do). I've also<br>
tried running<br>
> covariance analysis on a small ligand molecule, and I get the<br>
same error. I<br>
> can get g_covar to work with -nofit, but then I can't run g_anaeig.<br>
><br>
> I'm eagerly looking forward to your suggestions!<br>
><br>
> Have a great day,<br>
> Dayle Smith<br>
> Department of Physics<br>
> Whitman College<br>
><br>
<br></div></blockquote></blockquote></div><br>