| 41 |
|
int d; |
| 42 |
|
}; |
| 43 |
|
|
| 44 |
– |
// typedef struct{ |
| 45 |
– |
// int a, b; |
| 46 |
– |
// } bond_pair; |
| 44 |
|
|
| 48 |
– |
// typedef struct{ |
| 49 |
– |
// int a, b, c; |
| 50 |
– |
// } bend_set; |
| 45 |
|
|
| 52 |
– |
// typedef struct{ |
| 53 |
– |
// int a, b, c, d; |
| 54 |
– |
// } torsion_set; |
| 55 |
– |
|
| 56 |
– |
|
| 46 |
|
class ForceFields{ |
| 47 |
|
|
| 48 |
|
public: |
| 54 |
|
virtual void initializeBonds( bond_pair* the_bonds ) = 0; |
| 55 |
|
virtual void initializeBends( bend_set* the_bends ) = 0; |
| 56 |
|
virtual void initializeTorsions( torsion_set* the_torsions ) = 0; |
| 57 |
+ |
virtual void doForces( void ) = 0; |
| 58 |
|
|
| 59 |
|
protected: |
| 60 |
|
|
| 77 |
|
void initializeBonds( bond_pair* the_bonds ); |
| 78 |
|
void initializeBends( bend_set* the_bends ); |
| 79 |
|
void initializeTorsions( torsion_set* the_torsions ); |
| 80 |
+ |
void doForces( void ) {} |
| 81 |
|
}; |
| 82 |
|
|
| 83 |
|
|
| 91 |
|
void initializeBonds( bond_pair* the_bonds ); |
| 92 |
|
void initializeBends( bend_set* the_bends ); |
| 93 |
|
void initializeTorsions( torsion_set* the_torsions ); |
| 94 |
+ |
void doForces( void ) {} |
| 95 |
|
}; |
| 96 |
|
|
| 97 |
|
class TraPPE_ExFF : public ForceFields{ |
| 104 |
|
void initializeBonds( bond_pair* the_bonds ); |
| 105 |
|
void initializeBends( bend_set* the_bends ); |
| 106 |
|
void initializeTorsions( torsion_set* the_torsions ); |
| 107 |
+ |
void doForces( void ) {} |
| 108 |
|
}; |
| 109 |
|
|
| 110 |
|
class LJ_FF : public ForceFields{ |
| 123 |
|
short int* doPotentialCalc ) ){ |
| 124 |
|
doLJfortran = fortranSub; |
| 125 |
|
} |
| 126 |
+ |
void doForces( void ); |
| 127 |
|
|
| 134 |
– |
|
| 128 |
|
private: |
| 129 |
|
|
| 130 |
|
void fastForward( char* stopText, char* searchOwner ); |
| 135 |
|
double* forceArray, |
| 136 |
|
double* potentialEnergy, |
| 137 |
|
short int* doPotentialCalc ); |
| 138 |
< |
|
| 138 |
> |
void initFortran( void ); |
| 139 |
|
}; |
| 140 |
|
|
| 141 |
|
#endif |