nmrLDPSolver.h

Go to the documentation of this file.

/* -*- Mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-    */
/* ex: set filetype=cpp softtabstop=4 shiftwidth=4 tabstop=4 cindent expandtab: */

/*
  
  Author(s):    Ankur Kapoor
  Created on:   2004-10-30

  (C) Copyright 2004-2007 Johns Hopkins University (JHU), All Rights
  Reserved.

--- begin cisst license - do not edit ---

This software is provided "as is" under an open source license, with
no warranty.  The complete license can be found in license.txt and
http://www.cisst.org/cisst/license.txt.

--- end cisst license ---
*/


#ifndef _nmrLDPSolver_h
#define _nmrLDPSolver_h

#include <cisstVector/vctDynamicMatrix.h>
#include <cisstNumerical/nmrNetlib.h>
#include <cisstNumerical/nmrExport.h>


class nmrLDPSolver {
    
protected:
    CISSTNETLIB_INTEGER M;
    CISSTNETLIB_INTEGER N;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> E;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> CopyE;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> F;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> R;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> U;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> W;
    vctDynamicMatrix<CISSTNETLIB_INTEGER> Index;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> Zz;
    vctDynamicMatrix<CISSTNETLIB_DOUBLE> X;

public:
    nmrLDPSolver(void):
        M(0),
        N(0)
    {
        Allocate(M, N);
    }


    nmrLDPSolver(CISSTNETLIB_INTEGER m, CISSTNETLIB_INTEGER n) {
        Allocate(m, n);
    }


    nmrLDPSolver(vctDynamicMatrix<CISSTNETLIB_DOUBLE> &G, vctDynamicMatrix<CISSTNETLIB_DOUBLE> &h) {
        Allocate(G, h);
    }


    inline void Allocate(CISSTNETLIB_INTEGER m, CISSTNETLIB_INTEGER n) {
        M = m;
        N = n;
        E.SetSize(N + 1, M, VCT_COL_MAJOR);
        CopyE.SetSize(N + 1, M, VCT_COL_MAJOR);
        F.SetSize(N + 1, 1, VCT_COL_MAJOR);
        R.SetSize(N + 1, 1, VCT_COL_MAJOR);
        U.SetSize(M, 1, VCT_COL_MAJOR);
        W.SetSize(M, 1, VCT_COL_MAJOR);
        Index.SetSize(M, 1, VCT_COL_MAJOR);
        Zz.SetSize(N + 1, 1, VCT_COL_MAJOR);
        X.SetSize(N, 1, VCT_COL_MAJOR);
    }


    inline void Allocate(vctDynamicMatrix<CISSTNETLIB_DOUBLE> &G, vctDynamicMatrix<CISSTNETLIB_DOUBLE> & CMN_UNUSED(h)) {
        Allocate(G.rows(), G.cols());
    }
    

    CISST_EXPORT void Solve(vctDynamicMatrix<CISSTNETLIB_DOUBLE> &G, vctDynamicMatrix<CISSTNETLIB_DOUBLE> &h) throw (std::runtime_error);


    inline const vctDynamicMatrix<CISSTNETLIB_DOUBLE> &GetX(void) const {
        return X;
    }
};


#endif // _nmrLDPSolver_h