6 October 2025   6:18:58.486 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.737107          0.            0.            0.            0.      
    2:  0.121281E-01  0.636112          0.            0.            0.      
    3:  0.549154      0.948721      0.624258          0.            0.      
    4:  0.827766       2.48588      0.946426      0.833717          0.      
    5:   2.15890       5.30978       3.57634      0.178518      0.371458    
 
  The Hankel matrix H = L * L':
 
  Col          1             2             3             4             5      
  Row
 
    1:  0.543327      0.893972E-02  0.404785      0.610152       1.59134    
    2:  0.893972E-02  0.404785      0.610152       1.59134       3.40380    
    3:  0.404785      0.610152       1.59134       3.40380       8.45563    
    4:  0.610152       1.59134       3.40380       8.45563       18.5201    
    5:   1.59134       3.40380       8.45563       18.5201       45.8147    

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.

 6 October 2025   6:18:58.486 PM