12 |
|
#include "fortranWrappers.hpp" |
13 |
|
|
14 |
|
#ifdef IS_MPI |
15 |
– |
#include <mpi++.h> |
15 |
|
#include "mpiForceField.h" |
16 |
|
#endif // is_mpi |
17 |
|
|
98 |
|
return NULL; |
99 |
|
} |
100 |
|
|
101 |
+ |
void printMe( void ){ |
102 |
+ |
|
103 |
+ |
std::cerr << "LinkedAtype " << name << ": ident = " << ident << "\n"; |
104 |
+ |
if( next != NULL ) next->printMe(); |
105 |
+ |
|
106 |
+ |
} |
107 |
+ |
|
108 |
|
void add( atomStruct &info ){ |
109 |
|
|
110 |
|
// check for duplicates |
146 |
|
info.dipole = dipole; |
147 |
|
info.w0 = w0; |
148 |
|
info.v0 = v0; |
149 |
+ |
info.ident = ident; |
150 |
|
info.last = 0; |
151 |
|
} |
152 |
|
|
741 |
|
recieveFrcStruct( &atomInfo, mpiAtomStructType ); |
742 |
|
} |
743 |
|
} |
744 |
+ |
|
745 |
|
#endif // is_mpi |
746 |
+ |
|
747 |
+ |
|
748 |
|
|
749 |
|
// call new A_types in fortran |
750 |
|
|
852 |
|
"TraPPE_Ex bond structures read successfully." ); |
853 |
|
MPIcheckPoint(); |
854 |
|
|
855 |
< |
currentBondType = headBondType; |
855 |
> |
currentBondType = headBondType->next; |
856 |
|
while( currentBondType != NULL ){ |
857 |
|
currentBondType->duplicate( bondInfo ); |
858 |
|
sendFrcStruct( &bondInfo, mpiBondStructType ); |
867 |
|
|
868 |
|
// listen for node 0 to send out the force params |
869 |
|
|
870 |
< |
MPIcheckPoint(); |
870 |
> |
MPIcheckPoint(); |
871 |
|
|
872 |
|
headBondType = new LinkedBondType; |
873 |
|
recieveFrcStruct( &bondInfo, mpiBondStructType ); |
930 |
|
"TraPPE_Ex bend structures read successfully." ); |
931 |
|
MPIcheckPoint(); |
932 |
|
|
933 |
< |
currentBendType = headBendType; |
933 |
> |
currentBendType = headBendType->next; |
934 |
|
while( currentBendType != NULL ){ |
935 |
|
currentBendType->duplicate( bendInfo ); |
936 |
|
sendFrcStruct( &bendInfo, mpiBendStructType ); |
1010 |
|
"TraPPE_Ex torsion structures read successfully." ); |
1011 |
|
MPIcheckPoint(); |
1012 |
|
|
1013 |
< |
currentTorsionType = headTorsionType; |
1013 |
> |
currentTorsionType = headTorsionType->next; |
1014 |
|
while( currentTorsionType != NULL ){ |
1015 |
|
currentTorsionType->duplicate( torsionInfo ); |
1016 |
|
sendFrcStruct( &torsionInfo, mpiTorsionStructType ); |
1149 |
|
} |
1150 |
|
} |
1151 |
|
} |
1142 |
– |
|
1143 |
– |
|
1144 |
– |
#ifdef IS_MPI |
1145 |
– |
sprintf( checkPointMsg, "TraPPE_Ex atoms initialized succesfully" ); |
1146 |
– |
MPIcheckPoint(); |
1147 |
– |
#endif // is_mpi |
1148 |
– |
|
1152 |
|
} |
1153 |
|
|
1154 |
|
void TraPPE_ExFF::initializeBonds( int nBonds, Bond** bondArray, |
1187 |
|
entry_plug->n_constraints++; |
1188 |
|
} |
1189 |
|
} |
1187 |
– |
|
1188 |
– |
#ifdef IS_MPI |
1189 |
– |
sprintf( checkPointMsg, "TraPPE_Ex bonds initialized succesfully" ); |
1190 |
– |
MPIcheckPoint(); |
1191 |
– |
#endif // is_mpi |
1192 |
– |
|
1190 |
|
} |
1191 |
|
|
1192 |
|
void TraPPE_ExFF::initializeBends( int nBends, Bend** bendArray, |
1268 |
|
} |
1269 |
|
} |
1270 |
|
} |
1274 |
– |
|
1275 |
– |
#ifdef IS_MPI |
1276 |
– |
sprintf( checkPointMsg, "TraPPE_Ex bends initialized succesfully" ); |
1277 |
– |
MPIcheckPoint(); |
1278 |
– |
#endif // is_mpi |
1279 |
– |
|
1271 |
|
} |
1272 |
|
|
1273 |
|
void TraPPE_ExFF::initializeTorsions( int nTorsions, Torsion** torsionArray, |
1315 |
|
torsionArray[i] = cTors; |
1316 |
|
} |
1317 |
|
} |
1327 |
– |
|
1328 |
– |
#ifdef IS_MPI |
1329 |
– |
sprintf( checkPointMsg, "TraPPE_Ex torsions initialized succesfully" ); |
1330 |
– |
MPIcheckPoint(); |
1331 |
– |
#endif // is_mpi |
1332 |
– |
|
1318 |
|
} |
1319 |
|
|
1320 |
|
void TraPPE_ExFF::fastForward( char* stopText, char* searchOwner ){ |