ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/ForceFields.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/ForceFields.cpp (file contents):
Revision 423 by mmeineke, Thu Mar 27 20:12:15 2003 UTC vs.
Revision 428 by mmeineke, Thu Mar 27 21:07:14 2003 UTC

# Line 2 | Line 2
2  
3   #ifdef IS_MPI
4   #include <mpi.h>
5 + #include <mpi++.h>
6   #endif // is_mpi
7  
8  
# Line 20 | Line 21 | void ForceFields::calcRcut( void ){
21  
22    //calc rCut and rList
23  
24 <  entry_plug->rCut = 2.5 bigSigma;
24 >  entry_plug->rCut = 2.5 * bigSigma;
25    if(entry_plug->rCut > (entry_plug->box_x / 2.0))
26      entry_plug->rCut = entry_plug->box_x / 2.0;
27    if(entry_plug->rCut > (entry_plug->box_y / 2.0))
# Line 53 | Line 54 | void ForceFields::doForces( int calcPot, int calcStres
54    }
55  
56    for(i=0; i<entry_plug->n_mol; i++ ){
57 <    entry_plug->molecules[i]->calc_forces();
57 >    entry_plug->molecules[i].calcForces();
58    }
59  
60    frc = Atom::getFrcArray();
# Line 61 | Line 62 | void ForceFields::doForces( int calcPot, int calcStres
62    trq = Atom::getTrqArray();
63    A   = Atom::getAmatArray();
64    u_l = Atom::getUlArray();
64
65    tau = entry_plug->tau;
66 <    
66 >
67 >  
68    isError = 0;
69    entry_plug->lrPot = 0.0;
70  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines