30 void lm_lmpar(const int n, double* r, const int ldr, const int* Pivot,
33 void lm_qrfac(const int m, const int n, double* A, int* Pivot, double* Rdiag,
35 void lm_qrsolv(const int n, double* r, const int ldr, const int* Pivot,
221     int* Pivot = (int*)pws;  in lmmin()  local
304         lm_qrfac(m, n, fjac, Pivot, wa1, wa2, wa3);  in lmmin()
328             if (wa2[Pivot[j]] == 0)  in lmmin()
333             gnorm = MAX(gnorm, fabs(sum / wa2[Pivot[j]] / fnorm));  in lmmin()
389             lm_lmpar(n, fjac, m, Pivot, diag, qtf, delta, &lmpar,  in lmmin()
403                     wa3[i] -= fjac[j*m+i] * wa1[Pivot[j]];  in lmmin()
553 void lm_lmpar(const int n, double* r, const int ldr, const int* Pivot,  in lm_lmpar()  argument
648         x[Pivot[j]] = aux[j];  in lm_lmpar()
672             aux[j] = diag[Pivot[j]] * xdi[Pivot[j]] / dxnorm;  in lm_lmpar()
690         aux[j] = sum / diag[Pivot[j]];  in lm_lmpar()
716         lm_qrsolv(n, r, ldr, Pivot, aux, qtb, x, Sdiag, xdi);  in lm_lmpar()
740             aux[j] = diag[Pivot[j]] * xdi[Pivot[j]] / dxnorm;  in lm_lmpar()
766 void lm_qrfac(const int m, const int n, double* A, int* Pivot, double* Rdiag,  in lm_qrfac()  argument
818         Pivot[j] = j;  in lm_qrfac()
833             k = Pivot[j];  in lm_qrfac()
834             Pivot[j] = Pivot[kmax];  in lm_qrfac()
835             Pivot[kmax] = k;  in lm_qrfac()
900 void lm_qrsolv(const int n, double* r, const int ldr, const int* Pivot,  in lm_qrsolv()  argument
984         if (diag[Pivot[j]] != 0) {  in lm_qrsolv()
987             Sdiag[j] = diag[Pivot[j]];  in lm_qrsolv()
1054         x[Pivot[j]] = W[j];  in lm_qrsolv()