| 78 |  | return m; | 
| 79 |  | } | 
| 80 |  |  | 
| 81 | < | /** Retunrs  the inversion of this matrix. */ | 
| 81 | > | /** | 
| 82 | > | * Retunrs  the inversion of this matrix. | 
| 83 | > | * @todo need implementation | 
| 84 | > | */ | 
| 85 |  | SquareMatrix<Real, Dim>  inverse() { | 
| 86 |  | SquareMatrix<Real, Dim> result; | 
| 87 |  |  | 
| 88 |  | return result; | 
| 89 |  | } | 
| 90 |  |  | 
| 91 | < | /** Returns the determinant of this matrix. */ | 
| 91 | > | /** | 
| 92 | > | * Returns the determinant of this matrix. | 
| 93 | > | * @todo need implementation | 
| 94 | > | */ | 
| 95 |  | double determinant() const { | 
| 96 |  | double det; | 
| 97 |  | return det; | 
| 148 |  | return true; | 
| 149 |  | } | 
| 150 |  |  | 
| 151 | + | /** @todo need implementation */ | 
| 152 |  | void diagonalize() { | 
| 153 | < | jacobi(m, eigenValues, ortMat); | 
| 153 | > | //jacobi(m, eigenValues, ortMat); | 
| 154 |  | } | 
| 155 |  |  | 
| 156 |  | /** | 
| 165 |  | SquareMatrix<Real, Dim> ortMat; | 
| 166 |  |  | 
| 167 |  | if ( !isSymmetric()){ | 
| 168 | < | throw(); | 
| 168 | > | //throw(); | 
| 169 |  | } | 
| 170 |  |  | 
| 171 |  | SquareMatrix<Real, Dim> m(*this); | 
| 263 |  |  | 
| 264 |  | for (j=iq+1; j<N; j++) | 
| 265 |  | ROT(a,ip,j,iq,j); | 
| 266 | + |  | 
| 267 |  | for (j=0; j<N; j++) | 
| 268 |  | ROT(v,j,ip,j,iq); | 
| 269 |  | } |