crlibm-1.0beta4/maple/common-procedures.mpl

450 local n, deg, delta, deltap, i, S, P, Snorm, Smin, Smax, prec:
458 Smax:=abs(S):
459 Smin:=Smax:
465   Smin := abs(coeff(poly,x,i)) - xmax*Smax :
469   delta:= evalf(xmax*deltap + 2**(-prec)*xmax*Smax):
472     if abs(coeff(poly,x,i)) < xmax*Smax  # may be improved to xmax*Smax/2
474          printf("    coeff=%1.20e,  xmax*Smax=%1.20e"  ,  abs(coeff(poly,x,i)),  xmax*Smax ):
481   Smax := Snorm + deltap:
483 deltap, Smin, Smax:
492 local deg, p, rho, deltap, Smin, Smax:
498    deltap, Smin, Smax := compute_abs_rounding_error(poly/x,xmax, n-1):
499    rho :=  (2^(-p)*(Smax+deltap) +deltap ) / Smin :
501    deltap, Smin, Smax := compute_abs_rounding_error(poly,xmax, n):
522 local n, deg, delta, deltap, i, S, P, Snorm, Smin, Smax, prec:
530 Smax:=abs(S):
531 Smin:=Smax:
538   Smin := abs(coeff(poly,x,i)) - xmax*Smax :
542   delta:= evalf(xmax*deltap + 2**(-prec)*xmax*Smax):
545     if abs(coeff(poly,x,i)) < xmax*Smax  # may be improved to xmax*Smax/2
547          printf("    coeff=%1.20e,  xmax*Smax=%1.20e"  ,  abs(coeff(poly,x,i)),  xmax*Smax ):
553   Smax := Snorm + deltap:
555 deltap, Smin, Smax:
565 local deg, p, rho, deltap, Smin, Smax:
571    deltap, Smin, Smax := compute_abs_rounding_error_firstmult(poly/x,xmax, n-1):
572    rho :=  (2^(-p)*(Smax+deltap) +deltap ) / Smin :
574    deltap, Smin, Smax := compute_abs_rounding_error_firstmult(poly,xmax, n):
728 #    if abs(coeff(poly,x,k)) < xmax*maxSk  # may be improved to xmax*Smax/2
730 #         printf("    coeff=%1.20e,  xmax*Smax=%1.20e"  ,  abs(coeff(poly,x,i)),  xmax*maxSk ):