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)
MatrixID1 Muscat::Intersect1D(const MatrixID1 &A, const MatrixID1 &B)
MatrixID1 Muscat::Union1D(const MatrixID1 &A, const MatrixID1 &B)