C 代码 实现共轭梯度 （CG） 方法求解 一个正定稀疏线性系统 A x x=b.rar

19 January 2020 09:01:39 AM CG_RC_TEST: C version Test the CG_RC library. TEST01 Use CG_RC on the 1, -2, 1 matrix. Number of iterations was 9 Estimated error is 4.21596e-16 Loo error is 6.66134e-16 I X(I) X_EXACT(I) B(I) 0 7.80626e-17 0 0.309017 1 0.309017 0.309017 -0.0302487 2 0.587785 0.587785 -0.0575365 3 0.809017 0.809017 -0.0791922 4 0.951057 0.951057 -0.093096 5 1 1 -0.097887 6 0.951057 0.951057 -0.093096 7 0.809017 0.809017 -0.0791922 8 0.587785 0.587785 -0.0575365 9 0.309017 0.309017 -0.0302487 10 -4.79828e-17 1.22465e-16 3.88578e-16 11 -0.309017 -0.309017 0.0302487 12 -0.587785 -0.587785 0.0575365 13 -0.809017 -0.809017 0.0791922 14 -0.951057 -0.951057 0.093096 15 -1 -1 0.097887 16 -0.951057 -0.951057 0.093096 17 -0.809017 -0.809017 0.0791922 18 -0.587785 -0.587785 0.0575365 19 -0.309017 -0.309017 0.0302487 20 -4.10262e-16 -2.44929e-16 -0.309017 TEST02 Use CG_RC to solve a linear system involving the Wathen matrix. NX = 5 NY = 4 N = 79 Number of iterations was 18 Estimated error is 0.00730674 Loo error is 4.22526e-05 I X(I) X_EXACT(I) B(I) 0 0.218432 0.218418 -19.3113 1 0.956316 0.956318 88.7317 2 0.829512 0.829509 -36.3017 3 0.561714 0.561695 90.0773 4 0.415273 0.415307 -43.7437 5 0.066156 0.0661187 54.0011 6 0.25756 0.257578 -23.2089 7 0.109949 0.109957 32.2168 8 0.0438096 0.043829 -33.6528 9 0.633973 0.633966 65.2136 10 0.06171 0.0617272 -13.8091 11 0.449528 0.449539 70.363 12 0.401292 0.401306 154.809 13 0.754653 0.754673 168.305 14 0.797282 0.797287 118.754 15 0.00182286 0.00183837 72.6417 16 0.897471 0.897504 74.6033 17 0.350743 0.350752 -37.7129 18 0.0945624 0.0945448 107.05 19 0.0136232 0.0136169 -109.428 20 0.859112 0.859097 207.905 21 0.840875 0.840847 -64.218 22 0.123101 0.123104 99.2914 23 0.00751016 0.00751236 -67.2724 24 0.260308 0.260303 101.661 25 0.912477 0.912484 -42.2354 26 0.113677 0.113664 99.6324 27 0.351621 0.351629 -30.9739 28 0.822884 0.822887 78.3529 29 0.26713 0.267132 144.803 30 0.692067 0.692066 166.674 31 0.56167 0.561662 111.98 32 0.86121 0.861216 146.417 33 0.453772 0.453794 49.2886 34 0.911983 0.911977 -10.5401 35 0.597914 0.597917 114.666 36 0.18896 0.188955 -79.7521 37 0.761504 0.761492 176.378 38 0.396996 0.396988 -81.2853 39 0.185316 0.185314 108.356 40 0.574347 0.574366 -72.6881 41 0.367041 0.367027 156.156 42 0.61722 0.617205 -72.2069 43 0.361536 0.361529 124.845 44 0.212947 0.21293 -41.1752 45 0.714457 0.714471 46.6628 46 0.117699 0.117707 83.0738 47 0.299328 0.299329 131.474 48 0.824989 0.825003 189.243 49 0.824647 0.82466 191.901 50 0.0618704 0.0618618 62.9693 51 0.710738 0.710781 -1.73137 52 0.0882712 0.0882833 31.9862 53 0.77798 0.777994 -19.4569 54 0.745317 0.745303 131.412 55 0.308673 0.308675 -61.8462 56 0.899362 0.899373 151.142 57 0.763542 0.763537 -50.183 58 0.761737 0.761731 163.776 59 0.406984 0.40697 -88.3359 60 0.938776 0.938749 202.851 61 0.562085 0.562088 -47.2525 62 0.0178467 0.01782 -2.62116 63 0.501099 0.501103 66.5544 64 0.0418957 0.0419093 68.9155 65 0.368859 0.368851 96.42 66 0.271712 0.271724 138.43 67 0.858573 0.858573 114.962 68 0.0290067 0.0290366 -1.53808 69 0.0174254 0.0174423 -1.44348 70 0.152396 0.152384 -11.5233 71 0.114317 0.114319 32.6938 72 0.353917 0.353907 -24.8832 73 0.119309 0.119308 32.1821 74 0.206644 0.206653 -28.0064 75 0.212946 0.212924 38.958 76 0.612979 0.612948 -40.6021 77 0.809533 0.809519 107.951 78 0.587116 0.58709 -27.9974 CG_RC_TEST: Normal end of execution. 19 January 2020 09:01:39 AM