35     Matrix<Scalar,Dynamic,1>  wa(n);  in lmqrsolv()
45     wa = qtb;  in lmqrsolv()
72             temp = givens.c() * wa[k] + givens.s() * qtbpj;  in lmqrsolv()
73             qtbpj = -givens.s() * wa[k] + givens.c() * qtbpj;  in lmqrsolv()
74             wa[k] = temp;  in lmqrsolv()
90     wa.tail(n-nsing).setZero();  in lmqrsolv()
91 …s.topLeftCorner(nsing, nsing).transpose().template triangularView<Upper>().solveInPlace(wa.head(ns…  in lmqrsolv()
98     x = iPerm * wa;   in lmqrsolv()
115     Matrix<Scalar,Dynamic,1>  wa(n);  in lmqrsolv()
123     wa = qtb;  in lmqrsolv()
156         temp = givens.c() * wa(k) + givens.s() * qtbpj;   in lmqrsolv()
157         qtbpj = -givens.s() * wa(k) + givens.c() * qtbpj;  in lmqrsolv()
158         wa(k) = temp;  in lmqrsolv()
176     wa.tail(n-nsing).setZero();  in lmqrsolv()
178 …wa.head(nsing) = R.topLeftCorner(nsing,nsing).template triangularView<Upper>().solve/*InPlace*/(wa…  in lmqrsolv()
182     x = iPerm * wa;   in lmqrsolv()