Symbolic-Numeric Integration of Rational Functions

12/05/2017
by   Robert M. Corless, et al.
0

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.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/13/2022

Stability Problems in Symbolic Integration

This paper aims to initialize a dynamical aspect of symbolic integration...
research
12/31/2020

calculus: High Dimensional Numerical and Symbolic Calculus in R

The R package calculus implements C++ optimized functions for numerical ...
research
06/06/2017

Symbolic Multibody Methods for Real-Time Simulation of Railway Vehicles

In this work, recently developed state-of-the-art symbolic multibody met...
research
03/30/2021

Structural backward stability in rational eigenvalue problems solved via block Kronecker linearizations

We study the backward stability of running a backward stable eigenstruct...
research
06/27/2017

Symbolic Versus Numerical Computation and Visualization of Parameter Regions for Multistationarity of Biological Networks

We investigate models of the mitogenactivated protein kinases (MAPK) net...
research
05/09/2018

Generalized Hermite Reduction, Creative Telescoping and Definite Integration of D-Finite Functions

Hermite reduction is a classical algorithmic tool in symbolic integratio...
research
12/19/2017

Computing effectively stabilizing controllers for a class of nD systems

In this paper, we study the internal stabilizability and internal stabil...

Please sign up or login with your details

Forgot password? Click here to reset