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-2018 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) CISST_THROW(std::runtime_error);
 

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