02-Dec-2018 10:58:50

BLAS2_D_TEST
  MATLAB version
  Test the BLAS2_D library.

DGBMV_TEST
  DGBMV computes y := alpha * A * x + beta * y
  for a general band matrix A.

  Result vector y:

     1:           30
     2:           60
     3:           90
     4:          120
     5:          162

DGEMV_TEST
  DGEMV computes y := alpha * A * x + beta * y
  or             y := alpha * A' * x + beta * y,
  for a general matrix A.

  Matrix A:

  Col:       1             2             3             4       
  Row
      1 :          11            12            13            14  
      2 :          21            22            23            24  
      3 :          31            32            33            34  
      4 :          41            42            43            44  
      5 :          51            52            53            54  

  Vector X:

     1:            1
     2:            2
     3:            3
     4:            4

  Vector Y:

     1:           10
     2:           20
     3:           30
     4:           40
     5:           50

  Result Y = alpha * A  * x + beta * y

     1:          290
     2:          520
     3:          750
     4:          980
     5:         1210

  Matrix A:

  Col:       1             2             3             4       
  Row
      1 :          11            21            31            41  
      2 :          12            22            32            42  
      3 :          13            23            33            43  
      4 :          14            24            34            44  
      5 :          15            25            35            45  

  Vector X:

     1:            1
     2:            2
     3:            3
     4:            4
     5:            5

  Vector Y:

     1:           10
     2:           20
     3:           30
     4:           40

  Result Y = alpha * A' * x + beta * y

     1:          440
     2:          770
     3:         1100
     4:         1430

DGER_TEST
  DGER computes A := A + alpha * x * y'
  for a general matrix A.

  Matrix A:

  Col:       1             2             3             4       
  Row
      1 :          11            12            13            14  
      2 :          21            22            23            24  
      3 :          31            32            33            34  
      4 :          41            42            43            44  
      5 :          51            52            53            54  

  Vector X:

     1:            1
     2:            2
     3:            3
     4:            4
     5:            5

  Vector Y:

     1:           10
     2:           20
     3:           30
     4:           40

  Result A = A + alpha * x * y

  Col:       1             2             3             4       
  Row
      1 :          31            52            73            94  
      2 :          61           102           143           184  
      3 :          91           152           213           274  
      4 :         121           202           283           364  
      5 :         151           252           353           454  

DSBMV_TEST
  DSBMV computes y := alpha * A * x + beta * y
  for a symmetric band matrix A.

  Result Y:

     1:           30
     2:           60
     3:           90
     4:          120
     5:          162

DSPMV_TEST
  DSPMV computes y := alpha * A * x + beta * y
  for a symmetric packed band matrix A.

  Result Y:

     1:           30
     2:           60
     3:           90
     4:          120
     5:          162

DSPR_TEST
  DSPR computes A := A + alpha * x * x'
  for a symmetric packed band matrix A.

  A:

    2.00   -1.00    0.00    0.00    0.00
   -1.00    2.00   -1.00    0.00    0.00
    0.00   -1.00    2.00   -1.00    0.00
    0.00    0.00   -1.00    2.00   -1.00
    0.00    0.00    0.00   -1.00    2.00

  A = A + alpha*x*x':

    4.00    3.00    6.00    8.00   10.00
    3.00   10.00   11.00   16.00   20.00
    6.00   11.00   20.00   23.00   30.00
    8.00   16.00   23.00   34.00   39.00
   10.00   20.00   30.00   39.00   52.00

DTRMV_TEST
  DTRMV computes y := A * x or y := A' * x
  for a triangular matrix A.

  Result y = A * x

     1:           70
     2:           82
     3:           86
     4:           77
     5:           50

  Result y = A' * x

     1:            2
     2:           11
     3:           32
     4:           70
     5:          130

BLAS2_D_TEST
  Normal end of execution.

02-Dec-2018 10:58:50