simeq_accuracy.c running, big range of values A*X=Y, given X, generate random A, compute Y solve for computed Xc, compare accuracy to X n=6, xbig=5.000000e+01 Amin=-57.429023, Amax=58.586045, Aavg=4.140372 0 comp=0.000000 true=0.000000 err=8.975459e-15 1 comp=1.000000 true=1.000000 err=2.531308e-14 2 comp=2.000000 true=2.000000 err=1.088019e-14 3 comp=5.000000 true=5.000000 err=7.993606e-15 4 comp=10.000000 true=10.000000 err=7.105427e-15 5 comp=50.000000 true=50.000000 err=1.421085e-14 n=10, xbig=1.000000e+03 Amin=-1018.008496, Amax=986.364585, Aavg=-126.394554 0 comp=0.000000 true=0.000000 err=4.003464e-13 1 comp=1.000000 true=1.000000 err=1.255884e-12 2 comp=2.000000 true=2.000000 err=4.694023e-13 3 comp=5.000000 true=5.000000 err=7.789325e-13 4 comp=10.000000 true=10.000000 err=1.598721e-14 5 comp=50.000000 true=50.000000 err=6.110668e-13 6 comp=100.000000 true=100.000000 err=1.989520e-13 7 comp=200.000000 true=200.000000 err=3.126388e-13 8 comp=500.000000 true=500.000000 err=8.526513e-13 9 comp=1000.000000 true=1000.000000 err=3.410605e-13 n=19, xbig=1.000000e+06 Amin=-521194.670897, Amax=523203.445800, Aavg=-16576.555982 0 comp=-0.000000 true=0.000000 err=2.248726e-10 1 comp=1.000000 true=1.000000 err=1.465423e-10 2 comp=2.000000 true=2.000000 err=1.612079e-10 3 comp=5.000000 true=5.000000 err=2.747242e-10 4 comp=10.000000 true=10.000000 err=3.660716e-11 5 comp=50.000000 true=50.000000 err=1.478497e-10 6 comp=100.000000 true=100.000000 err=4.354206e-11 7 comp=200.000000 true=200.000000 err=3.377920e-10 8 comp=500.000000 true=500.000000 err=5.098286e-10 9 comp=1000.000000 true=1000.000000 err=1.375611e-11 10 comp=2000.000000 true=2000.000000 err=6.320988e-11 11 comp=5000.000000 true=5000.000000 err=1.828084e-10 12 comp=10000.000000 true=10000.000000 err=2.837623e-10 13 comp=20000.000000 true=20000.000000 err=1.418812e-10 14 comp=50000.000000 true=50000.000000 err=8.731149e-11 15 comp=100000.000000 true=100000.000000 err=1.309672e-10 16 comp=200000.000000 true=200000.000000 err=5.820766e-11 17 comp=500000.000000 true=500000.000000 err=1.164153e-10 18 comp=1000000.000000 true=1000000.000000 err=1.164153e-10 n=25, xbig=1.000000e+08 Amin=-33350494.093655, Amax=33321450.718642, Aavg=-319120.684824 0 comp=0.000000 true=0.000000 err=1.030494e-07 1 comp=1.000000 true=1.000000 err=9.500764e-08 2 comp=2.000000 true=2.000000 err=2.456227e-08 3 comp=5.000000 true=5.000000 err=4.662579e-07 4 comp=10.000000 true=10.000000 err=2.574525e-07 5 comp=50.000000 true=50.000000 err=2.034431e-07 6 comp=100.000001 true=100.000000 err=1.065200e-06 7 comp=200.000000 true=200.000000 err=3.996101e-07 8 comp=500.000000 true=500.000000 err=5.084075e-08 9 comp=999.999999 true=1000.000000 err=5.898473e-07 10 comp=2000.000000 true=2000.000000 err=1.719600e-07 11 comp=5000.000000 true=5000.000000 err=4.862668e-07 12 comp=10000.000000 true=10000.000000 err=2.211018e-07 13 comp=20000.000000 true=20000.000000 err=1.985100e-07 14 comp=50000.000000 true=50000.000000 err=4.325993e-07 15 comp=100000.000000 true=100000.000000 err=1.549779e-08 16 comp=200000.000000 true=200000.000000 err=1.569279e-07 17 comp=500000.000000 true=500000.000000 err=3.022142e-07 18 comp=1000000.000000 true=1000000.000000 err=3.697351e-07 19 comp=2000000.000000 true=2000000.000000 err=2.789311e-07 20 comp=5000000.000001 true=5000000.000000 err=5.746260e-07 21 comp=9999999.999999 true=10000000.000000 err=7.823110e-07 22 comp=20000000.000000 true=20000000.000000 err=1.192093e-07 23 comp=50000000.000001 true=50000000.000000 err=5.289912e-07 24 comp=100000000.000000 true=100000000.000000 err=3.129244e-07 n=31, xbig=1.000000e+10 Amin=-2146930723.999743, Amax=2143458001.998125, Aavg=-30591370.223402 0 comp=0.000008 true=0.000000 err=8.165836e-06 1 comp=0.999989 true=1.000000 err=1.070090e-05 2 comp=1.999998 true=2.000000 err=2.227258e-06 3 comp=4.999994 true=5.000000 err=5.573616e-06 4 comp=9.999997 true=10.000000 err=2.531800e-06 5 comp=49.999988 true=50.000000 err=1.218962e-05 6 comp=99.999992 true=100.000000 err=8.279458e-06 7 comp=200.000006 true=200.000000 err=5.915761e-06 8 comp=500.000003 true=500.000000 err=3.285415e-06 9 comp=1000.000014 true=1000.000000 err=1.411466e-05 10 comp=1999.999989 true=2000.000000 err=1.052860e-05 11 comp=4999.999996 true=5000.000000 err=4.136702e-06 12 comp=9999.999993 true=10000.000000 err=7.467519e-06 13 comp=19999.999995 true=20000.000000 err=5.138456e-06 14 comp=50000.000006 true=50000.000000 err=6.195623e-06 15 comp=99999.999998 true=100000.000000 err=2.425048e-06 16 comp=200000.000014 true=200000.000000 err=1.371745e-05 17 comp=499999.999993 true=500000.000000 err=6.523682e-06 18 comp=1000000.000003 true=1000000.000000 err=3.083958e-06 19 comp=2000000.000004 true=2000000.000000 err=4.045898e-06 20 comp=5000000.000005 true=5000000.000000 err=4.520640e-06 21 comp=9999999.999992 true=10000000.000000 err=7.884577e-06 22 comp=19999999.999994 true=20000000.000000 err=6.008893e-06 23 comp=50000000.000001 true=50000000.000000 err=1.281500e-06 24 comp=99999999.999991 true=100000000.000000 err=8.687377e-06 25 comp=200000000.000001 true=200000000.000000 err=5.960464e-07 26 comp=500000000.000003 true=500000000.000000 err=2.563000e-06 27 comp=999999999.999989 true=1000000000.000000 err=1.108646e-05 28 comp=1999999999.999992 true=2000000000.000000 err=7.867813e-06 29 comp=4999999999.999994 true=5000000000.000000 err=5.722046e-06 30 comp=10000000000.000021 true=10000000000.000000 err=2.098083e-05 note that some of the error comes from multiplying A*X end simeq_accuracy.c