/dports/science/py-scipy/scipy-1.7.1/scipy/integrate/quadpack/ |
H A D | dqk15i.f | 73 * dmax1,dmin1,d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth, local 152 fval1 = f(tabsc1) 153 if(inf.eq.2) fval1 = fval1+f(-tabsc1) 154 fc = (fval1/centr)/centr 168 fval1 = f(tabsc1) 170 if(inf.eq.2) fval1 = fval1+f(-tabsc1) 172 fval1 = (fval1/absc1)/absc1 174 fv1(j) = fval1 176 fsum = fval1+fval2 179 resabs = resabs+wgk(j)*(dabs(fval1)+dabs(fval2))
|
H A D | dqk15.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 140 fval1 = f(centr-absc) 142 fv1(jtw) = fval1 144 fsum = fval1+fval2 147 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 152 fval1 = f(centr-absc) 154 fv1(jtwm1) = fval1 156 fsum = fval1+fval2 158 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
H A D | dqk21.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 148 fval1 = f(centr-absc) 150 fv1(jtw) = fval1 152 fsum = fval1+fval2 155 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 160 fval1 = f(centr-absc) 162 fv1(jtwm1) = fval1 164 fsum = fval1+fval2 166 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
H A D | dqk31.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 157 fval1 = f(centr-absc) 159 fv1(jtw) = fval1 161 fsum = fval1+fval2 164 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 169 fval1 = f(centr-absc) 171 fv1(jtwm1) = fval1 173 fsum = fval1+fval2 175 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
H A D | dqk15w.f | 67 * dmax1,dmin1,d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth, local 144 fval1 = f(absc1)*w(absc1,p1,p2,p3,p4,kp) 146 fv1(jtw) = fval1 148 fsum = fval1+fval2 151 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 158 fval1 = f(absc1)*w(absc1,p1,p2,p3,p4,kp) 160 fv1(jtwm1) = fval1 162 fsum = fval1+fval2 164 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
H A D | dqk41.f | 55 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 173 fval1 = f(centr-absc) 175 fv1(jtw) = fval1 177 fsum = fval1+fval2 180 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 185 fval1 = f(centr-absc) 187 fv1(jtwm1) = fval1 189 fsum = fval1+fval2 191 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
H A D | dqk51.f | 54 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 186 fval1 = f(centr-absc) 188 fv1(jtw) = fval1 190 fsum = fval1+fval2 193 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 198 fval1 = f(centr-absc) 200 fv1(jtwm1) = fval1 202 fsum = fval1+fval2 204 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
H A D | dqk61.f | 55 * d1mach,epmach,f,fc,fsum,fval1,fval2,fv1,fv2,hlgth,resabs,resasc, local 197 fval1 = f(centr-dabsc) 199 fv1(jtw) = fval1 201 fsum = fval1+fval2 204 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 209 fval1 = f(centr-dabsc) 211 fv1(jtwm1) = fval1 213 fsum = fval1+fval2 215 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
/dports/math/octave/octave-6.4.0/liboctave/external/quadpack/ |
H A D | qk15i.f | 72 * dinf,r1mach,epmach,fc,fsum,fval1,fval2,fvalt,fv1, local 144 call f(tabsc1, ierr, fval1) 147 call f(-tabsc1, ierr, fval1) 149 fval1 = fval1 + fvalt 151 fc = (fval1/centr)/centr 165 call f(tabsc1, ierr, fval1) 172 fval1 = fval1 + fvalt 179 fval1 = (fval1/absc1)/absc1 181 fv1(j) = fval1 183 fsum = fval1+fval2 [all …]
|
H A D | qk21.f | 52 real a,absc,abserr,b,centr,dhlgth,epmach,fc,fsum,fval1,fval2, local 136 call f(centr-absc,ierr,fval1) 140 fv1(jtw) = fval1 142 fsum = fval1+fval2 145 resabs = resabs+wgk(jtw)*(abs(fval1)+abs(fval2)) 150 call f(centr-absc,ierr,fval1) 154 fv1(jtwm1) = fval1 156 fsum = fval1+fval2 158 resabs = resabs+wgk(jtwm1)*(abs(fval1)+abs(fval2))
|
/dports/audio/praat/praat-6.2.03/external/gsl/ |
H A D | gsl_integration__qk.c | 58 const double fval1 = GSL_FN_EVAL (f, center - abscissa); in gsl_integration_qk() local 60 const double fsum = fval1 + fval2; in gsl_integration_qk() 61 fv1[jtw] = fval1; in gsl_integration_qk() 65 result_abs += wgk[jtw] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qk() 72 const double fval1 = GSL_FN_EVAL (f, center - abscissa); in gsl_integration_qk() local 74 fv1[jtwm1] = fval1; in gsl_integration_qk() 76 result_kronrod += wgk[jtwm1] * (fval1 + fval2); in gsl_integration_qk() 77 result_abs += wgk[jtwm1] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qk()
|
H A D | gsl_integration__qng.c | 68 const double fval1 = GSL_FN_EVAL(f, center + abscissa); in gsl_integration_qng() local 70 const double fval = fval1 + fval2; in gsl_integration_qng() 73 resabs += w21a[k] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qng() 75 fv1[k] = fval1; in gsl_integration_qng() 82 const double fval1 = GSL_FN_EVAL(f, center + abscissa); in gsl_integration_qng() local 84 const double fval = fval1 + fval2; in gsl_integration_qng() 86 resabs += w21b[k] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qng() 88 fv3[k] = fval1; in gsl_integration_qng()
|
/dports/math/gsl/gsl-2.7/integration/ |
H A D | qk.c | 58 const double fval1 = GSL_FN_EVAL (f, center - abscissa); in gsl_integration_qk() local 60 const double fsum = fval1 + fval2; in gsl_integration_qk() 61 fv1[jtw] = fval1; in gsl_integration_qk() 65 result_abs += wgk[jtw] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qk() 72 const double fval1 = GSL_FN_EVAL (f, center - abscissa); in gsl_integration_qk() local 74 fv1[jtwm1] = fval1; in gsl_integration_qk() 76 result_kronrod += wgk[jtwm1] * (fval1 + fval2); in gsl_integration_qk() 77 result_abs += wgk[jtwm1] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qk()
|
H A D | qng.c | 68 const double fval1 = GSL_FN_EVAL(f, center + abscissa); in gsl_integration_qng() local 70 const double fval = fval1 + fval2; in gsl_integration_qng() 73 resabs += w21a[k] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qng() 75 fv1[k] = fval1; in gsl_integration_qng() 82 const double fval1 = GSL_FN_EVAL(f, center + abscissa); in gsl_integration_qng() local 84 const double fval = fval1 + fval2; in gsl_integration_qng() 86 resabs += w21b[k] * (fabs (fval1) + fabs (fval2)); in gsl_integration_qng() 88 fv3[k] = fval1; in gsl_integration_qng()
|
/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 150 fval1 = f(centr-absc) 154 fv1(jtw) = fval1 156 fsum = fval1+fval2 159 resabs = resabs+wgk(jtw)*(dabs(fval1)+dabs(fval2)) 164 fval1 = f(centr-absc) 168 fv1(jtwm1) = fval1 170 fsum = fval1+fval2 172 resabs = resabs+wgk(jtwm1)*(dabs(fval1)+dabs(fval2))
|
/dports/math/maxima/maxima-5.43.2/src/numerical/slatec/ |
H A D | dqk15i.lisp | 49 (fsum 0.0) (fval1 0.0) (fval2 0.0) (hlgth 0.0) (resg 0.0) (resk 0.0) 53 fval2 fval1 fsum fc epmach dinf centr absc2 67 (setf fval1 74 (if (= inf 2) (setf fval1 (+ fval1 (funcall f (- tabsc1))))) 75 (setf fc (/ (/ fval1 centr) centr)) 87 (setf fval1 101 (if (= inf 2) (setf fval1 (+ fval1 (funcall f (- tabsc1))))) 103 (setf fval1 (/ (/ fval1 absc1) absc1)) 105 (setf (f2cl-lib:fref fv1 (j) ((1 7))) fval1) 107 (setf fsum (+ fval1 fval2)) [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 73 (setf fval1 (funcall f (- centr absc))) 75 (setf (f2cl-lib:fref fv1 (jtw) ((1 7))) fval1) 77 (setf fsum (+ fval1 fval2)) 83 (+ (abs fval1) (abs fval2))))) 90 (setf fval1 (funcall f (- centr absc))) 92 (setf (f2cl-lib:fref fv1 (jtwm1) ((1 7))) fval1) 94 (setf fsum (+ fval1 fval2)) 99 (+ (abs fval1) (abs fval2)))))
|
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 92 (setf fval1 (funcall f (- centr absc))) 94 (setf (f2cl-lib:fref fv1 (jtw) ((1 20))) fval1) 96 (setf fsum (+ fval1 fval2)) 102 (+ (abs fval1) (abs fval2))))) 109 (setf fval1 (funcall f (- centr absc))) 111 (setf (f2cl-lib:fref fv1 (jtwm1) ((1 20))) fval1) 113 (setf fsum (+ fval1 fval2)) 118 (+ (abs fval1) (abs fval2)))))
|
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 101 (setf fval1 (funcall f (- centr absc))) 103 (setf (f2cl-lib:fref fv1 (jtw) ((1 25))) fval1) 105 (setf fsum (+ fval1 fval2)) 111 (+ (abs fval1) (abs fval2))))) 118 (setf fval1 (funcall f (- centr absc))) 120 (setf (f2cl-lib:fref fv1 (jtwm1) ((1 25))) fval1) 122 (setf fsum (+ fval1 fval2)) 127 (+ (abs fval1) (abs fval2)))))
|
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 78 (setf fval1 (funcall f (- centr absc))) 80 (setf (f2cl-lib:fref fv1 (jtw) ((1 10))) fval1) 82 (setf fsum (+ fval1 fval2)) 88 (+ (abs fval1) (abs fval2))))) 95 (setf fval1 (funcall f (- centr absc))) 97 (setf (f2cl-lib:fref fv1 (jtwm1) ((1 10))) fval1) 99 (setf fsum (+ fval1 fval2)) 104 (+ (abs fval1) (abs fval2)))))
|
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 106 (setf fval1 (funcall f (- centr dabsc))) 108 (setf (f2cl-lib:fref fv1 (jtw) ((1 30))) fval1) 110 (setf fsum (+ fval1 fval2)) 116 (+ (abs fval1) (abs fval2))))) 123 (setf fval1 (funcall f (- centr dabsc))) 125 (setf (f2cl-lib:fref fv1 (jtwm1) ((1 30))) fval1) 127 (setf fsum (+ fval1 fval2)) 132 (+ (abs fval1) (abs fval2)))))
|
/dports/cad/calculix-ccx/CalculiX/ccx_2.18/src/ |
H A D | dqag.f | 710 fv1(jtw) = fval1 712 fsum = fval1+fval2 722 fv1(jtwm1) = fval1 724 fsum = fval1+fval2 899 fv1(jtw) = fval1 901 fsum = fval1+fval2 911 fv1(jtwm1) = fval1 1097 fv1(jtw) = fval1 1310 fv1(jtw) = fval1 1536 fv1(jtw) = fval1 [all …]
|
/dports/math/teyjus/teyjus-2.1-7-ge63f40a/source/simulator/builtins/ |
H A D | evalexp.c | 170 float fval1, fval2; in BIEVAL_evalFloat() local 188 fval1 = BIEVAL_evalFloat(argVec); in BIEVAL_evalFloat() 191 return fval1 / fval2; in BIEVAL_evalFloat() 193 fval1 = BIEVAL_evalFloat(argVec); in BIEVAL_evalFloat() 194 if (fval1 < 0.0) EM_error(BI_ERROR_NEG_SQRT); in BIEVAL_evalFloat() 195 return sqrt(fval1); in BIEVAL_evalFloat() 203 fval1 = BIEVAL_evalFloat(argVec); in BIEVAL_evalFloat() 204 if (fval1 <= 0.0) EM_error(BI_ERROR_NEG_SQRT); in BIEVAL_evalFloat() 205 return log(fval1); in BIEVAL_evalFloat() 207 fval1 = BIEVAL_evalFloat(argVec); in BIEVAL_evalFloat() [all …]
|
/dports/math/kig/kig-21.12.3/misc/ |
H A D | kignumerics.cpp | 121 double fval1 = a; // double check... in calcCubicRoot() local 123 fval1 = b + xmin*fval1; in calcCubicRoot() 125 fval1 = c + xmin*fval1; in calcCubicRoot() 127 fval1 = d + xmin*fval1; in calcCubicRoot() 129 assert ( fval1 * fval2 <= 0 ); in calcCubicRoot() 198 double fval1, fval2, fpval1, fpval2, fppval1, fppval2; in calcCubicRootwithNewton() local 199 calcCubicDerivatives ( xmin, a, b, c, d, fval1, fpval1, fppval1 ); in calcCubicRootwithNewton() 201 assert ( fval1 * fval2 <= 0 ); in calcCubicRootwithNewton() 207 assert ( fval1 * fval2 <= 0 ); in calcCubicRootwithNewton() 212 if ( fval1*fval <= 0 ) in calcCubicRootwithNewton() [all …]
|
/dports/math/py-pandas/pandas-1.2.5/pandas/tests/frame/methods/ |
H A D | test_rank.py | 160 fval1 = float_frame.T 161 fval1 = fval1.fillna((fval1.min() - 1).to_dict()).T 162 fval1 = fval1.fillna(np.inf).values 165 exp1 = np.apply_along_axis(rankdata, 1, fval1) 191 fval1 = float_frame.T 192 fval1 = fval1.fillna((fval1.min() - 1).to_dict()).T 193 fval1 = fval1.fillna(np.inf).values 196 exp1 = np.apply_along_axis(rankdata, 1, -fval1)
|