20140404  143905.853
 
BLAS2_D_PRB
  FORTRAN77 version
  Test the BLAS library.
 
TEST01
  For a general matrix A,
  DGEMV computes y := alpha * A * x + beta * y
  or             y := alpha * A' * x + beta * y.
 
  Matrix A:
 
  Col         1             2             3             4       
  Row
 
    1:   11.0000       12.0000       13.0000       14.0000    
    2:   21.0000       22.0000       23.0000       24.0000    
    3:   31.0000       32.0000       33.0000       34.0000    
    4:   41.0000       42.0000       43.0000       44.0000    
    5:   51.0000       52.0000       53.0000       54.0000    
 
  Vector X:
 
         1:    1.0000000    
         2:    2.0000000    
         3:    3.0000000    
         4:    4.0000000    
 
  Vector Y:
 
         1:    10.000000    
         2:    20.000000    
         3:    30.000000    
         4:    40.000000    
         5:    50.000000    
 
  Result Y = alpha * A  * x + beta * y
 
         1:    290.00000    
         2:    520.00000    
         3:    750.00000    
         4:    980.00000    
         5:    1210.0000    
 
  Matrix A:
 
  Col         1             2             3             4       
  Row
 
    1:   11.0000       21.0000       31.0000       41.0000    
    2:   12.0000       22.0000       32.0000       42.0000    
    3:   13.0000       23.0000       33.0000       43.0000    
    4:   14.0000       24.0000       34.0000       44.0000    
    5:   15.0000       25.0000       35.0000       45.0000    
 
  Vector X:
 
         1:    1.0000000    
         2:    2.0000000    
         3:    3.0000000    
         4:    4.0000000    
         5:    5.0000000    
 
  Vector Y:
 
         1:    10.000000    
         2:    20.000000    
         3:    30.000000    
         4:    40.000000    
 
  Result Y = alpha * A' * x + beta * y
 
         1:    440.00000    
         2:    770.00000    
         3:    1100.0000    
         4:    1430.0000    
 
TEST02
  For a general band matrix A,
  DGBMV computes 
  y := alpha * A * x + beta * y
 
  Result vector Y
 
         1:    30.000000    
         2:    60.000000    
         3:    90.000000    
         4:    120.00000    
         5:    162.00000    
 
TEST03
  For a general symmetric matrix A,
  DSYMV computes 
  y := alpha * A * x + beta * y
 
  Result vector Y
 
         1:    30.000000    
         2:    60.000000    
         3:    90.000000    
         4:    120.00000    
         5:    162.00000    
 
TEST04
  For a symmetric band matrix A,
  DSBMV computes 
  y := alpha * A * x + beta * y
 
  Result vector Y
 
         1:    30.000000    
         2:    60.000000    
         3:    90.000000    
         4:    120.00000    
         5:    162.00000    
 
TEST05
  For a general matrix A,
  DGER computes A := A + alpha * x * y'
 
  Matrix A:
 
  Col         1             2             3             4       
  Row
 
    1:   11.0000       12.0000       13.0000       14.0000    
    2:   21.0000       22.0000       23.0000       24.0000    
    3:   31.0000       32.0000       33.0000       34.0000    
    4:   41.0000       42.0000       43.0000       44.0000    
    5:   51.0000       52.0000       53.0000       54.0000    
 
  Vector X:
 
         1:    1.0000000    
         2:    2.0000000    
         3:    3.0000000    
         4:    4.0000000    
         5:    5.0000000    
 
  Vector Y:
 
         1:    10.000000    
         2:    20.000000    
         3:    30.000000    
         4:    40.000000    
 
  Result A = A + alpha * x * y
 
  Col         1             2             3             4       
  Row
 
    1:   31.0000       52.0000       73.0000       94.0000    
    2:   61.0000       102.000       143.000       184.000    
    3:   91.0000       152.000       213.000       274.000    
    4:   121.000       202.000       283.000       364.000    
    5:   151.000       252.000       353.000       454.000    
 
TEST06
  For a triangular matrix A,
  DTRMV computes y := A * x or y := A' * x
 
  Result y = A * x
 
         1:    70.000000    
         2:    82.000000    
         3:    86.000000    
         4:    77.000000    
         5:    50.000000    
 
  Result y = A' * x
 
         1:    2.0000000    
         2:    11.000000    
         3:    32.000000    
         4:    70.000000    
         5:    130.00000    
 
BLAS2_D_PRB
  Normal end of execution.
 
20140404  143905.853