36     Matrix< Scalar, Dynamic, 1 >  wa1, wa2;  in lmpar()  local
41     wa1 = qtb;  in lmpar()
46             wa1[j] = 0.;  in lmpar()
49         wa1[j] /= r(j,j);  in lmpar()
50         temp = wa1[j];  in lmpar()
52             wa1[i] -= r(i,j) * temp;  in lmpar()
56         x[ipvt[j]] = wa1[j];  in lmpar()
77             wa1[j] = diag[l] * (wa2[l] / dxnorm);  in lmpar()
84                 sum += r(i,j) * wa1[i];  in lmpar()
85             wa1[j] = (wa1[j] - sum) / r(j,j);  in lmpar()
87         temp = wa1.blueNorm();  in lmpar()
93         wa1[j] = r.col(j).head(j+1).dot(qtb.head(j+1)) / diag[ipvt[j]];  in lmpar()
95     gnorm = wa1.stableNorm();  in lmpar()
114         wa1 = sqrt(par)* diag;  in lmpar()
117         qrsolv<Scalar>(r, ipvt, wa1, qtb, x, sdiag);  in lmpar()
133             wa1[j] = diag[l] * (wa2[l] / dxnorm);  in lmpar()
136             wa1[j] /= sdiag[j];  in lmpar()
137             temp = wa1[j];  in lmpar()
139                 wa1[i] -= r(i,j) * temp;  in lmpar()
141         temp = wa1.blueNorm();  in lmpar()
193     Matrix< Scalar, Dynamic, 1 >  wa1, wa2;  in lmpar2()  local
200     wa1 = qtb;  in lmpar2()
201     wa1.tail(n-rank).setZero();  in lmpar2()
202 …qr.matrixQR().topLeftCorner(rank, rank).template triangularView<Upper>().solveInPlace(wa1.head(ran…  in lmpar2()
204     x = qr.colsPermutation()*wa1;  in lmpar2()
223         wa1 = qr.colsPermutation().inverse() *  diag.cwiseProduct(wa2)/dxnorm;  in lmpar2()
224 …  qr.matrixQR().topLeftCorner(n, n).transpose().template triangularView<Lower>().solveInPlace(wa1);  in lmpar2()
225         temp = wa1.blueNorm();  in lmpar2()
231 …wa1[j] = qr.matrixQR().col(j).head(j+1).dot(qtb.head(j+1)) / diag[qr.colsPermutation().indices()(j…  in lmpar2()
233     gnorm = wa1.stableNorm();  in lmpar2()
253         wa1 = sqrt(par)* diag;  in lmpar2()
256         qrsolv<Scalar>(s, qr.colsPermutation().indices(), wa1, qtb, x, sdiag);  in lmpar2()
270         wa1 = qr.colsPermutation().inverse() * diag.cwiseProduct(wa2/dxnorm);  in lmpar2()
274             wa1[j] /= sdiag[j];  in lmpar2()
275             temp = wa1[j];  in lmpar2()
277                 wa1[i] -= s(i,j) * temp;  in lmpar2()
279         temp = wa1.blueNorm();  in lmpar2()