numlib/src/int.pas

36 Procedure int1fr(f: rfunc1r; a, b, ae: ArbFloat; Var integral, err: ArbFloat;
87  dhlgth := abs(hlgth);
91  resabs := abs(resk);
103        resabs := resabs+wgk[jtw]*(abs(fval1)+abs(fval2))
115        resabs := resabs+wgk[jtwm1]*(abs(fval1)+abs(fval2))
118  resasc := wgk[11]*abs(fc-reskh);
120      resasc := resasc+wgk[j]*(abs(fv1^[j]-reskh)+abs(fv2^[j]-reskh));
124    abserr := abs((resk-resg)*hlgth);
271         e0abs := abs(e0);
272         e1abs := abs(e1);
273         e2abs := abs(e2);
275         err2 := abs(delta2);
282         err3 := abs(delta3);
292            If abserr<5*macheps*abs(result)
293             Then abserr := 5*macheps*abs(result);
300         err1 := abs(delta1);
301         e3abs := abs(e3);
314            epsinf := abs(ss*e1);
322               error := err2+abs(res-e2)+err3;
358          abserr := abs(result-res3la[3]) +
359                    abs(result-res3la[2]) +
360                    abs(result-res3la[1]);
364          If abserr<5*macheps*abs(result)
365           Then abserr := 5*macheps*abs(result)
412  dres := abs(result);
466           If (abs(rlist^[maxerr]-area12) <= 1e-5*abs(area12))
476       errbnd := amax1(epsabs, epsrel*abs(area));
480       If amax1(abs(a1),abs(b2)) <= (1+100*macheps)*(abs(a2)+1000*midget)
516                       If abs(b1-a1) > small Then erlarg := erlarg+erro12;
518                          (abs(blist^[maxerr]-alist^[maxerr]) <= small) Then
533                                   If abs(blist^[maxerr]-alist^[maxerr]) > small
557                                   ertest := amax1(epsabs,epsrel*abs(reseps));
579                         small := abs(b-a)*0.375;
609       If abserr/abs(result) > errsum/abs(area) Then p0 := true
614       If (ksgn=-1) And (amax1(abs(result),abs(area)) <= defabs*0.01)
617       If (0.01 > result/area) Or (result/area > 100) Or (errsum>abs(area))
686       resabs := abs(resk);
705         resabs := resabs+wgk[j]*(abs(fval1)+abs(fval2))
709       resasc := wgk[8]*abs(fc-reskh);
713        resasc := resasc+wgk[j]*(abs(fv1[j]-reskh)+abs(fv2[j]-reskh));
718       abserr := abs((resk-resg)*hlgth);
781       dres := abs(result);
841            If (abs(rlist^[maxerr]-area12)<=1e-5*abs(area12)) And
850         errbnd := epsrel*abs(area);
856         h := abs(a1);
857        If h<abs(b2) Then h := abs(b2);
858         If h<=(1+100*macheps)*(abs(a2)+1000*midget) Then ier := 3;
893            If (abs(b1-a1)>small) Then erlarg := erlarg+erro12;
896             If (abs(blist^[maxerr]-alist^[maxerr])>small)
915                  If (abs(blist^[maxerr]-alist^[maxerr])>small)
936                   ertest := epsrel*abs(reseps);
975            h := abs(result);
976            If h<abs(area) Then h := abs(area);
978             If (0.01>result/area) Or (result/area>100) Or (errsum>abs(area))
987        If abserr/abs(result)>errsum/abs(area)
996            h := abs(result);
997            If h<abs(area) Then h := abs(area);
999             If (0.01>result/area) Or (result/area>100) Or (errsum>abs(area))
1014            h := abs(result);
1015            If h<abs(area) Then h := abs(area);
1017             If (0.01>result/area) Or (result/area>100) Or (errsum>abs(area))
1022 Procedure int1fr(f: rfunc1r; a, b, ae: ArbFloat; Var integral, err: ArbFloat;
1028  integral := NaN;
1029      If abs(a)=infinity
1030       Then If abs(b)=infinity
1035                     qagie(f, 0, 2, ae, epsrel, integral, err, ier);
1036                     If a=infinity Then integral := -integral
1039                 Then qagie(f, b, -1, ae, epsrel, integral, err, ier)
1042                     qagie(f, b, 1, ae, epsrel, integral, err, ier);
1043                     integral := -integral
1045      Else If abs(b)=infinity
1049                     qagie(f, a, -1, ae, epsrel, integral, err, ier);
1050                     integral := -integral
1052  Else qagie(f, a, 1, ae, epsrel, integral, err, ier)
1053           Else qagse(f, a, b, ae, epsrel, limit, integral, err, neval, ier, last);