11 #ifndef EIGEN_SPARSELU_UTILS_H 12 #define EIGEN_SPARSELU_UTILS_H 20 template <
typename Scalar,
typename Index>
24 nnzU = (glu.
xusub)(n);
25 Index nsuper = (glu.
supno)(n);
30 for (i = 0; i <= nsuper; i++)
35 for (j = fsupc; j < glu.
xsup(i+1); j++)
38 nnzU += j - fsupc + 1;
51 template <
typename Scalar,
typename Index>
54 Index fsupc, i, j, k, jstart;
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++)
80 #endif // EIGEN_SPARSELU_UTILS_H
iterative scaling algorithm to equilibrate rows and column norms in matrices
void fixupL(const Index n, const IndexVector &perm_r, GlobalLU_t &glu)
Fix up the data storage lsub for L-subscripts.
void countnz(const Index n, Index &nnzL, Index &nnzU, GlobalLU_t &glu)
Count Nonzero elements in the factors.