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