Wed Aug 17 12:17:59 2022

r83p_test():
  Python version: 3.6.9
  Test r83p().

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.867263     0.0363075                                          
    1:      0.719711       0.85351      0.126684                            
    2:                    0.486726      0.605199      0.260111              
    3:                                  0.241606      0.730423      0.120989
    4:                                                0.596943      0.664049
    5:                                                              0.803922

  Col:       5             6             7             8             9       
  Row
  ---
    4:      0.688299                                                        
    5:      0.526802      0.309067                                          
    6:      0.675371      0.661856      0.359801                            
    7:                    0.215644      0.111133      0.421912              
    8:                                 0.0930295      0.997044       0.28219
    9:                                               0.0498144      0.729357

  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

r83p_det_test():
  r83p_det() computes determinant of a tridiagonal periodic matrix.

  The periodic tridiagonal matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:           11            12                                        15  
    1:           21            22            23                              
    2:                         32            33            34                
    3:                                       43            44            45  
    4:           51                                        54            55  

  r83p_det() computes the determinant =  70772459.99999996
  r8ge_det() computes the determinant =  70772460.0
  np.linalg.det() computes            =  70772460.00000009

r83p_fa_test():
  r83p_fa() factors a tridiagonal periodic system
  which then can be solved by R83P_SL.

  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

  Solution to transposed system:

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

R83P_INDICATOR_TEST
  R83P_INDICATOR sets up an R83P indicator matrix.

  Matrix order N =  10

  The R83P indicator matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:          101           102                                            
    1:          201           202           203                              
    2:                        302           303           304                
    3:                                      403           404           405  
    4:                                                    504           505  
    5:                                                                  605  
    6:                                                                       
    7:                                                                       
    8:                                                                       
    9:         1001                                                          

  Col:       5             6             7             8             9       
  Row
  ---
    0:                                                                  110  
    1:                                                                       
    2:                                                                       
    3:                                                                       
    4:          506                                                          
    5:          606           607                                            
    6:          706           707           708                              
    7:                        807           808           809                
    8:                                      908           909           910  
    9:                                                   1009          1010  

r83p_ml_test()
  r83p_ml() computes A*x or A'*X
  where A has been factored by R83P_FA.

  Matrix order N =  10

  A*x and PLU*x

     0:        11.8711         11.8711
     1:        3.63858         3.63858
     2:        3.12429         3.12429
     3:        7.20759         7.20759
     4:        3.49942         3.49942
     5:        11.6599         11.6599
     6:        4.35901         4.35901
     7:        5.25943         5.25943
     8:        13.7402         13.7402
     9:        16.5588         16.5588

  A'*x and (PLU)'*x

     0:        5.88769         5.88769
     1:         1.6948          1.6948
     2:        4.25169         4.25169
     3:        5.29126         5.29126
     4:        10.3237         10.3237
     5:        12.7848         12.7848
     6:        11.3889         11.3889
     7:        6.07364         6.07364
     8:        16.0966         16.0966
     9:        5.22395         5.22395

R83P_MTV_TEST
  R83P_MTV computes A'*x=b for an R83P matrix.

  Matrix order N =  5

  The R83P matrix A:

  Col:       0             1             2             3             4       
  Row
  ---
    0:           11            12                                        15  
    1:           21            22            23                              
    2:                         32            33            34                
    3:                                       43            44            45  
    4:           51                                        54            55  

  The vector X:

     0:             1
     1:             2
     2:             3
     3:             4
     4:             5

  The product b = A'*x:

     0:           308
     1:           152
     2:           317
     3:           548
     4:           470

R83P_MV_TEST
  R83P_MV computes A*x=b for an R83P matrix.

  Matrix order N =  5

  The R83P matrix A:

  Col:       0             1             2             3             4       
  Row
  ---
    0:           11            12                                        15  
    1:           21            22            23                              
    2:                         32            33            34                
    3:                                       43            44            45  
    4:           51                                        54            55  

  The vector X:

     0:             1
     1:             2
     2:             3
     3:             4
     4:             5

  The product b = A*x:

     0:           110
     1:           134
     2:           299
     3:           530
     4:           542

R83P_PRINT_TEST
  R83P_PRINT prints an R83P matrix.

  Matrix order N =  5

  The R83P matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:           11            12                                        15  
    1:           21            22            23                              
    2:                         32            33            34                
    3:                                       43            44            45  
    4:           51                                        54            55  

R83P_PRINT_SOME_TEST
  R83P_PRINT_SOME prints some of an R83P matrix.

  Matrix order N =  10

  Rows 1:N, Cols 1:2:

  Col:       1             2       
  Row
  ---
    1:          202           203  
    2:          302           303  
    3:                        403  

R83P_RANDOM_TEST
  R83P_RANDOM sets up a random R83P matrix.

  Matrix order N =  5

  The R83P matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:     0.356607      0.688846                                   0.19011  
    1:      0.55198      0.578535      0.490332                              
    2:                   0.581784      0.344571      0.502409                
    3:                                 0.425894      0.560983     0.0757479  
    4:     0.606516                                  0.745573      0.702794  

R83P_SL_TEST
  R83P_SL solves a tridiagonal periodic system
  after it has been factored by R83P_FA.

  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

  Solution to transposed system:

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

R83P_TO_R8GE_TEST
  R83P_TO_R8GE converts an r83p matrix to R8GE format.

  Matrix order N =  5

  The R83P matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:     0.705881      0.275885                                  0.281176  
    1:     0.243801      0.920581     0.0828277                              
    2:                   0.981555       0.31577      0.212513                
    3:                                  0.54383       0.68663      0.455136  
    4:     0.801106                                  0.128439      0.880292  

  The R8GE matrix:

  Col:       0             1             2             3             4       
  Row
      0 :    0.705881      0.275885             0             0      0.281176  
      1 :    0.243801      0.920581     0.0828277             0             0  
      2 :           0      0.981555       0.31577      0.212513             0  
      3 :           0             0       0.54383       0.68663      0.455136  
      4 :    0.801106             0             0      0.128439      0.880292  

R83P_ZEROS_TEST
  R83P_ZEROS sets up a zero R83P matrix.

  Matrix order N =  5

  The R83P matrix:

  Col:       0             1             2             3             4       
  Row
  ---
    0:            0             0                                         0  
    1:            0             0             0                              
    2:                          0             0             0                
    3:                                        0             0             0  
    4:            0                                         0             0  

r83p_test():
  Normal end of execution.
Wed Aug 17 12:17:59 2022