14 January 2013 2:26:41.007 PM CG_RC_PRB: FORTRAN90 version Test the CG_RC library. TEST01 Use CG_RC on the 1, -2, 1 matrix. Number of iterations was 9 Estimated error is 0.421596E-15 Loo error is 0.666134E-15 I X(I) X_EXACT(I) B(I) 1 0.780626E-16 0.00000 0.309017 2 0.309017 0.309017 -0.302487E-01 3 0.587785 0.587785 -0.575365E-01 4 0.809017 0.809017 -0.791922E-01 5 0.951057 0.951057 -0.930960E-01 6 1.00000 1.00000 -0.978870E-01 7 0.951057 0.951057 -0.930960E-01 8 0.809017 0.809017 -0.791922E-01 9 0.587785 0.587785 -0.575365E-01 10 0.309017 0.309017 -0.302487E-01 11 -0.479828E-16 0.122465E-15 0.388578E-15 12 -0.309017 -0.309017 0.302487E-01 13 -0.587785 -0.587785 0.575365E-01 14 -0.809017 -0.809017 0.791922E-01 15 -0.951057 -0.951057 0.930960E-01 16 -1.00000 -1.00000 0.978870E-01 17 -0.951057 -0.951057 0.930960E-01 18 -0.809017 -0.809017 0.791922E-01 19 -0.587785 -0.587785 0.575365E-01 20 -0.309017 -0.309017 0.302487E-01 21 -0.410262E-15 -0.244929E-15 -0.309017 TEST02 Use CG_RC to solve a linear system involving the Wathen matrix. NX = 5 NY = 4 N = 79 Number of iterations was 18 Estimated error is 0.730674E-02 Loo error is 0.422526E-04 I X(I) X_EXACT(I) B(I) 1 0.218432 0.218418 -19.3113 2 0.956316 0.956318 88.7317 3 0.829512 0.829509 -36.3017 4 0.561714 0.561695 90.0773 5 0.415273 0.415307 -43.7437 6 0.661560E-01 0.661187E-01 54.0011 7 0.257560 0.257578 -23.2089 8 0.109949 0.109957 32.2168 9 0.438096E-01 0.438290E-01 -33.6528 10 0.633973 0.633966 65.2136 11 0.617100E-01 0.617272E-01 -13.8091 12 0.449528 0.449539 70.3630 13 0.401292 0.401306 154.809 14 0.754653 0.754673 168.305 15 0.797282 0.797287 118.754 16 0.182286E-02 0.183837E-02 72.6417 17 0.897471 0.897504 74.6033 18 0.350743 0.350752 -37.7129 19 0.945624E-01 0.945448E-01 107.050 20 0.136232E-01 0.136169E-01 -109.428 21 0.859112 0.859097 207.905 22 0.840875 0.840847 -64.2180 23 0.123101 0.123104 99.2914 24 0.751016E-02 0.751236E-02 -67.2724 25 0.260308 0.260303 101.661 26 0.912477 0.912484 -42.2354 27 0.113677 0.113664 99.6324 28 0.351621 0.351629 -30.9739 29 0.822884 0.822887 78.3529 30 0.267130 0.267132 144.803 31 0.692067 0.692066 166.674 32 0.561670 0.561662 111.980 33 0.861210 0.861216 146.417 34 0.453772 0.453794 49.2886 35 0.911983 0.911977 -10.5401 36 0.597914 0.597917 114.666 37 0.188960 0.188955 -79.7521 38 0.761504 0.761492 176.378 39 0.396996 0.396988 -81.2853 40 0.185316 0.185314 108.356 41 0.574347 0.574366 -72.6881 42 0.367041 0.367027 156.156 43 0.617220 0.617205 -72.2069 44 0.361536 0.361529 124.845 45 0.212947 0.212930 -41.1752 46 0.714457 0.714471 46.6628 47 0.117699 0.117707 83.0738 48 0.299328 0.299329 131.474 49 0.824989 0.825003 189.243 50 0.824647 0.824660 191.901 51 0.618704E-01 0.618618E-01 62.9693 52 0.710738 0.710781 -1.73137 53 0.882712E-01 0.882833E-01 31.9862 54 0.777980 0.777994 -19.4569 55 0.745317 0.745303 131.412 56 0.308673 0.308675 -61.8462 57 0.899362 0.899373 151.142 58 0.763542 0.763537 -50.1830 59 0.761737 0.761731 163.776 60 0.406984 0.406970 -88.3359 61 0.938776 0.938749 202.851 62 0.562085 0.562088 -47.2525 63 0.178467E-01 0.178200E-01 -2.62116 64 0.501099 0.501103 66.5544 65 0.418957E-01 0.419093E-01 68.9155 66 0.368859 0.368851 96.4200 67 0.271712 0.271724 138.430 68 0.858573 0.858573 114.962 69 0.290067E-01 0.290366E-01 -1.53808 70 0.174254E-01 0.174423E-01 -1.44348 71 0.152396 0.152384 -11.5233 72 0.114317 0.114319 32.6938 73 0.353917 0.353907 -24.8832 74 0.119309 0.119308 32.1821 75 0.206644 0.206653 -28.0064 76 0.212946 0.212924 38.9580 77 0.612979 0.612948 -40.6021 78 0.809533 0.809519 107.951 79 0.587116 0.587090 -27.9974 CG_RC_PRB: Normal end of execution. 14 January 2013 2:26:41.009 PM