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