ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Revert "change ordering of eigenvalues in jacobi3() function. update tests."

Browse Source 
This reverts commit c4eadd3a59.
This commit is contained in:
Axel Kohlmeyer
2024-04-02 23:02:04 -04:00
parent dfc68fb338
commit 21745538a7
99 changed files with 13271 additions and 13772 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/math_eigen.cpp
View File

@ -44,7 +44,7 @@ int MathEigen::jacobi3(double const mat[3][3], double *eval, double evec[3][3])
// create instance of generic Jacobi class and get eigenvalues and -vectors
Jacobi_v1 ecalc3(3, M, midx);
int ierror = ecalc3.Diagonalize(mat, eval, evec, Jacobi_v1::SORT_INCREASING_EVALS);
int ierror = ecalc3.Diagonalize(mat, eval, evec, Jacobi_v1::SORT_DECREASING_EVALS);
// transpose the evec matrix
@ -67,7 +67,7 @@ int MathEigen::jacobi3(double const *const *mat, double *eval, double **evec)
// create instance of generic Jacobi class and get eigenvalues and -vectors
Jacobi_v2 ecalc3(3, M, midx);
int ierror = ecalc3.Diagonalize(mat, eval, evec, Jacobi_v2::SORT_INCREASING_EVALS);
int ierror = ecalc3.Diagonalize(mat, eval, evec, Jacobi_v2::SORT_DECREASING_EVALS);
// transpose the evec matrix