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

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

Browse Source
This commit is contained in:
Axel Kohlmeyer
2024-03-26 01:52:07 -04:00
parent 45fbe73f2b
commit c4eadd3a59
99 changed files with 13772 additions and 13271 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_DECREASING_EVALS);
int ierror = ecalc3.Diagonalize(mat, eval, evec, Jacobi_v1::SORT_INCREASING_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_DECREASING_EVALS);
int ierror = ecalc3.Diagonalize(mat, eval, evec, Jacobi_v2::SORT_INCREASING_EVALS);
// transpose the evec matrix