28 #ifndef SPARSELU_RELAX_SNODE_H
29 #define SPARSELU_RELAX_SNODE_H
46 template <
typename Scalar,
typename Index>
47 void SparseLUImpl<Scalar,Index>::relax_snode (
const Index n, IndexVector& et,
const Index relax_columns, IndexVector& descendants, IndexVector& relax_end)
52 relax_end.setConstant(emptyIdxLU);
53 descendants.setZero();
54 for (j = 0; j < n; j++)
58 descendants(parent) += descendants(j) + 1;
66 while ( parent != n && descendants(parent) < relax_columns )
72 relax_end(snode_start) = j;
75 while (descendants(j) != 0 && j < n) j++;
Definition: Eigen_Colamd.h:50