25 #ifndef WLINEARALGEBRAFUNCTIONS_H
26 #define WLINEARALGEBRAFUNCTIONS_H
32 #include "linearAlgebra/WPosition.h"
34 template<
typename >
class WMatrix;
104 template<
typename T >
116 template<
typename T >
120 template<
typename T >
126 Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > eigenA( A );
127 Eigen::JacobiSVD< Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > > svd( eigenA, Eigen::ComputeFullU | Eigen::ComputeFullV );
133 template<
typename T >
140 computeSVD( input, U, V, Svec );
147 S( i, i ) = ( Svec[ i ] == 0.0 ) ? 0.0 : 1.0 / Svec[ i ];
153 #endif // WLINEARALGEBRAFUNCTIONS_H