06 October 2025 6:33:30.801 PM r8but_test(): Fortran90 version: Test r8but(). R8BUT_DET_TEST R8BUT_DET computes the determinant of an R8BUT matrix. Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.188664 0.615367 0.879385 0.251081 2 0.522449 0.942406 0.952775 0.192764 3 0.226079 0.126569 0.337081 4 0.862850 0.160726 5 0.859265 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.985304 0.299516 0.312571 0.917303 7 0.723570 0.473369 0.716847 0.304683 8 0.473224 0.783436 0.827496E-01 9 0.679290 0.459553 10 0.328646 Determinant = 0.124439E-02 R8BUT_INDICATOR_TEST R8BUT_INDICATOR sets up an R8BUT indicator matrix. Matrix order N = 6 Upper bandwidth MU = 2 The R8BUT indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 2 22.0000 23.0000 24.0000 3 33.0000 34.0000 35.0000 4 44.0000 45.0000 5 55.0000 6 Col: 6 Row --- 6 66.0000 R8BUT_MTV_TEST R8BUT_MTV computes b=A'*x, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.508323 0.203710 0.361804 0.532887 2 0.727039 0.317615 0.913331 0.599496 3 0.289964 0.592793 0.241227 4 0.796148 0.554503 5 0.350646 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A*x: 1 0.508323 2 1.65779 3 1.86693 4 7.32252 5 5.89391 R8BUT_MV_TEST R8BUT_MV computes b=A*x, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.525680E-01 0.275591 0.564405 0.101489E-01 2 0.782301 0.420080 0.610246 0.625325 3 0.264896 0.952019 0.932540 4 0.205491 0.356163 5 0.876340 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A*x: 1 2.33756 2 8.39245 3 9.26546 4 2.60278 5 4.38170 R8BUT_PRINT_TEST R8BUT_PRINT prints an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.510313 0.458115 0.974704 0.182337 2 0.907295 0.977329 0.880445 0.149698 3 0.295358 0.934124 0.266029 4 0.386011 0.555978 5 0.385895 R8BUT_PRINT_SOME_TEST R8BUT_PRINT_SOME prints some of an R8BUT matrix. Matrix order N = 10 Upper bandwidth MU = 3 Rows 1:4, Cols 2:4: Col: 2 3 4 Row --- 2 202.000 203.000 204.000 3 303.000 304.000 4 404.000 R8BUT_RANDOM_TEST R8BUT_RANDOM randomizes an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.450842E-01 0.551030 0.351933 0.752445 2 0.675654E-03 0.347731 0.327927 0.249294 3 0.192585 0.687915 0.447023 4 0.978413 0.807448 5 0.152758 R8BUT_SL_TEST R8BUT_SL solves A*x=b, where A is an R8BUT matrix. Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.868690 0.309387 0.592218 0.627443 2 0.983746E-01 0.596187 0.202975 0.152572 3 0.246400 0.521771 0.201240 4 0.699330 0.811579 5 0.210770 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.592368E-01 0.503818 0.804264 0.110018E-01 7 0.572953 0.871796 0.182394 0.417606 8 0.400094 0.449341 0.321607 9 0.129965 0.251039 10 0.692596 b: 1 5.77389 2 3.56007 3 9.57895 4 13.9048 5 11.8000 6 10.4153 7 16.8026 8 10.4609 9 3.68007 10 6.92596 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 9 9.00000 10 10.0000 R8BUT_SLT_TEST R8BUT_SLT solves A'*x=b, where A is an R8BUT matrix. Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.198207 0.838616 0.475795 0.220528 2 0.547056 0.632769 0.943823 0.473195 3 0.935864 0.831949 0.792573 4 0.224608 0.355230 5 0.675004 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.127114 0.152373 0.853424 0.375099E-01 7 0.594216 0.264225 0.138992 0.152629 8 0.323231E-01 0.788227 0.648037 9 0.642213 0.765299 10 0.719428 b: 1 0.198207 2 1.93273 3 4.54893 4 5.50245 5 8.12005 6 5.02335 7 10.9607 8 8.56227 9 13.2837 10 20.3347 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 9 9.00000 10 10.0000 R8BUT_TO_R8GE_TEST R8BUT_TO_R8GE converts a matrix from R8BUT to R8GE format. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.128188 0.683342 0.325726 0.729111 2 0.758278 0.756240 0.486434 0.278310 3 0.672651 0.967422 0.503119E-01 4 0.437179 0.569316 5 0.322081 The R8GE matrix Col: 1 2 3 4 5 Row --- 1 0.128188 0.683342 0.325726 0.729111 0.00000 2 0.00000 0.758278 0.756240 0.486434 0.278310 3 0.00000 0.00000 0.672651 0.967422 0.503119E-01 4 0.00000 0.00000 0.00000 0.437179 0.569316 5 0.00000 0.00000 0.00000 0.00000 0.322081 R8BUT_ZEROS_TEST R8BUT_ZEROS zeros an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 5 0.00000 r8but_test(): Normal end of execution. 06 October 2025 6:33:30.801 PM