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 |
|
|
751 |
|
int isError; |
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 ); |