ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/mdtools/md_code/Symplectic.cpp
(Generate patch)

Comparing trunk/mdtools/md_code/Symplectic.cpp (file contents):
Revision 11 by mmeineke, Tue Jul 9 18:40:59 2002 UTC vs.
Revision 25 by mmeineke, Tue Jul 9 21:46:32 2002 UTC

# Line 129 | Line 129 | void Symplectic::integrate( void ){
129    double angle;   // the angle through which to rotate the rotation matrix
130    double A[3][3]; // the rotation matrix
131  
132 +  int time;
133  
134    double dt          = entry_plug->dt;
135    double runTime     = entry_plug->run_time;
# Line 170 | Line 171 | void Symplectic::integrate( void ){
171      tStats->velocitize();
172    }
173    
174 +  dump_out->writeDump( 0.0 );
175 +  e_out->writeStat( 0.0 );
176 +
177    if( n_constrained ){
178  
179      double *Rx = new double[nAtoms];
# Line 358 | Line 362 | void Symplectic::integrate( void ){
362            dAtom->setJz( ji[2] );
363          }
364        }
365 <    
366 <      if( entry_plug->setTemp ){
367 <        if( !(tl % vel_n) ) tStats->velocitize();
365 >    
366 >      time = tl + 1;
367 >      
368 >      if( entry_plug->setTemp ){
369 >        if( !(time % vel_n) ) tStats->velocitize();
370        }
371 <      if( !(tl % sample_n) ) dump_out->writeDump( tl * dt );
372 <      if( !(tl % status_n) ) e_out->writeStat( tl * dt );
371 >      if( !(time % sample_n) ) dump_out->writeDump( time * dt );
372 >      if( !(time % status_n) ) e_out->writeStat( time * dt );
373      }
374    }
375    else{
# Line 521 | Line 527 | void Symplectic::integrate( void ){
527            dAtom->setJz( ji[2] );
528          }
529        }
530 <    
530 >      
531 >      time = tl + 1;
532 >      
533        if( entry_plug->setTemp ){
534 <        if( !(tl % vel_n) ) tStats->velocitize();
534 >        if( !(time % vel_n) ) tStats->velocitize();
535        }
536 <      if( !(tl % sample_n) ) dump_out->writeDump( tl * dt );
537 <      if( !(tl % status_n) ) e_out->writeStat( tl * dt );
536 >      if( !(time % sample_n) ) dump_out->writeDump( time * dt );
537 >      if( !(time % status_n) ) e_out->writeStat( time * dt );
538      }
539    }
540  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines