ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/branches/no-template-branch/OOPSE/libmdtools/Thermo.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/Thermo.cpp (file contents):
Revision 428 by mmeineke, Thu Mar 27 21:07:14 2003 UTC vs.
Revision 445 by gezelter, Thu Apr 3 19:58:24 2003 UTC

# Line 10 | Line 10 | using namespace std;
10   #include "Thermo.hpp"
11   #include "SRI.hpp"
12   #include "Integrator.hpp"
13 + #include "simError.h"
14  
15   #ifdef IS_MPI
16   #define __C
# Line 92 | Line 93 | double Thermo::getPotential(){
93    nSRI = entry_plug->n_SRI;
94  
95    potential_local = 0.0;
96 +  potential = 0.0;
97    potential_local += entry_plug->lrPot;
98  
99    for( el=0; el<entry_plug->n_mol; el++ ){    
100      potential_local += molecules[el].getPotential();
101    }
102  
103 + #ifdef IS_MPI
104 +  /*
105 +  std::cerr << "node " << worldRank << ": before LONG RANGE pot = " << entry_plug->lrPot
106 +            << "; pot_local = " << potential_local
107 +            << "; pot = " << potential << "\n";
108 +  */
109 + #endif
110 +
111    // Get total potential for entire system from MPI.
112   #ifdef IS_MPI
113    MPI::COMM_WORLD.Allreduce(&potential_local,&potential,1,MPI_DOUBLE,MPI_SUM);
# Line 105 | Line 115 | double Thermo::getPotential(){
115    potential = potential_local;
116   #endif // is_mpi
117  
118 + #ifdef IS_MPI
119 +  /*
120 +  std::cerr << "node " << worldRank << ": after pot = " << potential << "\n";
121 +  */
122 + #endif
123 +
124    return potential;
125   }
126  
# Line 138 | Line 154 | double Thermo::getPressure(){
154   }
155  
156   double Thermo::getPressure(){
157 +  // returns pressure in units amu*fs^-2*Ang^-1
158 +  // routine derived via viral theorem description in:
159 +  // Paci, E. and Marchi, M. J.Phys.Chem. 1996, 100, 4314-4322
160  
142 – //  const double conv_Pa_atm = 9.901E-6; // convert Pa -> atm
143 – // const double conv_internal_Pa = 1.661E-7; //convert amu/(fs^2 A) -> Pa
144 – //  const double conv_A_m = 1.0E-10; //convert A -> m
145 –
161    return 0.0;
162   }
163  
# Line 196 | Line 211 | void Thermo::velocitize() {
211  
212      av2 = 2.0 * kebar / atoms[vr]->getMass();
213      vbar = sqrt( av2 );
214 <
214 >
215   //     vbar = sqrt( 8.31451e-7 * temperature / atoms[vr]->getMass() );
216      
217      // picks random velocities from a gaussian distribution

Diff Legend

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