9 May 2025   9:14:44.613 PM

hankel_spd_test():
  FORTRAN90 version
  Test hankel_spd().

hankel_spd_cholesky_lower_test01():
  Test hankel_spd_cholesky_lower(),which computes a lower Cholesky
  matrix L such that the matrix H = L * L' is a
  symmetric positive definite (SPD) Hankel matrix.
 
  The Cholesky factor L:
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            0.            0.            0.            0.      
    2:      1.            1.            0.            0.            0.      
    3:      2.            1.            1.            0.            0.      
    4:      3.            3.            1.            1.            0.      
    5:      6.            4.            4.            1.            1.      
 
  The Hankel matrix H = L * L':
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            1.            2.            3.            6.      
    2:      1.            2.            3.            6.           10.      
    3:      2.            3.            6.           10.           20.      
    4:      3.            6.           10.           20.           35.      
    5:      6.           10.           20.           35.           70.      
 
  The Cholesky factor L:
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            0.            0.            0.            0.      
    2:      4.            2.            0.            0.            0.      
    3:     20.            3.            3.            0.            0.      
    4:     86.           37.            2.            4.            0.      
    5:    418.         82.5000       59.1667          1.            5.      
 
  The Hankel matrix H = L * L':
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            4.           20.           86.          418.      
    2:      4.           20.           86.          418.         1837.      
    3:     20.           86.          418.         1837.         8785.      
    4:     86.          418.         1837.         8785.         39122.8    
    5:    418.         1837.         8785.         39122.8       185057.    
 
  The Cholesky factor L:
 
  Col          1             2             3             4             5      
  Row
 
    1:  0.410323          0.            0.            0.            0.      
    2:  0.173841      0.929668          0.            0.            0.      
    3:   2.18000      0.109726      0.290234          0.            0.      
    4:   1.17220       4.99630      0.176810      0.136346          0.      
    5:   11.8167       1.18398       1.71070      0.232794      0.684076    
 
  The Hankel matrix H = L * L':
 
  Col          1             2             3             4             5      
  Row
 
    1:  0.168365      0.713310E-01  0.894504      0.480982       4.84868    
    2:  0.713310E-01  0.894504      0.480982       4.84868       3.15494    
    3:  0.894504      0.480982       4.84868       3.15494       26.3869    
    4:  0.480982       4.84868       3.15494       26.3869       20.1014    
    5:   4.84868       3.15494       26.3869       20.1014       144.486    

hankel_spd_cholesky_lower_test02():
  hankel_spd_cholesky_lower() computes a lower Cholesky
  matrix L such that the matrix H = L * L' is a
  symmetric positive definite (SPD) Hankel matrix.
 
  The Cholesky factor L:
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            0.            0.            0.            0.      
    2:      1.            1.            0.            0.            0.      
    3:      2.            1.            1.            0.            0.      
    4:      3.            3.            1.            1.            0.      
    5:      6.            4.            4.            1.            1.      
 
  The Hankel matrix H = L * L':
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            1.            2.            3.            6.      
    2:      1.            2.            3.            6.           10.      
    3:      2.            3.            6.           10.           20.      
    4:      3.            6.           10.           20.           35.      
    5:      6.           10.           20.           35.           70.      
 
  The Cholesky factor L2 of H:
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            0.            0.            0.            0.      
    2:      1.            1.            0.            0.            0.      
    3:      2.            1.            1.            0.            0.      
    4:      3.            3.            1.            1.            0.      
    5:      6.            4.            4.            1.            1.      
 
  The Hankel matrix H2 = L2 * L2':
 
  Col          1             2             3             4             5      
  Row
 
    1:      1.            1.            2.            3.            6.      
    2:      1.            2.            3.            6.           10.      
    3:      2.            3.            6.           10.           20.      
    4:      3.            6.           10.           20.           35.      
    5:      6.           10.           20.           35.           70.      

hankel_spd_test():
  Normal end of execution.

 9 May 2025   9:14:44.613 PM