Wed Sep 12 21:11:22 2018 CG_RC_TEST: Python version Test the CG_RC library. CG_RC_TEST01 Use CG_RC on the 1, -2, 1 matrix. Number of iterations was 10 Estimated error is 4.24555e-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_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.0844833 Loo error is 0.000103532 I X(I) X_EXACT(I) B(I) 0 0.218433 0.218418 -189.806895 1 0.956336 0.956318 872.128190 2 0.829519 0.829509 -1119.106460 3 0.561709 0.561695 3876.412887 4 0.415305 0.415307 -1786.402527 5 0.066121 0.066119 2015.747522 6 0.257558 0.257578 -781.844235 7 0.109946 0.109957 814.321482 8 0.043808 0.043829 -724.168582 9 0.633965 0.633966 1218.765971 10 0.061728 0.061727 -258.076313 11 0.449540 0.449539 691.585738 12 0.401290 0.401306 4342.166678 13 0.754668 0.754673 6849.368339 14 0.797289 0.797287 3797.946784 15 0.001837 0.001838 1538.553820 16 0.897479 0.897504 1394.248255 17 0.350748 0.350752 -250.319833 18 0.094557 0.094545 695.814806 19 0.013633 0.013617 -1978.320010 20 0.859111 0.859097 5520.089512 21 0.840845 0.840847 -1639.477418 22 0.123105 0.123104 2330.705419 23 0.007512 0.007512 -1099.725995 24 0.260291 0.260303 1139.930636 25 0.912479 0.912484 -707.954877 26 0.113676 0.113664 2332.043370 27 0.351625 0.351629 -731.368595 28 0.822859 0.822887 233.126686 29 0.267123 0.267132 1214.192570 30 0.692073 0.692066 1559.637022 31 0.561668 0.561662 365.154554 32 0.861204 0.861216 2096.301126 33 0.453787 0.453794 1406.126329 34 0.912002 0.911977 -31.396924 35 0.597922 0.597917 332.969167 36 0.188966 0.188955 -859.138889 37 0.761490 0.761492 2666.820882 38 0.396998 0.396988 -1153.483280 39 0.185307 0.185314 1363.607931 40 0.574347 0.574366 -1374.035684 41 0.367033 0.367027 3180.675408 42 0.617219 0.617205 -1914.550979 43 0.361535 0.361529 4077.493466 44 0.212932 0.212930 -1367.228249 45 0.714447 0.714471 129.616429 46 0.117710 0.117707 1287.341026 47 0.299324 0.299329 2517.164950 48 0.824998 0.825003 5041.185781 49 0.824652 0.824660 6691.890344 50 0.061872 0.061862 2259.207907 51 0.710722 0.710781 -14.097502 52 0.088302 0.088283 70.291761 53 0.777983 0.777994 -507.656585 54 0.745313 0.745303 3821.675629 55 0.308675 0.308675 -1385.275242 56 0.899369 0.899373 2587.955206 57 0.763526 0.763537 -1179.188092 58 0.761738 0.761731 3778.629569 59 0.406982 0.406970 -1662.028930 60 0.938761 0.938749 3345.140152 61 0.562083 0.562088 -736.548482 62 0.017761 0.017820 -0.216840 63 0.501101 0.501103 2066.912826 64 0.041910 0.041909 1915.133702 65 0.368874 0.368851 966.846379 66 0.271722 0.271724 248.370312 67 0.858557 0.858573 70.444082 68 0.029140 0.029037 -0.127240 69 0.017413 0.017442 -0.119415 70 0.152387 0.152384 -528.151063 71 0.114316 0.114319 1320.425788 72 0.353912 0.353907 -695.056570 73 0.119288 0.119308 507.958252 74 0.206684 0.206653 -274.373783 75 0.212896 0.212924 165.747232 76 0.613013 0.612948 -63.376707 77 0.809485 0.809519 66.147916 78 0.587182 0.587090 -17.155671 CG_RC_TEST02: Normal end of execution. R8_UNIFORM_01_TEST Python version: 3.6.5 R8_UNIFORM_01 produces a sequence of random values. Using random seed 123456789 SEED R8_UNIFORM_01(SEED) 469049721 0.218418 2053676357 0.956318 1781357515 0.829509 1206231778 0.561695 891865166 0.415307 141988902 0.066119 553144097 0.257578 236130416 0.109957 94122056 0.043829 1361431000 0.633966 Verify that the sequence can be restarted. Set the seed back to its original value, and see that we generate the same sequence. SEED R8_UNIFORM_01(SEED) 469049721 0.218418 2053676357 0.956318 1781357515 0.829509 1206231778 0.561695 891865166 0.415307 141988902 0.066119 553144097 0.257578 236130416 0.109957 94122056 0.043829 1361431000 0.633966 R8_UNIFORM_01_TEST Normal end of execution. R8VEC_UNIFORM_01_TEST Python version: 3.6.5 R8VEC_UNIFORM_01 computes a random R8VEC. Initial seed is 123456789 Random R8VEC: 0: 0.218418 1: 0.956318 2: 0.829509 3: 0.561695 4: 0.415307 5: 0.0661187 6: 0.257578 7: 0.109957 8: 0.043829 9: 0.633966 R8VEC_UNIFORM_01_TEST: Normal end of execution. WATHEN_TEST Python version: 3.6.5 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 6.5081e-15 WATHEN_TEST: Normal end of execution. WATHEN_ORDER_TEST: Python version: 3.6.5 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. Wed Sep 12 21:11:22 2018