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

---

MatSvdLapack.h

Go to the documentation of this file.

#ifndef Impala_Core_Matrix_MatSvd_h
#define Impala_Core_Matrix_MatSvd_h

#include "Core/Matrix/MatFunc.h"

namespace Impala
{
namespace Core
{
namespace Matrix
{


template<class ArrayT>
inline void
MatSvd(ArrayT* a, ArrayT*& u, ArrayT*& s, ArrayT*& v)
{
    char job = 'A'; // this specifies that all rows/columns should be computed
    int rows = a->H();
    int cols = a->W();
    double* aData;
    double* uData;
    double* sData;
    double* vData; // mind the transpose!
    int lwork = 1;
    double* workData;
    int info;
    int bla = dgesvd_(&job, &job, &rows, &cols, aData,
                      &rows, sData, uData, &rows, vData,
                      &cols, workData, &lwork, &info);

    //error checking...
    if(info < 0)
    {
        std::cerr << "[MatSvd] illegal value@" << info << ": " << aData[-info-1] << std::endl;
        return;
    }
    if(info < 0)
    {
        std::cerr << "[MatSvd] did not converge" << std::endl;
        return;
    }

    u = new ArrayT(uData);
    s = new ArrayT(sData);
    v = new ArrayT(vData); // mind the transpose!
}


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

#endif

---