DBLAS_PRB Tests for DBLAS, the Double Precision basic linear algebra subprograms. Today's date: 20020624 Today's time: 105138.584 TEST05 IDAMAX returns the index of maximum magnitude; The vector X: 1 2.0000 2 -2.0000 3 5.0000 4 1.0000 5 -3.0000 6 4.0000 7 0.0000 8 -4.0000 9 3.0000 10 -1.0000 11 -5.0000 The index of maximum magnitude = 3 TEST99 Demonstrate the use of SASUM which adds the absolute values of elements of a vector. X = 1, -2. 2, 4. 3, -6. 4, 8. 5, -10. 6, 12. 7, -14. 8, 16. 9, -18. 10, 20. DASUM ( NX, X, 1 ) = 11.6875 DASUM ( NX/2, X, 2 ) = 11.6875 DASUM ( 2, X, NX/2 ) = 2. Demonstrate with a matrix A: 11., -12., 13., -14. -21., 22., -23., 24. 31., -32., 33., -34. -41., 42., -43., 44. 51., -52., 53., -54. DASUM ( MA, A(1,2), 1 ) = 8.46875 DASUM ( NA, A(2,1), LDA ) = 11.9765625 TEST01 DDOT computes the dot product of vectors. Dot product of X and Y is -55. Product of row 2 of A and X is 85. Product of column 2 of A and X is 85. Matrix product computed with ddot: 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 TEST02 DNRM2 computes the Euclidean norm of a vector. The vector X: 1 1.0000 2 2.0000 3 3.0000 4 4.0000 5 5.0000 The 2-norm of X is 7.416198487095663 The 2-norm of row 2 of A is 11.618950038622252 The 2-norm of column 2 of A is 11.618950038622252 TEST03 Use IDAMAX, DAXPY, DDOT and DSCAL in a Gauss elimination routine. First five entries of solution: 1.00000 2.00000 3.00000 4.00000 5.00000 TEST04 Demonstrate Cholesky factor/solve routines built out of level 1, 2 and 3 blas First 10 entries of solution: 1.00000 2.00000 3.00000 4.00000 5.00000 6.00000 7.00000 8.00000 9.00000 10.0000 DBLAS_PRB Normal end of DBLAS tests.