12 July 2016 8:52:13.993 PM R8BB_PRB FORTRAN90 version: Test the R8BB library. R8BB_ADD_TEST R8BB_ADD adds a value to elements of an R8BB matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 0 Matrix before additions: Col: 1 2 3 4 5 Row --- 1 11.0000 0.00000 0.00000 14.0000 15.0000 2 21.0000 22.0000 0.00000 24.0000 25.0000 3 0.00000 32.0000 33.0000 34.0000 35.0000 4 41.0000 42.0000 43.0000 44.0000 45.0000 5 51.0000 52.0000 53.0000 54.0000 55.0000 The R8BB matrix after additions: Col: 1 2 3 4 5 Row --- 1 111.000 0.00000 0.00000 214.000 215.000 2 21.0000 122.000 0.00000 224.000 225.000 3 0.00000 32.0000 133.000 234.000 235.000 4 41.0000 42.0000 43.0000 44.0000 445.000 5 51.0000 52.0000 53.0000 454.000 55.0000 R8BB_DIF2_TEST R8BB_DIF2 sets up an R8BB second difference matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 6 0.00000 0.00000 0.00000 0.00000 -1.00000 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 Row --- 1 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 5 -1.00000 0.00000 0.00000 6 2.00000 -1.00000 0.00000 7 -1.00000 2.00000 -1.00000 8 0.00000 -1.00000 2.00000 R8BB_FA_TEST R8BB_FA factors an R8BB matrix. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The border-banded matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.829509 0.00000 0.00000 0.00000 2 0.956318 0.561695 0.661187E-01 0.00000 0.00000 3 0.00000 0.415307 0.257578 0.438290E-01 0.00000 4 0.00000 0.00000 0.109957 0.633966 0.449539 5 0.00000 0.00000 0.00000 0.617272E-01 0.401306 6 0.00000 0.00000 0.00000 0.00000 0.754673 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.396988 0.574366 0.617205 0.212930 0.117707 10 0.185314 0.367027 0.361529 0.714471 0.299329 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.123104 0.692066 2 0.00000 0.00000 0.00000 0.751236E-02 0.561662 3 0.00000 0.00000 0.00000 0.260303 0.861216 4 0.00000 0.00000 0.00000 0.912484 0.453794 5 0.797287 0.00000 0.00000 0.113664 0.911977 6 0.183837E-02 0.350752 0.00000 0.351629 0.597917 7 0.897504 0.945448E-01 0.859097 0.822887 0.188955 8 0.00000 0.136169E-01 0.840847 0.267132 0.761492 9 0.825003 0.618618E-01 0.882833E-01 0.745303 0.899373 10 0.824660 0.710781 0.777994 0.308675 0.763537 The right hand side vector: 1 9.90604 2 7.96230 3 12.7335 4 17.8637 5 17.1799 6 15.3835 7 22.2151 8 16.8412 9 26.6284 10 32.9193 The FACTORED border-banded matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.561695 0.661187E-01 0.00000 0.00000 2 -0.228395 0.701221 -0.151012E-01 0.00000 0.00000 3 0.00000 -0.592263 0.266522 0.438290E-01 0.00000 4 0.00000 0.00000 -0.412562 0.615884 0.449539 5 0.00000 0.00000 0.00000 -0.100225 0.754673 6 0.00000 0.00000 0.00000 0.00000 -0.472060 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.396988 0.574366 0.617205 0.212930 0.117707 10 0.185314 0.367027 0.361529 0.714471 0.299329 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.125967 0.710759E-01 2 0.00000 0.00000 0.00000 -0.181574 -0.853023 3 0.00000 0.00000 0.00000 -0.393040 -2.27612 4 0.00000 0.00000 0.00000 -1.90868 1.80997 5 0.183837E-02 0.350752 0.00000 0.758046 -3.00525 6 0.897504 0.945448E-01 0.859097 -0.376345 0.228679 7 -0.887371 -0.249472 -0.762338 -2.63153 4.76017 8 0.00000 0.545828E-01 0.799237 -0.275079 -0.982712 9 0.825003 0.618618E-01 0.882833E-01 -3.40831 2.84183 10 0.824660 0.710781 0.777994 -0.107477 -0.844588 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8BB_GET_TEST R8BB_GET gets a value of an element of an R8BB matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 0 The matrix to be queried: Col: 1 2 3 4 5 Row --- 1 11.0000 0.00000 0.00000 14.0000 15.0000 2 21.0000 22.0000 0.00000 24.0000 25.0000 3 0.00000 32.0000 33.0000 34.0000 35.0000 4 41.0000 42.0000 43.0000 44.0000 45.0000 5 51.0000 52.0000 53.0000 54.0000 55.0000 A( 2, 5) = 25.0000 A( 5, 3) = 53.0000 A( 3, 1) = 0.00000 A( 2, 1) = 21.0000 A( 1, 4) = 14.0000 A( 1, 3) = 0.00000 A( 3, 4) = 34.0000 A( 4, 1) = 41.0000 A( 5, 2) = 52.0000 A( 1, 1) = 11.0000 R8BB_INDICATOR_TEST R8BB_INDICATOR sets up an R8BB indicator matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 0.00000 32.0000 33.0000 34.0000 0.00000 4 0.00000 0.00000 43.0000 44.0000 45.0000 5 0.00000 0.00000 0.00000 54.0000 55.0000 6 0.00000 0.00000 0.00000 0.00000 65.0000 7 71.0000 72.0000 73.0000 74.0000 75.0000 8 81.0000 82.0000 83.0000 84.0000 85.0000 Col: 6 7 8 Row --- 1 0.00000 17.0000 18.0000 2 0.00000 27.0000 28.0000 3 0.00000 37.0000 38.0000 4 0.00000 47.0000 48.0000 5 56.0000 57.0000 58.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 R8BB_MTV_TEST R8BB_MTV computes b=A'*x, where A is an R8BB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix A: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 0.00000 32.0000 33.0000 34.0000 0.00000 4 0.00000 0.00000 43.0000 44.0000 45.0000 5 0.00000 0.00000 0.00000 54.0000 55.0000 6 0.00000 0.00000 0.00000 0.00000 65.0000 7 71.0000 72.0000 73.0000 74.0000 75.0000 8 81.0000 82.0000 83.0000 84.0000 85.0000 Col: 6 7 8 Row --- 1 0.00000 17.0000 18.0000 2 0.00000 27.0000 28.0000 3 0.00000 37.0000 38.0000 4 0.00000 47.0000 48.0000 5 56.0000 57.0000 58.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 The vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 The product b=A'*x: 1 1198.00 2 1312.00 3 1492.00 4 1738.00 5 2050.00 6 1896.00 7 2292.00 8 2328.00 R8BB_MV_TEST R8BB_MV computes b=A*x, where A is an R8BB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix A: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 0.00000 32.0000 33.0000 34.0000 0.00000 4 0.00000 0.00000 43.0000 44.0000 45.0000 5 0.00000 0.00000 0.00000 54.0000 55.0000 6 0.00000 0.00000 0.00000 0.00000 65.0000 7 71.0000 72.0000 73.0000 74.0000 75.0000 8 81.0000 82.0000 83.0000 84.0000 85.0000 Col: 6 7 8 Row --- 1 0.00000 17.0000 18.0000 2 0.00000 27.0000 28.0000 3 0.00000 37.0000 38.0000 4 0.00000 47.0000 48.0000 5 56.0000 57.0000 58.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 The vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 The product b=A*x: 1 298.000 2 547.000 3 862.000 4 1243.00 5 1690.00 6 1734.00 7 2724.00 8 3084.00 R8BB_PRINT_TEST R8BB_PRINT prints an R8BB indicator matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 0.00000 32.0000 33.0000 34.0000 0.00000 4 0.00000 0.00000 43.0000 44.0000 45.0000 5 0.00000 0.00000 0.00000 54.0000 55.0000 6 0.00000 0.00000 0.00000 0.00000 65.0000 7 71.0000 72.0000 73.0000 74.0000 75.0000 8 81.0000 82.0000 83.0000 84.0000 85.0000 Col: 6 7 8 Row --- 1 0.00000 17.0000 18.0000 2 0.00000 27.0000 28.0000 3 0.00000 37.0000 38.0000 4 0.00000 47.0000 48.0000 5 56.0000 57.0000 58.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 R8BB_PRINT_SOME_TEST R8BB_PRINT_SOME prints some of an R8BB indicator matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The Lower Right Block: Col: 7 8 Row --- 7 77.0000 78.0000 8 87.0000 88.0000 R8BB_RANDOM_TEST R8BB_RANDOM randomizes an R8BB matrix; R8BB_SL solves. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The border-banded matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.829509 0.00000 0.00000 0.00000 2 0.956318 0.561695 0.661187E-01 0.00000 0.00000 3 0.00000 0.415307 0.257578 0.438290E-01 0.00000 4 0.00000 0.00000 0.109957 0.633966 0.449539 5 0.00000 0.00000 0.00000 0.617272E-01 0.401306 6 0.00000 0.00000 0.00000 0.00000 0.754673 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.396988 0.574366 0.617205 0.212930 0.117707 10 0.185314 0.367027 0.361529 0.714471 0.299329 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.123104 0.692066 2 0.00000 0.00000 0.00000 0.751236E-02 0.561662 3 0.00000 0.00000 0.00000 0.260303 0.861216 4 0.00000 0.00000 0.00000 0.912484 0.453794 5 0.797287 0.00000 0.00000 0.113664 0.911977 6 0.183837E-02 0.350752 0.00000 0.351629 0.597917 7 0.897504 0.945448E-01 0.859097 0.822887 0.188955 8 0.00000 0.136169E-01 0.840847 0.267132 0.761492 9 0.825003 0.618618E-01 0.882833E-01 0.745303 0.899373 10 0.824660 0.710781 0.777994 0.308675 0.763537 R8BB_SET_TEST R8BB_SET sets elements of an R8BB matrix. Matrix order N = 5 Matrix suborder N1 = 4 Matrix suborder N2 = 1 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The R8BB matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 15.0000 2 21.0000 22.0000 23.0000 0.00000 25.0000 3 31.0000 32.0000 33.0000 34.0000 35.0000 4 0.00000 42.0000 43.0000 44.0000 45.0000 5 51.0000 52.0000 53.0000 54.0000 55.0000 R8BB_SL_TEST R8BB_SL solves a linear system factored by R8BB_FA. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The border-banded matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.829509 0.00000 0.00000 0.00000 2 0.956318 0.561695 0.661187E-01 0.00000 0.00000 3 0.00000 0.415307 0.257578 0.438290E-01 0.00000 4 0.00000 0.00000 0.109957 0.633966 0.449539 5 0.00000 0.00000 0.00000 0.617272E-01 0.401306 6 0.00000 0.00000 0.00000 0.00000 0.754673 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.396988 0.574366 0.617205 0.212930 0.117707 10 0.185314 0.367027 0.361529 0.714471 0.299329 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.123104 0.692066 2 0.00000 0.00000 0.00000 0.751236E-02 0.561662 3 0.00000 0.00000 0.00000 0.260303 0.861216 4 0.00000 0.00000 0.00000 0.912484 0.453794 5 0.797287 0.00000 0.00000 0.113664 0.911977 6 0.183837E-02 0.350752 0.00000 0.351629 0.597917 7 0.897504 0.945448E-01 0.859097 0.822887 0.188955 8 0.00000 0.136169E-01 0.840847 0.267132 0.761492 9 0.825003 0.618618E-01 0.882833E-01 0.745303 0.899373 10 0.824660 0.710781 0.777994 0.308675 0.763537 The right hand side vector: 1 9.90604 2 7.96230 3 12.7335 4 17.8637 5 17.1799 6 15.3835 7 22.2151 8 16.8412 9 26.6284 10 32.9193 The FACTORED border-banded matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.561695 0.661187E-01 0.00000 0.00000 2 -0.228395 0.701221 -0.151012E-01 0.00000 0.00000 3 0.00000 -0.592263 0.266522 0.438290E-01 0.00000 4 0.00000 0.00000 -0.412562 0.615884 0.449539 5 0.00000 0.00000 0.00000 -0.100225 0.754673 6 0.00000 0.00000 0.00000 0.00000 -0.472060 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 9 0.396988 0.574366 0.617205 0.212930 0.117707 10 0.185314 0.367027 0.361529 0.714471 0.299329 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.125967 0.710759E-01 2 0.00000 0.00000 0.00000 -0.181574 -0.853023 3 0.00000 0.00000 0.00000 -0.393040 -2.27612 4 0.00000 0.00000 0.00000 -1.90868 1.80997 5 0.183837E-02 0.350752 0.00000 0.758046 -3.00525 6 0.897504 0.945448E-01 0.859097 -0.376345 0.228679 7 -0.887371 -0.249472 -0.762338 -2.63153 4.76017 8 0.00000 0.545828E-01 0.799237 -0.275079 -0.982712 9 0.825003 0.618618E-01 0.882833E-01 -3.40831 2.84183 10 0.824660 0.710781 0.777994 -0.107477 -0.844588 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8BB_TO_R8GE_TEST R8BB_TO_R8GE converts an R8BB matrix to R8GE format. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 0.00000 32.0000 33.0000 34.0000 0.00000 4 0.00000 0.00000 43.0000 44.0000 45.0000 5 0.00000 0.00000 0.00000 54.0000 55.0000 6 0.00000 0.00000 0.00000 0.00000 65.0000 7 71.0000 72.0000 73.0000 74.0000 75.0000 8 81.0000 82.0000 83.0000 84.0000 85.0000 Col: 6 7 8 Row --- 1 0.00000 17.0000 18.0000 2 0.00000 27.0000 28.0000 3 0.00000 37.0000 38.0000 4 0.00000 47.0000 48.0000 5 56.0000 57.0000 58.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 0.00000 0.00000 0.00000 2 21.0000 22.0000 23.0000 0.00000 0.00000 3 0.00000 32.0000 33.0000 34.0000 0.00000 4 0.00000 0.00000 43.0000 44.0000 45.0000 5 0.00000 0.00000 0.00000 54.0000 55.0000 6 0.00000 0.00000 0.00000 0.00000 65.0000 7 71.0000 72.0000 73.0000 74.0000 75.0000 8 81.0000 82.0000 83.0000 84.0000 85.0000 Col: 6 7 8 Row --- 1 0.00000 17.0000 18.0000 2 0.00000 27.0000 28.0000 3 0.00000 37.0000 38.0000 4 0.00000 47.0000 48.0000 5 56.0000 57.0000 58.0000 6 66.0000 67.0000 68.0000 7 76.0000 77.0000 78.0000 8 86.0000 87.0000 88.0000 R8BB_ZEROS_TEST R8BB_ZEROS zeros an R8BB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB zero matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 0.00000 6 0.00000 0.00000 0.00000 0.00000 0.00000 7 0.00000 0.00000 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 Row --- 1 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 6 0.00000 0.00000 0.00000 7 0.00000 0.00000 0.00000 8 0.00000 0.00000 0.00000 R8BB_PRB Normal end of execution. 12 July 2016 8:52:13.997 PM