Hi, I am trying to understand normal mode analysis in gromacs. I just took one amino acid(Glycine), energy minimized then calculated normal mode eigenvectors as described in Gromacs manual(grompp,mdrun,g_nmeig). Then I calculated the inner product between them using 'g_anaeig'. example: g_anaeig -v eigenvec.trr -v2 eigenvec.trr -inpr I would expect all diagonal elements are 1 and the off-diagonal elements are zero. But I didn't get that. Could anyone explain to me what went wrong?
PS: I did everything in double precision Thanks in advance.. ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ Forschungszentrum Juelich GmbH 52425 Juelich Sitz der Gesellschaft: Juelich Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498 Vorsitzender des Aufsichtsrats: MinDirig Dr. Karl Eugen Huthmacher Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender), Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt, Prof. Dr. Sebastian M. Schmidt ------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------ -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users 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 gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists