Changes in src/molecule_geometry.cpp [1883f9:4bb63c]

File:

• src/molecule_geometry.cpp (modified) (19 diffs)

src/molecule_geometry.cpp

@@ -48 +48 @@
 
   // go through all atoms
-  BOOST_FOREACH(atom* iter, atoms){
-    *iter -= *Center;
-    *iter -= *CenterBox;
-  }
+  ActOnAllVectors( &Vector::SubtractVector, *Center);
+  ActOnAllVectors( &Vector::SubtractVector, *CenterBox);
   atoms.transformNodes(boost::bind(&Box::WrapPeriodically,domain,_1));
 
@@ -68 +66 @@
   Box &domain = World::getInstance().getDomain();
 
-  // go through all atoms
   atoms.transformNodes(boost::bind(&Box::WrapPeriodically,domain,_1));
 
@@ -86 +83 @@
   if (iter != end()) { //list not empty?
     for (int i=NDIM;i--;) {
-      max->at(i) = (*iter)->at(i);
-      min->at(i) = (*iter)->at(i);
+      max->at(i) = (*iter)->x[i];
+      min->at(i) = (*iter)->x[i];
     }
     for (; iter != end(); ++iter) {// continue with second if present
       //(*iter)->Output(1,1,out);
       for (int i=NDIM;i--;) {
-        max->at(i) = (max->at(i) < (*iter)->at(i)) ? (*iter)->at(i) : max->at(i);
-        min->at(i) = (min->at(i) > (*iter)->at(i)) ? (*iter)->at(i) : min->at(i);
+        max->at(i) = (max->at(i) < (*iter)->x[i]) ? (*iter)->x[i] : max->at(i);
+        min->at(i) = (min->at(i) > (*iter)->x[i]) ? (*iter)->x[i] : min->at(i);
       }
     }
@@ -104 +101 @@
     (*max) += (*min);
     Translate(min);
+    Center.Zero();
   }
   delete(min);
@@ -117 +115 @@
   int Num = 0;
   molecule::const_iterator iter = begin();  // start at first in list
-  Vector Center;
 
   Center.Zero();
+
   if (iter != end()) {   //list not empty?
     for (; iter != end(); ++iter) {  // continue with second if present
       Num++;
-      Center += (*iter)->getPosition();
+      Center += (*iter)->x;
     }
     Center.Scale(-1./(double)Num); // divide through total number (and sign for direction)
     Translate(&Center);
+    Center.Zero();
   }
 };
@@ -144 +143 @@
     for (; iter != end(); ++iter) {  // continue with second if present
       Num++;
-      (*a) += (*iter)->getPosition();
+      (*a) += (*iter)->x;
     }
     a->Scale(1./(double)Num); // divide through total mass (and sign for direction)
@@ -177 +176 @@
   if (iter != end()) {   //list not empty?
     for (; iter != end(); ++iter) {  // continue with second if present
-      Num += (*iter)->getType()->mass;
-      tmp = (*iter)->getType()->mass * (*iter)->getPosition();
+      Num += (*iter)->type->mass;
+      tmp = (*iter)->type->mass * (*iter)->x;
       (*a) += tmp;
     }
@@ -195 +194 @@
 void molecule::CenterPeriodic()
 {
-  Vector NewCenter;
-  DeterminePeriodicCenter(NewCenter);
-  // go through all atoms
-  BOOST_FOREACH(atom* iter, atoms){
-    *iter -= NewCenter;
-  }
+  DeterminePeriodicCenter(Center);
 };
 
@@ -210 +204 @@
 void molecule::CenterAtVector(Vector *newcenter)
 {
-  // go through all atoms
-  BOOST_FOREACH(atom* iter, atoms){
-    *iter -= *newcenter;
-  }
+  Center = *newcenter;
 };
 
@@ -230 +221 @@
     for (int j=0;j<MDSteps;j++)
       (*iter)->Trajectory.R.at(j).ScaleAll(*factor);
-    (*iter)->ScaleAll(*factor);
+    (*iter)->x.ScaleAll(*factor);
   }
 };
@@ -242 +233 @@
     for (int j=0;j<MDSteps;j++)
       (*iter)->Trajectory.R.at(j) += (*trans);
-    *(*iter) += (*trans);
+    (*iter)->x += (*trans);
   }
 };
@@ -255 +246 @@
 
   // go through all atoms
-  BOOST_FOREACH(atom* iter, atoms){
-    *iter += *trans;
-  }
+  ActOnAllVectors( &Vector::AddVector, *trans);
   atoms.transformNodes(boost::bind(&Box::WrapPeriodically,domain,_1));
 
@@ -283 +272 @@
   bool flag;
   Vector Testvector, Translationvector;
-  Vector Center;
 
   do {
@@ -290 +278 @@
     for (molecule::const_iterator iter = begin(); iter != end(); ++iter) {
 #ifdef ADDHYDROGEN
-      if ((*iter)->getType()->Z != 1) {
+      if ((*iter)->type->Z != 1) {
 #endif
-        Testvector = inversematrix * (*iter)->getPosition();
+        Testvector = inversematrix * (*iter)->x;
         Translationvector.Zero();
         for (BondList::const_iterator Runner = (*iter)->ListOfBonds.begin(); Runner != (*iter)->ListOfBonds.end(); (++Runner)) {
          if ((*iter)->nr < (*Runner)->GetOtherAtom((*iter))->nr) // otherwise we shift one to, the other fro and gain nothing
             for (int j=0;j<NDIM;j++) {
-              tmp = (*iter)->at(j) - (*Runner)->GetOtherAtom(*iter)->at(j);
+              tmp = (*iter)->x[j] - (*Runner)->GetOtherAtom(*iter)->x[j];
               if ((fabs(tmp)) > BondDistance) {
                 flag = false;
@@ -315 +303 @@
         // now also change all hydrogens
         for (BondList::const_iterator Runner = (*iter)->ListOfBonds.begin(); Runner != (*iter)->ListOfBonds.end(); (++Runner)) {
-          if ((*Runner)->GetOtherAtom((*iter))->getType()->Z == 1) {
-            Testvector = inversematrix * (*Runner)->GetOtherAtom((*iter))->getPosition();
+          if ((*Runner)->GetOtherAtom((*iter))->type->Z == 1) {
+            Testvector = inversematrix * (*Runner)->GetOtherAtom((*iter))->x;
             Testvector += Translationvector;
             Testvector *= matrix;
@@ -329 +317 @@
 
   Center.Scale(1./static_cast<double>(getAtomCount()));
-  CenterAtVector(&Center);
 };
 
@@ -348 +335 @@
   DoLog(1) && (Log() << Verbose(1) << "Z-X-angle: " << alpha << " ... ");
   for (molecule::const_iterator iter = begin(); iter != end(); ++iter) {
-    tmp = (*iter)->at(0);
-    (*iter)->set(0,  cos(alpha) * tmp + sin(alpha) * (*iter)->at(2));
-    (*iter)->set(2, -sin(alpha) * tmp + cos(alpha) * (*iter)->at(2));
+    tmp = (*iter)->x[0];
+    (*iter)->x[0] =  cos(alpha) * tmp + sin(alpha) * (*iter)->x[2];
+    (*iter)->x[2] = -sin(alpha) * tmp + cos(alpha) * (*iter)->x[2];
     for (int j=0;j<MDSteps;j++) {
       tmp = (*iter)->Trajectory.R.at(j)[0];
@@ -367 +354 @@
   DoLog(1) && (Log() << Verbose(1) << "Z-Y-angle: " << alpha << " ... ");
   for (molecule::const_iterator iter = begin(); iter != end(); ++iter) {
-    tmp = (*iter)->at(1);
-    (*iter)->set(1,  cos(alpha) * tmp + sin(alpha) * (*iter)->at(2));
-    (*iter)->set(2, -sin(alpha) * tmp + cos(alpha) * (*iter)->at(2));
+    tmp = (*iter)->x[1];
+    (*iter)->x[1] =  cos(alpha) * tmp + sin(alpha) * (*iter)->x[2];
+    (*iter)->x[2] = -sin(alpha) * tmp + cos(alpha) * (*iter)->x[2];
     for (int j=0;j<MDSteps;j++) {
       tmp = (*iter)->Trajectory.R.at(j)[1];
@@ -407 +394 @@
   // go through all atoms
   for (molecule::const_iterator iter = par->mol->begin(); iter != par->mol->end(); ++iter) {
-    if ((*iter)->getType() == ((struct lsq_params *)params)->type) { // for specific type
-      c = (*iter)->getPosition() - a;
+    if ((*iter)->type == ((struct lsq_params *)params)->type) { // for specific type
+      c = (*iter)->x - a;
       t = c.ScalarProduct(b);           // get direction parameter
       d = t*b;       // and create vector