Lines Matching refs:glu

47  * \param glu Global LU data. 
53 Index SparseLUImpl<Scalar,Index>::column_bmod(const Index jcol, const Index nseg, BlockScalarVector dense, ScalarVector& tempv, BlockIndexVector segrep, BlockIndexVector repfnz, Index fpanelc, GlobalLU_t& glu)
67 jsupno = glu.supno(jcol);
77 ksupno = glu.supno(krep);
81 fsupc = glu.xsup(ksupno);
88 luptr = glu.xlusup(fst_col) + d_fsupc;
89 lptr = glu.xlsub(fsupc) + d_fsupc;
96 nsupr = glu.xlsub(fsupc+1) - glu.xlsub(fsupc);
98 Index lda = glu.xlusup(fst_col+1) - glu.xlusup(fst_col);
105 LU_kernel_bmod<1>::run(segsize, dense, tempv, glu.lusup, luptr, lda, nrow, glu.lsub, lptr, no_zeros);
107 LU_kernel_bmod<Dynamic>::run(segsize, dense, tempv, glu.lusup, luptr, lda, nrow, glu.lsub, lptr, no_zeros);
112 nextlu = glu.xlusup(jcol);
113 fsupc = glu.xsup(jsupno);
117 new_next = nextlu + glu.xlsub(fsupc + 1) - glu.xlsub(fsupc);
121 while (new_next > glu.nzlumax )
123 mem = memXpand<ScalarVector>(glu.lusup, glu.nzlumax, nextlu, LUSUP, glu.num_expansions);
127 for (isub = glu.xlsub(fsupc); isub < glu.xlsub(fsupc+1); isub++)
129 irow = glu.lsub(isub);
130 glu.lusup(nextlu) = dense(irow);
137 glu.lusup.segment(nextlu,offset).setZero();
140 glu.xlusup(jcol + 1) = nextlu; // close L\U(*,jcol);
156 lptr = glu.xlsub(fsupc) + d_fsupc;
157 luptr = glu.xlusup(fst_col) + d_fsupc;
158 nsupr = glu.xlsub(fsupc+1) - glu.xlsub(fsupc); // leading dimension
163 ufirst = glu.xlusup(jcol) + d_fsupc;
164 Index lda = glu.xlusup(jcol+1) - glu.xlusup(jcol);
165 Map<Matrix<Scalar,Dynamic,Dynamic>, 0, OuterStride<> > A( &(glu.lusup.data()[luptr]), nsupc, nsupc, OuterStride<>(lda) );
166 VectorBlock<ScalarVector> u(glu.lusup, ufirst, nsupc);
169 new (&A) Map<Matrix<Scalar,Dynamic,Dynamic>, 0, OuterStride<> > ( &(glu.lusup.data()[luptr+nsupc]), nrow, nsupc, OuterStride<>(lda) );
170 VectorBlock<ScalarVector> l(glu.lusup, ufirst+nsupc, nrow);