| 92 |
|
* Returns the determinant of this matrix. |
| 93 |
|
* @todo need implementation |
| 94 |
|
*/ |
| 95 |
< |
double determinant() const { |
| 96 |
< |
double det; |
| 95 |
> |
Real determinant() const { |
| 96 |
> |
Real det; |
| 97 |
|
return det; |
| 98 |
|
} |
| 99 |
|
|
| 100 |
|
/** Returns the trace of this matrix. */ |
| 101 |
< |
double trace() const { |
| 102 |
< |
double tmp = 0; |
| 101 |
> |
Real trace() const { |
| 102 |
> |
Real tmp = 0; |
| 103 |
|
|
| 104 |
|
for (unsigned int i = 0; i < Dim ; i++) |
| 105 |
|
tmp += data_[i][i]; |
| 182 |
|
* @param w output eigenvalues |
| 183 |
|
* @param v output eigenvectors |
| 184 |
|
*/ |
| 185 |
< |
bool jacobi(const SquareMatrix<Real, Dim>& a, Vector<Real, Dim>& w, |
| 185 |
> |
bool jacobi(SquareMatrix<Real, Dim>& a, Vector<Real, Dim>& w, |
| 186 |
|
SquareMatrix<Real, Dim>& v); |
| 187 |
|
};//end SquareMatrix |
| 188 |
|
|
| 191 |
|
#define MAX_ROTATIONS 60 |
| 192 |
|
|
| 193 |
|
template<typename Real, int Dim> |
| 194 |
< |
bool SquareMatrix<Real, Dim>::jacobi(const SquareMatrix<Real, Dim>& a, Vector<Real, Dim>& w, |
| 194 |
> |
bool SquareMatrix<Real, Dim>::jacobi(SquareMatrix<Real, Dim>& a, Vector<Real, Dim>& w, |
| 195 |
|
SquareMatrix<Real, Dim>& v) { |
| 196 |
|
const int N = Dim; |
| 197 |
|
int i, j, k, iq, ip; |
| 198 |
< |
double tresh, theta, tau, t, sm, s, h, g, c; |
| 199 |
< |
double tmp; |
| 198 |
> |
Real tresh, theta, tau, t, sm, s, h, g, c; |
| 199 |
> |
Real tmp; |
| 200 |
|
Vector<Real, Dim> b, z; |
| 201 |
|
|
| 202 |
|
// initialize |
