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

Comparing trunk/OOPSE/libmdtools/Molecule.cpp (file contents):
Revision 1097 by gezelter, Mon Apr 12 20:32:20 2004 UTC vs.
Revision 1104 by gezelter, Tue Apr 13 16:26:03 2004 UTC

# Line 12 | Line 12 | Molecule::Molecule( void ){
12    myBonds = NULL;
13    myBends = NULL;
14    myTorsions = NULL;
15 –  myRigidBodies = NULL;
15  
16   }
17  
19 –
20 –
18   Molecule::~Molecule( void ){
19    int i;
20    
# Line 41 | Line 38 | Molecule::~Molecule( void ){
38      delete[] myTorsions;
39    }
40  
44 –  if( myRigidBodies != NULL ){
45 –    for(i=0; i<nRigidBodies; i++) if(myRigidBodies[i] != NULL )
46 –      delete myRigidBodies[i];
47 –    delete[] myRigidBodies;
48 –  }
49 –  
41   }
42  
43  
# Line 65 | Line 56 | void Molecule::initialize( molInit &theInit ){
56    myBends = theInit.myBends;
57    myTorsions = theInit.myTorsions;
58    myRigidBodies = theInit.myRigidBodies;
59 +
60 +  myIntegrableObjects = theInit.myIntegrableObjects;
61      
62   }
63  
# Line 72 | Line 65 | void Molecule::calcForces( void ){
65    
66    int i;
67  
68 <  for(i=0; i<nRigidBodies; i++) {
68 >  for(i=0; i<myRigidBodies.size(); i++) {
69      myRigidBodies[i]->updateAtoms();
70    }
71  
# Line 147 | Line 140 | void Molecule::moveCOM(double delta[3]){
140      }
141    }
142  
143 <  for(i=0; i<nRigidBodies; i++) {
143 >  for(i=0; i<myRigidBodies.size(); i++) {
144  
152 –    if (myRigidBodies[i] != NULL) {
153 –      
145        myRigidBodies[i]->getPos( aPos );
146  
147        for (j=0; j< 3; j++)
# Line 158 | Line 149 | void Molecule::moveCOM(double delta[3]){
149        
150        myRigidBodies[i]->setPos( aPos );
151      }
161 –  }  
152   }
153  
154   void Molecule::atoms2rigidBodies( void ) {
155    int i;
156 <  for (i = 0; i < nRigidBodies; i++) {
157 <    if (myRigidBodies[i] != NULL) {
168 <      myRigidBodies[i]->calcForcesAndTorques();  
169 <    }
156 >  for (i = 0; i < myRigidBodies.size(); i++) {
157 >    myRigidBodies[i]->calcForcesAndTorques();  
158    }
159   }
160  

Diff Legend

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