| 55 |
|
stringToEnumMap_["MAW"] = MAW; |
| 56 |
|
stringToEnumMap_["ShiftedMorse"] = ShiftedMorse; |
| 57 |
|
stringToEnumMap_["LennardJones"] = LennardJones; |
| 58 |
< |
stringToEnumMap_["RepulsiveMorse"] = LennardJones; |
| 58 |
> |
stringToEnumMap_["RepulsiveMorse"] = RepulsiveMorse; |
| 59 |
|
|
| 60 |
|
} |
| 61 |
|
|
| 73 |
|
|
| 74 |
|
std::string at1 = tokenizer.nextToken(); |
| 75 |
|
std::string at2 = tokenizer.nextToken(); |
| 76 |
< |
MetalNonMetalInteractionTypeEnum nbit = getMetalNonMetalInteractionTypeEnum(tokenizer.nextToken()); |
| 76 |
> |
std::string itype = tokenizer.nextToken(); |
| 77 |
> |
|
| 78 |
> |
MetalNonMetalInteractionTypeEnum nbit = getMetalNonMetalInteractionTypeEnum(itype); |
| 79 |
|
nTokens -= 3; |
| 80 |
|
NonBondedInteractionType* interactionType; |
| 81 |
< |
|
| 81 |
> |
|
| 82 |
|
//switch is a nightmare to maintain |
| 83 |
|
switch(nbit) { |
| 84 |
|
case MAW : |
| 85 |
< |
if (nTokens < 6) { |
| 85 |
> |
if (nTokens < 4) { |
| 86 |
|
sprintf(painCave.errMsg, "MetalNonMetalInteractionsSectionParser Error: Not enough tokens at line %d\n", |
| 87 |
|
lineNo); |
| 88 |
|
painCave.isFatal = 1; |
| 89 |
|
simError(); |
| 90 |
|
} else { |
| 91 |
< |
RealType r0 = tokenizer.nextTokenAsDouble(); |
| 92 |
< |
RealType D0 = tokenizer.nextTokenAsDouble(); |
| 93 |
< |
RealType beta0 = tokenizer.nextTokenAsDouble(); |
| 92 |
< |
RealType betaH = tokenizer.nextTokenAsDouble(); |
| 93 |
< |
RealType gamma = tokenizer.nextTokenAsDouble(); |
| 91 |
> |
RealType r_e = tokenizer.nextTokenAsDouble(); |
| 92 |
> |
RealType D_e = tokenizer.nextTokenAsDouble(); |
| 93 |
> |
RealType beta = tokenizer.nextTokenAsDouble(); |
| 94 |
|
RealType alpha = tokenizer.nextTokenAsDouble(); |
| 95 |
< |
interactionType = new MAWInteractionType(D0, beta0, r0, betaH, gamma, alpha); |
| 95 |
> |
interactionType = new MAWInteractionType(D_e, beta, r_e, alpha); |
| 96 |
|
} |
| 97 |
|
break; |
| 98 |
|
|