| /dports/math/cminpack/cminpack-1.3.6/examples/ |
| H A D | ssqjac.c | 59 fjac -= 1 + ldfjac; in ssqjac()
130 fjac[2 + 2 * ldfjac] = 0.; in ssqjac()
145 fjac[2 + 3 * ldfjac] = 0.; in ssqjac()
162 fjac[2 + 4 * ldfjac] = -fjac[2 + 3 * ldfjac]; in ssqjac()
164 fjac[3 + 3 * ldfjac] = -2. * fjac[3 + 2 * ldfjac]; in ssqjac()
166 fjac[4 + 4 * ldfjac] = -fjac[4 + 1 * ldfjac]; in ssqjac()
205 fjac[i + 3 * ldfjac] = fjac[i + 1 * ldfjac] * fjac[i + 2 * ldfjac]; in ssqjac()
206 fjac[i + 4 * ldfjac] = fjac[i + 3 * ldfjac] / v[i - 1]; in ssqjac()
219 fjac[i + 3 * ldfjac] = -tmp1 * fjac[i + 2 * ldfjac]; in ssqjac()
282 fjac[i + 2 * ldfjac] = temp * fjac[i + 1 * ldfjac]; in ssqjac()
[all …]
|
| H A D | errjac.c | 72 fjac_offset = 1 + ldfjac; in errjac()
99 fjac[1 * ldfjac + 1] = 1.; in errjac()
100 fjac[2 * ldfjac + 1] = 0.; in errjac()
116 fjac[1 * ldfjac + 1] = 1.; in errjac()
119 fjac[4 * ldfjac + 2] = -fjac[3 * ldfjac + 2]; in errjac()
121 fjac[3 * ldfjac + 3] = 2. * fjac[2 * ldfjac + 3]; in errjac()
123 fjac[4 * ldfjac + 4] = -fjac[1 * ldfjac + 4]; in errjac()
237 fjac[1 * ldfjac + 1] = fjac[1 * ldfjac + 1] + 6. * (d__1 * d__1) - 2. * x[2] + 3.; in errjac()
242 fjac[j + k * ldfjac] = fjac[k + j * ldfjac]; in errjac()
283 fjac[j + j * ldfjac] = 2.; in errjac()
[all …]
|
| H A D | vecjac.c | 69 fjac_offset = 1 + ldfjac; in vecjac()
97 fjac[2 * ldfjac + 1] = 0.; in vecjac()
112 fjac[1 * ldfjac + 1] = 1.; in vecjac()
115 fjac[4 * ldfjac + 2] = -fjac[3 * ldfjac + 2]; in vecjac()
117 fjac[3 * ldfjac + 3] = -2. * fjac[2 * ldfjac + 3]; in vecjac()
119 fjac[4 * ldfjac + 4] = -fjac[1 * ldfjac + 4]; in vecjac()
230 fjac[1 * ldfjac + 1] = fjac[1 * ldfjac + 1] + 6. * (d__1 * d__1) - 2. * x[2] + 3.; in vecjac()
235 fjac[j + k * ldfjac] = fjac[k + j * ldfjac]; in vecjac()
276 fjac[j + j * ldfjac] = 2.; in vecjac()
339 fjac[k + k * ldfjac] += 1.; in vecjac()
[all …]
|
| H A D | tlmderc.c | 27 int ldfjac, int iflag);
31 int i, j, ldfjac, maxfev, mode, nprint, info, nfev, njev; in main() local
67 ldfjac = 15; in main()
115 fjac1[i*ldfjac+j] *= jacfac[i] * jacfac[j]; in main()
128 fjac[i*ldfjac+j] *= jacfac[i] * jacfac[j]; in main()
140 if (fjac[i*ldfjac+j] != fjac1[i*ldfjac+j]*covfac) { in main()
152 int ldfjac, int iflag) in fcn() argument
203 fjac[i + ldfjac*0] = -1.; in fcn()
204 fjac[i + ldfjac*1] = tmp1*tmp2/tmp4; in fcn()
205 fjac[i + ldfjac*2] = tmp1*tmp3/tmp4; in fcn()
[all …]
|
| H A D | tchkder_.c | 10 real *fjac, const int *ldfjac, int *iflag);
14 int i, m, n, ldfjac, mode; in main() local
29 ldfjac = 15; in main()
32 __minpack_func__(chkder)(&m, &n, x, NULL, NULL, &ldfjac, xp, NULL, &mode, NULL); in main()
34 fcn(&m, &n, x, fvec, NULL, &ldfjac, &one); in main()
35 fcn(&m, &n, x, NULL, fjac, &ldfjac, &two); in main()
36 fcn(&m, &n, xp, fvecp, NULL, &ldfjac, &one); in main()
59 real *fjac, const int *ldfjac, int *iflag) in fcn() argument
103 fjac[i-1+ *ldfjac*(1-1)] = -1.; in fcn()
104 fjac[i-1+ *ldfjac*(2-1)] = tmp1*tmp2/tmp4; in fcn()
[all …]
|
| H A D | tchkderc.c | 20 real *fjac, int ldfjac, int iflag);
24 int i, ldfjac; in main() local
51 ldfjac = 15; in main()
56 fcn(&data, m, n, x, fvec, NULL, ldfjac, 1); in main()
58 fcn(&data, m, n, x, NULL, fjac, ldfjac, 2); in main()
60 fcn(&data, m, n, xp, fvecp, NULL, ldfjac, 1); in main()
89 real *fjac, int ldfjac, int iflag) in fcn() argument
143 fjac[i + ldfjac*0] = -1.; in fcn()
144 fjac[i + ldfjac*1] = tmp1*tmp2/tmp4; in fcn()
145 fjac[i + ldfjac*2] = tmp1*tmp3/tmp4; in fcn()
[all …]
|
| H A D | tfdjac2_.c | 13 void fcnjac(int m, int n, const real *x, real *fjac, int ldfjac);
17 int i, m, n, ldfjac, mode; in main() local
33 ldfjac = 15; in main()
38 __minpack_func__(chkder)(&m, &n, x, NULL, NULL, &ldfjac, xp, NULL, &mode, NULL); in main()
41 __minpack_func__(fdjac2)(fcn, &m, &n, x, fvec, fdjac, &ldfjac, &iflag, &epsfcn, wa); in main()
42 fcnjac(m, n, x, fjac, ldfjac); in main()
45 __minpack_func__(chkder)(&m, &n, x, fvec, fjac, &ldfjac, NULL, fvecp, &mode, err); in main()
101 real *fjac, int ldfjac) in fcnjac() argument
118 fjac[i-1+ ldfjac*(1-1)] = -1.; in fcnjac()
119 fjac[i-1+ ldfjac*(2-1)] = tmp1*tmp2/tmp4; in fcnjac()
[all …]
|
| H A D | tfdjac2c.c | 18 void fcnjac(int m, int n, const real *x, real *fjac, int ldfjac);
22 int i, ldfjac; in main() local
42 ldfjac = 15; in main()
47 __cminpack_func__(chkder)(m, n, x, NULL, NULL, ldfjac, xp, NULL, 1, NULL); in main()
53 fcnjac(m, n, x, fjac, ldfjac); in main()
57 __cminpack_func__(chkder)(m, n, x, fvec, fdjac, ldfjac, NULL, fvecp, 2, errd); in main()
59 __cminpack_func__(chkder)(m, n, x, fvec, fjac, ldfjac, NULL, fvecp, 2, err); in main()
123 real *fjac, int ldfjac) in fcnjac() argument
138 fjac[i-1+ ldfjac*(1-1)] = -1.; in fcnjac()
139 fjac[i-1+ ldfjac*(2-1)] = tmp1*tmp2/tmp4; in fcnjac()
[all …]
|
| H A D | thybrj1_.c | 9 void fcn(const int *n, const real *x, real *fvec, real *fjac, const int *ldfjac,
14 int j, n, ldfjac, info, lwa; in main() local
27 ldfjac = 9; in main()
36 __minpack_func__(hybrj1)(&fcn, &n, x, fvec, fjac, &ldfjac, &tol, &info, wa, &lwa); in main()
51 void fcn(const int *n, const real *x, real *fvec, real *fjac, const int *ldfjac, in fcn() argument
83 fjac[k-1 + *ldfjac*(j-1)] = 0; in fcn()
85 fjac[k-1 + *ldfjac*(k-1)] = 3 - 4*x[k-1]; in fcn()
87 fjac[k-1 + *ldfjac*(k-1-1)] = -1; in fcn()
90 fjac[k-1 + *ldfjac*(k+1-1)] = -2; in fcn()
|
| H A D | thybrj1c.c | 9 int fcn(void *p, int n, const real *x, real *fvec, real *fjac, int ldfjac,
14 int j, n, ldfjac, info, lwa; in main() local
26 ldfjac = 9; in main()
35 info = __cminpack_func__(hybrj1)(fcn, 0, n, x, fvec, fjac, ldfjac, tol, wa, lwa); in main()
50 int fcn(void *p, int n, const real *x, real *fvec, real *fjac, int ldfjac, in fcn() argument
83 fjac[k-1 + ldfjac*(j-1)] = 0; in fcn()
85 fjac[k-1 + ldfjac*(k-1)] = 3 - 4*x[k-1]; in fcn()
87 fjac[k-1 + ldfjac*(k-1-1)] = -1; in fcn()
90 fjac[k-1 + ldfjac*(k+1-1)] = -2; in fcn()
|
| H A D | thybrj_.c | 9 void fcn(const int *n, const real *x, real *fvec, real *fjac, const int *ldfjac,
15 int j, n, ldfjac, maxfev, mode, nprint, info, nfev, njev, lr; in main() local
29 ldfjac = 9; in main()
46 __minpack_func__(hybrj)(&fcn, &n, x, fvec, fjac, &ldfjac, &xtol, &maxfev, diag, in main()
63 void fcn(const int *n, const real *x, real *fvec, real *fjac, const int *ldfjac, in fcn() argument
100 fjac[k-1 + *ldfjac*(j-1)] = 0; in fcn()
102 fjac[k-1 + *ldfjac*(k-1)] = 3 - 4*x[k-1]; in fcn()
104 fjac[k-1 + *ldfjac*(k-1-1)] = -1; in fcn()
107 fjac[k-1 + *ldfjac*(k+1-1)] = -2; in fcn()
|
| H A D | tlmder_.c | 11 const int *ldfjac, int *iflag);
15 int i, j, m, n, ldfjac, maxfev, mode, nprint, info, nfev, njev; in main() local
32 ldfjac = 15; in main()
47 __minpack_func__(lmder)(&fcn, &m, &n, x, fvec, fjac, &ldfjac, &ftol, &xtol, >ol, in main()
62 __minpack_func__(covar)(&n, fjac, &ldfjac, ipvt, &ftol, wa1); in main()
66 printf("%s%15.7g", j%3==1?"\n ":"", (double)(fjac[(i-1)*ldfjac+j-1]*covfac)); in main()
74 const int *ldfjac, int *iflag) in fcn() argument
116 fjac[i-1 + *ldfjac*(1-1)] = -1.; in fcn()
117 fjac[i-1 + *ldfjac*(2-1)] = tmp1*tmp2/tmp4; in fcn()
118 fjac[i-1 + *ldfjac*(3-1)] = tmp1*tmp3/tmp4; in fcn()
|
| H A D | tlmdifc.c | 22 int i, j, maxfev, mode, nprint, info, nfev, ldfjac; in main() local
47 ldfjac = 15; in main()
64 diag, mode, factor, nprint, &nfev, fjac, ldfjac, in main()
82 __cminpack_func__(covar)(n, fjac1, ldfjac, ipvt, ftol, wa1); in main()
93 k = __cminpack_func__(covar1)(m, n, fnorm*fnorm, fjac, ldfjac, ipvt, ftol, wa1); in main()
97 printf("%s%15.7g", j%3==0?"\n ":"", (double)fjac[i*ldfjac+j]); in main()
106 if (fjac[i*ldfjac+j] != fjac1[i*ldfjac+j]*covfac) { in main()
107 …covar and covar1 differ: %g != %g\n", i, j, (double)fjac[i*ldfjac+j], (double)(fjac1[i*ldfjac+j]*c… in main()
|
| H A D | tlmder1_.c | 11 const int *ldfjac, int *iflag);
15 int j, m, n, ldfjac, info, lwa; in main() local
30 ldfjac = 15; in main()
39 __minpack_func__(lmder1)(&fcn, &m, &n, x, fvec, fjac, &ldfjac, &tol, in main()
54 const int *ldfjac, int *iflag) in fcn() argument
95 fjac[i-1 + *ldfjac*(1-1)] = -1.; in fcn()
96 fjac[i-1 + *ldfjac*(2-1)] = tmp1*tmp2/tmp4; in fcn()
97 fjac[i-1 + *ldfjac*(3-1)] = tmp1*tmp3/tmp4; in fcn()
|
| H A D | tlmder1c.c | 17 int ldfjac, int iflag);
21 int j, ldfjac, info, lwa; in main() local
40 ldfjac = 15; in main()
49 info = __cminpack_func__(lmder1)(fcn, &data, m, n, x, fvec, fjac, ldfjac, tol, in main()
64 int ldfjac, int iflag) in fcn() argument
98 fjac[i + ldfjac*0] = -1.; in fcn()
99 fjac[i + ldfjac*1] = tmp1*tmp2/tmp4; in fcn()
100 fjac[i + ldfjac*2] = tmp1*tmp3/tmp4; in fcn()
|
| H A D | thybrjc.c | 16 int fcn(void *p, int n, const real *x, real *fvec, real *fjac, int ldfjac,
22 int j, n, ldfjac, maxfev, mode, nprint, info, nfev, njev, lr; in main() local
44 ldfjac = 9; in main()
61 info = __cminpack_func__(hybrj)(fcn, p, n, x, fvec, fjac, ldfjac, xtol, maxfev, diag, in main()
87 int fcn(void *p, int n, const real *x, real *fvec, real *fjac, int ldfjac, in fcn() argument
136 fjac[k + ldfjac*j] = 0; in fcn()
138 fjac[k + ldfjac*k] = 3 - 4*x[k]; in fcn()
140 fjac[k + ldfjac*(k-1)] = -1; in fcn()
143 fjac[k + ldfjac*(k+1)] = -2; in fcn()
147 fjac[k + ldfjac*j] *= jacfac[j]; in fcn()
|
| /dports/math/cminpack/cminpack-1.3.6/fortran/ |
| H A D | minpack.f90 | 18 integer, intent(in) :: m,n,ldfjac,mode local
60 integer :: n,ldfjac,nprob local
61 real(r8) :: x(n),fjac(ldfjac,n)
66 integer :: n,ldfjac,iflag,ml,mu local
80 integer :: m,n,ldfjac,iflag local
153 integer :: n,ldfjac,iflag local
161 integer :: n,ldfjac,info,lwa local
167 integer :: n,ldfjac,iflag local
193 integer :: m,n,ldfjac,info,lwa local
337 integer :: m,n,ldfjac,nprob local
[all …]
|
| /dports/misc/visp/visp-3.4.0/modules/vision/src/pose-estimation/ |
| H A D | vpPoseLowe.cpp | 222 *MIJ(jac, 0, i, ldfjac) = 1 / rz; in fcn()
223 *MIJ(jac, 1, i, ldfjac) = 0.0; in fcn()
231 *MIJ(jac, 3, i, ldfjac) = 0.0; in fcn()
232 *MIJ(jac, 4, i, ldfjac) = 0.0; in fcn()
233 *MIJ(jac, 5, i, ldfjac) = 0.0; in fcn()
235 *MIJ(jac, 0, npt + i, ldfjac) = 0.0; in fcn()
239 *MIJ(jac, 3, npt + i, ldfjac) = in fcn()
241 *MIJ(jac, 4, npt + i, ldfjac) = in fcn()
243 *MIJ(jac, 5, npt + i, ldfjac) = in fcn()
246 *MIJ(jac, 3, npt + i, ldfjac) = 0.0; in fcn()
[all …]
|
| /dports/math/gretl/gretl-2021d/minpack/ |
| H A D | lmder.c | 235 (*fcn)(m, n, x, fvec, fjac, ldfjac, &iflag, p); in lmder_()
251 (*fcn)(m, n, x, fvec, fjac, ldfjac, &iflag, p); in lmder_()
301 if (fjac[j + j * ldfjac] != 0.0) { in lmder_()
304 sum += fjac[i + j * ldfjac] * wa4[i]; in lmder_()
306 temp = -sum / fjac[j + j * ldfjac]; in lmder_()
308 wa4[i] += fjac[i + j * ldfjac] * temp; in lmder_()
311 fjac[j + j * ldfjac] = wa1[j]; in lmder_()
351 lmpar_(n, fjac, ldfjac, ipvt, diag, qtf, delta, in lmder_()
369 (*fcn)(m, n, wa2, wa4, fjac, ldfjac, &iflag, p); in lmder_()
391 wa3[i] += fjac[i + j * ldfjac] * temp; in lmder_()
[all …]
|
| H A D | lmdif.c | 199 int ldfjac, int *ipvt, double *qtf, double *wa1, in lmdif_() argument
222 if (n <= 0 || m < n || ldfjac < m || in lmdif_()
277 qrfac_(m, n, fjac, ldfjac, ipvt, wa1, wa2, wa3); in lmdif_()
309 if (fjac[j + j * ldfjac] != 0.0) { in lmdif_()
312 sum += fjac[i + j * ldfjac] * wa4[i]; in lmdif_()
314 temp = -sum / fjac[j + j * ldfjac]; in lmdif_()
316 wa4[i] += fjac[i + j * ldfjac] * temp; in lmdif_()
319 fjac[j + j * ldfjac] = wa1[j]; in lmdif_()
332 sum += fjac[i + j * ldfjac] * (qtf[i] / fnorm); in lmdif_()
359 lmpar_(n, fjac, ldfjac, ipvt, diag, qtf, delta, &par, in lmdif_()
[all …]
|
| /dports/astro/p5-Starlink-AST/Starlink-AST-1.05/ast/cminpack/ |
| H A D | lmder.c | 7 real *fvec, real *fjac, int ldfjac, real ftol, in __cminpack_func__()
240 iflag = fcnder_mn(p, m, n, x, fvec, fjac, ldfjac, 1); in __cminpack_func__()
258 iflag = fcnder_mn(p, m, n, x, fvec, fjac, ldfjac, 2); in __cminpack_func__()
314 if (fjac[j + j * ldfjac] != 0.) { in __cminpack_func__()
317 sum += fjac[i + j * ldfjac] * wa4[i]; in __cminpack_func__()
319 temp = -sum / fjac[j + j * ldfjac]; in __cminpack_func__()
321 wa4[i] += fjac[i + j * ldfjac] * temp; in __cminpack_func__()
324 fjac[j + j * ldfjac] = wa1[j]; in __cminpack_func__()
391 iflag = fcnder_mn(p, m, n, wa2, wa4, fjac, ldfjac, 1); in __cminpack_func__()
415 wa3[i] += fjac[i + j * ldfjac] * temp; in __cminpack_func__()
[all …]
|
| /dports/math/cminpack/cminpack-1.3.6/ |
| H A D | lmder.c | 7 real *fvec, real *fjac, int ldfjac, real ftol, in __cminpack_func__()
240 iflag = fcnder_mn(p, m, n, x, fvec, fjac, ldfjac, 1); in __cminpack_func__()
258 iflag = fcnder_mn(p, m, n, x, fvec, fjac, ldfjac, 2); in __cminpack_func__()
314 if (fjac[j + j * ldfjac] != 0.) { in __cminpack_func__()
317 sum += fjac[i + j * ldfjac] * wa4[i]; in __cminpack_func__()
319 temp = -sum / fjac[j + j * ldfjac]; in __cminpack_func__()
321 wa4[i] += fjac[i + j * ldfjac] * temp; in __cminpack_func__()
324 fjac[j + j * ldfjac] = wa1[j]; in __cminpack_func__()
391 iflag = fcnder_mn(p, m, n, wa2, wa4, fjac, ldfjac, 1); in __cminpack_func__()
415 wa3[i] += fjac[i + j * ldfjac] * temp; in __cminpack_func__()
[all …]
|
| H A D | lmstr.c | 7 real *fvec, real *fjac, int ldfjac, real ftol, in __cminpack_func__()
225 if (n <= 0 || m < n || ldfjac < n || ftol < 0. || xtol < 0. || in __cminpack_func__()
276 fjac[i + j * ldfjac] = 0.; in __cminpack_func__()
296 if (fjac[j + j * ldfjac] == 0.) { in __cminpack_func__()
300 wa2[j] = __cminpack_enorm__(j+1, &fjac[j * ldfjac + 0]); in __cminpack_func__()
306 if (fjac[j + j * ldfjac] != 0.) { in __cminpack_func__()
309 sum += fjac[i + j * ldfjac] * qtf[i]; in __cminpack_func__()
311 temp = -sum / fjac[j + j * ldfjac]; in __cminpack_func__()
313 qtf[i] += fjac[i + j * ldfjac] * temp; in __cminpack_func__()
316 fjac[j + j * ldfjac] = wa1[j]; in __cminpack_func__()
[all …]
|
| /dports/math/maxima/maxima-5.43.2/share/minpack/ |
| H A D | minpack-interface.lisp | 17 (ldfjac m)
39 (declare (type f2cl-lib:integer4 m n ldfjac iflag)
60 (setf (aref fjac (+ row-index (* ldfjac col-index)))
64 (values m n nil nil nil ldfjac iflag)))
72 ldfjac
91 (declare (type f2cl-lib:integer4 m n ldfjac iflag)
127 (ldfjac n)
148 (declare (type f2cl-lib:integer4 n ldfjac iflag)
173 (values n nil nil nil ldfjac iflag)))
180 ldfjac
[all …]
|
| /dports/science/py-scipy/scipy-1.7.1/scipy/optimize/minpack/ |
| H A D | hybrj.f | 2 *subroutine hybrj(fcn,n,x,fvec,fjac,ldfjac,xtol,maxfev,diag,mode,
5 integer n,ldfjac,maxfev,mode,nprint,info,nfev,njev,lr
7 double precision x(n),fvec(n),fjac(ldfjac,n),diag(n),r(lr),
173 if (n .le. 0 .or. ldfjac .lt. n .or. xtol .lt. zero
186 call fcn(n,x,fvec,fjac,ldfjac,iflag)
207 call fcn(n,x,fvec,fjac,ldfjac,iflag)
213 call qrfac(n,n,fjac,ldfjac,.false.,iwa,1,wa1,wa2,wa3)
273 call qform(n,n,fjac,ldfjac,wa1)
292 * call fcn(n,x,fvec,fjac,ldfjac,iflag)
316 call fcn(n,wa2,wa4,fjac,ldfjac,iflag)
[all …]
|