--- trunk/src/io/ForceFieldOptions.cpp 2005/12/23 23:04:25 820 +++ trunk/src/io/ForceFieldOptions.cpp 2006/07/16 02:48:00 1005 @@ -40,12 +40,12 @@ */ #define __C #include "io/ForceFieldOptions.hpp" -#include "UseTheForce/fForceOptions.h" -#include "UseTheForce/DarkSide/fForceOptions_interface.h" + namespace oopse { ForceFieldOptions::ForceFieldOptions() { DefineOptionalParameter(Name, "Name"); + DefineOptionalParameter(CutoffPolicy, "cutoffPolicy"); DefineOptionalParameterWithDefaultValue(vdWtype, "vdWtype", "Lennard-Jones"); DefineOptionalParameterWithDefaultValue(DistanceMixingRule, "DistanceMixingRule", "arithmetic"); DefineOptionalParameterWithDefaultValue(DistanceType, "DistanceType", "sigma"); @@ -53,20 +53,47 @@ namespace oopse { DefineOptionalParameterWithDefaultValue(EnergyUnitScaling, "EnergyUnitScaling", 1.0); DefineOptionalParameterWithDefaultValue(DistanceUnitScaling, "DistanceUnitScaling", 1.0); DefineOptionalParameterWithDefaultValue(AngleUnitScaling, "AngleUnitScaling", 1.0); - DefineOptionalParameterWithDefaultValue(TorsionAngleConvention, "TorsionAngleConvention", "180 is trans"); + DefineOptionalParameterWithDefaultValue(TorsionAngleConvention, "TorsionAngleConvention", "180_is_trans"); DefineOptionalParameterWithDefaultValue(vdw14scale, "vdW-14-scale", 0.0); - DefineOptionalParameterWithDefaultValue(electrostatic14scale, "electrostatic-14-scale", 0.0); + DefineOptionalParameterWithDefaultValue(electrostatic14scale, "electrostatic-14-scale", 0.0) DefineOptionalParameterWithDefaultValue(dielectric, "dielectric", 1.0); + DefineOptionalParameterWithDefaultValue(GayBerneMu, "GayBerneMu", 2.0); + DefineOptionalParameterWithDefaultValue(GayBerneNu, "GayBerneNu", 1.0); } - void ForceFieldOptions::tellFortranOptions(){ - ForceOptions fortranForceOptions; + void ForceFieldOptions::makeFortranOptions(ForceOptions& fortranForceOptions){ fortranForceOptions.vdw14scale = this->getvdw14scale(); - fortranForceOptions.vdw14scale = this->getvdw14scale(); - + fortranForceOptions.electrostatic14scale = this->getelectrostatic14scale(); + fortranForceOptions.GayBerneMu = this->getGayBerneMu(); + fortranForceOptions.GayBerneNu = this->getGayBerneNu(); + + std::string DistanceMix = this->getDistanceMixingRule(); + toUpper(DistanceMix); + if(DistanceMix == "ARITHMETIC"){ + fortranForceOptions.DistanceMixingRule = ARITHMETIC_MIXING_RULE; + }else if(DistanceMix == "GEOMETRIC"){ + fortranForceOptions.DistanceMixingRule = GEOMETRIC_MIXING_RULE; + }else if(DistanceMix == "CUBIC"){ + fortranForceOptions.DistanceMixingRule = CUBIC_MIXING_RULE; + }else{ + fortranForceOptions.DistanceMixingRule = ARITHMETIC_MIXING_RULE; + } + std::string EnergyMix = this->getEnergyMixingRule(); + toUpper(EnergyMix); + if(EnergyMix == "ARITHMETIC"){ + fortranForceOptions.EnergyMixingRule = ARITHMETIC_MIXING_RULE; + }else if(EnergyMix == "GEOMETRIC"){ + fortranForceOptions.EnergyMixingRule = GEOMETRIC_MIXING_RULE; + }else if(EnergyMix == "HHG"){ + fortranForceOptions.EnergyMixingRule = HHG_MIXING_RULE; + }else{ + fortranForceOptions.EnergyMixingRule = GEOMETRIC_MIXING_RULE; + } + + } }