Abstract:One of the three main processes associated with polynomials is evaluation; the two other ones being interpolation and root finding. Higham [1, chap. 5] devotes an entire chapter to polynomials and more especially to polynomial evaluation. The small backward error the Horner scheme introduce when evaluated in floating point arithmetic justifies its practical interest. It is well known that the computed evaluation of p(x) is the exact value at x of a polynomial obtained by making relative perturbations of at most size 2nu to the coefficients of p where n denotes the polynomial degree and u the finite precision of the computation [1]. Nevertheless, the computed result can be arbitrary less accurate than the working precision u when evaluating p(x) is ill-conditioned. This is the case for example in the neighborhood of multiple roots where all the digits or even the order of the computed value of p(x) can be false. Numerous multiprecision libraries are available when the computing precision u is not sufficient to guarantee a prescribed accuracy. Fixed-length expansions such as in “double-double” libraries [2] are actual and effective solutions to simulate twice the IEEE-754 double precision. These fixed-length expansions are currently embedded in major developments such as the new extended and mixed precision BLAS [2]. On the other hand, an approach similar to Kahan’s compensated summation (see [1, p.83–88]) can be used, such as in [3] where efficient algorithms are proposed for the computation of summations and dot products in a specified precision. In this paper we present an accurate and fast compensated algorithm to evaluate univariate polynomials in floating point arithmetic. By accurate, we mean that the accuracy of our compensated Horner scheme is similar to the one given by the Horner scheme computed in doubled working precision, for example using fixed-length expansions. Of course the accuracy of the result still depends on the condition number. By fast, we mean that the algorithm should run at least twice as fast as the fixed-length expansion counterparts that produce the same output accuracy. We constraint the computation to be performed in the same precision as the data, for example the IEEE-754 double precision. Thanks to an a priori error analysis, we provide a sufficiant criterium on the condition number of the polynomial evaluation to ensure that our compensated Horner scheme provides a faithful rounding of the exact evaluation. By faithful rounding we mean that the computed result is one of the two floating point neighbours of the exact result. We also derive a validated running error analysis of our algorithm, which leads to a sharper enclosure of the computed result. In particular, our numerical experiments show that the computed evaluation can be prooved to be faithfully rounded at the running time, as long as the condition number is smaller than the inverse of the working precision.

Crystallization and preliminary X-ray characterization of VanA from Enterococcus faecium BM4147: towards the molecular basis of bacterial resistance to the glycopeptide antibiotic vancomycin.

Strategy for compression of medical image in real time with bit-plane runlength coding

Male circumcision and HIV status among Latino immigrant MSM in New York City.

CompaCT: Fractal-Based Heuristic Pixel Segmentation for Lossless Compression of High-Color DICOM Medical Images

A Lossless Compression Method for Multi-Component Medical Images Based on Big Data Mining

An optimized JPEG-XT-based algorithm for the lossy and lossless compression of 16-bit depth medical image

Batch Images Compression Algorithm Based on the Common Features

Machine Learning Methods in Medical Image Compression

Development of a Lossless Compression Algorithm Using Content-Based Classification of Image Block for Wireless Capsule Endoscopic Image

Low-Complexity and High-Efficiency Image Compression Algorithm for Wireless Endoscopy System

Study on Hybrid Compression Algorithm of CT Images

Joint Reconstruction of Dynamic PET Activity and Kinetic Parametric Images Using Total Variation Constrained Dictionary Sparse Coding

Accurate Polynomial Evaluation in Floating Point Arithmetic

Ueber einen Fall von retrograder Amnesie nach traumatischer Epilepsie

Lossy Image Compression in a Preclinical Multimodal Imaging Study

Dynamic PET cardiac and parametric image reconstruction: a fixed-point proximity gradient approach using patch-based DCT and tensor SVD regularization

Regularized Compression of MRI Data: Modular Optimization of Joint Reconstruction and Coding

Extreme MRI: Large-Scale Volumetric Dynamic Imaging from Continuous Non-Gated Acquisitions

Learning Lossless Compression for High Bit-Depth Volumetric Medical Image

Nearly-lossless-to-lossy medical image compression by the optimized JPEGXT and JPEG algorithms through the anatomical regions of interest

Content-based retrieval of dynamic PET functional images