| 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 ); |