| 45 |
|
#include "types/MAWInteractionType.hpp" |
| 46 |
|
#include "types/LennardJonesInteractionType.hpp" |
| 47 |
|
#include "types/RepulsiveMorseInteractionType.hpp" |
| 48 |
+ |
#include "types/RepulsivePowerInteractionType.hpp" |
| 49 |
|
#include "UseTheForce/ForceField.hpp" |
| 50 |
|
#include "utils/simError.h" |
| 51 |
|
namespace OpenMD { |
| 54 |
|
setSectionName("NonBondedInteractions"); |
| 55 |
|
|
| 56 |
|
stringToEnumMap_["MAW"] = MAW; |
| 57 |
< |
stringToEnumMap_["ShiftedMorse"] = ShiftedMorse; |
| 57 |
> |
stringToEnumMap_["ShiftedMorse"] = ShiftedMorse; |
| 58 |
|
stringToEnumMap_["LennardJones"] = LennardJones; |
| 59 |
|
stringToEnumMap_["RepulsiveMorse"] = RepulsiveMorse; |
| 60 |
+ |
stringToEnumMap_["RepulsivePower"] = RepulsivePower; |
| 61 |
|
|
| 62 |
|
} |
| 63 |
|
|
| 139 |
|
interactionType = new LennardJonesInteractionType(sigma, epsilon); |
| 140 |
|
} |
| 141 |
|
break; |
| 142 |
+ |
|
| 143 |
+ |
case RepulsivePower : |
| 144 |
+ |
if (nTokens < 3) { |
| 145 |
+ |
sprintf(painCave.errMsg, "NonBondedInteractionsSectionParser Error: Not enough tokens at line %d\n", |
| 146 |
+ |
lineNo); |
| 147 |
+ |
painCave.isFatal = 1; |
| 148 |
+ |
simError(); |
| 149 |
+ |
} else { |
| 150 |
+ |
RealType sigma = tokenizer.nextTokenAsDouble(); |
| 151 |
+ |
RealType epsilon = tokenizer.nextTokenAsDouble(); |
| 152 |
+ |
int nRep = tokenizer.nextTokenAsInt(); |
| 153 |
+ |
interactionType = new RepulsivePowerInteractionType(sigma, epsilon, nRep); |
| 154 |
+ |
} |
| 155 |
+ |
break; |
| 156 |
|
|
| 157 |
|
case Unknown : |
| 158 |
|
default: |