11#ifndef EIGEN_SPARSELU_UTILS_H
12#define EIGEN_SPARSELU_UTILS_H
20template <
typename Scalar,
typename StorageIndex>
21void SparseLUImpl<Scalar,StorageIndex>::countnz(
const Index n,
Index& nnzL,
Index& nnzU, GlobalLU_t& glu)
24 nnzU = (glu.xusub)(n);
25 Index nsuper = (glu.supno)(n);
30 for (i = 0; i <= nsuper; i++)
33 jlen = glu.xlsub(fsupc+1) - glu.xlsub(fsupc);
35 for (j = fsupc; j < glu.xsup(i+1); j++)
38 nnzU += j - fsupc + 1;
51template <
typename Scalar,
typename StorageIndex>
52void SparseLUImpl<Scalar,StorageIndex>::fixupL(
const Index n,
const IndexVector& perm_r, GlobalLU_t& glu)
54 Index fsupc, i, j, k, jstart;
56 StorageIndex nextl = 0;
57 Index nsuper = (glu.supno)(n);
60 for (i = 0; i <= nsuper; i++)
63 jstart = glu.xlsub(fsupc);
64 glu.xlsub(fsupc) = nextl;
65 for (j = jstart; j < glu.xlsub(fsupc + 1); j++)
67 glu.lsub(nextl) = perm_r(glu.lsub(j));
70 for (k = fsupc+1; k < glu.xsup(i+1); k++)
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