Exact arithmetic as a tool for convergence assessment of the IRM-CG method

06/14/2019
by   J. Dvornik, et al.
0

Using exact computer arithmetic, it is possible to determine the (exact) solution of a numerical model without rounding error. For such purposes, a corresponding system of equations should be exactly defined, either directly or by rationalisation of numerically given input data. In the latter case there is an initial round off error, but this does not propagate during the solution process. If this system is first exactly solved, then by the floating-point arithmetic, convergence of the numerical method is easily followed. As one example, IRM-CG, a special case of the more general Iterated Ritz method and interesting replacement for a standard or preconditioned CG, is verified. Further, because the computer demands and execution time grow enourmously with the number of unknowns using this strategy, the possibilities for larger systems are also provided.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
09/16/2021

Computing the exact sign of sums of products with floating point arithmetic

IIn computational geometry, the construction of essential primitives lik...
research
02/19/2018

Multithreading for the expression-dag-based number type Real_algebraic

Many algorithms, especially in the field of computational geometry, are ...
research
06/03/2021

When does the Lanczos algorithm compute exactly?

In theory, the Lanczos algorithm generates an orthogonal basis of the co...
research
05/20/2021

Indirect predicates for geometric constructions

Geometric predicates are a basic ingredient to implement a vast range of...
research
01/26/2000

Numerical Replication of Computer Simulations: Some Pitfalls and How To Avoid Them

A computer simulation, such as a genetic algorithm, that uses IEEE stand...
research
10/30/2017

VLSI Computational Architectures for the Arithmetic Cosine Transform

The discrete cosine transform (DCT) is a widely-used and important signa...
research
04/28/2020

A posteriori error estimation for the non-self-consistent Kohn-Sham equations

We address the problem of bounding rigorously the errors in the numerica...

Please sign up or login with your details

Forgot password? Click here to reset