Your first, really simple, huge error is your
"The backward error remains and those powers of 10 depend on that error critically."
You are totally, 100% WRONG. All you had to do was just READ what the heck I wrote. So far you have let your arrogance and resistance at believing something new to you keep you from just elementary READING.
So, return and just READ the post.
Or again, yet one more time, just for you this time, to repeat, to pound the head of the nail a half inch below the surface of the wood, for some positive integer n, consider a system of n linear equations in n unknowns. To this end, let A be an n x n real matrix, x be an n x 1 real vector, and let b be an n x 1 real vector.
So, the system to be solved is just Ax = b. Here we seek x.
Matrix A is assumed to be non-singular but otherwise may be horribly 'ill-conditioned'.
Assume that the numerical data is given in base 10.
Then for some positive integer m, all the components of
(10^m)A and (10^m)b
are just whole numbers. It is also permitted to pick a different m for each of the n equations, that is, each row of Ax = b.
This is the second time I've made this very elementary point. The point is correct, 100% rock solidly correct. All that is needed for the technique is that the components of
(10^m)A and (10^m)b
be whole numbers, and clearly, with total crystal clarity, the 'condition number' of A is not, Not, NOT
N O T
involved.
Twice now you've missed the VERY simple point. TWICE. And it's DIRT SIMPLE.
You continued:
"The backward error remains and those powers of 10 depend on that error critically."
Wrong. Totally wrong. Uninformed, misinformed, just plain wrong, badly wrong, totally off the subject wrong.
There is no such "depend" at all. None. Zip, zilch, zero.
If you are sober and paying any attention at all, which so far is doubtful, then you have just ASSUMED that you understand the technique and drawn totally wrong conclusions.
Maybe you just bitterly resent the idea of learning something new to you.
With 99% probability, you have missed the technique. Just MISSED it. You just didn't 'get it'. It didn't register with you. You totally missed it.
Once again, yet again, over again, just for you because everyone else in this thread long since 'got it', given the system of equations where all the given data is integers, we pick a list of single precision integer prime numbers. There are a LOT of such prime numbers, far, far, far more than there is any danger of our needing.
Then for each such prime on the list, consider the 'algebraic field' of the integers modulo that prime. In that field, the GCD algorithm is JUST how to find a multiplicative inverse.
Now, so armed with the GCD algorithm, use Gaussian elimination, essentially as usual, to solve the system of linear equations in the field of the integers modulo that prime. In this work, will need essentially only single precision integer arithmetic. The result will be just single precision integers.
A large condition number does NOT make anything about this use of Gaussian elimination more difficult.
Do this work for each prime on the list.
Then using the Chinese remainder theorem, we can assemble the exact, rational answer from the results from the primes. Of course the solution will be in the rationals, that is, a ratio of integers. In this case, the integers will be multiple precision. Finding these integers is nearly the only multiple precision arithmetic needed by the technique and is routine.
This is NOT a joke. It's rock solid.
It is true that we need an estimate of the number of prime numbers to have on our list so that the result of the Chinese remainder theorem will be correct and not too small, but this estimate is easy and, for anything like reasonable input data, is comfortably small, typically no more than a few dozen.
As I warned you, I will provide rock solid support for my claims at which time your uninformed objections will make you look like a fool before HN.
Actually what I am describing was worked out by M. Newman in the 1960s at the old National Bureau of Standards. He wrote a paper on the work:
Morris Newman, "Solving Equations Exactly," 'Journal of Research of the National Bureau of Standards -- B. Mathematics and Mathematical 'Physics', volume 71B, number 4, October-December, 1967, pages 171-179.
That is peer-reviewed mathematics from a quite good and fully competent mathematician.
He also wrote software to implement the technique and ran it on a Univac 1108 computer. It worked GREAT.
His paper shows how to find a sufficient list of prime numbers.
You have shown all of HN that you think like a FOOL.
Apparently you don't mind acting like a fool in public.
It is difficult to fathom how someone sober could have learned about condition number and backward error analysis yet fail to see the middle school point that
(10^m)A and (10^m)b
can be all integers, presented, now, TWICE.