| 188 |
|
q[3] = (this->data_[0][1] - this->data_[1][0]) * s; |
| 189 |
|
} else { |
| 190 |
|
|
| 191 |
< |
ad1 = fabs( this->data_[0][0] ); |
| 192 |
< |
ad2 = fabs( this->data_[1][1] ); |
| 193 |
< |
ad3 = fabs( this->data_[2][2] ); |
| 191 |
> |
ad1 = this->data_[0][0]; |
| 192 |
> |
ad2 = this->data_[1][1]; |
| 193 |
> |
ad3 = this->data_[2][2]; |
| 194 |
|
|
| 195 |
|
if( ad1 >= ad2 && ad1 >= ad3 ){ |
| 196 |
|
|
| 316 |
|
|
| 317 |
|
m /= det; |
| 318 |
|
return m; |
| 319 |
+ |
} |
| 320 |
+ |
|
| 321 |
+ |
SquareMatrix3<Real> transpose() const{ |
| 322 |
+ |
SquareMatrix3<Real> result; |
| 323 |
+ |
|
| 324 |
+ |
for (unsigned int i = 0; i < 3; i++) |
| 325 |
+ |
for (unsigned int j = 0; j < 3; j++) |
| 326 |
+ |
result(j, i) = this->data_[i][j]; |
| 327 |
+ |
|
| 328 |
+ |
return result; |
| 329 |
|
} |
| 330 |
|
/** |
| 331 |
|
* Extract the eigenvalues and eigenvectors from a 3x3 matrix. |
