Robert H. C. Moir

Robert M. Corless, C. Yalcin Kaya, Robert H. C. Moir

We show how to compute the \emph{optimal relative backward error} for the numerical solution of the Dahlquist test problem by one-step methods. This is an example of a general approach that uses results from optimal control theory to compute optimal residuals, but elementary methods can also be used here because the problem is so simple. This analysis produces new insight into the numerical solution of stiff problems.

Robert M. Corless, Robert H. C. Moir, Marc Moreno Maza et al.

We consider the problem of symbolic-numeric integration of symbolic functions, focusing on rational functions. Using a hybrid method allows the stable yet efficient computation of symbolic antiderivatives while avoiding issues of ill-conditioning to which numerical methods are susceptible. We propose two alternative methods for exact input that compute the rational part of the integral using Hermite reduction and then compute the transcendental part two different ways using a combination of exact integration and efficient numerical computation of roots. The symbolic computation is done within BPAS, or Basic Polynomial Algebra Subprograms, which is a highly optimized environment for polynomial computation on parallel architectures, while the numerical computation is done using the highly optimized multiprecision rootfinding package MPSolve. We show that both methods are forward and backward stable in a structured sense and away from singularities tolerance proportionality is achieved by adjusting the precision of the rootfinding tasks.