| 39 |
|
void getFrc( double theF[3] ); |
| 40 |
|
void addFrc( double theF[3] ); |
| 41 |
|
void zeroForces(); |
| 42 |
+ |
|
| 43 |
+ |
virtual bool isLinear() {return is_linear;} |
| 44 |
+ |
virtual int linearAxis() {return linear_axis;} |
| 45 |
|
|
| 46 |
|
double getMass( void ) { return mass; } |
| 47 |
|
|
| 48 |
|
void printAmatIndex( void ); |
| 49 |
< |
void setEulerAngles( double phi, double theta, double psi ); |
| 49 |
> |
void setEuler( double phi, double theta, double psi ); |
| 50 |
|
void getQ( double the_q[4] ); // get the quanternions |
| 51 |
|
void setQ( double the_q[4] ); |
| 52 |
|
|
| 56 |
|
void getJ( double theJ[3] ); |
| 57 |
|
void setJ( double theJ[3] ); |
| 58 |
|
|
| 59 |
+ |
virtual void setType(char* type) {strcpy(rbName, type);} |
| 60 |
+ |
virtual char* getType() { return rbName;} |
| 61 |
+ |
|
| 62 |
|
void getTrq( double theT[3] ); |
| 63 |
|
void addTrq( double theT[3] ); |
| 64 |
|
|
| 86 |
|
// utility routines |
| 87 |
|
|
| 88 |
|
void findCOM( void ); |
| 89 |
< |
void findOrient( void ); |
| 89 |
> |
|
| 90 |
> |
virtual void accept(BaseVisitor* v); |
| 91 |
> |
|
| 92 |
> |
vector<Atom*> getAtoms() { return myAtoms;} |
| 93 |
|
|
| 94 |
|
protected: |
| 95 |
|
|
| 103 |
|
double I[3][3]; // the inertial tensor (body fixed) |
| 104 |
|
double sU[3][3]; // the standard unit vectors (body fixed) |
| 105 |
|
|
| 106 |
+ |
bool is_linear; |
| 107 |
+ |
int linear_axis; |
| 108 |
+ |
const static double momIntTol = 1e-6; |
| 109 |
+ |
|
| 110 |
|
vector<Atom*> myAtoms; // the vector of atoms |
| 111 |
|
vector<vec3> refCoords; |
| 112 |
|
vector<mat3x3> refOrients; |
| 113 |
|
|
| 114 |
< |
bool com_good; |
| 102 |
< |
bool forces_good; |
| 103 |
< |
bool precalc_done; |
| 104 |
< |
bool orient_good; |
| 114 |
> |
char rbName[100]; //it will eventually be converted into string |
| 115 |
|
}; |
| 116 |
|
|
| 117 |
|
#endif |
