ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/brains/SimInfo.cpp
(Generate patch)

Comparing trunk/src/brains/SimInfo.cpp (file contents):
Revision 682 by chrisfen, Tue Oct 18 15:01:42 2005 UTC vs.
Revision 734 by chuckv, Tue Nov 15 16:05:38 2005 UTC

# Line 54 | Line 54
54   #include "primitives/Molecule.hpp"
55   #include "UseTheForce/fCutoffPolicy.h"
56   #include "UseTheForce/DarkSide/fElectrostaticSummationMethod.h"
57 + #include "UseTheForce/DarkSide/fElectrostaticScreeningMethod.h"
58 + #include "UseTheForce/DarkSide/fSwitchingFunctionType.h"
59   #include "UseTheForce/doForces_interface.h"
60   #include "UseTheForce/DarkSide/electrostatic_interface.h"
61   #include "UseTheForce/notifyCutoffs_interface.h"
62 + #include "UseTheForce/DarkSide/switcheroo_interface.h"
63   #include "utils/MemoryUtils.hpp"
64   #include "utils/simError.h"
65   #include "selection/SelectionManager.hpp"
# Line 471 | Line 474 | namespace oopse {
474      int isError = 0;
475      
476      setupElectrostaticSummationMethod( isError );
477 +    setupSwitchingFunction();
478  
479      if(isError){
480        sprintf( painCave.errMsg,
# Line 515 | Line 519 | namespace oopse {
519      int useLennardJones = 0;
520      int useElectrostatic = 0;
521      int useEAM = 0;
522 +    int useSC = 0;
523      int useCharge = 0;
524      int useDirectional = 0;
525      int useDipole = 0;
# Line 528 | Line 533 | namespace oopse {
533      //usePBC and useRF are from simParams
534      int usePBC = simParams_->getUsePeriodicBoundaryConditions();
535      int useRF;
536 +    int useSF;
537      std::string myMethod;
538  
539      // set the useRF logical
540      useRF = 0;
541 +    useSF = 0;
542 +
543 +
544      if (simParams_->haveElectrostaticSummationMethod()) {
545 <        myMethod = simParams_->getElectrostaticSummationMethod();
546 <        if (myMethod == "REACTION_FIELD")
547 <             useRF = 1;
545 >      std::string myMethod = simParams_->getElectrostaticSummationMethod();
546 >      toUpper(myMethod);
547 >      if (myMethod == "REACTION_FIELD") {
548 >        useRF=1;
549 >      } else {
550 >        if (myMethod == "SHIFTED_FORCE") {
551 >          useSF = 1;
552 >        }
553 >      }
554      }
555  
556      //loop over all of the atom types
# Line 543 | Line 558 | namespace oopse {
558        useLennardJones |= (*i)->isLennardJones();
559        useElectrostatic |= (*i)->isElectrostatic();
560        useEAM |= (*i)->isEAM();
561 +      useSC |= (*i)->isSC();
562        useCharge |= (*i)->isCharge();
563        useDirectional |= (*i)->isDirectional();
564        useDipole |= (*i)->isDipole();
# Line 593 | Line 609 | namespace oopse {
609      temp = useEAM;
610      MPI_Allreduce(&temp, &useEAM, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);    
611  
612 +    temp = useSC;
613 +    MPI_Allreduce(&temp, &useSC, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);
614 +    
615      temp = useShape;
616      MPI_Allreduce(&temp, &useShape, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);  
617  
# Line 602 | Line 621 | namespace oopse {
621      temp = useRF;
622      MPI_Allreduce(&temp, &useRF, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);    
623  
624 +    temp = useSF;
625 +    MPI_Allreduce(&temp, &useSF, 1, MPI_INT, MPI_LOR, MPI_COMM_WORLD);    
626 +
627   #endif
628  
629      fInfo_.SIM_uses_PBC = usePBC;    
# Line 614 | Line 636 | namespace oopse {
636      fInfo_.SIM_uses_StickyPower = useStickyPower;
637      fInfo_.SIM_uses_GayBerne = useGayBerne;
638      fInfo_.SIM_uses_EAM = useEAM;
639 +    fInfo_.SIM_uses_SC = useSC;
640      fInfo_.SIM_uses_Shapes = useShape;
641      fInfo_.SIM_uses_FLARB = useFLARB;
642      fInfo_.SIM_uses_RF = useRF;
643 +    fInfo_.SIM_uses_SF = useSF;
644  
645 <    if( fInfo_.SIM_uses_Dipoles && myMethod == "REACTION_FIELD") {
646 <
645 >    if( myMethod == "REACTION_FIELD") {
646 >      
647        if (simParams_->haveDielectric()) {
648          fInfo_.dielect = simParams_->getDielectric();
649        } else {
# Line 629 | Line 653 | namespace oopse {
653                  "\tsetting a dielectric constant!\n");
654          painCave.isFatal = 1;
655          simError();
656 <      }
633 <        
634 <    } else {
635 <      fInfo_.dielect = 0.0;
656 >      }      
657      }
658  
659    }
# Line 824 | Line 845 | namespace oopse {
845          sprintf(painCave.errMsg,
846                  "SimCreator Warning: No value was set for switchingRadius.\n"
847                  "\tOOPSE will use a default value of\n"
848 <                "\t0.95 * cutoffRadius for the switchingRadius\n");
848 >                "\t0.85 * cutoffRadius for the switchingRadius\n");
849          painCave.isFatal = 0;
850          simError();
851 <        rsw = 0.95 * rcut;
851 >        rsw = 0.85 * rcut;
852        } else{
853          rsw = simParams_->getSwitchingRadius();
854        }
# Line 895 | Line 916 | namespace oopse {
916      
917      int errorOut;
918      int esm =  NONE;
919 +    int sm = UNDAMPED;
920      double alphaVal;
921      double dielectric;
922  
# Line 908 | Line 930 | namespace oopse {
930        if (myMethod == "NONE") {
931          esm = NONE;
932        } else {
933 <        if (myMethod == "UNDAMPED_WOLF") {
934 <          esm = UNDAMPED_WOLF;
933 >        if (myMethod == "SWITCHING_FUNCTION") {
934 >          esm = SWITCHING_FUNCTION;
935          } else {
936 <          if (myMethod == "DAMPED_WOLF") {            
937 <            esm = DAMPED_WOLF;
938 <            if (!simParams_->haveDampingAlpha()) {
939 <              //throw error
940 <              sprintf( painCave.errMsg,
919 <                       "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);
920 <              painCave.isFatal = 0;
921 <              simError();
922 <            }
923 <          } else {
924 <            if (myMethod == "REACTION_FIELD") {
925 <              esm = REACTION_FIELD;
936 >          if (myMethod == "SHIFTED_POTENTIAL") {
937 >            esm = SHIFTED_POTENTIAL;
938 >          } else {
939 >            if (myMethod == "SHIFTED_FORCE") {            
940 >              esm = SHIFTED_FORCE;
941              } else {
942 <              // throw error        
943 <              sprintf( painCave.errMsg,
944 <                       "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() );
945 <              painCave.isFatal = 1;
946 <              simError();
947 <            }    
948 <          }          
942 >              if (myMethod == "REACTION_FIELD") {            
943 >                esm = REACTION_FIELD;
944 >              } else {
945 >                // throw error        
946 >                sprintf( painCave.errMsg,
947 >                         "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() );
948 >                painCave.isFatal = 1;
949 >                simError();
950 >              }    
951 >            }          
952 >          }
953 >        }
954 >      }
955 >    }
956 >    
957 >    if (simParams_->haveElectrostaticScreeningMethod()) {
958 >      std::string myScreen = simParams_->getElectrostaticScreeningMethod();
959 >      toUpper(myScreen);
960 >      if (myScreen == "UNDAMPED") {
961 >        sm = UNDAMPED;
962 >      } else {
963 >        if (myScreen == "DAMPED") {
964 >          sm = DAMPED;
965 >          if (!simParams_->haveDampingAlpha()) {
966 >            //throw error
967 >            sprintf( painCave.errMsg,
968 >                     "SimInfo warning: dampingAlpha was not specified in the input file. A default value of %f (1/ang) will be used.", alphaVal);
969 >            painCave.isFatal = 0;
970 >            simError();
971 >          }
972 >        } else {
973 >          // throw error        
974 >          sprintf( painCave.errMsg,
975 >                   "SimInfo error: Unknown electrostaticScreeningMethod. (Input file specified %s .)\n\telectrostaticScreeningMethod must be one of: \"undamped\" or \"damped\".", myScreen.c_str() );
976 >          painCave.isFatal = 1;
977 >          simError();
978          }
979        }
980      }
981 +    
982      // let's pass some summation method variables to fortran
983      setElectrostaticSummationMethod( &esm );
984 <    setDampedWolfAlpha( &alphaVal );
984 >    setScreeningMethod( &sm );
985 >    setDampingAlpha( &alphaVal );
986      setReactionFieldDielectric( &dielectric );
987      initFortranFF( &esm, &errorOut );
988    }
989  
990 +  void SimInfo::setupSwitchingFunction() {    
991 +    int ft = CUBIC;
992 +
993 +    if (simParams_->haveSwitchingFunctionType()) {
994 +      std::string funcType = simParams_->getSwitchingFunctionType();
995 +      toUpper(funcType);
996 +      if (funcType == "CUBIC") {
997 +        ft = CUBIC;
998 +      } else {
999 +        if (funcType == "FIFTH_ORDER_POLYNOMIAL") {
1000 +          ft = FIFTH_ORDER_POLY;
1001 +        } else {
1002 +          // throw error        
1003 +          sprintf( painCave.errMsg,
1004 +                   "SimInfo error: Unknown switchingFunctionType. (Input file specified %s .)\n\tswitchingFunctionType must be one of: \"cubic\" or \"fifth_order_polynomial\".", funcType.c_str() );
1005 +          painCave.isFatal = 1;
1006 +          simError();
1007 +        }          
1008 +      }
1009 +    }
1010 +
1011 +    // send switching function notification to switcheroo
1012 +    setFunctionType(&ft);
1013 +
1014 +  }
1015 +
1016    void SimInfo::addProperty(GenericData* genData) {
1017      properties_.addProperty(genData);  
1018    }

Diff Legend

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