| 57 |
|
hmat_ = m; |
| 58 |
|
invHmat_ = hmat_.inverse(); |
| 59 |
|
|
| 60 |
+ |
|
| 61 |
|
//prepare fortran Hmat |
| 62 |
|
RealType fortranHmat[9]; |
| 63 |
|
RealType fortranInvHmat[9]; |
| 156 |
|
|
| 157 |
|
} |
| 158 |
|
|
| 159 |
+ |
Vector3d Snapshot::getCOM() { |
| 160 |
+ |
if( !hasCOM_ ) { |
| 161 |
+ |
sprintf( painCave.errMsg, "COM was requested before COM was computed!\n"); |
| 162 |
+ |
painCave.severity = OOPSE_ERROR; |
| 163 |
+ |
simError(); |
| 164 |
+ |
} |
| 165 |
+ |
return COM_; |
| 166 |
+ |
} |
| 167 |
+ |
|
| 168 |
+ |
Vector3d Snapshot::getCOMvel() { |
| 169 |
+ |
if( !hasCOM_ ) { |
| 170 |
+ |
sprintf( painCave.errMsg, "COMvel was requested before COM was computed!\n"); |
| 171 |
+ |
painCave.severity = OOPSE_ERROR; |
| 172 |
+ |
simError(); |
| 173 |
+ |
} |
| 174 |
+ |
return COMvel_; |
| 175 |
+ |
} |
| 176 |
+ |
|
| 177 |
+ |
Vector3d Snapshot::getCOMw() { |
| 178 |
+ |
if( !hasCOM_ ) { |
| 179 |
+ |
sprintf( painCave.errMsg, "COMw was requested before COM was computed!\n"); |
| 180 |
+ |
painCave.severity = OOPSE_ERROR; |
| 181 |
+ |
simError(); |
| 182 |
+ |
} |
| 183 |
+ |
return COMw_; |
| 184 |
+ |
} |
| 185 |
+ |
|
| 186 |
|
} |
| 187 |
|
|