Wed Oct 8 07:27:28 2025 cg_rc_test(): python version: 3.10.12 numpy version: 1.26.4 Test cg_rc(). cg_rc_test01(): Use cg_rc() on the 1, -2, 1 matrix. Number of iterations was 10 Estimated error is 4.6391e-16 Loo error is 4.44089e-16 I X(I) X_EXACT(I) B(I) 0 -0.000000 0.000000 0.309017 1 0.309017 0.309017 -0.030249 2 0.587785 0.587785 -0.057537 3 0.809017 0.809017 -0.079192 4 0.951057 0.951057 -0.093096 5 1.000000 1.000000 -0.097887 6 0.951057 0.951057 -0.093096 7 0.809017 0.809017 -0.079192 8 0.587785 0.587785 -0.057537 9 0.309017 0.309017 -0.030249 10 -0.000000 0.000000 0.000000 11 -0.309017 -0.309017 0.030249 12 -0.587785 -0.587785 0.057537 13 -0.809017 -0.809017 0.079192 14 -0.951057 -0.951057 0.093096 15 -1.000000 -1.000000 0.097887 16 -0.951057 -0.951057 0.093096 17 -0.809017 -0.809017 0.079192 18 -0.587785 -0.587785 0.057537 19 -0.309017 -0.309017 0.030249 20 -0.000000 -0.000000 -0.309017 cg_rc_test02(): cg_rc() solves a linear system involving the Wathen matrix. NX = 5 NY = 4 N = 79 Number of iterations was 19 Estimated error is 0.136727 Loo error is 0.000133136 I X(I) X_EXACT(I) B(I) 0 0.554623 0.554638 -223.616122 1 0.470487 0.470474 1017.608031 2 0.168500 0.168478 -942.097882 3 0.355896 0.355905 1947.682869 4 0.057881 0.057909 -1496.835360 5 0.224651 0.224654 1816.267736 6 0.091790 0.091779 -989.739657 7 0.454300 0.454307 1218.061966 8 0.136371 0.136363 -553.127600 9 0.558703 0.558706 596.616450 10 0.737018 0.737013 -133.611302 11 0.033544 0.033559 833.184608 12 0.243513 0.243528 3116.690118 13 0.049463 0.049458 3696.343859 14 0.751062 0.751068 3952.722108 15 0.266173 0.266156 1665.304736 16 0.089447 0.089400 544.946244 17 0.031021 0.031062 -659.156787 18 0.563010 0.563011 2189.304347 19 0.202878 0.202893 -1238.539995 20 0.997305 0.997295 3183.316512 21 0.584509 0.584519 -1782.532005 22 0.751720 0.751732 5278.185222 23 0.566177 0.566192 -1909.253370 24 0.417464 0.417476 3353.226442 25 0.324581 0.324576 -1869.813636 26 0.827646 0.827651 3343.699905 27 0.784941 0.784929 -652.744722 28 0.311312 0.311332 1113.040921 29 0.259252 0.259244 1360.992813 30 0.253654 0.253659 2499.901644 31 0.783767 0.783743 5308.850990 32 0.492053 0.492046 4472.740494 33 0.771515 0.771508 2652.369926 34 0.578758 0.578754 -279.925791 35 0.595546 0.595561 2365.961407 36 0.510600 0.510600 -698.320048 37 0.254948 0.254937 1232.088636 38 0.317356 0.317346 -1561.885415 39 0.609851 0.609869 4378.911627 40 0.621420 0.621427 -2631.106289 41 0.540732 0.540756 5620.613744 42 0.922461 0.922452 -1929.621248 43 0.530543 0.530549 2765.906427 44 0.763127 0.763116 -709.787095 45 0.358164 0.358144 912.198062 46 0.099203 0.099194 1409.133723 47 0.647264 0.647247 3088.714559 48 0.416407 0.416395 4925.680069 49 0.983064 0.983046 4310.185183 50 0.931071 0.931032 389.625457 51 0.603227 0.603227 -556.855027 52 0.073831 0.073846 2135.259668 53 0.440472 0.440462 -1349.463479 54 0.387623 0.387627 2372.784013 55 0.141099 0.141075 -1259.023327 56 0.496050 0.496056 1807.524432 57 0.964945 0.964953 -1775.879116 58 0.977010 0.977022 6599.368977 59 0.600117 0.600107 -1992.923054 60 0.422777 0.422775 1874.215019 61 0.582595 0.582602 -413.203764 62 0.620501 0.620497 2146.684816 63 0.536250 0.536237 3601.422579 64 0.172106 0.172084 1264.667590 65 0.803990 0.803971 2680.227130 66 0.475839 0.475842 4065.048707 67 0.206282 0.206277 1263.636603 68 0.275088 0.275100 -514.312664 69 0.561298 0.561306 2221.467050 70 0.523969 0.523969 -768.112946 71 0.775693 0.775711 1655.873727 72 0.576862 0.576874 -348.444046 73 0.552796 0.552663 35.120173 74 0.176665 0.176702 -780.503094 75 0.773764 0.773770 2495.220015 76 0.766442 0.766433 -890.678233 77 0.488608 0.488607 1544.326641 78 0.025177 0.025165 -474.586347 wathen_test(): Assemble, factor and solve a Wathen system. Elements in X direction NX = 4 Elements in Y direction NY = 4 Number of elements = 16 Number of nodes N = 65 Norm of solution error is 1.74838e-15 wathen_order_test(): wathen_order() returns N, the order of a Wathen finite element matrix given NX and NY, the number of rows and columns of nodes in the underlying grid. NX / NY: 1 2 3 4 5 6 1 8 13 18 23 28 33 2 13 21 29 37 45 53 3 18 29 40 51 62 73 4 23 37 51 65 79 93 5 28 45 62 79 96 113 6 33 53 73 93 113 133 7 38 61 84 107 130 153 8 43 69 95 121 147 173 9 48 77 106 135 164 193 10 53 85 117 149 181 213 cg_rc_test(): Normal end of execution. Wed Oct 8 07:27:28 2025