211     double* fjac = (double*)pws;  in lmmin()  local
268                 fjac[j*m+i] = (wf[i] - fvec[i]) / step;  in lmmin()
277                     printf("%.5e ", fjac[j*m+i]);  in lmmin()
304         lm_qrfac(m, n, fjac, Pivot, wa1, wa2, wa3);  in lmmin()
312             temp3 = fjac[j*m+j];  in lmmin()
316                     sum += fjac[j*m+i] * wf[i];  in lmmin()
319                     wf[i] += fjac[j*m+i] * temp;  in lmmin()
321             fjac[j*m+j] = wa1[j];  in lmmin()
332                 sum += fjac[j*m+i] * qtf[i];  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()