28 March 2017 8:41:03.890 PM BLAS2_S_PRB FORTRAN90 version Test the BLAS2_S library. SGBMV_TEST SGBMV computes y := alpha * A * x + beta * y For a general band matrix A. Result vector y 1: 30.000000 2: 60.000000 3: 90.000000 4: 120.00000 5: 162.00000 SGEMV_TEST SGEMV 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.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. 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.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 SGER_TEST SGER 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.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. 52. 73. 94. 2: 61. 102. 143. 184. 3: 91. 152. 213. 274. 4: 121. 202. 283. 364. 5: 151. 252. 353. 454. SSBMV_TEST SSBMV computes y := alpha * A * x + beta * y for a symmetric band matrix A. Result vector y: 1: 30.000000 2: 60.000000 3: 90.000000 4: 120.00000 5: 162.00000 SSPMV_TEST SSPMV computes y := alpha * A * x + beta * y for a symmetric packed band matrix A. Result Y: 1: 30.000000 2: 60.000000 3: 90.000000 4: 120.00000 5: 162.00000 SSPR_TEST SSPR 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 + 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 SSPR2_TEST SSPR2 computes A := A + alpha*x*y'+alpha*y*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*y'+alpha*y*x': 42.00 79.00 120.00 160.00 200.00 79.00 162.00 239.00 320.00 400.00 120.00 239.00 362.00 479.00 600.00 160.00 320.00 479.00 642.00 799.00 200.00 400.00 600.00 799.00 1002.00 SSYMV_TEST SSYMV computes y := alpha * A * x + beta * y for a general symmetric matrix A. Result vector y: 1: 30.000000 2: 60.000000 3: 90.000000 4: 120.00000 5: 162.00000 SSYR_TEST SSYR computes A := A + alpha * x * x' for a symmetric 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 Vector X: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 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 SSYR2_TEST SSYR2 computes A := A + alpha*x*y' + alpha*y*x' for a symmetric 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 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 5: 50.000000 A = A + alpha*x*y' + alpha*y*x': 42.00 79.00 120.00 160.00 200.00 79.00 162.00 239.00 320.00 400.00 120.00 239.00 362.00 479.00 600.00 160.00 320.00 479.00 642.00 799.00 200.00 400.00 600.00 799.00 1002.00 STBMV_TEST STBMV computes x=A*x or x=A'*x for a triangular band matrix A. x := A*x 1: 0.0000000 2: 1.0000000 3: 2.0000000 4: 3.0000000 5: 10.000000 x := A'*x 1: 2.0000000 2: 3.0000000 3: 4.0000000 4: 5.0000000 5: 6.0000000 STRSV_TEST STRSV solves A * x = b or A' * x = b for a triangular matrix A. Right hand side b: 1: 70.000000 2: 82.000000 3: 86.000000 4: 77.000000 5: 50.000000 Solution x of A * x = b 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 Right hand side b: 1: 2.0000000 2: 11.000000 3: 32.000000 4: 70.000000 5: 130.00000 Solution x of A' * x = b 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 STPMV_TEST STPMV computes x := A * x or x := A' * x for a packed triangular matrix A. 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 STPSV_TEST STPSV computes solves A*x=b or A'*x=b for a packed triangular matrix A. b: 1: 70.000000 2: 82.000000 3: 86.000000 4: 77.000000 5: 50.000000 Solution x so that A*x=b 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 b: 1: 2.0000000 2: 11.000000 3: 32.000000 4: 70.000000 5: 130.00000 Solution x so that A'*x=b 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 STRMV_TEST STRMV computes y := A * x or y := A' * x for a triangular matrix A. 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 STRSV_TEST STRSV solves A * x = b or A' * x = b for a triangular matrix A. Right hand side b: 1: 70.000000 2: 82.000000 3: 86.000000 4: 77.000000 5: 50.000000 Solution x of A * x = b 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 Right hand side b: 1: 2.0000000 2: 11.000000 3: 32.000000 4: 70.000000 5: 130.00000 Solution x of A' * x = b 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 BLAS2_S_PRB Normal end of execution. 28 March 2017 8:41:03.892 PM