Tue May 20 21:19:20 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.115878 Loo error is 4.16372e-05 I X(I) X_EXACT(I) B(I) 0 0.707047 0.707036 -736.028374 1 0.734649 0.734643 3145.292352 2 0.434955 0.434958 -852.125755 3 0.378464 0.378433 312.563935 4 0.605704 0.605716 -847.727816 5 0.692459 0.692458 3035.233745 6 0.369084 0.369085 -1553.942408 7 0.516637 0.516638 2457.103194 8 0.032550 0.032543 -1367.300936 9 0.951370 0.951349 1772.699946 10 0.326560 0.326565 -429.016360 11 0.408164 0.408158 2788.790856 12 0.809011 0.809010 3626.192858 13 0.951500 0.951493 3667.271443 14 0.443130 0.443121 5242.760254 15 0.237192 0.237193 3962.624333 16 0.035276 0.035296 1377.409975 17 0.783714 0.783709 -830.278429 18 0.556714 0.556718 3314.449031 19 0.790330 0.790319 -1277.119832 20 0.387470 0.387488 2391.785948 21 0.698177 0.698172 -1831.188252 22 0.541046 0.541053 5256.443179 23 0.043112 0.043099 -3036.483301 24 0.900692 0.900705 4774.320424 25 0.743758 0.743743 -1409.249283 26 0.766863 0.766851 1753.096775 27 0.030790 0.030801 -550.326656 28 0.294838 0.294812 299.566046 29 0.502216 0.502197 2486.122468 30 0.677427 0.677417 4915.318923 31 0.465007 0.465005 3978.295412 32 0.268184 0.268177 1708.109185 33 0.189682 0.189640 59.483966 34 0.984588 0.984576 63.892974 35 0.630236 0.630253 1728.029497 36 0.747727 0.747728 -456.477184 37 0.522054 0.522073 2356.592319 38 0.509405 0.509414 -1570.268348 39 0.526549 0.526543 3802.672971 40 0.607014 0.607031 -1612.539408 41 0.811257 0.811258 3756.813287 42 0.969436 0.969433 -630.946066 43 0.697606 0.697622 1566.679032 44 0.099786 0.099786 -322.157980 45 0.367084 0.367075 960.555764 46 0.433243 0.433234 1088.789264 47 0.271757 0.271762 1388.773744 48 0.737271 0.737272 3507.199018 49 0.098077 0.098071 2093.179817 50 0.273754 0.273729 1130.998846 51 0.951721 0.951734 -353.745378 52 0.041504 0.041506 1930.173453 53 0.539827 0.539832 -793.236301 54 0.607103 0.607114 838.406224 55 0.743116 0.743110 -795.235114 56 0.857044 0.857051 3186.563826 57 0.986256 0.986251 -752.441112 58 0.242875 0.242879 1789.542734 59 0.469130 0.469128 -1270.135155 60 0.332532 0.332538 2613.516619 61 0.776145 0.776157 -400.110637 62 0.535954 0.535956 1852.264119 63 0.598701 0.598704 2564.373796 64 0.408005 0.407992 1947.988757 65 0.126419 0.126411 1810.365884 66 0.345480 0.345486 2298.869620 67 0.464048 0.464047 1827.685515 68 0.115626 0.115633 -498.603595 69 0.593429 0.593424 1965.864254 70 0.299366 0.299356 -585.845372 71 0.559587 0.559605 637.662824 72 0.202121 0.202128 -562.643125 73 0.769514 0.769520 1394.567209 74 0.750013 0.750006 -356.235058 75 0.724880 0.724872 1004.733588 76 0.501079 0.501077 -401.924128 77 0.575457 0.575462 1899.429386 78 0.792439 0.792451 -251.917795 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 3.77751e-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. Tue May 20 21:19:20 2025