Tue Oct 19 11:25:19 2021 cg_rc_test(): Python version Test cg_rc(). cg_rc_test01(): Use cg_rc() on the 1, -2, 1 matrix. Number of iterations was 10 Estimated error is 3.79195e-16 Loo error is 5.55112e-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_test01(): Normal end of execution. cg_rc_test02(): Use cg_rc() to solve a linear system involving the Wathen matrix. NX = 5 NY = 4 N = 79 Number of iterations was 19 Estimated error is 0.138994 Loo error is 4.22478e-05 I X(I) X_EXACT(I) B(I) 0 0.908716 0.908740 -255.491602 1 0.470417 0.470423 1150.656251 2 0.464610 0.464592 -1424.642943 3 0.907689 0.907679 4435.971711 4 0.067365 0.067356 -1804.155069 5 0.248088 0.248079 1336.691789 6 0.862560 0.862518 -573.815356 7 0.687204 0.687201 1204.476262 8 0.294817 0.294796 -1315.429133 9 0.369745 0.369739 3093.417896 10 0.010340 0.010343 -1093.971949 11 0.649622 0.649611 1283.308700 12 0.597681 0.597684 5583.417784 13 0.297511 0.297517 5264.008011 14 0.171333 0.171350 2476.673029 15 0.614860 0.614855 4472.210855 16 0.959135 0.959142 3687.106118 17 0.874050 0.874049 -398.511922 18 0.704528 0.704536 2262.322165 19 0.822112 0.822115 -1464.349459 20 0.686569 0.686560 4257.236419 21 0.410631 0.410635 -2252.075465 22 0.956259 0.956241 4710.267568 23 0.488381 0.488390 -2172.055030 24 0.490822 0.490817 4147.444689 25 0.332448 0.332472 -2590.296769 26 0.792053 0.792056 5248.660191 27 0.761087 0.761071 -979.797072 28 0.883049 0.883047 932.254599 29 0.503225 0.503202 799.191743 30 0.938223 0.938225 2572.748605 31 0.950245 0.950257 6273.456342 32 0.207562 0.207566 3795.847440 33 0.533158 0.533162 1428.362418 34 0.409388 0.409373 -887.000807 35 0.101244 0.101241 2658.531635 36 0.474956 0.474938 -1107.751206 37 0.967824 0.967832 923.876743 38 0.602872 0.602858 -1007.028560 39 0.486476 0.486473 3265.492874 40 0.426383 0.426381 -2986.126060 41 0.722941 0.722935 7965.562521 42 0.653556 0.653560 -2274.173164 43 0.163575 0.163582 1556.964677 44 0.702590 0.702600 -405.663846 45 0.926336 0.926344 2671.260705 46 0.346416 0.346419 2814.549627 47 0.191669 0.191657 1276.176745 48 0.681453 0.681471 5356.603319 49 0.975409 0.975409 5361.314342 50 0.462201 0.462178 570.676607 51 0.063843 0.063866 -1832.927605 52 0.837982 0.837968 6945.587826 53 0.214492 0.214485 -2849.582672 54 0.066022 0.066027 3020.457080 55 0.058002 0.058014 -1738.922518 56 0.176814 0.176791 1278.305818 57 0.908410 0.908399 -1762.027204 58 0.327878 0.327880 5890.347415 59 0.102260 0.102255 -2818.832055 60 0.189229 0.189222 2412.809855 61 0.534625 0.534620 -743.997266 62 0.564180 0.564193 3920.003030 63 0.689313 0.689327 7191.031395 64 0.635341 0.635349 4046.568740 65 0.298280 0.298284 2723.187357 66 0.815097 0.815103 4683.378502 67 0.472730 0.472734 2042.395661 68 0.920980 0.921007 -972.570390 69 0.171373 0.171373 3421.642801 70 0.144120 0.144144 -2168.045341 71 0.780318 0.780297 3548.260049 72 0.747636 0.747627 -704.378727 73 0.249917 0.249926 702.408992 74 0.939847 0.939872 -539.786032 75 0.834696 0.834692 2516.963795 76 0.167205 0.167218 -1153.976532 77 0.584108 0.584110 2218.370883 78 0.372072 0.372075 -591.230742 cg_rc_test02(): Normal end of execution. wathen_test(): Python version: 3.6.9 Assemble, factor and solve a Wathen system defined by wathen(). 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.1648e-14 wathen_test(): Normal end of execution. wathen_order_test(): Python version: 3.6.9 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 wathen_order_test(): Normal end of execution. cg_rc_test(): Normal end of execution. Tue Oct 19 11:25:19 2021