18 September 2023  12:19:45.789 PM
 
blas_test():
  FORTRAN77 version
  Test blas().
 
TEST01
  CAXPY adds a multiple of one
  single precision complex vector to another.
 
  X = 
 
       1     2.0    -1.0
       2    -4.0    -2.0
       3     3.0     1.0
       4     2.0     2.0
       5    -1.0    -1.0
 
  Y = 
 
       1    -1.0     0.0
       2     0.0    -3.0
       3     4.0     0.0
       4    -3.0     4.0
       5    -2.0     0.0
 
  The scalar multiplier is:   0.500000      -1.00000    
 
  A * X + Y = 
 
       1 -1.000000 -2.500000
       2 -4.000000  0.000000
       3  6.500000 -2.500000
       4  0.000000  3.000000
       5 -3.500000  0.500000
 
test02
  DASUM adds the absolute values of 
  elements of a double precision vector.
 
  X = 
 
       1  -2.00000    
       2   4.00000    
       3  -6.00000    
       4   8.00000    
       5  -10.0000    
       6   12.0000    
       7  -14.0000    
       8   16.0000    
       9  -18.0000    
      10   20.0000    
 
  DASUM ( NX, X, 1 ) =      110.000    
  DASUM ( NX/2, X, 2 ) =    50.0000    
  DASUM ( 2, X, NX/2 ) =    14.0000    
 
  Demonstrate with a matrix A:
 
     11.0000      -12.0000       13.0000      -14.0000    
    -21.0000       22.0000      -23.0000       24.0000    
     31.0000      -32.0000       33.0000      -34.0000    
    -41.0000       42.0000      -43.0000       44.0000    
     51.0000      -52.0000       53.0000      -54.0000    
 
  DASUM(MA,A(1,2),1) =      160.000    
  DASUM(NA,A(2,1),LDA) =    90.0000    
 
TEST03
  SDOT computes the dot product of vectors.
 
  Dot product of X and Y is   -55.0000    
 
  Product of row 2 of A and X is    85.0000    
 
  Product of column 2 of A and X is    85.0000    
 
  Matrix product computed with SDOT:
 
     50.0000       30.0000       10.0000      -10.0000      -30.0000    
     60.0000       35.0000       10.0000      -15.0000      -40.0000    
     70.0000       40.0000       10.0000      -20.0000      -50.0000    
     80.0000       45.0000       10.0000      -25.0000      -60.0000    
     90.0000       50.0000       10.0000      -30.0000      -70.0000    
 
TEST04
  ZSCAL multiplies a double complex scalar 
  times a vector.
 
  X = 
 
       1    10.0     1.0
       2    20.0     2.0
       3    30.0     3.0
       4    40.0     4.0
       5    50.0     5.0
       6    60.0     6.0
 
  ZSCAL ( N, (  5.0000  0.0000), X, 1 )
 
       1    50.0     5.0
       2   100.0    10.0
       3   150.0    15.0
       4   200.0    20.0
       5   250.0    25.0
       6   300.0    30.0
 
  ZSCAL ( 3, ( -2.0000  1.0000), X, 2 )
 
       1   -21.0     8.0
       2    20.0     2.0
       3   -63.0    24.0
       4    40.0     4.0
       5  -105.0    40.0
       6    60.0     6.0
 
TEST05
  For a general matrix A,
  CGEMV computes y := alpha * A * x + beta * y
 
  Result vector Y = 
 
     283.000       72.0000    
     300.000       3.00000    
     300.000       3.00000    
     300.000       3.00000    
     355.000      -42.0000    
 
TEST06
  For a general symmetric matrix A,
  DSYMV computes 
  y := alpha * A * x + beta * y

  Result vector y:

   30.0000    
   60.0000    
   90.0000    
   120.000    
   162.000    
 
TEST07
  For a general band matrix A,
  SGBMV computes 
  y := alpha * A * x + beta * y

  Result vector y:

   30.0000    
   60.0000    
   90.0000    
   120.000    
   162.000    
 
TEST08
  For a general matrix A,
  ZGEMV computes y := alpha * A * x + beta * y
 
  Result vector Y = 
 
     283.000       72.0000    
     300.000       3.00000    
     300.000       3.00000    
     300.000       3.00000    
     355.000      -42.0000    
 
blas_test():
  Normal end of execution.
 
18 September 2023  12:19:45.790 PM