/dports/math/gsl/gsl-2.7/integration/ |
H A D | test.c | 715 double result = 0, abserr=0; in main() local 786 double result = 0, abserr=0; in main() local 861 double result = 0, abserr=0; in main() local 905 double result = 0, abserr=0; in main() local 943 double result = 0, abserr=0; in main() local 1019 double result = 0, abserr=0; in main() local 1091 double result = 0, abserr=0; in main() local 1189 double result = 0, abserr=0; in main() local 1279 double result = 0, abserr=0; in main() local 2282 prev_abserr = abserr; in main() [all …]
|
H A D | gsl_integration.h | 128 double *result, double *abserr, 132 double *result, double *abserr, 136 double *result, double *abserr, 140 double *result, double *abserr, 144 double *result, double *abserr, 148 double *result, double *abserr, 172 double * result, double * abserr, 179 double *result, double *abserr, 187 double *result, double *abserr); 192 double *result, double *abserr); [all …]
|
/dports/math/octave/octave-6.4.0/liboctave/numeric/ |
H A D | Quad.h | 56 double abserr; in integrate() local 63 float abserr; in float_integrate() local 70 double abserr; in integrate() local 77 float abserr; in float_integrate() local 83 double abserr; in integrate() local 89 float abserr; in float_integrate() local 94 double& abserr) in integrate() argument 141 double& abserr); 171 double& abserr); 209 float& abserr); [all …]
|
/dports/math/R/R-4.1.2/src/appl/ |
H A D | integrate.c | 238 *abserr = 0.; in Rdqagi() 502 *abserr = 0.; in rdqagie() 535 if (*ier != 0 || (*abserr <= errbnd && *abserr != resabs) in rdqagie() 709 *abserr = abseps; in rdqagie() 958 *abserr = 0.; in Rdqags() 1216 *abserr = 0.; in rdqagse() 1242 if (*abserr <= epmach * 100. * defabs && *abserr > errbnd) in rdqagse() 1246 if (*ier != 0 || (*abserr <= errbnd && *abserr != resabs) in rdqagse() 1705 *abserr = fmax2(epmach * 50. * *resabs, *abserr); in rdqk15i() 1907 *abserr = fmax2(*abserr, epmach * 5. * fabs(*result)); in rdqelg() [all …]
|
/dports/math/libRmath/R-4.1.1/src/appl/ |
H A D | integrate.c | 238 *abserr = 0.; in Rdqagi() 502 *abserr = 0.; in rdqagie() 535 if (*ier != 0 || (*abserr <= errbnd && *abserr != resabs) in rdqagie() 709 *abserr = abseps; in rdqagie() 958 *abserr = 0.; in Rdqags() 1216 *abserr = 0.; in rdqagse() 1242 if (*abserr <= epmach * 100. * defabs && *abserr > errbnd) in rdqagse() 1246 if (*ier != 0 || (*abserr <= errbnd && *abserr != resabs) in rdqagse() 1705 *abserr = fmax2(epmach * 50. * *resabs, *abserr); in rdqk15i() 1907 *abserr = fmax2(*abserr, epmach * 5. * fabs(*result)); in rdqelg() [all …]
|
/dports/audio/praat/praat-6.2.03/external/gsl/ |
H A D | gsl_integration.h | 128 double *result, double *abserr, 132 double *result, double *abserr, 136 double *result, double *abserr, 140 double *result, double *abserr, 144 double *result, double *abserr, 148 double *result, double *abserr, 172 double * result, double * abserr, 179 double *result, double *abserr, 187 double *result, double *abserr); 192 double *result, double *abserr); [all …]
|
/dports/math/PDL/PDL-2.019/t/ |
H A D | gsl_integ.t | 28 my ($res,$abserr,$neval,$ierr) = gslinteg_qng(\&f1,0,1,0,1e-9); 30 ($res,$abserr,$neval,$ierr) = gslinteg_qng(\&f1,0,1,0,1e-9,{Warn => 'y'}); 33 ($res,$abserr,$ierr) = gslinteg_qag(\&f1,0,1,0,1e-10,1000,1); 38 ($res,$abserr,$ierr) = gslinteg_qags(\&f1,0,1,0,1e-10,1000); 44 ($res,$abserr,$ierr) = gslinteg_qagp(\&f454,$points,0,1e-3,1000); 49 ($res,$abserr,$ierr) = gslinteg_qagi(\&myfn1,1e-7,0,1000); 51 ($res,$abserr,$ierr) = gslinteg_qagi(\&myfn1,1e-7,0,1000,{Warn => 'y'}); 55 ($res,$abserr,$ierr) = gslinteg_qagiu(\&f16,99.9,1e-7,0,1000); 60 ($res,$abserr,$ierr) = gslinteg_qagil(\&myfn2,1.0,1e-7,0,1000); 66 ($res,$abserr,$ierr) = gslinteg_qawc(\&f459,-1,5,0,0,1e-3,1000); [all …]
|
/dports/math/p5-Math-GSL/Math-GSL-0.43/t/ |
H A D | Integration.t | 38 my ($status, $result, $abserr) = gsl_integration_qag ( 47 sprintf('gsl_integration_qag: res=%.18f, abserr=%.18f',$res,$abserr), 48 $abserr 82 my ($status, $result, $abserr, @params); 96 abs($result-$actual)/abs($actual),$abserr 113 sprintf('gsl_integration_qags: res=%.18f, abserr=%.18f',$res,$abserr), 114 $abserr 131 sprintf('gsl_integration_qags: res=%.18f, abserr=%.18f',$res,$abserr), 132 $abserr 152 sprintf('gsl_integration_qng: res=%.18f, abserr=%.18f',$res,$abserr), [all …]
|
/dports/math/PDL/PDL-2.019/GENERATED/PDL/GSL/ |
H A D | INTEG.pm | 144 return ($res,$abserr,$ierr,$neval); 192 return ($res,$abserr,$ierr); 239 return ($res,$abserr,$ierr); 286 return ($res,$abserr,$ierr); 333 return ($res,$abserr,$ierr); 380 return ($res,$abserr,$ierr); 427 return ($res,$abserr,$ierr); 474 return ($res,$abserr,$ierr); 521 return ($res,$abserr,$ierr); 576 return ($res,$abserr,$ierr); [all …]
|
/dports/math/fpc-numlib/fpc-3.2.2/packages/numlib/src/ |
H A D | int.pas | 126 Then abserr := resasc*amin1(1,exp(1.5*ln(200*abserr/resasc))); 128 Then abserr := amax1((50*macheps)*resabs, abserr); 395 abserr := 0; 420 If (ier <> 0) Or ((abserr <= errbnd) And (abserr <> resabs)) Or (abserr=0) 432 errmax := abserr; 435 errsum := abserr; 436 abserr := giant; 598 If ierro=3 Then abserr := abserr+correc; 733 If abserr<reskh Then abserr := reskh 789 If ((abserr<=errbnd) And (abserr<>resabs)) Or (abserr=0) Then exit; [all …]
|
/dports/lang/fpc-source/fpc-3.2.2/packages/numlib/src/ |
H A D | int.pas | 126 Then abserr := resasc*amin1(1,exp(1.5*ln(200*abserr/resasc))); 128 Then abserr := amax1((50*macheps)*resabs, abserr); 395 abserr := 0; 420 If (ier <> 0) Or ((abserr <= errbnd) And (abserr <> resabs)) Or (abserr=0) 432 errmax := abserr; 435 errsum := abserr; 436 abserr := giant; 598 If ierro=3 Then abserr := abserr+correc; 733 If abserr<reskh Then abserr := reskh 789 If ((abserr<=errbnd) And (abserr<>resabs)) Or (abserr=0) Then exit; [all …]
|
/dports/science/py-scipy/scipy-1.7.1/scipy/integrate/quadpack/ |
H A D | dqawfe.f | 2 * result,abserr,neval,ier,rslst,erlst,ierlst,lst,alist,blist, argument 244 abserr = 0.0d+00 259 erlst(1) = abserr 284 abserr = 0.0d+00 331 if(abseps.gt.abserr.and.lst.ne.3) go to 30 332 abserr = abseps 340 if((abserr+0.1d+02*correc).le.epsabs.or. 351 60 abserr = abserr+0.1d+02*correc 354 if(abserr.gt.errsum) go to 80 358 if(ier.ge.1.and.ier.ne.7) abserr = abserr+drl [all …]
|
H A D | dqagse.f | 225 abserr = 0.0d+00 248 elist(1) = abserr 250 if(abserr.le.1.0d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 252 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs).or. 259 errmax = abserr 262 errsum = abserr 263 abserr = oflow 389 abserr = abseps 416 100 if(abserr.eq.oflow) go to 115 418 if(ierro.eq.3) abserr = abserr+correc [all …]
|
H A D | dqagie.f | 227 abserr = 0.0d+00 255 elist(1) = abserr 259 if(abserr.le.1.0d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 261 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs).or. 270 errmax = abserr 273 errsum = abserr 274 abserr = oflow 397 abserr = abseps 423 100 if(abserr.eq.oflow) go to 115 425 if(ierro.eq.3) abserr = abserr+correc [all …]
|
H A D | dqage.f | 1 subroutine dqage(f,a,b,epsabs,epsrel,key,limit,result,abserr, argument 199 abserr = 0.0d+00 216 if(keyf.eq.1) call dqk15(f,a,b,result,abserr,defabs,resabs) 217 if(keyf.eq.2) call dqk21(f,a,b,result,abserr,defabs,resabs) 224 elist(1) = abserr 230 if(abserr.le.0.5d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 232 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs) 233 * .or.abserr.eq.0.0d+00) go to 60 239 errmax = abserr 242 errsum = abserr [all …]
|
H A D | dqawoe.f | 281 abserr = 0.0d+00 308 elist(1) = abserr 310 if(abserr.le.0.1d+03*epmach*defabs.and.abserr.gt.errbnd) ier = 2 319 errmax = abserr 322 errsum = abserr 323 abserr = oflow 473 if(abseps.ge.abserr) go to 100 475 abserr = abseps 505 if(ierro.eq.3) abserr = abserr+correc 508 if(abserr.gt.errsum) go to 170 [all …]
|
/dports/science/libxc/libxc-5.0.0/src/ |
H A D | integrate.c | 302 *abserr = 0.; in xc_rdqagse() 328 if (*abserr <= epmach * 100. * defabs && *abserr > errbnd) in xc_rdqagse() 332 if (*ier != 0 || (*abserr <= errbnd && *abserr != resabs) in xc_rdqagse() 339 errmax = *abserr; in xc_rdqagse() 342 errsum = *abserr; in xc_rdqagse() 343 *abserr = oflow; in xc_rdqagse() 571 *abserr = errsum; in xc_rdqagse() 671 *abserr = oflow; in rdqelg() 778 *abserr = max(*abserr, epmach * 5. * fabs(*result)); in rdqelg() 983 *abserr = *resasc * min(1., pow(*abserr * 200. / *resasc, 1.5)); in rdqk21() [all …]
|
/dports/math/octave/octave-6.4.0/liboctave/external/quadpack/ |
H A D | qagie.f | 227 abserr = 0.0e+00 256 elist(1) = abserr 260 if(abserr.le.1.0e+02*epmach*defabs.and.abserr.gt. 263 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs).or. 272 errmax = abserr 275 errsum = abserr 276 abserr = oflow 405 abserr = abseps 431 100 if(abserr.eq.oflow) go to 115 433 if(ierro.eq.3) abserr = abserr+correc [all …]
|
/dports/math/scilab/scilab-6.1.1/scilab/modules/differential_equations/src/fortran/ |
H A D | dqagse.f | 224 abserr = 0.0d+00 247 elist(1) = abserr 249 if(abserr.le.1.0d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 251 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs).or. 258 errmax = abserr 261 errsum = abserr 262 abserr = oflow 388 abserr = abseps 415 100 if(abserr.eq.oflow) go to 115 417 if(ierro.eq.3) abserr = abserr+correc [all …]
|
/dports/lang/gcc12-devel/gcc-12-20211205/gcc/testsuite/gfortran.dg/g77/ |
H A D | 980310-4.f | 19 subroutine dqage(f,a,b,epsabs,epsrel,limit,result,abserr, argument 175 double precision a,abserr,alist,area,area1,area12,area2,a1,a2,b, 226 abserr = 0.0d+00 240 call dqk15(f,a,b,result,abserr,defabs,resabs) 243 elist(1) = abserr 249 if(abserr.le.0.5d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 251 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs) 252 * .or.abserr.eq.0.0d+00) go to 60 258 errmax = abserr 261 errsum = abserr [all …]
|
/dports/lang/gcc10/gcc-10.3.0/gcc/testsuite/gfortran.dg/g77/ |
H A D | 980310-4.f | 19 subroutine dqage(f,a,b,epsabs,epsrel,limit,result,abserr, argument 175 double precision a,abserr,alist,area,area1,area12,area2,a1,a2,b, 226 abserr = 0.0d+00 240 call dqk15(f,a,b,result,abserr,defabs,resabs) 243 elist(1) = abserr 249 if(abserr.le.0.5d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 251 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs) 252 * .or.abserr.eq.0.0d+00) go to 60 258 errmax = abserr 261 errsum = abserr [all …]
|
/dports/devel/riscv64-none-elf-gcc/gcc-8.4.0/gcc/testsuite/gfortran.dg/g77/ |
H A D | 980310-4.f | 19 subroutine dqage(f,a,b,epsabs,epsrel,limit,result,abserr, argument 175 double precision a,abserr,alist,area,area1,area12,area2,a1,a2,b, 226 abserr = 0.0d+00 240 call dqk15(f,a,b,result,abserr,defabs,resabs) 243 elist(1) = abserr 249 if(abserr.le.0.5d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 251 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs) 252 * .or.abserr.eq.0.0d+00) go to 60 258 errmax = abserr 261 errsum = abserr [all …]
|
/dports/lang/gcc11-devel/gcc-11-20211009/gcc/testsuite/gfortran.dg/g77/ |
H A D | 980310-4.f | 19 subroutine dqage(f,a,b,epsabs,epsrel,limit,result,abserr, argument 175 double precision a,abserr,alist,area,area1,area12,area2,a1,a2,b, 226 abserr = 0.0d+00 240 call dqk15(f,a,b,result,abserr,defabs,resabs) 243 elist(1) = abserr 249 if(abserr.le.0.5d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 251 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs) 252 * .or.abserr.eq.0.0d+00) go to 60 258 errmax = abserr 261 errsum = abserr [all …]
|
/dports/lang/gcc11/gcc-11.2.0/gcc/testsuite/gfortran.dg/g77/ |
H A D | 980310-4.f | 19 subroutine dqage(f,a,b,epsabs,epsrel,limit,result,abserr, argument 175 double precision a,abserr,alist,area,area1,area12,area2,a1,a2,b, 226 abserr = 0.0d+00 240 call dqk15(f,a,b,result,abserr,defabs,resabs) 243 elist(1) = abserr 249 if(abserr.le.0.5d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 251 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs) 252 * .or.abserr.eq.0.0d+00) go to 60 258 errmax = abserr 261 errsum = abserr [all …]
|
/dports/devel/arm-none-eabi-gcc/gcc-8.4.0/gcc/testsuite/gfortran.dg/g77/ |
H A D | 980310-4.f | 19 subroutine dqage(f,a,b,epsabs,epsrel,limit,result,abserr, argument 175 double precision a,abserr,alist,area,area1,area12,area2,a1,a2,b, 226 abserr = 0.0d+00 240 call dqk15(f,a,b,result,abserr,defabs,resabs) 243 elist(1) = abserr 249 if(abserr.le.0.5d+02*epmach*defabs.and.abserr.gt.errbnd) ier = 2 251 if(ier.ne.0.or.(abserr.le.errbnd.and.abserr.ne.resabs) 252 * .or.abserr.eq.0.0d+00) go to 60 258 errmax = abserr 261 errsum = abserr [all …]
|