09 May 2025 9:29:09.460 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.111673 0.104672 0.255714 0.926115 2 0.962941 0.581018 0.333466 0.104342 3 0.999617 0.297442 0.658724 4 0.884015 0.241723 5 0.934900 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.955942 0.460859 0.219477E-01 0.641485 7 0.261372 0.135537 0.774612 0.887840 8 0.232689 0.665635 0.888340 9 0.131335 0.684241 10 0.795144E-01 Determinant = 0.539387E-04 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.748508 0.870291 0.151134 0.855678 2 0.307270 0.869312 0.195586 0.967871 3 0.598819 0.495847 0.134951 4 0.675378 0.941043 5 0.325241 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A*x: 1 0.748508 2 1.48483 3 3.68622 4 5.43590 5 7.73097 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.752663 0.238951 0.184401 0.494848 2 0.223552 0.530574 0.311915 0.690494 3 0.212034 0.379469 0.620931 4 0.432157E-01 0.373607 5 0.349297 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A*x: 1 3.76316 2 6.73896 3 5.25864 4 2.04090 5 1.74649 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.751250 0.517006 0.492435 0.956631 2 0.667200 0.763021 0.151288 0.484598 3 0.199588 0.268466 0.902078 4 0.236874 0.308943 5 0.805242 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.541531E-01 0.877790E-01 0.253246 0.141094 2 0.335793 0.640135 0.583745 0.969073 3 0.608223 0.501339 0.801876 4 0.678653 0.642257 5 0.130105 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.743176 0.417328 0.358434 0.727512 2 0.772454 0.785703 0.362951 0.972314 3 0.785944 0.832144 0.196476 4 0.298885 0.147830 5 0.289748 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.651489 0.813560 0.822073 0.873779 7 0.378553 0.267838 0.271075 0.525659 8 0.764381 0.275927 0.945327 9 0.585105 0.587728 10 0.185277 b: 1 5.56318 2 10.2154 3 9.96345 4 4.91850 5 9.20069 6 24.0444 7 12.4888 8 18.0517 9 11.1432 10 1.85277 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.396614 0.603108 0.672639 0.166354 2 0.698350 0.536745 0.818765 0.579919 3 0.580941 0.970004 0.793159 4 0.222763 0.995919 5 0.105873 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.935469 0.680999 0.671875 0.782471 7 0.238548 0.930483 0.660275 0.221724 8 0.277897E-01 0.558661 0.437259 9 0.425546 0.636326 10 0.264716 b: 1 0.396614 2 1.99981 3 3.48895 4 5.60495 5 8.05236 6 13.2319 7 10.4756 8 12.3346 9 17.6159 10 13.4242 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.153968 0.695819 0.305918 0.737510 2 0.245661E-01 0.444355 0.459150 0.364940 3 0.797407 0.718703 0.841591 4 0.827690 0.476303 5 0.276903 The R8GE matrix Col: 1 2 3 4 5 Row --- 1 0.153968 0.695819 0.305918 0.737510 0.00000 2 0.00000 0.245661E-01 0.444355 0.459150 0.364940 3 0.00000 0.00000 0.797407 0.718703 0.841591 4 0.00000 0.00000 0.00000 0.827690 0.476303 5 0.00000 0.00000 0.00000 0.00000 0.276903 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. 09 May 2025 9:29:09.460 PM