37 #include <blitz/array.h> 39 #define SORT_eigenValueBy_ascendingAlgebra 0 40 #define SORT_eigenValueBy_descendingAlgebra 1 41 #define SORT_eigenValueBy_descendingABS 2 46 double A[][3],
double V[][3],
double d[],
int sortMode);
48 template<
typename DataT,
int Dim>
50 blitz::Array<blitz::TinyVector<DataT,6>,3> &hessianMatrix,
51 blitz::Array<blitz::TinyVector<DataT,Dim>,Dim> &hessianEigValues,
52 blitz::Array<blitz::TinyVector<DataT,Dim>,Dim> &v1,
53 blitz::Array<blitz::TinyVector<DataT,Dim>,Dim> &v2,
54 blitz::Array<blitz::TinyVector<DataT,Dim>,Dim> &v3,
int sortMode);
void eigen_decomposition(double A[][3], double V[][3], double d[], int sortMode)
void eig3DSym(blitz::Array< blitz::TinyVector< DataT, 6 >, 3 > &hessianMatrix, blitz::Array< blitz::TinyVector< DataT, Dim >, Dim > &hessianEigValues, blitz::Array< blitz::TinyVector< DataT, Dim >, Dim > &v1, blitz::Array< blitz::TinyVector< DataT, Dim >, Dim > &v2, blitz::Array< blitz::TinyVector< DataT, Dim >, Dim > &v3, int sortMode)