Home
last modified time | relevance | path

Searched refs:ewt (Results 1 – 25 of 492) sorted by relevance

12345678910>>...20

/dports/science/py-scipy/scipy-1.7.1/scipy/integrate/odepack/
H A Dewset.f1 subroutine ewset (n, itol, rtol, atol, ycur, ewt) argument
11 double precision rtol, atol, ycur, ewt local
12 dimension rtol(1), atol(1), ycur(n), ewt(n)
17 15 ewt(i) = rtol(1)*dabs(ycur(i)) + atol(1)
21 25 ewt(i) = rtol(1)*dabs(ycur(i)) + atol(i)
25 35 ewt(i) = rtol(i)*dabs(ycur(i)) + atol(1)
29 45 ewt(i) = rtol(i)*dabs(ycur(i)) + atol(i)
H A Dprja.f1 subroutine prja (neq, y, yh, nyh, ewt, ftem, savf, wm, iwm, argument
12 double precision y, yh, ewt, ftem, savf, wm, rsav local
18 dimension neq(1), y(1), yh(nyh,*), ewt(1), ftem(1), savf(1),
80 200 fac = vmnorm (n, savf, ewt)
87 r = dmax1(srur*dabs(yj),r0/ewt(j))
103 pdnorm = fnorm (n, wm(3), ewt)/dabs(hl0)
144 fac = vmnorm (n, savf, ewt)
150 r = dmax1(srur*dabs(yi),r0/ewt(i))
160 r = dmax1(srur*dabs(yjj),r0/ewt(jj))
172 pdnorm = bnorm (n, wm(3), meband, ml, mu, ewt)/dabs(hl0)
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/kinsol/serial/
H A DkinRoberts_fp.c318 N_Vector ewt; /* error weight vector */ in check_ans() local
323 ewt = N_VClone(u); in check_ans()
331 N_VAbs(ref, ewt); in check_ans()
332 N_VScale(rtol, ewt, ewt); in check_ans()
333 N_VAddConst(ewt, atol, ewt); in check_ans()
334 if (N_VMin(ewt) <= ZERO) { in check_ans()
338 N_VInv(ewt, ewt); in check_ans()
342 err = N_VWrmsNorm(ref, ewt); in check_ans()
353 N_VDestroy(ewt); in check_ans()
/dports/math/sundials/sundials-5.7.0/examples/kinsol/serial/
H A DkinRoberts_fp.c328 N_Vector ewt; /* error weight vector */ in check_ans() local
333 ewt = N_VClone(u); in check_ans()
341 N_VAbs(ref, ewt); in check_ans()
342 N_VScale(rtol, ewt, ewt); in check_ans()
343 N_VAddConst(ewt, atol, ewt); in check_ans()
344 if (N_VMin(ewt) <= ZERO) { in check_ans()
348 N_VInv(ewt, ewt); in check_ans()
352 err = N_VWrmsNorm(ref, ewt); in check_ans()
363 N_VDestroy(ewt); in check_ans()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/arkode/C_serial/
H A Dark_robertson.c304 N_Vector ewt; /* error weight vector */ in check_ans() local
311 ewt = N_VClone(y); in check_ans()
319 N_VAbs(ref, ewt); in check_ans()
320 N_VScale(rtol, ewt, ewt); in check_ans()
321 N_VAddConst(ewt, atol, ewt); in check_ans()
322 if (N_VMin(ewt) <= ZERO) { in check_ans()
326 N_VInv(ewt, ewt); in check_ans()
330 err = N_VWrmsNorm(ref, ewt); in check_ans()
341 N_VDestroy(ewt); in check_ans()
H A Dark_robertson_constraints.c319 N_Vector ewt; /* error weight vector */ in check_ans() local
326 ewt = N_VClone(y); in check_ans()
334 N_VAbs(ref, ewt); in check_ans()
335 N_VScale(rtol, ewt, ewt); in check_ans()
336 N_VAddConst(ewt, atol, ewt); in check_ans()
337 if (N_VMin(ewt) <= ZERO) { in check_ans()
341 N_VInv(ewt, ewt); in check_ans()
345 err = N_VWrmsNorm(ref, ewt); in check_ans()
356 N_VDestroy(ewt); in check_ans()
/dports/math/sundials/sundials-5.7.0/examples/arkode/C_serial/
H A Dark_robertson_constraints.c319 N_Vector ewt; /* error weight vector */ in check_ans() local
326 ewt = N_VClone(y); in check_ans()
334 N_VAbs(ref, ewt); in check_ans()
335 N_VScale(rtol, ewt, ewt); in check_ans()
336 N_VAddConst(ewt, atol, ewt); in check_ans()
337 if (N_VMin(ewt) <= ZERO) { in check_ans()
341 N_VInv(ewt, ewt); in check_ans()
345 err = N_VWrmsNorm(ref, ewt); in check_ans()
356 N_VDestroy(ewt); in check_ans()
H A Dark_robertson.c304 N_Vector ewt; /* error weight vector */ in check_ans() local
311 ewt = N_VClone(y); in check_ans()
319 N_VAbs(ref, ewt); in check_ans()
320 N_VScale(rtol, ewt, ewt); in check_ans()
321 N_VAddConst(ewt, atol, ewt); in check_ans()
322 if (N_VMin(ewt) <= ZERO) { in check_ans()
326 N_VInv(ewt, ewt); in check_ans()
330 err = N_VWrmsNorm(ref, ewt); in check_ans()
341 N_VDestroy(ewt); in check_ans()
/dports/math/scilab/scilab-6.1.1/scilab/modules/differential_equations/src/fortran/
H A Dprja.f2 subroutine prja (neq, y, yh, nyh, ewt, ftem, savf, wm, iwm, argument
14 double precision y, yh, ewt, ftem, savf, wm local
20 dimension neq(*), y(*), yh(nyh,*), ewt(*), ftem(*), savf(*),
84 200 fac = vmnorm (n, savf, ewt)
91 r = max(srur*abs(yj),r0/ewt(j))
104 pdnorm = fnorm (n, wm(3), ewt)/abs(hl0)
139 fac = vmnorm (n, savf, ewt)
145 r = max(srur*abs(yi),r0/ewt(i))
152 r = max(srur*abs(yjj),r0/ewt(jj))
164 pdnorm = bnorm (n, wm(3), meband, ml, mu, ewt)/abs(hl0)
H A Dewset.f2 subroutine ewset (n, itol, rtol, atol, ycur, ewt) argument
14 double precision rtol, atol, ycur, ewt local
16 dimension rtol(*), atol(*), ycur(n), ewt(n)
22 ewt(i) = rtoli*abs(ycur(i)) + atoli
H A Dprepj.f2 subroutine prepj (neq, y, yh, nyh, ewt, ftem, savf, wm, iwm, argument
13 double precision y, yh, ewt, ftem, savf, wm local
18 dimension neq(*), y(*), yh(nyh,*), ewt(*), ftem(*), savf(*),
80 200 fac = vnorm (n, savf, ewt)
87 r = max(srur*abs(yj),r0/ewt(j))
119 if (abs(r0) .lt. uround/ewt(i)) go to 320
152 fac = vnorm (n, savf, ewt)
158 r = max(srur*abs(yi),r0/ewt(i))
165 r = max(srur*abs(yjj),r0/ewt(jj))
/dports/math/maxima/maxima-5.43.2/share/odepack/src/
H A Ddewset.lisp20 (defun dewset (n itol rtol atol ycur ewt)
21 (declare (type (array double-float (*)) ewt ycur atol rtol)
31 (setf (f2cl-lib:fref ewt (i) ((1 n)))
42 (setf (f2cl-lib:fref ewt (i) ((1 n)))
53 (setf (f2cl-lib:fref ewt (i) ((1 n)))
64 (setf (f2cl-lib:fref ewt (i) ((1 n)))
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/cvodes/serial/
H A DcvsRoberts_dns_constraints.c425 N_Vector ewt; /* error weight vector */ in check_ans() local
430 ewt = N_VClone(y); in check_ans()
438 N_VAbs(ref, ewt); in check_ans()
439 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
440 if (N_VMin(ewt) <= ZERO) { in check_ans()
444 N_VInv(ewt, ewt); in check_ans()
448 err = N_VWrmsNorm(ref, ewt); in check_ans()
459 N_VDestroy(ewt); in check_ans()
H A DcvsRoberts_dns.c409 N_Vector ewt; /* error weight vector */ in check_ans() local
415 ewt = N_VClone(y); in check_ans()
423 N_VAbs(ref, ewt); in check_ans()
424 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
425 if (N_VMin(ewt) <= ZERO) { in check_ans()
429 N_VInv(ewt, ewt); in check_ans()
433 err = N_VWrmsNorm(ref, ewt); in check_ans()
444 N_VDestroy(ewt); in check_ans()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/cvode/serial/
H A DcvRoberts_dns.c409 N_Vector ewt; /* error weight vector */ in check_ans() local
415 ewt = N_VClone(y); in check_ans()
423 N_VAbs(ref, ewt); in check_ans()
424 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
425 if (N_VMin(ewt) <= ZERO) { in check_ans()
429 N_VInv(ewt, ewt); in check_ans()
433 err = N_VWrmsNorm(ref, ewt); in check_ans()
444 N_VDestroy(ewt); in check_ans()
H A DcvRoberts_dns_constraints.c422 N_Vector ewt; /* error weight vector */ in check_ans() local
427 ewt = N_VClone(y); in check_ans()
435 N_VAbs(ref, ewt); in check_ans()
436 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
437 if (N_VMin(ewt) <= ZERO) { in check_ans()
441 N_VInv(ewt, ewt); in check_ans()
445 err = N_VWrmsNorm(ref, ewt); in check_ans()
456 N_VDestroy(ewt); in check_ans()
/dports/math/sundials/sundials-5.7.0/examples/cvode/serial/
H A DcvRoberts_dns.c409 N_Vector ewt; /* error weight vector */ in check_ans() local
415 ewt = N_VClone(y); in check_ans()
423 N_VAbs(ref, ewt); in check_ans()
424 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
425 if (N_VMin(ewt) <= ZERO) { in check_ans()
429 N_VInv(ewt, ewt); in check_ans()
433 err = N_VWrmsNorm(ref, ewt); in check_ans()
444 N_VDestroy(ewt); in check_ans()
H A DcvRoberts_dns_constraints.c422 N_Vector ewt; /* error weight vector */ in check_ans() local
427 ewt = N_VClone(y); in check_ans()
435 N_VAbs(ref, ewt); in check_ans()
436 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
437 if (N_VMin(ewt) <= ZERO) { in check_ans()
441 N_VInv(ewt, ewt); in check_ans()
445 err = N_VWrmsNorm(ref, ewt); in check_ans()
456 N_VDestroy(ewt); in check_ans()
/dports/math/sundials/sundials-5.7.0/examples/cvodes/serial/
H A DcvsRoberts_dns.c409 N_Vector ewt; /* error weight vector */ in check_ans() local
415 ewt = N_VClone(y); in check_ans()
423 N_VAbs(ref, ewt); in check_ans()
424 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
425 if (N_VMin(ewt) <= ZERO) { in check_ans()
429 N_VInv(ewt, ewt); in check_ans()
433 err = N_VWrmsNorm(ref, ewt); in check_ans()
444 N_VDestroy(ewt); in check_ans()
H A DcvsRoberts_dns_constraints.c425 N_Vector ewt; /* error weight vector */ in check_ans() local
430 ewt = N_VClone(y); in check_ans()
438 N_VAbs(ref, ewt); in check_ans()
439 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
440 if (N_VMin(ewt) <= ZERO) { in check_ans()
444 N_VInv(ewt, ewt); in check_ans()
448 err = N_VWrmsNorm(ref, ewt); in check_ans()
459 N_VDestroy(ewt); in check_ans()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/idas/serial/
H A DidasRoberts_dns.c445 N_Vector ewt; /* error weight vector */ in check_ans() local
450 ewt = N_VClone(y); in check_ans()
458 N_VAbs(ref, ewt); in check_ans()
459 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
460 if (N_VMin(ewt) <= ZERO) { in check_ans()
464 N_VInv(ewt, ewt); in check_ans()
468 err = N_VWrmsNorm(ref, ewt); in check_ans()
479 N_VDestroy(ewt); in check_ans()
/dports/science/cantera/cantera-2.5.1-611-gc4d6ecc15/ext/sundials/examples/ida/serial/
H A DidaRoberts_dns.c460 N_Vector ewt; /* error weight vector */ in check_ans() local
465 ewt = N_VClone(y); in check_ans()
473 N_VAbs(ref, ewt); in check_ans()
474 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
475 if (N_VMin(ewt) <= ZERO) { in check_ans()
479 N_VInv(ewt, ewt); in check_ans()
483 err = N_VWrmsNorm(ref, ewt); in check_ans()
494 N_VDestroy(ewt); in check_ans()
/dports/math/sundials/sundials-5.7.0/examples/ida/serial/
H A DidaRoberts_dns.c460 N_Vector ewt; /* error weight vector */ in check_ans() local
465 ewt = N_VClone(y); in check_ans()
473 N_VAbs(ref, ewt); in check_ans()
474 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
475 if (N_VMin(ewt) <= ZERO) { in check_ans()
479 N_VInv(ewt, ewt); in check_ans()
483 err = N_VWrmsNorm(ref, ewt); in check_ans()
494 N_VDestroy(ewt); in check_ans()
/dports/math/sundials/sundials-5.7.0/examples/idas/serial/
H A DidasRoberts_dns.c445 N_Vector ewt; /* error weight vector */ in check_ans() local
450 ewt = N_VClone(y); in check_ans()
458 N_VAbs(ref, ewt); in check_ans()
459 N_VLinearSum(rtol, ewt, RCONST(10.0), atol, ewt); in check_ans()
460 if (N_VMin(ewt) <= ZERO) { in check_ans()
464 N_VInv(ewt, ewt); in check_ans()
468 err = N_VWrmsNorm(ref, ewt); in check_ans()
479 N_VDestroy(ewt); in check_ans()
/dports/science/simbody/simbody-Simbody-3.7/SimTKmath/Integrators/src/CPodes/sundials/src/cpodes/
H A Dcpodes_nls.c64 #define ewt (cp_mem->cp_ewt) macro
267 del = N_VWrmsNorm(acor, ewt); in cpNlsFunctionalExpl()
277 acnrm = (m == 0) ? del : N_VWrmsNorm(acor, ewt); in cpNlsFunctionalExpl()
462 retval = lsolve(cp_mem, b, ewt, y, NULL, ftemp); in cpNewtonIterationExpl()
476 del = N_VWrmsNorm(b, ewt); in cpNewtonIterationExpl()
497 acnrm = (m==0) ? del : N_VWrmsNorm(acor, ewt); in cpNewtonIterationExpl()
600 del = N_VWrmsNorm(tempv, ewt); in cpNlsFunctionalImpl()
610 acnrm = (m == 0) ? del : N_VWrmsNorm(acor, ewt); in cpNlsFunctionalImpl()
767 retval = lsolve(cp_mem, tempv, ewt, y, yp, ftemp); in cpNewtonIterationImpl()
773 del = N_VWrmsNorm(tempv, ewt); in cpNewtonIterationImpl()
[all …]

12345678910>>...20