Tue Aug 16 12:06:57 2022

r83_np_test():
  Python version: 3.6.9
  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.543081      0.885338                                          
    1:      0.299569      0.342346      0.568441                            
    2:                    0.292904      0.293687      0.649803              
    3:                                 0.0724743      0.563777     0.0547356
    4:                                                0.854588      0.404391
    5:                                                              0.970195

  Col:       5             6             7             8             9       
  Row
  ---
    4:      0.851008                                                        
    5:      0.956989      0.977807                                          
    6:      0.776726      0.957145      0.200744                            
    7:                    0.383979      0.488641      0.518432              
    8:                                  0.344702      0.849382     0.0555662
    9:                                                0.545662      0.440529

  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.37674         1.37674
     1:        2.84778         2.84778
     2:         2.8591          2.8591
     3:         6.8071          6.8071
     4:        4.87787         4.87787
     5:        8.13162         8.13162
     6:        12.9133         12.9133
     7:        5.86067         5.86067
     8:        15.2796         15.2796
     9:        3.75188         3.75188

  A'*x and (PLU)'*x

     0:       0.713395        0.713395
     1:        2.83836         2.83836
     2:        4.31231         4.31231
     3:        2.69593         2.69593
     4:        9.64095         9.64095
     5:        8.08702         8.08702
     6:        6.72158         6.72158
     7:        9.50925         9.50925
     8:        17.9782         17.9782
     9:        7.06252         7.06252

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.799796      0.127965                                          
    1:      0.639924      0.154573      0.602818                            
    2:                    0.849302       0.34825      0.891448              
    3:                                  0.358966      0.853392      0.794225
    4:                                                0.281341        0.6686
    5:                                                              0.537398

  Col:       5             6             7             8             9       
  Row
  ---
    4:      0.321289                                                        
    5:      0.897704     0.0481679                                          
    6:      0.955134      0.569751      0.965453                            
    7:                    0.285188      0.309461      0.574781              
    8:                                  0.634949      0.948259      0.952792
    9:                                                0.894094      0.818907

  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 Aug 16 12:06:57 2022