--- trunk/mdtools/md_code/Symplectic.cpp 2003/01/17 21:53:33 237 +++ trunk/mdtools/md_code/Symplectic.cpp 2003/01/17 21:53:36 238 @@ -4,6 +4,7 @@ #include "Integrator.hpp" #include "Thermo.hpp" #include "ReadWrite.hpp" +#include "ForceField.hpp" #include "simError.h" extern "C"{ @@ -34,7 +35,6 @@ Symplectic::Symplectic( SimInfo* the_entry_plug ){ isFirst = 1; srInteractions = entry_plug->sr_interactions; - longRange = entry_plug->longRange; nSRI = entry_plug->n_SRI; // give a little love back to the SimInfo object @@ -142,6 +142,8 @@ void Symplectic::integrate( void ){ int status_n = (int)( statusTime / dt ); int vel_n = (int)( thermalTime / dt ); + ForceFields *ff = entry_plug-> + Thermo *tStats = new Thermo( entry_plug ); StatWriter* e_out = new StatWriter( entry_plug ); @@ -154,6 +156,7 @@ void Symplectic::integrate( void ){ // initialize the forces the before the first step + for(i = 0; i < nAtoms; i++){ atoms[i]->zeroForces(); }