10#define SPARSELU_IMPL_H
19template <
typename Scalar,
typename StorageIndex>
23 typedef Matrix<Scalar,Dynamic,1> ScalarVector;
24 typedef Matrix<StorageIndex,Dynamic,1> IndexVector;
25 typedef Matrix<Scalar,Dynamic,Dynamic,ColMajor> ScalarMatrix;
26 typedef Map<ScalarMatrix, 0, OuterStride<> > MappedMatrixBlock;
27 typedef typename ScalarVector::RealScalar RealScalar;
28 typedef Ref<Matrix<Scalar,Dynamic,1> > BlockScalarVector;
29 typedef Ref<Matrix<StorageIndex,Dynamic,1> > BlockIndexVector;
30 typedef LU_GlobalLU_t<IndexVector, ScalarVector> GlobalLU_t;
31 typedef SparseMatrix<Scalar,ColMajor,StorageIndex> MatrixType;
34 template <
typename VectorType>
37 template <
typename VectorType>
38 Index memXpand(VectorType& vec,
Index& maxlen,
Index nbElts, MemType memtype,
Index& num_expansions);
39 void heap_relax_snode (
const Index n, IndexVector& et,
const Index relax_columns, IndexVector& descendants, IndexVector& relax_end);
40 void relax_snode (
const Index n, IndexVector& et,
const Index relax_columns, IndexVector& descendants, IndexVector& relax_end);
41 Index snode_dfs(
const Index jcol,
const Index kcol,
const MatrixType& mat, IndexVector& xprune, IndexVector& marker, GlobalLU_t& glu);
42 Index snode_bmod (
const Index jcol,
const Index fsupc, ScalarVector& dense, GlobalLU_t& glu);
43 Index pivotL(
const Index jcol,
const RealScalar& diagpivotthresh, IndexVector& perm_r, IndexVector& iperm_c,
Index& pivrow, GlobalLU_t& glu);
44 template <
typename Traits>
45 void dfs_kernel(
const StorageIndex jj, IndexVector& perm_r,
46 Index& nseg, IndexVector& panel_lsub, IndexVector& segrep,
47 Ref<IndexVector> repfnz_col, IndexVector& xprune, Ref<IndexVector> marker, IndexVector& parent,
48 IndexVector& xplore, GlobalLU_t& glu,
Index& nextl_col,
Index krow, Traits& traits);
49 void panel_dfs(
const Index m,
const Index w,
const Index jcol, MatrixType& A, IndexVector& perm_r,
Index& nseg, ScalarVector& dense, IndexVector& panel_lsub, IndexVector& segrep, IndexVector& repfnz, IndexVector& xprune, IndexVector& marker, IndexVector& parent, IndexVector& xplore, GlobalLU_t& glu);
51 void panel_bmod(
const Index m,
const Index w,
const Index jcol,
const Index nseg, ScalarVector& dense, ScalarVector& tempv, IndexVector& segrep, IndexVector& repfnz, GlobalLU_t& glu);
52 Index column_dfs(
const Index m,
const Index jcol, IndexVector& perm_r,
Index maxsuper,
Index& nseg, BlockIndexVector lsub_col, IndexVector& segrep, BlockIndexVector repfnz, IndexVector& xprune, IndexVector& marker, IndexVector& parent, IndexVector& xplore, GlobalLU_t& glu);
53 Index column_bmod(
const Index jcol,
const Index nseg, BlockScalarVector dense, ScalarVector& tempv, BlockIndexVector segrep, BlockIndexVector repfnz,
Index fpanelc, GlobalLU_t& glu);
54 Index copy_to_ucol(
const Index jcol,
const Index nseg, IndexVector& segrep, BlockIndexVector repfnz ,IndexVector& perm_r, BlockScalarVector dense, GlobalLU_t& glu);
55 void pruneL(
const Index jcol,
const IndexVector& perm_r,
const Index pivrow,
const Index nseg,
const IndexVector& segrep, BlockIndexVector repfnz, IndexVector& xprune, GlobalLU_t& glu);
57 void fixupL(
const Index n,
const IndexVector& perm_r, GlobalLU_t& glu);
59 template<
typename ,
typename >
60 friend struct column_dfs_traits;
Namespace containing all symbols from the Eigen library.
Definition: Core:141
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74