33 void lm_qrfac(const int m, const int n, double* A, int* Pivot, double* Rdiag,
766 void lm_qrfac(const int m, const int n, double* A, int* Pivot, double* Rdiag,  in lm_qrfac()  argument
817         W[j] = Rdiag[j] = Acnorm[j] = lm_enorm(m, &A[j*m]);  in lm_qrfac()
828             if (Rdiag[k] > Rdiag[kmax])  in lm_qrfac()
842             Rdiag[kmax] = Rdiag[j];  in lm_qrfac()
850             Rdiag[j] = 0;  in lm_qrfac()
878             if (Rdiag[k] != 0) {  in lm_qrfac()
879                 temp = A[m*k+j] / Rdiag[k];  in lm_qrfac()
881                     Rdiag[k] *= sqrt(1 - SQR(temp));  in lm_qrfac()
882                     temp = Rdiag[k] / W[k];  in lm_qrfac()
886                     Rdiag[k] = lm_enorm(m-j-1, &A[m*k+j+1]);  in lm_qrfac()
887                     W[k] = Rdiag[k];  in lm_qrfac()
892         Rdiag[j] = -ajnorm;  in lm_qrfac()