--- trunk/src/brains/SimInfo.cpp	2005/11/01 19:14:27	705
+++ trunk/src/brains/SimInfo.cpp	2005/11/11 15:22:11	726
@@ -54,9 +54,12 @@
 #include "primitives/Molecule.hpp"
 #include "UseTheForce/fCutoffPolicy.h"
 #include "UseTheForce/DarkSide/fElectrostaticSummationMethod.h"
+#include "UseTheForce/DarkSide/fElectrostaticScreeningMethod.h"
+#include "UseTheForce/DarkSide/fSwitchingFunctionType.h"
 #include "UseTheForce/doForces_interface.h"
 #include "UseTheForce/DarkSide/electrostatic_interface.h"
 #include "UseTheForce/notifyCutoffs_interface.h"
+#include "UseTheForce/DarkSide/switcheroo_interface.h"
 #include "utils/MemoryUtils.hpp"
 #include "utils/simError.h"
 #include "selection/SelectionManager.hpp"
@@ -471,6 +474,7 @@ namespace oopse {
     int isError = 0;
     
     setupElectrostaticSummationMethod( isError );
+    setupSwitchingFunction();
 
     if(isError){
       sprintf( painCave.errMsg,
@@ -528,12 +532,12 @@ namespace oopse {
     //usePBC and useRF are from simParams
     int usePBC = simParams_->getUsePeriodicBoundaryConditions();
     int useRF;
-    int useDW;
+    int useSF;
     std::string myMethod;
 
     // set the useRF logical
     useRF = 0;
-    useDW = 0;
+    useSF = 0;
 
 
     if (simParams_->haveElectrostaticSummationMethod()) {
@@ -542,8 +546,8 @@ namespace oopse {
       if (myMethod == "REACTION_FIELD") {
         useRF=1;
       } else {
-	if (myMethod == "DAMPED_WOLF") {
-	  useDW = 1;
+	if (myMethod == "SHIFTED_FORCE") {
+	  useSF = 1;
 	}
       }
     }
@@ -612,8 +616,8 @@ namespace oopse {
     temp = useRF;
     MPI_Allreduce(&temp, &useRF, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);    
 
-    temp = useDW;
-    MPI_Allreduce(&temp, &useDW, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);    
+    temp = useSF;
+    MPI_Allreduce(&temp, &useSF, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);    
 
 #endif
 
@@ -630,7 +634,7 @@ namespace oopse {
     fInfo_.SIM_uses_Shapes = useShape;
     fInfo_.SIM_uses_FLARB = useFLARB;
     fInfo_.SIM_uses_RF = useRF;
-    fInfo_.SIM_uses_DampedWolf = useDW;
+    fInfo_.SIM_uses_SF = useSF;
 
     if( myMethod == "REACTION_FIELD") {
       
@@ -906,6 +910,7 @@ namespace oopse {
      
     int errorOut;
     int esm =  NONE;
+    int sm = UNDAMPED;
     double alphaVal;
     double dielectric;
 
@@ -919,39 +924,89 @@ namespace oopse {
       if (myMethod == "NONE") {
         esm = NONE;
       } else {
-        if (myMethod == "UNDAMPED_WOLF") {
-          esm = UNDAMPED_WOLF;
+        if (myMethod == "SWITCHING_FUNCTION") {
+          esm = SWITCHING_FUNCTION;
         } else {
-          if (myMethod == "DAMPED_WOLF") {            
-            esm = DAMPED_WOLF;
-	    if (!simParams_->haveDampingAlpha()) {
-	      //throw error
-	      sprintf( painCave.errMsg,
-		       "SimInfo warning: dampingAlpha was not specified in the input file. A default value of %f (1/ang) will be used for the Damped Wolf Method.", alphaVal);
-	      painCave.isFatal = 0;
-	      simError();
-	    }
-          } else {
-	    if (myMethod == "REACTION_FIELD") {	      
-	      esm = REACTION_FIELD;
+	  if (myMethod == "SHIFTED_POTENTIAL") {
+	    esm = SHIFTED_POTENTIAL;
+	  } else {
+	    if (myMethod == "SHIFTED_FORCE") {            
+	      esm = SHIFTED_FORCE;
 	    } else {
-	      // throw error        
-	      sprintf( painCave.errMsg,
-		       "SimInfo error: Unknown electrostaticSummationMethod. (Input file specified %s .)\n\telectrostaticSummationMethod must be one of: \"none\", \"undamped_wolf\", \"damped_wolf\", or \"reaction_field\".", myMethod.c_str() );
-	      painCave.isFatal = 1;
-	      simError();
-	    }     
-	  }           
+	      if (myMethod == "REACTION_FIELD") {	      
+		esm = REACTION_FIELD;
+	      } else {
+		// throw error        
+		sprintf( painCave.errMsg,
+			 "SimInfo error: Unknown electrostaticSummationMethod. (Input file specified %s .)\n\telectrostaticSummationMethod must be one of: \"none\", \"shifted_potential\", \"shifted_force\", or \"reaction_field\".", myMethod.c_str() );
+		painCave.isFatal = 1;
+		simError();
+	      }     
+	    }           
+	  }
 	}
       }
     }
+    
+    if (simParams_->haveElectrostaticScreeningMethod()) {
+      std::string myScreen = simParams_->getElectrostaticScreeningMethod();
+      toUpper(myScreen);
+      if (myScreen == "UNDAMPED") {
+	sm = UNDAMPED;
+      } else {
+	if (myScreen == "DAMPED") {
+	  sm = DAMPED;
+	  if (!simParams_->haveDampingAlpha()) {
+	    //throw error
+	    sprintf( painCave.errMsg,
+		     "SimInfo warning: dampingAlpha was not specified in the input file. A default value of %f (1/ang) will be used.", alphaVal);
+	    painCave.isFatal = 0;
+	    simError();
+	  }
+	} else {
+	  // throw error        
+	  sprintf( painCave.errMsg,
+		   "SimInfo error: Unknown electrostaticScreeningMethod. (Input file specified %s .)\n\telectrostaticScreeningMethod must be one of: \"undamped\" or \"damped\".", myScreen.c_str() );
+	  painCave.isFatal = 1;
+	  simError();
+	}
+      }
+    }
+    
     // let's pass some summation method variables to fortran
     setElectrostaticSummationMethod( &esm );
-    setDampedWolfAlpha( &alphaVal );
+    setScreeningMethod( &sm );
+    setDampingAlpha( &alphaVal );
     setReactionFieldDielectric( &dielectric );
     initFortranFF( &esm, &errorOut );
   }
 
+  void SimInfo::setupSwitchingFunction() {    
+    int ft = CUBIC;
+
+    if (simParams_->haveSwitchingFunctionType()) {
+      std::string funcType = simParams_->getSwitchingFunctionType();
+      toUpper(funcType);
+      if (funcType == "CUBIC") {
+        ft = CUBIC;
+      } else {
+        if (funcType == "FIFTH_ORDER_POLYNOMIAL") {
+          ft = FIFTH_ORDER_POLY;
+	} else {
+	  // throw error        
+	  sprintf( painCave.errMsg,
+		   "SimInfo error: Unknown switchingFunctionType. (Input file specified %s .)\n\tswitchingFunctionType must be one of: \"cubic\" or \"fifth_order_polynomial\".", funcType.c_str() );
+	  painCave.isFatal = 1;
+	  simError();
+        }           
+      }
+    }
+
+    // send switching function notification to switcheroo
+    setFunctionType(&ft);
+
+  }
+
   void SimInfo::addProperty(GenericData* genData) {
     properties_.addProperty(genData);  
   }