07-Jan-2022 17:59:22 CG_RC_TEST: MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2. Test CG_RC. CG_RC_TEST01 Use CG_RC on the 1, -2, 1 matrix. Number of iterations was 9 Estimated error is 2.33739e-16 Loo error is 6.66134e-16 I X(I) X_EXACT(I) B(I) 1 -0.000000 0.000000 0.309017 2 0.309017 0.309017 -0.030249 3 0.587785 0.587785 -0.057537 4 0.809017 0.809017 -0.079192 5 0.951057 0.951057 -0.093096 6 1.000000 1.000000 -0.097887 7 0.951057 0.951057 -0.093096 8 0.809017 0.809017 -0.079192 9 0.587785 0.587785 -0.057537 10 0.309017 0.309017 -0.030249 11 -0.000000 0.000000 0.000000 12 -0.309017 -0.309017 0.030249 13 -0.587785 -0.587785 0.057537 14 -0.809017 -0.809017 0.079192 15 -0.951057 -0.951057 0.093096 16 -1.000000 -1.000000 0.097887 17 -0.951057 -0.951057 0.093096 18 -0.809017 -0.809017 0.079192 19 -0.587785 -0.587785 0.057537 20 -0.309017 -0.309017 0.030249 21 -0.000000 -0.000000 -0.309017 CG_RC_TEST02 Use CG_RC to solve a linear system involving the Wathen matrix. NX = 5 NY = 4 N = 79 ans = 79 1 ans = 79 79 Number of iterations was 19 Estimated error is 0.00411276 Loo error is 4.02795e-05 I X(I) X_EXACT(I) B(I) 1 0.814748 0.814724 -34.224530 2 0.905785 0.905792 132.363927 3 0.127005 0.126987 -76.715606 4 0.913372 0.913376 124.562270 5 0.632366 0.632359 -52.813244 6 0.097535 0.097540 64.883778 7 0.278498 0.278498 -40.161932 8 0.546867 0.546882 58.340911 9 0.957495 0.957507 -27.894679 10 0.964880 0.964889 84.775663 11 0.157653 0.157613 -21.070406 12 0.970579 0.970593 136.131192 13 0.957177 0.957167 261.490009 14 0.485394 0.485376 184.645131 15 0.800286 0.800280 157.963175 16 0.141896 0.141886 116.943402 17 0.421766 0.421761 73.124325 18 0.915721 0.915736 -51.128413 19 0.792224 0.792207 215.745886 20 0.959493 0.959492 -87.842624 21 0.655730 0.655741 170.828654 22 0.035722 0.035712 -87.912020 23 0.849122 0.849129 148.648245 24 0.933996 0.933993 -25.132363 25 0.678722 0.678735 96.643583 26 0.757732 0.757740 -45.397322 27 0.743135 0.743132 140.260117 28 0.392226 0.392227 -34.579202 29 0.655480 0.655478 81.590690 30 0.171181 0.171187 108.151506 31 0.706056 0.706046 110.035332 32 0.031838 0.031833 44.657357 33 0.276930 0.276923 82.333061 34 0.046169 0.046171 46.698326 35 0.097140 0.097132 -50.585530 36 0.823459 0.823458 177.454725 37 0.694841 0.694829 -60.228583 38 0.317101 0.317099 86.850607 39 0.950223 0.950222 -34.031423 40 0.034452 0.034446 93.072050 41 0.438735 0.438744 -67.786593 42 0.381547 0.381558 105.780783 43 0.765518 0.765517 -58.908015 44 0.795193 0.795200 172.043650 45 0.186875 0.186873 -49.199583 46 0.489752 0.489764 90.783171 47 0.445581 0.445586 126.433703 48 0.646315 0.646313 128.829770 49 0.709367 0.709365 173.739532 50 0.754685 0.754687 205.980716 51 0.276038 0.276025 91.712650 52 0.679703 0.679703 -39.409263 53 0.655111 0.655098 172.590824 54 0.162607 0.162612 -79.372126 55 0.119002 0.118998 80.140912 56 0.498359 0.498364 -68.781008 57 0.959733 0.959744 210.807743 58 0.340379 0.340386 -99.207154 59 0.585259 0.585268 184.803664 60 0.223814 0.223812 -98.745921 61 0.751269 0.751267 176.190105 62 0.255075 0.255095 -48.533400 63 0.505952 0.505957 72.505259 64 0.699073 0.699077 141.411745 65 0.890910 0.890903 180.326802 66 0.959302 0.959291 217.943822 67 0.547218 0.547216 153.709549 68 0.138635 0.138624 54.774117 69 0.149298 0.149294 -24.421668 70 0.257507 0.257508 66.986500 71 0.840729 0.840717 -28.059776 72 0.254282 0.254282 59.161719 73 0.814276 0.814285 -41.690515 74 0.243497 0.243525 95.514015 75 0.929249 0.929264 -49.414664 76 0.349984 0.349984 87.611487 77 0.196602 0.196595 -47.815884 78 0.251070 0.251084 55.699244 79 0.616075 0.616045 -16.566724 CG_RC_TEST: Normal end of execution. 07-Jan-2022 17:59:23