Extrapolation | EPFL Graph Search

In mathematics, extrapolation is a type of estimation, beyond the original observation range, of the value of a variable on the basis of its relationship with another variable. It is similar to interpolation, which produces estimates between known observations, but extrapolation is subject to greater uncertainty and a higher risk of producing meaningless results. Extrapolation may also mean extension of a method, assuming similar methods will be applicable. Extrapolation may also apply to human experience to project, extend, or expand known experience into an area not known or previously experienced so as to arrive at a (usually conjectural) knowledge of the unknown (e.g. a driver extrapolates road conditions beyond his sight while driving). The extrapolation method can be applied in the interior reconstruction problem. Methods A sound choice of which extrapolation method to apply relies on a priori knowledge of the process that created the existing data points. Some ex

Numerical analysis of optimization-constrained differential equations

Chantal Landry

The modeling of a system composed by a gas phase and organic aerosol particles, and its numerical resolution are studied. The gas-aerosol system is modeled by ordinary differential equations coupled with a mixed-constrained optimization problem. This coupling induces discontinuities when inequality constraints are activated or deactivated. Two approaches for the solution of the optimization-constrained differential equations are presented. The first approach is a time splitting scheme together with a fixed-point method that alternates between the differential and optimization parts. The ordinary differential equations are approximated by the Crank-Nicolson scheme and a primal-dual interior-point method combined with a warm-start strategy is used to solve the minimization problem. The second approach considers the set of equations as a system of differential algebraic equations after replacing the minimization problem by its first order optimality conditions. An implicit 5th-order Runge-Kutta method (RADAU5) is then used. Both approaches are completed by numerical techniques for the detection and computation of the events (activation and deactivation of inequality constraints) when the system evolves in time. The computation of the events is based on continuation techniques and geometric arguments. Moreover the first approach completes the computation with extrapolation polynomials and sensitivity analysis, whereas the second approach uses dense output formulas. Numerical results for gas-aerosol system made of several chemical species are proposed for both approaches. These examples show the efficiency and accuracy of each method. They also indicate that the second approach is more efficient than the first one. Furthermore theoretical examples show that the method for the computation of the activation is of second order for the first approach and exact for the second one.

EPFL2009

A higher-order numerical framework for stochastic simulation of chemical reaction systems

Background: In this paper, we present a framework for improving the accuracy of fixed-step methods for Monte Carlo simulation of discrete stochastic chemical kinetics. Stochasticity is ubiquitous in many areas of cell biology, for example in gene regulation, biochemical cascades and cell-cell interaction. However most discrete stochastic simulation techniques are slow. We apply Richardson extrapolation to the moments of three fixed-step methods, the Euler, midpoint and theta-trapezoidal tau-leap methods, to demonstrate the power of stochastic extrapolation. The extrapolation framework can increase the order of convergence of any fixed-step discrete stochastic solver and is very easy to implement; the only condition for its use is knowledge of the appropriate terms of the global error expansion of the solver in terms of its stepsize. In practical terms, a higher-order method with a larger stepsize can achieve the same level of accuracy as a lower-order method with a smaller one, potentially reducing the computational time of the system. Results: By obtaining a global error expansion for a general weak first-order method, we prove that extrapolation can increase the weak order of convergence for the moments of the Euler and the midpoint tau-leap methods, from one to two. This is supported by numerical simulations of several chemical systems of biological importance using the Euler, midpoint and theta-trapezoidal tau-leap methods. In almost all cases, extrapolation results in an improvement of accuracy. As in the case of ordinary and stochastic differential equations, extrapolation can be repeated to obtain even higher-order approximations. Conclusions: Extrapolation is a general framework for increasing the order of accuracy of any fixed-step stochastic solver. This enables the simulation of complicated systems in less time, allowing for more realistic biochemical problems to be solved.

Biomed Central Ltd2012

Efficient algorithms for computing Sommerfeld integral tails

Ruzica Golubovic, Juan Ramon Mosig

Sommerfeld-integrals (SIs) are ubiquitous in the analysis of problems involving antennas and scatterers embedded in planar multilayered media. It is well known that the oscillating and slowly decaying nature of their integrands makes the numerical evaluation of the SI real-axis tail segment a very time consuming and computationally expensive task. Therefore, SI tails have to be specially treated. In this paper we compare two recently developed techniques for their efficient numerical evaluation. First, a partition-extrapolation method, in which the integration-then-summation procedure is combined with a new version of the weighted averages (WA) extrapolation technique, is summarized. The previous variants of WA technique are also discussed. Then, a review of double-exponential (DE) quadrature formulas for direct integration of the SI tails is presented. The efficient way of implementing the algorithms, their pros and cons, as well as comparisons of their performance are discussed in detail.

Institute of Electrical and Electronics Engineers2012

A higher-order numerical framework for stochastic simulation of chemical reaction systems

Biomed Central Ltd2012