Linear Algebra Tools¶

template<typename TM, typename TIM, typename S> void Muscat::inv22(TM &m, TIM &minv, S &det)¶

template<typename TM, typename TIM, typename S> void Muscat::inv33(TM &m, TIM &minv, S &det)¶

void Muscat::GetInv_Jacobian_Det(MapMatrixDDD &valdphidxi, const Eigen::Ref<const MatrixDDD> &xcoor, int &Dimensionality, MatrixDDD &Jacobian, double &Jdet, Eigen::ColPivHouseholderQR<MatrixDDD> &Jinv)¶

void Muscat::GetNormal(const int &SpaceDim, const int &elementDim, MatrixDDD &Jacobian, MatrixD31 &Normal)¶

struct EigenSolvers¶

Public Types

typedef Eigen::ConjugateGradient<SpMatD, Upper | Lower> EigenSpCG¶

typedef Eigen::SparseLU<SpMatD> EigenSpLU¶

typedef Eigen::SparseQR<SpMatD, Eigen::COLAMDOrdering<int>> EigenSpQR¶

typedef Eigen::BiCGSTAB<SpMatDR> EigenSpBiCGSTAB¶

Public Functions

inline EigenSolvers()¶

inline void Clean()¶

inline void ForceNumberOfThreads(int n)¶

inline void SetSolverType(int i)¶

inline void SetOp(const CMuscatIndex &sizem, const CMuscatIndex &sizen, const CMuscatIndex &ev_size, CMuscatFloat *ev, CMuscatIndex *ei, CMuscatIndex *ej, const double &tolerance)¶

inline void Solve(CMuscatIndex size, CMuscatFloat *_rhs, CMuscatFloat *_sol)¶

inline CMuscatIndex GetSPQRRank()¶

inline CMuscatIndex GetSPQR_R_nonZeros()¶

inline void GetSPQR_R(CMuscatIndex *sizei, CMuscatIndex *sizej, CMuscatFloat *ev, CMuscatIndex *ei, CMuscatIndex *ej)¶

inline void GetSPQR_P(CMuscatIndex *p)¶

inline CMuscatIndex GetSPQR_Q_nonZeros()¶

inline void GetSPQR_Q(CMuscatIndex *sizei, CMuscatIndex *sizej, CMuscatFloat *ev, CMuscatIndex *ei, CMuscatIndex *ej)¶

Public Members

std::unique_ptr<EigenSpCG> cgSolver¶

std::unique_ptr<EigenSpLU> luSolver¶

std::unique_ptr<EigenSpQR> spqrSolver¶

std::unique_ptr<EigenSpBiCGSTAB> bicgstabSolver¶

std::unique_ptr<SpMatD> A¶

std::unique_ptr<SpMatDR> A_rowmajor¶

int solverType¶

Eigen::SparseMatrix<double> Q¶

Eigen::SparseMatrix<double> R¶

Numpy like operators¶

MatrixID1 Muscat::NonZero(const MatrixBD1 &input)¶

template<typename TypeA, typename TypeB> MatrixID1 Muscat::Intersect1D(const TypeA &A, const TypeB &B)¶

template<typename TypeA, typename TypeB> MatrixID1 Muscat::Union1D(const TypeA &A, const TypeB &B)¶