26 |
|
#define NPTxyz_ENS 4 |
27 |
|
|
28 |
|
|
29 |
< |
#define FF_DUFF 0 |
30 |
< |
#define FF_LJ 1 |
31 |
< |
#define FF_EAM 2 |
32 |
< |
#define FF_H2O 3 |
29 |
> |
#define FF_DUFF 0 |
30 |
> |
#define FF_LJ 1 |
31 |
> |
#define FF_EAM 2 |
32 |
> |
#define FF_H2O 3 |
33 |
> |
#define FF_SHAPES 4 |
34 |
|
|
35 |
|
using namespace std; |
36 |
|
using namespace oopse; |
715 |
|
} |
716 |
|
else if (!strcasecmp(force_field, "WATER")){ |
717 |
|
ffCase = FF_H2O; |
718 |
+ |
} |
719 |
+ |
else if (!strcasecmp(force_field, "SHAPES")){ |
720 |
+ |
ffCase = FF_SHAPES; |
721 |
|
} |
722 |
|
else{ |
723 |
|
sprintf(painCave.errMsg, "SimSetup Error. Unrecognized force field -> %s\n", |
724 |
|
force_field); |
725 |
< |
painCave.isFatal = 1; |
726 |
< |
simError(); |
725 |
> |
painCave.isFatal = 1; |
726 |
> |
simError(); |
727 |
|
} |
728 |
|
if (globals->haveForceFieldVariant()) { |
729 |
< |
strcpy(forcefield_variant, globals->getForceFieldVariant()); |
729 |
> |
forcefield_variant = globals->getForceFieldVariant(); |
730 |
|
has_forcefield_variant = 1; |
731 |
|
} |
732 |
|
|
1261 |
|
|
1262 |
|
case FF_H2O: |
1263 |
|
the_ff = new WATER(); |
1264 |
+ |
break; |
1265 |
+ |
|
1266 |
+ |
case FF_SHAPES: |
1267 |
+ |
the_ff = new Shapes_FF(); |
1268 |
|
break; |
1269 |
|
|
1270 |
|
default: |