09 July 2016 3:07:21.282 PM LINPLUS_C8_PRB FORTRAN90 version: Test the LINPLUS_C8 library. C83_CR_FA_TEST C83_CR_FA factors a complex tridiagonal matrix; Matrix order N = 10 Solution: Col: 1 Row --- 1 1.00 10.0 2 2.00 20.0 3 3.00 30.0 4 4.00 40.0 5 5.00 50.0 6 6.00 60.0 7 7.00 70.0 8 8.00 80.0 9 9.00 90.0 10 10.0 100. C83_CR_SLS_TEST C83_CR_SLS solves one or more linear systems that were factored by C83_CR_FA. Matrix order N = 10 Solution: Col: 1 Row --- 1 1.00 10.0 2 2.00 20.0 3 3.00 30.0 4 4.00 40.0 5 5.00 50.0 6 6.00 60.0 7 7.00 70.0 8 8.00 80.0 9 9.00 90.0 10 10.0 100. C83_NP_FA_TEST C83_NP_FA factors a C83 matrix with no pivoting; Matrix order N = 10 The tridiagonal matrix Col: 1 1 2 2 3 3 Row --- 1 0.30636 0.26275E-01 0.24909 0.57874 2 0.44986 -0.12667 0.50080 -0.77993 -0.50839 0.56206 3 -0.84320 -0.34428 0.35047 0.16555E-01 4 0.58963 0.26009 Col: 4 4 5 5 6 6 Row --- 3 -0.50657 0.60050 4 0.43499 -0.26662 -0.10216 -0.44999 5 0.39114 0.32340 -0.20095 0.27071 -0.10464 0.32674 6 -0.13947 -0.15614 -0.97460E-01 0.90188 7 -0.23607 0.77459E-01 Col: 7 7 8 8 9 9 Row --- 6 0.41063 -0.81018 7 -0.77024 -0.31431 -0.60672E-01-0.24121 8 0.18599E-01-0.63321 -0.88918 0.26566 0.51993E-01-0.29254 9 0.89285 0.10314E-01-0.77987 -0.55116 10 -0.56047 0.76380 Col: 10 10 Row --- 9 -0.31111 0.80530 10 0.31360E-01-0.43356 The right hand side 1 3.47667 11.8337 2 18.8607 -14.8029 3 -19.7483 26.1456 4 -21.3440 24.1115 5 11.5822 -23.8848 6 -78.6323 -20.0172 7 118.490 -39.9184 8 -4.79867 5.50455 9 -15.5794 -121.738 10 -31.6080 -21.9519 The solution 1 1.00000 10.0000 2 2.00000 20.0000 3 3.00000 30.0000 4 4.00000 40.0000 5 5.00000 50.0000 6 6.00000 60.0000 7 7.00000 70.0000 8 8.00000 80.0000 9 9.00000 90.0000 10 10.0000 100.000 The second right hand side 1 12.2878 -1.06452 2 -10.7750 -21.4185 3 21.7173 24.5346 4 19.4075 25.6959 5 -21.1183 -16.0825 6 -35.1915 73.1114 7 -15.6646 -124.048 8 4.44532 5.79366 9 -122.413 -8.83572 10 -27.7762 26.6352 The second solution 1 10.0000 1.00000 2 20.0000 2.00000 3 30.0000 3.00000 4 40.0000 4.00000 5 50.0000 5.00000 6 60.0000 6.00000 7 70.0000 7.00000 8 80.0000 8.00000 9 90.0000 9.00000 10 100.000 10.0000 C83_NP_ML_TEST C83_NP_ML multiplies A*X after A has been factored by C83_NP_FA. Matrix order N = 10 The tridiagonal matrix Col: 1 1 2 2 3 3 Row --- 1 0.30636 0.26275E-01 0.24909 0.57874 2 0.44986 -0.12667 0.50080 -0.77993 -0.50839 0.56206 3 -0.84320 -0.34428 0.35047 0.16555E-01 4 0.58963 0.26009 Col: 4 4 5 5 6 6 Row --- 3 -0.50657 0.60050 4 0.43499 -0.26662 -0.10216 -0.44999 5 0.39114 0.32340 -0.20095 0.27071 -0.10464 0.32674 6 -0.13947 -0.15614 -0.97460E-01 0.90188 7 -0.23607 0.77459E-01 Col: 7 7 8 8 9 9 Row --- 6 0.41063 -0.81018 7 -0.77024 -0.31431 -0.60672E-01-0.24121 8 0.18599E-01-0.63321 -0.88918 0.26566 0.51993E-01-0.29254 9 0.89285 0.10314E-01-0.77987 -0.55116 10 -0.56047 0.76380 Col: 10 10 Row --- 9 -0.31111 0.80530 10 0.31360E-01-0.43356 The right hand side 1 3.47667 11.8337 2 18.8607 -14.8029 3 -19.7483 26.1456 4 -21.3440 24.1115 5 11.5822 -23.8848 6 -78.6323 -20.0172 7 118.490 -39.9184 8 -4.79867 5.50455 9 -15.5794 -121.738 10 -31.6080 -21.9519 The solution 1 1.00000 10.0000 2 2.00000 20.0000 3 3.00000 30.0000 4 4.00000 40.0000 5 5.00000 50.0000 6 6.00000 60.0000 7 7.00000 70.0000 8 8.00000 80.0000 9 9.00000 90.0000 10 10.0000 100.000 The second right hand side 1 12.2878 -1.06452 2 -10.7750 -21.4185 3 21.7173 24.5346 4 19.4075 25.6959 5 -21.1183 -16.0825 6 -35.1915 73.1114 7 -15.6646 -124.048 8 4.44532 5.79366 9 -122.413 -8.83572 10 -27.7762 26.6352 The second solution 1 10.0000 1.00000 2 20.0000 2.00000 3 30.0000 3.00000 4 40.0000 4.00000 5 50.0000 5.00000 6 60.0000 6.00000 7 70.0000 7.00000 8 80.0000 8.00000 9 90.0000 9.00000 10 100.000 10.0000 C83_NP_SL_TEST C83_NP_SL solves a linear system factored by C83_NP_FA. We will look at the TRANSPOSED linear system. Matrix order N = 10 The tridiagonal matrix Col: 1 1 2 2 3 3 Row --- 1 0.30636 0.26275E-01 0.24909 0.57874 2 0.44986 -0.12667 0.50080 -0.77993 -0.50839 0.56206 3 -0.84320 -0.34428 0.35047 0.16555E-01 4 0.58963 0.26009 Col: 4 4 5 5 6 6 Row --- 3 -0.50657 0.60050 4 0.43499 -0.26662 -0.10216 -0.44999 5 0.39114 0.32340 -0.20095 0.27071 -0.10464 0.32674 6 -0.13947 -0.15614 -0.97460E-01 0.90188 7 -0.23607 0.77459E-01 Col: 7 7 8 8 9 9 Row --- 6 0.41063 -0.81018 7 -0.77024 -0.31431 -0.60672E-01-0.24121 8 0.18599E-01-0.63321 -0.88918 0.26566 0.51993E-01-0.29254 9 0.89285 0.10314E-01-0.77987 -0.55116 10 -0.56047 0.76380 Col: 10 10 Row --- 9 -0.31111 0.80530 10 0.31360E-01-0.43356 The right hand side B1 1 -11.0331 9.22911 2 -0.701468E-01 -0.737258 3 -20.2921 -24.8494 4 28.3599 27.4444 5 -46.1438 3.92724 6 11.9986 14.8825 7 29.3570 -76.5999 8 42.8859 -70.0933 9 -34.7381 -26.6965 10 -30.1165 -44.7672 The solution to At * X1 = B1 1 1.00000 10.0000 2 2.00000 20.0000 3 3.00000 30.0000 4 4.00000 40.0000 5 5.00000 50.0000 6 6.00000 60.0000 7 7.00000 70.0000 8 8.00000 80.0000 9 9.00000 90.0000 10 10.0000 100.000 C8CI_SL_TEST C8CI_SL solves a complex circulant system. Matrix order N = 10 The circulant matrix: Col: 1 2 3 4 Row --- 1 0.450 -0.127 -0.843 -0.344 0.590 0.260 0.391 0.323 2 0.306 0.263E-01 0.450 -0.127 -0.843 -0.344 0.590 0.260 3 -0.560 0.764 0.306 0.263E-01 0.450 -0.127 -0.843 -0.344 4 0.893 0.103E-01-0.560 0.764 0.306 0.263E-01 0.450 -0.127 5 0.186E-01-0.633 0.893 0.103E-01-0.560 0.764 0.306 0.263E-01 6 -0.236 0.775E-01 0.186E-01-0.633 0.893 0.103E-01-0.560 0.764 7 -0.139 -0.156 -0.236 0.775E-01 0.186E-01-0.633 0.893 0.103E-01 8 0.391 0.323 -0.139 -0.156 -0.236 0.775E-01 0.186E-01-0.633 9 0.590 0.260 0.391 0.323 -0.139 -0.156 -0.236 0.775E-01 10 -0.843 -0.344 0.590 0.260 0.391 0.323 -0.139 -0.156 Col: 5 6 7 8 Row --- 1 -0.139 -0.156 -0.236 0.775E-01 0.186E-01-0.633 0.893 0.103E-01 2 0.391 0.323 -0.139 -0.156 -0.236 0.775E-01 0.186E-01-0.633 3 0.590 0.260 0.391 0.323 -0.139 -0.156 -0.236 0.775E-01 4 -0.843 -0.344 0.590 0.260 0.391 0.323 -0.139 -0.156 5 0.450 -0.127 -0.843 -0.344 0.590 0.260 0.391 0.323 6 0.306 0.263E-01 0.450 -0.127 -0.843 -0.344 0.590 0.260 7 -0.560 0.764 0.306 0.263E-01 0.450 -0.127 -0.843 -0.344 8 0.893 0.103E-01-0.560 0.764 0.306 0.263E-01 0.450 -0.127 9 0.186E-01-0.633 0.893 0.103E-01-0.560 0.764 0.306 0.263E-01 10 -0.236 0.775E-01 0.186E-01-0.633 0.893 0.103E-01-0.560 0.764 Col: 9 10 Row --- 1 -0.560 0.764 0.306 0.263E-01 2 0.893 0.103E-01-0.560 0.764 3 0.186E-01-0.633 0.893 0.103E-01 4 -0.236 0.775E-01 0.186E-01-0.633 5 -0.139 -0.156 -0.236 0.775E-01 6 0.391 0.323 -0.139 -0.156 7 0.590 0.260 0.391 0.323 8 -0.843 -0.344 0.590 0.260 9 0.450 -0.127 -0.843 -0.344 10 0.306 0.263E-01 0.450 -0.127 Solution: 1 1.00000 10.0000 2 2.00000 20.0000 3 3.00000 30.0000 4 4.00000 40.0000 5 5.00000 50.0000 6 6.00000 60.0000 7 7.00000 70.0000 8 8.00000 80.0000 9 9.00000 90.0000 10 10.0000 100.000 Solution to transposed system: 1 1.00000 10.0000 2 2.00000 20.0000 3 3.00000 30.0000 4 4.00000 40.0000 5 5.00000 50.0000 6 6.00000 60.0000 7 7.00000 70.0000 8 8.00000 80.0000 9 9.00000 90.0000 10 10.0000 100.000 C8TO_SL_TEST C8TO_SL solves a complex Toeplitz system. Matrix order N = 4 The Toeplitz matrix: Col: 1 2 3 4 Row --- 1 0.450 -0.127 -0.843 -0.344 0.590 0.260 0.391 0.323 2 -0.139 -0.156 0.450 -0.127 -0.843 -0.344 0.590 0.260 3 -0.236 0.775E-01-0.139 -0.156 0.450 -0.127 -0.843 -0.344 4 0.186E-01-0.633 -0.236 0.775E-01-0.139 -0.156 0.450 -0.127 Desired solution: 1 1.00000 -1.00000 2 2.00000 -2.00000 3 3.00000 -3.00000 4 4.00000 -4.00000 Right hand side: 1 3.35555 -0.838263 2 0.187223 -0.991121 3 -4.53014 0.546271 4 -0.525864 -2.38088 Computed solution: 1 1.00000 -1.00000 2 2.00000 -2.00000 3 3.00000 -3.00000 4 4.00000 -4.00000 Desired solution: 1 1.00000 -1.00000 2 2.00000 -2.00000 3 3.00000 -3.00000 4 4.00000 -4.00000 Right hand side: 1 -3.20229 -2.27654 2 -2.06232 0.549956 3 -1.73806 -1.12796 4 0.144314 -1.53617 Computed solution: 1 1.00000 -1.00000 2 2.00000 -2.00000 3 3.00000 -3.00000 4 4.00000 -4.00000 C8VEC_UNITY_TEST C8VEC_UNITY returns the N complex roots of unity. Roots of unity: 1 1.00000 0.00000 Roots of unity: 1 1.00000 0.00000 2 -1.00000 0.122465E-15 Roots of unity: 1 1.00000 0.00000 2 -0.500000 0.866025 3 -0.500000 -0.866025 Roots of unity: 1 1.00000 0.00000 2 0.612323E-16 1.00000 3 -1.00000 0.122465E-15 4 -0.183697E-15 -1.00000 Roots of unity: 1 1.00000 0.00000 2 0.309017 0.951057 3 -0.809017 0.587785 4 -0.809017 -0.587785 5 0.309017 -0.951057 LINPLUS_C8_PRB Normal end of execution. 09 July 2016 3:07:21.284 PM