Lines Matching refs:glu
145 * \param glu persistent data to facilitate multiple factors : will be deleted later ??
152 Index SparseLUImpl<Scalar,Index>::memInit(Index m, Index n, Index annz, Index lwork, Index fillratio, Index panel_size, GlobalLU_t& glu)
154 Index& num_expansions = glu.num_expansions; //No memory expansions so far
156 glu.nzumax = glu.nzlumax = (std::min)(fillratio * annz / n, m) * n; // estimated number of nonzeros in U
157 glu.nzlmax = (std::max)(Index(4), fillratio) * annz / 4; // estimated nnz in L factor
165 + (glu.nzlmax + glu.nzumax) * sizeof(Index) + (glu.nzlumax+glu.nzumax) * sizeof(Scalar) + n;
172 glu.xsup.resize(n+1);
173 glu.supno.resize(n+1);
174 glu.xlsub.resize(n+1);
175 glu.xlusup.resize(n+1);
176 glu.xusub.resize(n+1);
181 if( (expand<ScalarVector>(glu.lusup, glu.nzlumax, 0, 0, num_expansions)<0)
182 || (expand<ScalarVector>(glu.ucol, glu.nzumax, 0, 0, num_expansions)<0)
183 || (expand<IndexVector> (glu.lsub, glu.nzlmax, 0, 0, num_expansions)<0)
184 || (expand<IndexVector> (glu.usub, glu.nzumax, 0, 1, num_expansions)<0) )
187 glu.nzlumax /= 2;
188 glu.nzumax /= 2;
189 glu.nzlmax /= 2;
190 if (glu.nzlumax < annz ) return glu.nzlumax;
192 } while (!glu.lusup.size() || !glu.ucol.size() || !glu.lsub.size() || !glu.usub.size());