|| || || || || || || ||

---

MatSVDLpp.h

Go to the documentation of this file.

#ifndef Impala_Core_Matrix_MatSVD_LAPACK++_h
#define Impala_Core_Matrix_MatSVD_LAPACK++_h

#include "Core/Matrix/MatFunc.h"
#include "Core/Matrix/MatMul.h"
//#include "lapackpp/lapackpp.h"
//
//
//#include "lapackpp/lafnames.h"
#include "Link/lapackpp/gmd.h"
#include "Link/lapackpp/lavd.h"
#include "Link/lapackpp/lasvd.h"

namespace Impala
{
namespace Core
{
namespace Matrix
{


template<class ArrayT>
inline void
MatSvdLpp(ArrayT* a, ArrayT*& u, ArrayT*& s, ArrayT*& v)
{
    int m = a->H();
    int n = a->W();

    /* 
     * la::mat(double*v, int m, int n, bool row_ordering=false);->row_ordering is
     * true for c type arrays 
     * */
    
    LaGenMatDouble A(a->PB(),m,n,true);
    LaGenMatDouble ATmp(A);
    LaVectorDouble S(n);
    LaGenMatDouble U(m,m);
    LaGenMatDouble VT(n,n);
    
    /*
    * void LaSVD_IP(LaGenMatDouble& A, LaVectorDouble &Sigma,
    *  LaGenMatDouble& U, LaGenMatDouble& VT );
    */
    LaSVD_IP(ATmp,S,U,VT);

    //Copy the Results
    u=MatCreate<Mat>(m,m);
    for(int i=0;i<m;i++)
        for(int j=0;j<m;j++)
            u->SetValue(U(i,j),j,i);

    s=MatCreate<Mat>(m,n);
    for(int i=0;i<m;i++){
        for(int j=0;j<n;j++)
            s->SetValue(0,j,i);
        s->SetValue(S(i),i,i);
    }

    v=MatCreate<Mat>(n,n);
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            v->SetValue(VT(i,j),j,i);

}


} // namespace Matrix
} // namespace Core
} // namespace Impala

#endif

---