Tue May 20 22:27:31 2025

r83_np_test():
  python version: 3.10.12
  numpy version:  1.26.4
  Test r83_np().

r83_np_det_test():
  r83_np_det() computes the determinant of a tridiagonal
  matrix factored by r83_np_fa.

  Matrix order N =  10

  The factored R83 matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:             2            -1                                          
    1:          -0.5           1.5            -1                            
    2:                   -0.666667       1.33333            -1              
    3:                                     -0.75          1.25            -1
    4:                                                    -0.8           1.2
    5:                                                             -0.833333

  Col:       5             6             7             8             9       
  Row
  ---
    4:            -1                                                        
    5:       1.16667            -1                                          
    6:     -0.857143       1.14286            -1                            
    7:                      -0.875         1.125            -1              
    8:                                 -0.888889       1.11111            -1
    9:                                                    -0.9           1.1

  r83_np_det() computes determinant =  10.999999999999998
  Exact determinant is =  11

r83_np_fa_test():
  r83_np_fa() factors a tridiagonal matrix with no pivoting

  Matrix order N =  10

  The tridiagonal matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:      0.982958      0.797004                                          
    1:      0.520802      0.114438      0.836984                            
    2:                    0.131832      0.717529     0.0415045              
    3:                                  0.769138     0.0731872      0.962241
    4:                                                0.907622       0.25375
    5:                                                             0.0156681

  Col:       5             6             7             8             9       
  Row
  ---
    4:       0.59729                                                        
    5:      0.387392      0.317376                                          
    6:       0.43372      0.893855      0.199122                            
    7:                    0.472497      0.150687      0.594193              
    8:                                  0.589586      0.581731      0.217997
    9:                                                0.100327      0.256817

  Solution:

     0:             1
     1:             2
     2:             3
     3:             4
     4:             5
     5:             6
     6:             7
     7:             8
     8:             9
     9:            10

  Solution to tranposed system:

     0:             1
     1:             2
     2:             3
     3:             4
     4:             5
     5:             6
     6:             7
     7:             8
     8:             9
     9:            10

r83_np_fs_test():
  r83_np_fs() factors and solves a tridiagonal linear system.

  Matrix order N =  10

  Solution:

     0:             1
     1:             2
     2:             3
     3:             4
     4:             5
     5:             6
     6:             7
     7:             8
     8:             9
     9:            10

r83_np_ml_test():
  r83_np_ml() computes A*x or A'*x
  where A has been factored by r83_fa().

  Matrix order N =  10

  A*x and PLU*x:

     0:        1.93047         1.93047
     1:        4.25772         4.25772
     2:        4.21906         4.21906
     3:        4.85608         4.85608
     4:        4.86612         4.86612
     5:         10.727          10.727
     6:        11.8251         11.8251
     7:        11.3579         11.3579
     8:        18.0066         18.0066
     9:        9.57299         9.57299

  A'*x and (PLU)'*x

     0:        2.20674         2.20674
     1:        3.20924         3.20924
     2:        1.22748         1.22748
     3:        7.00161         7.00161
     4:        10.4903         10.4903
     5:        13.5317         13.5317
     6:        12.8859         12.8859
     7:        11.7683         11.7683
     8:        14.8257         14.8257
     9:        8.08427         8.08427

r83_np_sl_test():
  r83_np_sl() solves a linear system after the tridiagonal
  matrix has been factored by r83_np_fa().

  Matrix order N =  10

  The tridiagonal matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:      0.669602      0.774832                                          
    1:      0.582059      0.461572      0.512386                            
    2:                    0.231519     0.0810168      0.738669              
    3:                                  0.412376      0.505719      0.802806
    4:                                                0.472111      0.936768
    5:                                                             0.0708405

  Col:       5             6             7             8             9       
  Row
  ---
    4:      0.897983                                                        
    5:      0.324646      0.332334                                          
    6:      0.729679       0.64428      0.561154                            
    7:                    0.612538     0.0252873      0.618898              
    8:                                  0.594636      0.937497      0.676942
    9:                                                0.925259      0.125711

  Solution:

     0:             1
     1:             2
     2:             3
     3:             4
     4:             5
     5:             6
     6:             7
     7:             8
     8:             9
     9:            10

  Solution to tranposed system:

     0:             1
     1:             2
     2:             3
     3:             4
     4:             5
     5:             6
     6:             7
     7:             8
     8:             9
     9:            10

r83_np_test():
  Normal end of execution.
Tue May 20 22:27:32 2025