Mon Aug 22 19:37:09 2022

r85_test():
  Python version: 3.6.9
  Test R85().

R85_DIF2_TEST
  R85_DIF2 sets up a R85 second difference matrix.

  Matrix order N =  5

  The R85 second difference matrix:

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

R85_INDICATOR_TEST
  R85_INDICATOR sets up a R85 indicator matrix.

  Matrix order N =  5

  The R85 indicator matrix:

  Col:       0             1             2             3             4       
  Row
      0 :          11            12            13                              
      1 :          21            22            23            24                
      2 :          31            32            33            34            35  
      3 :                        42            43            44            45  
      4 :                                      53            54            55  

R85_MTV_TEST
  R85_MTV computes b=A'*x, where A is an R85 matrix.

  Matrix order N =  5

  The R85 matrix:

  Col:       0             1             2             3             4       
  Row
      0 :          11            12            13                              
      1 :          21            22            23            24                
      2 :          31            32            33            34            35  
      3 :                        42            43            44            45  
      4 :                                      53            54            55  

  The vector x:

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

  The product b=A'*x:

     0:            74
     1:           230
     2:           505
     3:           614
     4:           650

R85_MV_TEST
  R85_MV computes b=A*x, where A is an R85 matrix.

  Matrix order N =  5

  The R85 matrix:

  Col:       0             1             2             3             4       
  Row
      0 :          11            12            13                              
      1 :          21            22            23            24                
      2 :          31            32            33            34            35  
      3 :                        42            43            44            45  
      4 :                                      53            54            55  

  The vector x:

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

  The product b=A*x:

     0:           146
     1:           320
     2:           595
     3:           596
     4:           560

R85_NP_FS_TEST
  R85_NP_FS factors and solves a pentadiagonal
  linear system, with no pivoting.

  Matrix order N =  10

  The pentadiagonal matrix:

  Col:       0             1             2             3             4       
  Row
      0 :    0.073011      0.413359      0.963101                              
      1 :    0.365398      0.593823      0.815742      0.861717                
      2 :    0.664827      0.865385      0.940526     0.0403024      0.881022  
      3 :                   0.46526     0.0978217      0.725564     0.0897018  
      4 :                                0.894938      0.304446      0.341678  
      5 :                                              0.583939     0.0797615  
      6 :                                                            0.150225  

  Col:       5             6             7             8             9       
  Row
      3 :    0.257548                                                          
      4 :    0.802664      0.181949                                            
      5 :    0.103542      0.622108      0.787917                              
      6 :    0.650034       0.98223      0.829524      0.958816                
      7 :    0.132616      0.055362      0.554453     0.0123034      0.169606  
      8 :                  0.146515      0.984925      0.313077      0.701836  
      9 :                                0.608923      0.863459      0.245265  

  The right hand side:

     0:       2.79829
     1:        6.0582
     2:       10.2821
     3:       9.77246
     4:       6.24041
     5:       11.2759
     6:       13.2795
     7:       29.9234
     8:       18.2624
     9:        10.126

  Solution:

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

R85_PRINT_TEST
  R85_PRINT prints an R85 matrix.

  Matrix order N =  5

  The R85 matrix:

  Col:       0             1             2             3             4       
  Row
      0 :          11            12            13                              
      1 :          21            22            23            24                
      2 :          31            32            33            34            35  
      3 :                        42            43            44            45  
      4 :                                      53            54            55  

R85_PRINT_SOME_TEST
  R85_PRINT_SOME prints some of an R85 matrix.

  Matrix order N =  8

  Rows 1-5, Cols 2-4:

  Col:       2             3             4       
  Row
      1 :          23            24                
      2 :          33            34            35  
      3 :          43            44            45  
      4 :          53            54            55  
      5 :                        64            65  

R85_RANDOM_TEST
  R85_RANDOM sets up a random R85 matrix.

  Matrix order N =  5

  The R85 random matrix:

  Col:       0             1             2             3             4       
  Row
      0 :     0.65429      0.762833        0.8834                              
      1 :    0.348314      0.879689      0.880899      0.763473                
      2 :    0.179312      0.876165       0.95641      0.640945      0.497574  
      3 :                  0.461735      0.678476      0.591083       0.35201  
      4 :                                0.887339      0.830445      0.514588  

R85_TO_R8GE_TEST
  R85_TO_R8GE converts an R85 matrix to R8GE form.

  Matrix order N =  5

  The R85 matrix:

  Col:       0             1             2             3             4       
  Row
      0 :          11            12            13                              
      1 :          21            22            23            24                
      2 :          31            32            33            34            35  
      3 :                        42            43            44            45  
      4 :                                      53            54            55  

  The R8GE matrix:
[[11. 21. 31.  0.  0.]
 [12. 22. 32. 42.  0.]
 [13. 23. 33. 43. 53.]
 [ 0. 24. 34. 44. 54.]
 [ 0.  0. 35. 45. 55.]]

R85_ZEROS_TEST
  R85_ZEROS sets up a R85 zero matrix.

  Matrix order N =  5

  The R85 zero matrix:

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

r85_test():
  Normal end of execution.
Mon Aug 22 19:37:10 2022