54 …                  BlockIndexVector segrep, BlockIndexVector repfnz, Index fpanelc, GlobalLU_t& glu)  in column_bmod()  argument
68   jsupno = glu.supno(jcol);  in column_bmod()
78     ksupno = glu.supno(krep);   in column_bmod()
82       fsupc = glu.xsup(ksupno);   in column_bmod()
89       luptr = glu.xlusup(fst_col) + d_fsupc;   in column_bmod()
90       lptr = glu.xlsub(fsupc) + d_fsupc;   in column_bmod()
97       nsupr = glu.xlsub(fsupc+1) - glu.xlsub(fsupc);   in column_bmod()
99       Index lda = glu.xlusup(fst_col+1) - glu.xlusup(fst_col);  in column_bmod()
106 …LU_kernel_bmod<1>::run(segsize, dense, tempv, glu.lusup, luptr, lda, nrow, glu.lsub, lptr, no_zero…  in column_bmod()
108 …LU_kernel_bmod<Dynamic>::run(segsize, dense, tempv, glu.lusup, luptr, lda, nrow, glu.lsub, lptr, n…  in column_bmod()
113   nextlu = glu.xlusup(jcol);   in column_bmod()
114   fsupc = glu.xsup(jsupno);  in column_bmod()
118   new_next = nextlu + glu.xlsub(fsupc + 1) - glu.xlsub(fsupc);   in column_bmod()
122   while (new_next > glu.nzlumax )  in column_bmod()
124     mem = memXpand<ScalarVector>(glu.lusup, glu.nzlumax, nextlu, LUSUP, glu.num_expansions);    in column_bmod()
128   for (isub = glu.xlsub(fsupc); isub < glu.xlsub(fsupc+1); isub++)  in column_bmod()
130     irow = glu.lsub(isub);  in column_bmod()
131     glu.lusup(nextlu) = dense(irow);  in column_bmod()
138     glu.lusup.segment(nextlu,offset).setZero();  in column_bmod()
141   glu.xlusup(jcol + 1) = StorageIndex(nextlu);  // close L\U(*,jcol);   in column_bmod()
157     lptr = glu.xlsub(fsupc) + d_fsupc;   in column_bmod()
158     luptr = glu.xlusup(fst_col) + d_fsupc;   in column_bmod()
159     nsupr = glu.xlsub(fsupc+1) - glu.xlsub(fsupc); // leading dimension  in column_bmod()
164     ufirst = glu.xlusup(jcol) + d_fsupc;   in column_bmod()
165     Index lda = glu.xlusup(jcol+1) - glu.xlusup(jcol);  in column_bmod()
166     MappedMatrixBlock A( &(glu.lusup.data()[luptr]), nsupc, nsupc, OuterStride<>(lda) );  in column_bmod()
167     VectorBlock<ScalarVector> u(glu.lusup, ufirst, nsupc);   in column_bmod()
170 …  new (&A) MappedMatrixBlock ( &(glu.lusup.data()[luptr+nsupc]), nrow, nsupc, OuterStride<>(lda) );  in column_bmod()
171     VectorBlock<ScalarVector> l(glu.lusup, ufirst+nsupc, nrow);   in column_bmod()