Hi Steven,
> Thanks for your reply. I just thought I'd point out that a while ago I
> tried obtaining the forces just by using the debugf options, but they
> don't play nicely with QM/MM:
>
> http://amber.ch.ic.ac.uk/archive/200702/0386.html
Ah ha, now I remember. However, I think that was referring to the dumpfrc
option rather than the actual debugf option - that is testing against
numerical derivatives definitely works (although perhaps only in serial)
since I have used it several times to check that the analytical gradients
are coming out correctly so what I suggested about trying do_debugf=1 and
atomn=1,2,3... should work as a way for you to access the force array via
two different methods, with your write statement and with the info the
debugf prints to the output file, so that you can check that they match
okay.
If you are having problems here please send me your input files and I can
try it out myself.
All the best
Ross
/\
\/
|\oss Walker
| Assistant Research Professor |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
|
http://www.rosswalker.co.uk | PGP Key available on request |
Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Wed Jun 25 2008 - 06:07:50 PDT
