/dports/science/py-scipy/scipy-1.7.1/scipy/integrate/quadpack/ |
H A D | dqc25s.f | 148 10 hlgth = 0.5d+00*(br-bl) 151 fval(1) = 0.5d+00*f(hlgth+centr)*(fix-hlgth)**beta 153 fval(25) = 0.5d+00*f(centr-hlgth)*(fix+hlgth)**beta 155 u = hlgth*x(i-1) 160 factor = hlgth**(alfa+0.1d+01) 203 u = hlgth*x(i-1) 247 140 hlgth = 0.5d+00*(br-bl) 250 fval(1) = 0.5d+00*f(hlgth+centr)*(fix+hlgth)**alfa 252 fval(25) = 0.5d+00*f(centr-hlgth)*(fix-hlgth)**alfa 254 u = hlgth*x(i-1) [all …]
|
H A D | dqng.f | 276 hlgth = 0.5d+00*(b-a) 277 dhlgth = dabs(hlgth) 291 absc = hlgth*x1(k) 305 absc = hlgth*x2(k) 318 result = res21*hlgth 326 abserr = dabs((res21-res10)*hlgth) 339 absc = hlgth*x3(k) 347 result = res43*hlgth 348 abserr = dabs((res43-res21)*hlgth) 359 absc = hlgth*x4(k) [all …]
|
H A D | dqk15.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 127 hlgth = 0.5d+00*(b-a) 128 dhlgth = dabs(hlgth) 139 absc = hlgth*xgk(jtw) 151 absc = hlgth*xgk(jtwm1) 165 result = resk*hlgth 168 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqk15i.f | 73 * dmax1,dmin1,d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth, local 150 hlgth = 0.5d+00*(b-a) 163 absc = hlgth*xgk(j) 186 result = resk*hlgth 187 resasc = resasc*hlgth 188 resabs = resabs*hlgth 189 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqk21.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 135 hlgth = 0.5d+00*(b-a) 136 dhlgth = dabs(hlgth) 147 absc = hlgth*xgk(jtw) 159 absc = hlgth*xgk(jtwm1) 173 result = resk*hlgth 176 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqk31.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 144 hlgth = 0.5d+00*(b-a) 145 dhlgth = dabs(hlgth) 156 absc = hlgth*xgk(jtw) 168 absc = hlgth*xgk(jtwm1) 182 result = resk*hlgth 185 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqk15w.f | 67 * dmax1,dmin1,d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth, local 129 hlgth = 0.5d+00*(b-a) 130 dhlgth = dabs(hlgth) 141 absc = hlgth*xgk(jtw) 155 absc = hlgth*xgk(jtwm1) 171 result = resk*hlgth 174 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqk41.f | 55 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 160 hlgth = 0.5d+00*(b-a) 161 dhlgth = dabs(hlgth) 172 absc = hlgth*xgk(jtw) 184 absc = hlgth*xgk(jtwm1) 198 result = resk*hlgth 201 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqc25f.f | 100 * d1mach,d2,estc,ests,f,fval,hlgth,oflow,omega,parint,par2,par22, local 158 hlgth = 0.5d+00*(b-a) 159 parint = omega*hlgth 174 10 conc = hlgth*dcos(centr*omega) 175 cons = hlgth*dsin(centr*omega) 315 fval(1) = 0.5d+00*f(centr+hlgth) 317 fval(25) = 0.5d+00*f(centr-hlgth) 320 fval(i) = f(hlgth*x(i-1)+centr) 321 fval(isym) = f(centr-hlgth*x(i-1)) 347 resabs = resabs*dabs(hlgth)
|
H A D | dqk51.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 173 hlgth = 0.5d+00*(b-a) 174 dhlgth = dabs(hlgth) 185 absc = hlgth*xgk(jtw) 197 absc = hlgth*xgk(jtwm1) 211 result = resk*hlgth 214 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqk61.f | 55 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 183 hlgth = 0.5d+00*(b-a) 184 dhlgth = dabs(hlgth) 196 dabsc = hlgth*xgk(jtw) 208 dabsc = hlgth*xgk(jtwm1) 222 result = resk*hlgth 225 abserr = dabs((resk-resg)*hlgth)
|
H A D | dqc25c.f | 58 * cheb12,cheb24,dabs,dlog,dqwgtc,f,fval,hlgth,p2,p3,p4,resabs, local 117 10 hlgth = 0.5d+00*(b-a) 120 fval(1) = 0.5d+00*f(hlgth+centr) 122 fval(25) = 0.5d+00*f(centr-hlgth) 124 u = hlgth*x(i-1)
|
/dports/math/octave/octave-6.4.0/liboctave/external/quadpack/ |
H A D | qk21.f | 53 * fv1,fv2,hlgth,resabs,resg,resk,reskh,result,r1mach,uflow,wg,wgk, local 122 hlgth = 0.5e+00*(b-a) 123 dhlgth = abs(hlgth) 135 absc = hlgth*xgk(jtw) 149 absc = hlgth*xgk(jtwm1) 165 result = resk*hlgth 168 abserr = abs((resk-resg)*hlgth)
|
H A D | qk15i.f | 73 * fv2,hlgth,resabs,resasc,resg,resk,reskh,result,tabsc1,tabsc2, local 142 hlgth = 0.5e+00*(b-a) 160 absc = hlgth*xgk(j) 193 result = resk*hlgth 194 resasc = resasc*hlgth 195 resabs = resabs*hlgth 196 abserr = abs((resk-resg)*hlgth)
|
/dports/math/scilab/scilab-6.1.1/scilab/modules/differential_equations/src/fortran/ |
H A D | dqk21.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 136 hlgth = 0.5d+00*(b-a) 137 dhlgth = dabs(hlgth) 149 absc = hlgth*xgk(jtw) 163 absc = hlgth*xgk(jtwm1) 179 result = resk*hlgth 182 abserr = dabs((resk-resg)*hlgth)
|
/dports/cad/calculix-ccx/CalculiX/ccx_2.18/src/ |
H A D | dqag.f | 696 dhlgth = dabs(hlgth) 733 result = resk*hlgth 885 dhlgth = dabs(hlgth) 922 result = resk*hlgth 1083 dhlgth = dabs(hlgth) 1120 result = resk*hlgth 1296 dhlgth = dabs(hlgth) 1333 result = resk*hlgth 1522 dhlgth = dabs(hlgth) 1559 result = resk*hlgth [all …]
|
/dports/math/maxima/maxima-5.43.2/src/numerical/slatec/ |
H A D | dqc25s.lisp | 73 (setf hlgth (* 0.5 (- br bl))) 77 (* 0.5 (funcall f (+ hlgth centr)) (expt (- fix hlgth) beta))) 88 (* 0.5 (funcall f (- centr hlgth)) (expt (+ fix hlgth) beta))) 93 (* hlgth 101 (setf factor (expt hlgth (+ alfa 1.0))) 168 (* hlgth 230 (setf hlgth (* 0.5 (- br bl))) 234 (* 0.5 (funcall f (+ hlgth centr)) (expt (+ fix hlgth) alfa))) 245 (* 0.5 (funcall f (- centr hlgth)) (expt (- fix hlgth) alfa))) 250 (* hlgth [all …]
|
H A D | dqng.lisp | 167 (setf hlgth (* 0.5 (- b a))) 168 (setf dhlgth (abs hlgth)) 190 (setf absc (* hlgth (f2cl-lib:fref x1 (k) ((1 5))))) 209 (setf absc (* hlgth (f2cl-lib:fref x2 (k) ((1 5))))) 222 (setf result (* res21 hlgth)) 241 (setf abserr (abs (* (- res21 res10) hlgth))) 259 (setf absc (* hlgth (f2cl-lib:fref x3 (k) ((1 11))))) 266 (setf result (* res43 hlgth)) 267 (setf abserr (abs (* (- res43 res21) hlgth))) 290 (setf result (* res87 hlgth)) [all …]
|
H A D | dqk15.lisp | 47 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 50 (type (double-float) uflow reskh resk resg hlgth fval2 fval1 56 (setf hlgth (* 0.5 (- b a))) 57 (setf dhlgth (abs hlgth)) 72 (setf absc (* hlgth (f2cl-lib:fref xgk (jtw) ((1 8))))) 89 (setf absc (* hlgth (f2cl-lib:fref xgk (jtwm1) ((1 8))))) 112 (setf result (* resk hlgth)) 115 (setf abserr (abs (* (- resk resg) hlgth)))
|
H A D | dqk15i.lisp | 49 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 52 (type (double-float) uflow tabsc2 tabsc1 reskh resk resg hlgth 65 (setf hlgth (* 0.5 (- b a))) 82 (setf absc (* hlgth (f2cl-lib:fref xgk (j) ((1 8))))) 126 (setf result (* resk hlgth)) 127 (setf resasc (* resasc hlgth)) 128 (setf resabs (* resabs hlgth)) 129 (setf abserr (abs (* (- resk resg) hlgth)))
|
H A D | dqk41.lisp | 66 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 69 (type (double-float) uflow reskh resk resg hlgth fval2 fval1 75 (setf hlgth (* 0.5 (- b a))) 76 (setf dhlgth (abs hlgth)) 91 (setf absc (* hlgth (f2cl-lib:fref xgk (jtw) ((1 21))))) 108 (setf absc (* hlgth (f2cl-lib:fref xgk (jtwm1) ((1 21))))) 131 (setf result (* resk hlgth)) 134 (setf abserr (abs (* (- resk resg) hlgth)))
|
H A D | dqk51.lisp | 75 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 78 (type (double-float) uflow reskh resk resg hlgth fval2 fval1 84 (setf hlgth (* 0.5 (- b a))) 85 (setf dhlgth (abs hlgth)) 100 (setf absc (* hlgth (f2cl-lib:fref xgk (jtw) ((1 26))))) 117 (setf absc (* hlgth (f2cl-lib:fref xgk (jtwm1) ((1 26))))) 140 (setf result (* resk hlgth)) 143 (setf abserr (abs (* (- resk resg) hlgth)))
|
H A D | dqk21.lisp | 52 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 55 (type (double-float) uflow reskh resk resg hlgth fval2 fval1 61 (setf hlgth (* 0.5 (- b a))) 62 (setf dhlgth (abs hlgth)) 77 (setf absc (* hlgth (f2cl-lib:fref xgk (jtw) ((1 11))))) 94 (setf absc (* hlgth (f2cl-lib:fref xgk (jtwm1) ((1 11))))) 117 (setf result (* resk hlgth)) 120 (setf abserr (abs (* (- resk resg) hlgth)))
|
H A D | dqk61.lisp | 80 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 83 (type (double-float) uflow reskh resk resg hlgth fval2 fval1 89 (setf hlgth (* 0.5 (- b a))) 90 (setf dhlgth (abs hlgth)) 105 (setf dabsc (* hlgth (f2cl-lib:fref xgk (jtw) ((1 31))))) 122 (setf dabsc (* hlgth (f2cl-lib:fref xgk (jtwm1) ((1 31))))) 145 (setf result (* resk hlgth)) 148 (setf abserr (abs (* (- resk resg) hlgth)))
|
H A D | dqk31.lisp | 58 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 61 (type (double-float) uflow reskh resk resg hlgth fval2 fval1 67 (setf hlgth (* 0.5 (- b a))) 68 (setf dhlgth (abs hlgth)) 83 (setf absc (* hlgth (f2cl-lib:fref xgk (jtw) ((1 16))))) 100 (setf absc (* hlgth (f2cl-lib:fref xgk (jtwm1) ((1 16))))) 123 (setf result (* resk hlgth)) 126 (setf abserr (abs (* (- resk resg) hlgth)))
|