Sun Aug 21 20:55:02 2022 r8cbb_test(): Python version: 3.6.9 Test r8cbb(). R8CBB_ADD_TEST R8CBB_ADD adds a value to elements of an R8CBB 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: 0 1 2 3 4 Row 0 : 0 0 0 0 15 1 : 21 22 0 0 25 2 : 0 32 33 0 35 3 : 0 0 43 14 41 4 : 51 52 53 54 55 The matrix after additions: Col: 0 1 2 3 4 Row 0 : 100 0 0 0 215 1 : 21 122 0 0 225 2 : 0 32 133 0 235 3 : 0 0 43 214 41 4 : 51 52 53 454 55 R8CBB_DIF2_TEST R8CBB_DIF2 sets up an R8CBB second difference matrix. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB second difference matrix: Col: 0 1 2 3 4 Row 0 : 2 -1 0 0 0 1 : -1 2 -1 0 0 2 : 0 -1 2 -1 0 3 : 0 0 -1 2 -1 4 : 0 0 0 -1 2 5 : 0 0 0 0 -1 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : -1 0 0 0 0 5 : 2 -1 0 0 0 6 : -1 2 -1 0 0 7 : 0 -1 2 0 0 8 : 0 0 0 0 0 9 : 0 0 0 -1 2 R8CBB_FA_TEST R8CBB_FA factors an R8CBB matrix, with no pivoting. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB matrix: Col: 0 1 2 3 4 Row 0 : 0.165512 0.407578 0 0 0 1 : 0.874734 0.610803 0.146375 0 0 2 : 0 0.028462 0.111448 0.641631 0 3 : 0 0 0.04227 0.554477 0.353925 4 : 0 0 0 0.076505 0.473532 5 : 0 0 0 0 0.206923 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0.0434475 0.428982 0.200522 0.35036 0.394299 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.814567 1 : 0 0 0 0 0.0145098 2 : 0 0 0 0 0.273021 3 : 0 0 0 0 0.485526 4 : 0.344834 0 0 0 0.210331 5 : 0.282506 0.338827 0 0 0.556461 6 : 0.542698 0.162269 0.73677 0 0.167214 7 : 0 0.589235 0.468655 0.903827 0.781859 8 : 0 0 0 0.266018 0.164467 9 : 0.845244 0.285127 0.461612 0.360216 0.780593 The factored matrix: Col: 0 1 2 3 4 Row 0 : 0.165512 0.407578 0 0 0 1 : -5.28501 -1.54325 0.146375 0 0 2 : 0 0.0184429 0.114148 0.641631 0 3 : 0 0 -0.37031 0.316875 0.353925 4 : 0 0 0 -0.241436 0.388082 5 : 0 0 0 0 -0.533195 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0.0434475 0.428982 0.200522 0.35036 0.394299 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.922257 1 : 0 0 0 0 -2.37307 2 : 0 0 0 0 4.29198 3 : 0 0 0 0 -1.06574 4 : 0.344834 0 0 0 -0.214792 5 : 0.0986425 0.338827 0 0 -0.0785489 6 : -5.50166 -1.70184 0.73677 0 -1.44565 7 : 0 0.346234 0.72375 3.57394 0.149297 8 : 0 0 0 -3.66888 0.164467 9 : 0.845244 0.285127 0.461612 -0.36197 0.275457 The right hand side vector: 0: 17.2608 1: 5.07473 2: 8.08916 3: 10.3297 4: 8.16347 5: 12.8022 6: 15.079 7: 24.3384 8: 24.5567 9: 26.684 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R8CBB_GET_TEST R8CBB_GET gets a value of an element of an R8CBB matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 0 The R8CBB matrix to be queried: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 15 1 : 21 22 0 0 25 2 : 0 32 33 0 35 3 : 0 0 43 14 41 4 : 51 52 53 54 55 A( [4] , [2] ) = [53.] A( [4] , [4] ) = [55.] A( [4] , [1] ) = [52.] A( [0] , [0] ) = [0.] A( [0] , [0] ) = [0.] A( [4] , [4] ) = [55.] A( [1] , [2] ) = 0.0 A( [0] , [3] ) = 0.0 A( [0] , [1] ) = 0.0 A( [4] , [4] ) = [55.] R8CBB_INDICATOR_TEST R8CBB_INDICATOR sets up an R8CBB indicator matrix. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The compact border-banded matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 201 202 203 0 0 2 : 0 302 303 304 0 3 : 0 0 403 404 405 4 : 0 0 0 504 505 5 : 0 0 0 0 605 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 1001 1002 1003 1004 1005 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 110 1 : 0 0 0 0 210 2 : 0 0 0 0 310 3 : 0 0 0 0 410 4 : 506 0 0 0 510 5 : 606 607 0 0 610 6 : 706 707 708 0 710 7 : 0 807 808 109 810 8 : 0 0 908 209 901 9 : 1006 1007 1008 1009 1010 R8CBB_MTV_TEST R8CBB_MTV computes b=A'*x, where A is an R8CBB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB matrix A: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 21 22 23 0 0 2 : 0 32 33 34 0 3 : 0 0 43 44 45 4 : 0 0 0 54 55 5 : 0 0 0 0 65 6 : 0 0 0 0 0 7 : 81 82 83 84 85 Col: 5 6 7 Row 0 : 0 0 18 1 : 0 0 28 2 : 0 0 38 3 : 0 0 48 4 : 56 0 58 5 : 66 17 68 6 : 76 27 71 7 : 86 87 88 The vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 The product b=A'*x: 0: 1187 1: 1300 2: 1492 3: 1738 4: 2050 5: 1896 6: 2292 7: 2328 R8CBB_MV_TEST R8CBB_MV computes b=A*x, where A is an R8CBB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB matrix A: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 21 22 23 0 0 2 : 0 32 33 34 0 3 : 0 0 43 44 45 4 : 0 0 0 54 55 5 : 0 0 0 0 65 6 : 0 0 0 0 0 7 : 81 82 83 84 85 Col: 5 6 7 Row 0 : 0 0 18 1 : 0 0 28 2 : 0 0 38 3 : 0 0 48 4 : 56 0 58 5 : 66 17 68 6 : 76 27 71 7 : 86 87 88 The vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 The product b=A*x: 0: 263 1: 547 2: 862 3: 1243 4: 1690 5: 1734 6: 2724 7: 3084 R8CBB_PRINT_TEST R8CBB_PRINT prints an R8CBB matrix Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB matrix: Col: 0 1 2 3 4 Row 0 : 0.592524 0.85277 0 0 0 1 : 0.600942 0.626002 0.933746 0 0 2 : 0 0.384237 0.82163 0.89027 0 3 : 0 0 0.535288 0.516338 0.974479 4 : 0 0 0 0.686014 0.398557 5 : 0 0 0 0 0.583092 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0.0910836 0.109657 0.836676 0.0397571 0.87174 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.911487 1 : 0 0 0 0 0.408304 2 : 0 0 0 0 0.943705 3 : 0 0 0 0 0.327025 4 : 0.521204 0 0 0 0.537426 5 : 0.0308943 0.332065 0 0 0.863374 6 : 0.697468 0.217756 0.0705907 0 0.961828 7 : 0 0.949902 0.749121 0.558154 0.333542 8 : 0 0 0 0.512339 0.721881 9 : 0.791373 0.627904 0.962444 0.0610695 0.0402064 R8CBB_PRINT_SOME_TEST R8CBB_PRINT_SOME prints some of an R8CBB matrix Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 Rows 1-10, Cols 9-10 Col: 9 Row 1 : 0.0317999 2 : 0.291067 3 : 0.0160248 4 : 0.865154 5 : 0.311839 6 : 0.852541 7 : 0.570552 8 : 0.0665591 9 : 0.93194 R8CBB_RANDOM_TEST R8CBB_RANDOM generates a random R8CBB matrix Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB matrix: Col: 0 1 2 3 4 Row 0 : 0.466004 0.0577937 0 0 0 1 : 0.00558318 0.82967 0.49611 0 0 2 : 0 0.979914 0.334839 0.19626 0 3 : 0 0 0.0628345 0.590166 0.501345 4 : 0 0 0 0.380084 0.891012 5 : 0 0 0 0 0.187791 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0.294381 0.570454 0.233806 0.00580827 0.0412103 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.435036 1 : 0 0 0 0 0.2737 2 : 0 0 0 0 0.0136107 3 : 0 0 0 0 0.202899 4 : 0.85826 0 0 0 0.978101 5 : 0.823066 0.210805 0 0 0.152642 6 : 0.264008 0.335545 0.507384 0 0.374974 7 : 0 0.00964202 0.135558 0.581709 0.260126 8 : 0 0 0 0.137028 0.168963 9 : 0.264625 0.570943 0.342477 0.614416 0.511269 R8CBB_SET_TEST R8CBB_SET sets the value of an element of an R8CBB matrix. Matrix order N = 5 Matrix suborder N1 = 4 Matrix suborder N2 = 1 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The matrix after additions: Col: 0 1 2 3 4 Row 0 : 11 12 0 0 0 1 : 21 22 23 0 0 2 : 31 32 33 34 0 3 : 0 42 43 44 15 4 : 0 0 0 0 25 R8CBB_SL_TEST R8CBB_SL solves a linear system factored by R8CBB_FA Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB matrix: Col: 0 1 2 3 4 Row 0 : 0.534096 0.731356 0 0 0 1 : 0.710446 0.0410714 0.236789 0 0 2 : 0 0.659148 0.0715156 0.688364 0 3 : 0 0 0.365094 0.574558 0.530478 4 : 0 0 0 0.137671 0.962985 5 : 0 0 0 0 0.0561374 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0.305328 0.300227 0.516647 0.687435 0.0302034 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0.310492 1 : 0 0 0 0 0.372682 2 : 0 0 0 0 0.599007 3 : 0 0 0 0 0.896912 4 : 0.503064 0 0 0 0.295764 5 : 0.195933 0.2317 0 0 0.0155111 6 : 0.704881 0.759678 0.338365 0 0.318406 7 : 0 0.436597 0.201918 0.784014 0.447396 8 : 0 0 0 0.280751 0.116086 9 : 0.370243 0.977503 0.375368 0.362124 0.410322 The factored matrix: Col: 0 1 2 3 4 Row 0 : 0.534096 0.731356 0 0 0 1 : -1.33019 -0.931767 0.236789 0 0 2 : 0 0.707417 0.239024 0.688364 0 3 : 0 0 -1.52744 -0.476876 0.530478 4 : 0 0 0 0.288694 1.11613 5 : 0 0 0 0 -0.0502964 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0.305328 0.300227 0.516647 0.687435 0.0302034 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 -0.257403 1 : 0 0 0 0 -0.236567 2 : 0 0 0 0 -0.760572 3 : 0 0 0 0 -0.564645 4 : 0.503064 0 0 0 -0.555744 5 : 0.17063 0.2317 0 0 0.630425 6 : -4.13104 -0.197487 0.338365 0 -0.465402 7 : 0 2.21076 0.949963 -1.39554 -1.20942 8 : 0 0 0 -0.0528655 0.116086 9 : 0.370243 0.977503 0.375368 0.0596684 -1.28927 The right hand side vector: 0: 12.1579 1: 7.75654 2: 15.8963 3: 17.4415 4: 20.1599 5: 5.45939 6: 17.145 7: 13.2971 8: 32.995 9: 24.7857 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R8CBB_TO_R8GE_TEST R8CBB_TO_R8GE converts an R8CBB matrix to R8GE format. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 201 202 203 0 0 2 : 0 302 303 304 0 3 : 0 0 403 404 405 4 : 0 0 0 504 505 5 : 0 0 0 0 605 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 1001 1002 1003 1004 1005 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 110 1 : 0 0 0 0 210 2 : 0 0 0 0 310 3 : 0 0 0 0 410 4 : 506 0 0 0 510 5 : 606 607 0 0 610 6 : 706 707 708 0 710 7 : 0 807 808 109 810 8 : 0 0 908 209 901 9 : 1006 1007 1008 1009 1010 The R8GE matrix: [[ 0. 0. 0. 0. 0. 0. 0. 0. 109. 110.] [ 201. 202. 203. 0. 0. 0. 0. 0. 209. 210.] [ 0. 302. 303. 304. 0. 0. 0. 0. 309. 310.] [ 0. 0. 403. 404. 405. 0. 0. 0. 409. 410.] [ 0. 0. 0. 504. 505. 506. 0. 0. 509. 510.] [ 0. 0. 0. 0. 605. 606. 607. 0. 609. 610.] [ 0. 0. 0. 0. 0. 706. 707. 708. 709. 710.] [ 0. 0. 0. 0. 0. 0. 807. 808. 809. 810.] [ 901. 902. 903. 904. 905. 906. 907. 908. 909. 910.] [1001. 1002. 1003. 1004. 1005. 1006. 1007. 1008. 1009. 1010.]] R8CBB_ZEROS_TEST R8CBB_ZEROS zeros an R8CBB matrix. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8CBB zero matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0 0 0 0 0 5 : 0 0 0 0 0 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 Col: 5 6 7 8 9 Row 0 : 0 0 0 0 0 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0 0 0 0 0 5 : 0 0 0 0 0 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0 0 0 0 0 9 : 0 0 0 0 0 r8cbb_test(): Normal end of execution. Sun Aug 21 20:55:02 2022