--- trunk/OOPSE/libmdtools/ZConstraint.cpp	2003/08/26 20:12:51	723
+++ trunk/OOPSE/libmdtools/ZConstraint.cpp	2003/08/26 20:37:30	726
@@ -32,11 +32,11 @@ template<typename T> ZConstraint<T>::ZConstraint(SimIn
   if(!data){
     sprintf( painCave.errMsg,
                "ZConstraint Warning: User does not set force substraction policy, "
-               "average force substraction policy is used\n");
+               "PolicyByMass is used\n");
     painCave.isFatal = 0;
     simError();      
 
-    forcePolicy = (ForceSubstractionPolicy*) new PolicyByNumber(this);
+    forcePolicy = (ForceSubstractionPolicy*) new PolicyByMass(this);
   }
   else{
     policy = dynamic_cast<StringData*>(data);
@@ -44,11 +44,11 @@ template<typename T> ZConstraint<T>::ZConstraint(SimIn
     if(!policy){
       sprintf( painCave.errMsg,
                  "ZConstraint Error: Convertion from GenericData to StringData failure, "
-                 "average force substraction policy is used\n");
+                 "PolicyByMass is used\n");
       painCave.isFatal = 0;
       simError();      
 
-      forcePolicy = (ForceSubstractionPolicy*) new PolicyByNumber(this);
+      forcePolicy = (ForceSubstractionPolicy*) new PolicyByMass(this);
     }
     else{
       if(policy->getData() == "BYNUMBER")
@@ -825,10 +825,11 @@ template<typename T> void ZConstraint<T>::doZconstrain
       } 
       totalFZ_local += fz[i];
 
-//       cout << "Fixed Molecule\tindex: " << indexOfZConsMols[i] 
-//              <<"\tcurrent zpos: " << COM[whichDirection] 
-//              << "\tcurrent fz: " <<fz[i] << endl;
+      //cout << "Fixed Molecule\tindex: " << indexOfZConsMols[i] 
+      //      <<"\tcurrent zpos: " << COM[whichDirection] 
+      //       << "\tcurrent fz: " <<fz[i] << endl;
 
+
     }
     
   } 
@@ -854,8 +855,8 @@ template<typename T> void ZConstraint<T>::doZconstrain
       zconsAtoms = zconsMols[i]->getMyAtoms();  
     
       for(int j =0; j < nAtomOfCurZConsMol; j++) {
-      force[whichDirection] = -fz[i]/ nAtomOfCurZConsMol;
-        //force[whichDirection] = - forcePolicy->getZFOfFixedZMols(zconsMols[i], zconsAtoms[j], fz[i]);
+        //force[whichDirection] = -fz[i]/ nAtomOfCurZConsMol;
+        force[whichDirection] = - forcePolicy->getZFOfFixedZMols(zconsMols[i], zconsAtoms[j], fz[i]);
         zconsAtoms[j]->addFrc(force);
       }
     
@@ -892,8 +893,8 @@ template<typename T> void ZConstraint<T>::doZconstrain
      Atom** unconsAtoms = unconsMols[i]->getMyAtoms();
      
      for(int j = 0; j < unconsMols[i]->getNAtoms(); j++){          
-       force[whichDirection] = totalFZ / (totNumOfUnconsAtoms + nMovingZAtoms);
-       //force[whichDirection] = forcePolicy->getZFOfMovingMols(unconsAtoms[j],totalFZ);
+       //force[whichDirection] = totalFZ / (totNumOfUnconsAtoms + nMovingZAtoms);
+       force[whichDirection] = forcePolicy->getZFOfMovingMols(unconsAtoms[j],totalFZ);
        unconsAtoms[j]->addFrc(force);
      }
      
@@ -906,8 +907,8 @@ template<typename T> void ZConstraint<T>::doZconstrain
       Atom** movingZAtoms = zconsMols[i]->getMyAtoms();     
 
       for(int j = 0; j < zconsMols[i]->getNAtoms(); j++){
-        force[whichDirection] = totalFZ / (totNumOfUnconsAtoms + nMovingZAtoms);
-        //force[whichDirection] = forcePolicy->getZFOfMovingMols(movingZAtoms[j],totalFZ);
+        //force[whichDirection] = totalFZ / (totNumOfUnconsAtoms + nMovingZAtoms);
+        force[whichDirection] = forcePolicy->getZFOfMovingMols(movingZAtoms[j],totalFZ);
         movingZAtoms[j]->addFrc(force);
       }
     }
@@ -944,11 +945,11 @@ template<typename T> void ZConstraint<T>::doHarmonic()
       zconsMols[i]->getCOM(COM);
 //       cout << "Moving Molecule\tindex: " << indexOfZConsMols[i] 
 // 	   << "\tcurrent zpos: " << COM[whichDirection] << endl;
+
+      diff = COM[whichDirection] -zPos[i];
     
-    diff = COM[whichDirection] -zPos[i];
-    
       harmonicU = 0.5 * kz[i] * diff * diff;  
-    info->lrPot += harmonicU;
+      info->lrPot += harmonicU;
 
       harmonicF =  - kz[i] * diff;
       totalFZ_local += harmonicF;
@@ -958,8 +959,8 @@ template<typename T> void ZConstraint<T>::doHarmonic()
       Atom** movingZAtoms = zconsMols[i]->getMyAtoms();     
 
        for(int j = 0; j < zconsMols[i]->getNAtoms(); j++){          
-        force[whichDirection] = harmonicF / zconsMols[i]->getNAtoms();
-         //force[whichDirection] = forcePolicy->getHFOfFixedZMols(zconsMols[i], movingZAtoms[j], harmonicF);
+        //force[whichDirection] = harmonicF / zconsMols[i]->getNAtoms();
+        force[whichDirection] = forcePolicy->getHFOfFixedZMols(zconsMols[i], movingZAtoms[j], harmonicF);
          movingZAtoms[j]->addFrc(force);
        }
     }
@@ -982,8 +983,8 @@ template<typename T> void ZConstraint<T>::doHarmonic()
      Atom** unconsAtoms = unconsMols[i]->getMyAtoms();
      
      for(int j = 0; j < unconsMols[i]->getNAtoms(); j++){          
-       force[whichDirection] = - totalFZ /totNumOfUnconsAtoms;
-       //force[whichDirection] = - forcePolicy->getHFOfUnconsMols(unconsAtoms[j], totalFZ);
+       //force[whichDirection] = - totalFZ /totNumOfUnconsAtoms;
+       force[whichDirection] = - forcePolicy->getHFOfUnconsMols(unconsAtoms[j], totalFZ);
        unconsAtoms[j]->addFrc(force);     
      }
   }   
@@ -1023,8 +1024,9 @@ template<typename T> bool ZConstraint<T>::checkZConsSt
 #else
   MPI_Allreduce(&changed_local, &changed, 1, MPI_INT,MPI_SUM, MPI_COMM_WORLD);
 #endif
-  
+
   return (changed > 0);
+
 }
 
 template<typename T> bool ZConstraint<T>::haveFixedZMols(){
@@ -1046,7 +1048,7 @@ template<typename T> bool ZConstraint<T>::haveFixedZMo
   MPI_Allreduce(&havingFixed_local, &havingFixed, 1, MPI_INT,MPI_SUM, MPI_COMM_WORLD);
 #endif
 
-  return havingFixed > 0 ? true : false;
+  return (havingFixed > 0);
 }
 
 
@@ -1072,7 +1074,7 @@ template<typename T> bool ZConstraint<T>::haveMovingZM
   MPI_Allreduce(&havingMoving_local, &havingMoving, 1, MPI_INT,MPI_SUM, MPI_COMM_WORLD);
 #endif
 
-  return havingMoving > 0 ? true : false;
+  return (havingMoving > 0);
   
 }