This article is a case study in the implementation of a portable, proven and efficient correctly rounded elementary function in double-precision. We describe the methodology used to achieve these goals in the crlibm library. There are two novel aspects to this approach. The first is the proof framework, and in general the techniques used to balance performance and provability. The second is the introduction of processor-specific optimization to get performance equivalent to the best current...

The Treatise on Quadratureof Fermat (c. 1659), besides containing the first known proof of the computation of the area under a higher parabola, $\int x^{+m/n}dx$, or under a higher hyperbola, $\int x^{-m/n}dx$—with the appropriate limits of integration in each case—has a second part which was mostly unnoticed by Fermat’s contemporaries. This second part of theTreatise is obscure and difficult to read. In it Fermat reduced the quadrature of a great number of algebraic curves in implicit form to the quadrature of known curves: the...

In the paper we give an analogue of the Filippov Lemma for the second order differential inclusions with the initial conditions y(0) = 0, y′(0) = 0, where the matrix A ∈ ℝd×d and multifunction is Lipschitz continuous in y with a t-independent constant l. The main result is the following: Assume that F is measurable in t and integrably bounded. Let y 0 ∈ W 2,1 be an arbitrary function fulfilling the above initial conditions and such that where p 0 ∈ L 1[0, 1]. Then there exists a solution y ∈ W 2,1...

We first prove that given any analytic filter ℱ on ω the set of all functions f on $2^{\omega }$ which can be represented as the pointwise limit relative to ℱ of some sequence ${\left(fₙ\right)}_{n\in \omega }$ of continuous functions ($f=li{m}_{ℱ}fₙ$), is exactly the set of all Borel functions of class ξ for some countable ordinal ξ that we call the rank of ℱ. We discuss several structural properties of this rank. For example, we prove that any free Π⁰₄ filter is of rank 1.

A new theorem in the theory of first return representations of Baire class one functions of a real variable is presented which has as immediate consequences several known characterizations of standard subclasses of the Baire one functions. Further, this theorem yields new insights into how finely Baire one functions can be recovered and yields a characterization of another subclass of Baire one functions.

In this paper we study a model problem describing the movement of a glacier under Glen’s flow law and investigated by Colinge and Rappaz [Colinge and Rappaz, ESAIM: M2AN 33 (1999) 395–406]. We establish error estimates for finite element approximation using the results of Chow [Chow, SIAM J. Numer. Analysis 29 (1992) 769–780] and Liu and Barrett [Liu and Barrett, SIAM J. Numer. Analysis 33 (1996) 98–106] and give an analysis of the convergence of the successive approximations used in [Colinge and...

In this paper we study a model problem describing the movement of a glacier under Glen's flow law and investigated by Colinge and Rappaz [Colinge and Rappaz, ESAIM: M2AN33 (1999) 395–406]. We establish error estimates for finite element approximation using the results of Chow [Chow, SIAM J. Numer. Analysis29 (1992) 769–780] and Liu and Barrett [Liu and Barrett, SIAM J. Numer. Analysis33 (1996) 98–106] and give an analysis of the convergence of the successive approximations used in [Colinge and...

We discuss two ways to construct standard probability measures, called push-down measures, from internal probability measures. We show that the Wasserstein distance between an internal probability measure and its push-down measure is infinitesimal. As an application to standard probability theory, we show that every finitely-additive Borel probability measure $P$ on a separable metric space is a limit of a sequence of countably-additive Borel probability measures ${\left\{P_n\right\}}_{n\in \mathbb{N}}$ in the sense that $\int f\mathrm{d}P={lim}_{n\to \infty }\int f\mathrm{d}{P}_n$ for all bounded...

Let $T_{\gamma }f\left(x\right)={ʃ}_0^{x}k{(x,y)}^{\gamma }f\left(y\right)dy$, where k is a nonnegative kernel increasing in x, decreasing in y, and satisfying a triangle inequality. An nth-order Opial inequality has the form ${ʃ}_0^{\infty }({\prod }_{i=1}^n|T_{{\gamma }_i}{f\left(x\right)|}^{q_i}{\left|\right)\left|f\left(x\right)\right|}^{q_0}w\left(x\right)dx\le C({ʃ}_0^{\infty }{\left|f\left(x\right)\right|}^p{v\left(x\right)dx)}^{(q_0+\dots +q_n)/p}$. Such inequalities can always be simplified to nth-order reduced inequalities, where the exponent $q_0=0$. When n = 1, the reduced inequality is a standard weighted norm inequality, and characterizing the weights is easy. We also find necessary and sufficient conditions on the weights for second-order reduced Opial inequalities to hold.

First order characterizations of pseudoconvex functions are investigated in terms of generalized directional derivatives. A connection with the invexity is analysed. Well-known first order characterizations of the solution sets of pseudolinear programs are generalized to the case of pseudoconvex programs. The concepts of pseudoconvexity and invexity do not depend on a single definition of the generalized directional derivative.