| Project Detail |
Transforming science with faster, more reliable solvers of differential equations
Differential equations are ubiquitous in science and technology, often requiring approximate numeric or symbolic solutions. The ERC-funded ODELIX project aims to enhance the speed and reliability of solving these equations. To this end, researchers will analyse computational costs and develop new theories using advanced computer algebra techniques like numerical homotopy continuation and sparse interpolation. The team will create open-source implementations and high-performance libraries, such as medium-precision arithmetics and reliable homotopy continuation, to validate their approach. These innovations should allow solving complex problems currently beyond reach.
Being a language of nature, differential equations are ubiquitous in science and technology. Solving them is a fundamental computational task with a long and rich history. Applications usually require approximate solutions, which can be computed using numerical methods such as Runge-Kutta schemes. Alternatively, one may search for symbolic solutions, which have the advantage of presenting the solutions in an exact and more intelligible way. However, such solutions do not always exist and may be hard to compute.
The present proposal aims at making the resolution of differential equations both faster and more reliable. We will undertake a systematic analysis of the cost to compute both numeric and symbolic solutions, as a function of the required precision, special properties of the equation and its solutions, and hardware specifics of the computer. This includes the cost to certify approximate numeric solutions, e.g. through the computation of provable error bounds. In order to compute symbolic solutions more efficiently, we will develop a new theory that relies on two techniques from computer algebra that were improved significantly in the past decade: numerical homotopy continuation and sparse interpolation.
Theoretical progress on the above problems will be accompanied by open source implementations. For this purpose, we will also implement several high performance libraries of independent interest: non-conventional medium precision arithmetic, reliable homotopy continuation, sparse interpolation, faster-than-just-in-time compilation, etc. Altogether, these implementations will validate the correctness and efficiency of our approach. They should also allow us to tackle problems from applications that are currently out of reach. |
