4 |
|
#include <vector> |
5 |
|
#include "forcerCmd.h" |
6 |
|
#include "PDBReader.hpp" |
7 |
+ |
#include "RigidBody.hpp" |
8 |
|
|
9 |
|
#define MK_STR(s) # s |
10 |
|
#define STR_DEFINE(t, s) t = MK_STR(s) |
48 |
|
vector<Atype*> vdwAtypes; |
49 |
|
vector<Atype*>::iterator i; |
50 |
|
Atype* at; |
51 |
+ |
RigidBody* rb; |
52 |
|
vector<VDWAtom*> theAtoms; |
53 |
|
vector<VDWAtom*>::iterator j; |
54 |
|
VDWAtom* atom; |
60 |
|
char temp[200]; |
61 |
|
char readLine[500]; |
62 |
|
FILE *vdwFile, *structureFile; |
63 |
< |
char* ffPath_env = "FORCE_PARAM_PATH"; |
63 |
> |
char* ffPath_env = "VDW_PATH"; |
64 |
|
char* ffPath; |
65 |
|
char* eof_test; |
66 |
|
char* foo; |
129 |
|
if( vdwFile == NULL ){ |
130 |
|
|
131 |
|
printf( "Error opening the vanDerWaals parameter file: %s\n" |
132 |
< |
"Have you tried setting the FORCE_PARAM_PATH environment " |
132 |
> |
"Have you tried setting the VDW_PARAM_DIR environment " |
133 |
|
"variable?\n", |
134 |
|
vdwFileName ); |
135 |
|
exit(-1); |
224 |
|
} |
225 |
|
} |
226 |
|
|
227 |
+ |
rb = new RigidBody(); |
228 |
+ |
for( j = theAtoms.begin(); j != theAtoms.end(); ++j){ |
229 |
+ |
rb->addAtom(*j); |
230 |
+ |
} |
231 |
+ |
|
232 |
|
//GridBuilder gb = new GridBuilder(); |
233 |
|
//gb->findAxesAndOrigin(theAtoms); |
234 |
|
//gb->launchProbe(FF, sigmaGrid, sGrid, epsGrid); |