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