/dports/science/py-scipy/scipy-1.7.1/scipy/integrate/odepack/ |
H A D | ewset.f | 1 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 D | prja.f | 1 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 D | kinRoberts_fp.c | 318 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 D | kinRoberts_fp.c | 328 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 D | ark_robertson.c | 304 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 D | ark_robertson_constraints.c | 319 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 D | ark_robertson_constraints.c | 319 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 D | ark_robertson.c | 304 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 D | prja.f | 2 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 D | ewset.f | 2 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 D | prepj.f | 2 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 D | dewset.lisp | 20 (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 D | cvsRoberts_dns_constraints.c | 425 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 D | cvsRoberts_dns.c | 409 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 D | cvRoberts_dns.c | 409 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 D | cvRoberts_dns_constraints.c | 422 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 D | cvRoberts_dns.c | 409 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 D | cvRoberts_dns_constraints.c | 422 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 D | cvsRoberts_dns.c | 409 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 D | cvsRoberts_dns_constraints.c | 425 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 D | idasRoberts_dns.c | 445 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 D | idaRoberts_dns.c | 460 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 D | idaRoberts_dns.c | 460 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 D | idasRoberts_dns.c | 445 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 D | cpodes_nls.c | 64 #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 …]
|