Algorithm 5.4.1 (Preconditioned CG-Method for solving LSE with p.d. matrix)
A= [[ 10. 1. 0. ..., 0. 0. 0.] [ 1. 10. 1. ..., 0. 0. 0.] [ 0. 1. 10. ..., 0. 0. 0.] ..., [ 0. 0. 0. ..., 10. 1. 0.] [ 0. 0. 0. ..., 1. 10. 1.] [ 0. 0. 0. ..., 0. 1. 10.]] b= [ 11. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 11.] x= [ 1.00000001 0.99999998 1. 1. 1.00000001 1. 1. 1.00000001 1. 1. 1.00000001 1. 1. 0.99999998 1. 1. 1. 1. 1.00000001 0.99999998 1.00000001 1. 1. 1. 1. 1.00000001 1. 1.00000001 1.00000001 1. 1. 1. 0.99999999 1. 1.00000001 1. 1.00000001 0.99999999 1.00000001 0.99999999 1. 1.00000001 1. 1.00000001 1. 1.00000002 1. 1. 1.00000001 1. 1.00000001 1. 1. 1. 1. 0.99999999 1.00000001 0.99999999 1. 1.00000001 1. 0.99999999 1. 0.99999999 1. 1. 1. 0.99999999 1.00000001 1.00000001 0.99999999 1. 1. 1. 1.00000001 1.00000001 1. 1. 1. 1. 1. 0.99999999 1. 1. 1. 1. 1.00000001 1.00000002 0.99999999 1.00000001 0.99999998 1. 1. 1. 1.00000001 0.99999998 0.99999999 0.99999999 1. 0.99999999] ||r||= 1.91648814166e-07 9 A= [[ 10. 1. 0. ..., 0. 0. 0.] [ 1. 10. 1. ..., 0. 0. 0.] [ 0. 1. 10. ..., 0. 0. 0.] ..., [ 0. 0. 0. ..., 10. 1. 0.] [ 0. 0. 0. ..., 1. 10. 1.] [ 0. 0. 0. ..., 0. 1. 10.]] b= [ 11. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 12. 11.] x= [ 1.00000001 0.99999998 1. 1. 1.00000001 1. 1. 1.00000001 1. 1. 1.00000001 1. 1. 0.99999998 1. 1. 1. 1. 1.00000001 0.99999998 1.00000001 1. 1. 1. 1. 1.00000001 1. 1.00000001 1.00000001 1. 1. 1. 0.99999999 1. 1.00000001 1. 1.00000001 0.99999999 1.00000001 0.99999999 1. 1.00000001 1. 1.00000001 1. 1.00000002 1. 1. 1.00000001 1. 1.00000001 1. 1. 1. 1. 0.99999999 1.00000001 0.99999999 1. 1.00000001 1. 0.99999999 1. 0.99999999 1. 1. 1. 0.99999999 1.00000001 1.00000001 0.99999999 1. 1. 1. 1.00000001 1.00000001 1. 1. 1. 1. 1. 0.99999999 1. 1. 1. 1. 1.00000001 1.00000002 0.99999999 1.00000001 0.99999998 1. 1. 1. 1.00000001 0.99999998 0.99999999 0.99999999 1. 0.99999999] ||r||= 1.91648814166e-07 9 |
|