iRoCS Toolbox  1.1.0
lmath.hh
Go to the documentation of this file.
1 /**************************************************************************
2  *
3  * Copyright (C) 2015 Kun Liu, Thorsten Falk
4  *
5  * Image Analysis Lab, University of Freiburg, Germany
6  *
7  * This program is free software; you can redistribute it and/or modify
8  * it under the terms of the GNU General Public License as published by
9  * the Free Software Foundation; either version 3 of the License, or
10  * (at your option) any later version.
11  *
12  * This program is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15  * GNU General Public License for more details.
16  *
17  * You should have received a copy of the GNU General Public License
18  * along with this program; if not, write to the Free Software Foundation,
19  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  *
21  **************************************************************************/
22 
23 /*
24  * lMath.hh
25  *
26  * Created on: Nov 1, 2010
27  * Author: liu
28  */
29 
30 #ifndef LMATH_HH_
31 #define LMATH_HH_
32 
33 #ifdef HAVE_CONFIG_H
34 #include <config.hh>
35 #endif
36 
37 #include <blitz/array.h>
38 
39 #define SORT_eigenValueBy_ascendingAlgebra 0
40 #define SORT_eigenValueBy_descendingAlgebra 1
41 #define SORT_eigenValueBy_descendingABS 2
42 
43 /* Calculation of eigenvalues (D) and eigenvectors (V) of 3x3 matrix A
44  */
46  double A[][3], double V[][3], double d[], int sortMode);
47 
48 template<typename DataT, int Dim>
49 void eig3DSym(
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);
55 
56 #include "lmath.icc"
57 
58 #endif /* LMATH_HH_ */
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)