superlu.h

Go to the documentation of this file.

#ifndef TBCI_SOLVER_SUPERLU_H
#define TBCI_SOLVER_SUPERLU_H

#include "tbci/vector.h"
#include "tbci/cscmatrix.h"

#if !defined(COLPERM_T_DECLARED) //and !defined(__SUPERLU_UTIL)

typedef enum {NATURAL, MMD_ATA, MMD_AT_PLUS_A, COLAMD, MY_PERMC} colperm_t;  
#define COLPERM_T_DECLARED
#endif

NAMESPACE_TBCI

int lu_solve (CSCMatrix<double>& M, Vector<double>& x, 
              const Vector<double>& rhs, 
              colperm_t permc_spec = NATURAL,
              bool verbose = false, bool symm = false);

inline
TVector<double> lu_solve(CSCMatrix<double>& M, 
                         const Vector<double>& rhs,
                         colperm_t permc_spec = NATURAL,
                         bool verbose = false, bool symm = false)
{
    TVector<double> x(rhs.size());
    const int info = lu_solve(M, (Vector<double>&)x, rhs, permc_spec, verbose, symm);
    return x;
}

#ifdef TBCI_CPLX_H
int lu_solve (CSCMatrix< cplx<double> >& M, Vector< cplx<double> >& x, 
              const Vector< cplx<double> >& rhs, 
              colperm_t permc_spec = NATURAL,
              bool verbose = false, bool symm = false);

inline
TVector< cplx<double> > lu_solve(CSCMatrix< cplx<double> >& M, 
                                 const Vector< cplx<double> >& rhs, 
                                 colperm_t permc_spec = NATURAL,
                                 bool verbose = false, bool symm = false)
{
    TVector< cplx<double> > x(rhs.size());
    const int info = lu_solve(M, (Vector< cplx<double> >&)x,
                              rhs, permc_spec, verbose, symm);
    return x;
}

#endif
NAMESPACE_END

#endif /* TBCI_SOLVER_SUPERLU_H */