5 October 2018   9:08:50.857 AM
 
BLAS3_D_TEST
  FORTRAN90 version
  Test the BLAS library.

DGEMM_TEST
  DGEMM multiplies two matrices A and B.

  1: C = alpha * A  * B  + beta * C;
  2: C = alpha * A' * B  + beta * C;
  3: C = alpha * A  * B' + beta * C;
  4: C = alpha * A' * B' + beta * C;

  We carry out all four calculations, but in each case,
  we choose our input matrices so that we get the same result.
 
  C = alpha * A * B + beta * C:
 
  Col          1             2             3             4             5      
  Row
 
    1:   1585.         1660.         1735.         1810.         1885.      
    2:   2875.         3010.         3145.         3280.         3415.      
    3:   4165.         4360.         4555.         4750.         4945.      
    4:   5455.         5710.         5965.         6220.         6475.      
 
  C = alpha * A' * B + beta * C:
 
  Col          1             2             3             4             5      
  Row
 
    1:   1585.         1660.         1735.         1810.         1885.      
    2:   2875.         3010.         3145.         3280.         3415.      
    3:   4165.         4360.         4555.         4750.         4945.      
    4:   5455.         5710.         5965.         6220.         6475.      
 
  C = alpha * A * B' + beta * C:
 
  Col          1             2             3             4             5      
  Row
 
    1:   1585.         1660.         1735.         1810.         1885.      
    2:   2875.         3010.         3145.         3280.         3415.      
    3:   4165.         4360.         4555.         4750.         4945.      
    4:   5455.         5710.         5965.         6220.         6475.      
 
  C = alpha * A' * B' + beta * C:
 
  Col          1             2             3             4             5      
  Row
 
    1:   1585.         1660.         1735.         1810.         1885.      
    2:   2875.         3010.         3145.         3280.         3415.      
    3:   4165.         4360.         4555.         4750.         4945.      
    4:   5455.         5710.         5965.         6220.         6475.      

DTRMM_TEST
  DTRMM multiplies a triangular matrix A and a
  rectangular matrix B

  1: B = alpha * A  * B;
  2: B = alpha * A' * B;
 
  B = alpha * A * B:
 
  Col          1             2             3             4             5      
  Row
 
    1:    828.          856.          884.          912.          940.      
    2:    970.         1000.         1030.         1060.         1090.      
    3:    946.          972.          998.         1024.         1050.      
    4:    656.          672.          688.          704.          720.      
 
  B = alpha * A * B:
 
  Col          1             2             3             4             5      
  Row
 
    1:     44.           48.           52.           56.           60.      
    2:    234.          248.          262.          276.          290.      
    3:    670.          700.          730.          760.          790.      
    4:   1452.         1504.         1556.         1608.         1660.      

DTRSM_TEST
  DTRSM solves a linear system involving a triangular
  matrix A and a rectangular matrix B.

  1: Solve A  * X  = alpha * B;
  2: Solve A' * X  = alpha * B;
  3: Solve X  * A  = alpha * B;
  4: Solve X  * A' = alpha * B;
 
  X = inv ( A ) * alpha * B:
 
  Col          1             2             3             4             5      
  Row
 
    1:  -7.10938      -6.92708      -6.74479      -6.56250      -6.38021    
    2:  -2.84375      -2.77083      -2.69792      -2.62500      -2.55208    
    3:  -1.62500      -1.58333      -1.54167      -1.50000      -1.45833    
    4:   10.2500       10.5000       10.7500         11.         11.2500    
 
  X = inv ( A' ) * alpha * B:
 
  Col          1             2             3             4             5      
  Row
 
    1:     11.           12.           13.           14.           15.      
    2:   2.25000          2.         1.75000       1.50000       1.25000    
    3:   1.12500          1.        0.875000      0.750000      0.625000    
    4:  0.703125      0.625000      0.546875      0.468750      0.390625    
 
  X = alpha * B * inv ( A ):
 
  Col          1             2             3             4             5      
  Row
 
    1:     11.        -2.25000      -1.12500     -0.703125     -0.492188    
    2:     21.        -4.75000      -2.37500      -1.48438      -1.03906    
    3:     31.        -7.25000      -3.62500      -2.26562      -1.58594    
    4:     41.        -9.75000      -4.87500      -3.04688      -2.13281    
 
  X = alpha * B * inv ( A' ):
 
  Col          1             2             3             4             5      
  Row
 
    1:  0.820312      0.328125      0.187500      0.125000          3.      
    2:   2.46094      0.984375      0.562500      0.375000          5.      
    3:   4.10156       1.64062      0.937500      0.625000          7.      
    4:   5.74219       2.29688       1.31250      0.875000          9.      
 
BLAS3_D_TEST
  Normal end of execution.
 
 5 October 2018   9:08:50.864 AM