Re: AMBER: restrt

From: Bill Ross <ross.cgl.ucsf.edu>
Date: Thu, 22 Jan 2004 09:15:58 -0800 (PST)

> The trajectory in mdcrd only contains coordinates. If you need the
> velocities, you need to save them by setting ntvel(?)=(e.g.)500. This will
> trigger saving of velocities to mdvel everything 500 steps. You then need
> to combine the velocities and coordinates to restrt.
>
> HOWEVER, if you start the simulations from the restrt files combined this
> way you will produce really the same trajectory delayed by (e.g.) 500
> steps (they will eventually diverge if the simulations are done on
> parallel machines).

The trajectories will also diverge on a single cpu because
the precision is highly truncated in the mdcrd/mdvel format.
This will happen even with the restrt format - only saving
the info in binary (machine-specific) form would overcome
this.

Bill Ross

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Jan 22 2004 - 17:53:01 PST
Custom Search