superlu_stdcplx.h

Go to the documentation of this file.

#ifndef TBCI_SOLVER_SUPERLU_STDCPLX_H
#define TBCI_SOLVER_SUPERLU_STDCPLX_H

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

#ifndef COLPERM_T_DECLARED

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< CPLX__ complex<double> >& M, 
              Vector< CPLX__ complex<double> >& x,
              const Vector< CPLX__ complex<double> >& rhs, 
              colperm_t permc_spec = NATURAL,
              bool verbose = false, bool symm = false);

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

NAMESPACE_END

#endif /* TBCI_SOLVER_SUPERLU_STDCPLX_H */