02-Dec-2018 09:14:55 BLAS1_D_TEST MATLAB version Test the BLAS1_D library. DASUM_TEST DASUM adds the absolute values of elements of a vector. X = 1 -2.000000 2 4.000000 3 -6.000000 4 8.000000 5 -10.000000 6 12.000000 7 -14.000000 8 16.000000 9 -18.000000 10 20.000000 DASUM ( NX, X, 1 ) = 110.000000 DASUM ( NX/2, X, 2 ) = 50.000000 DASUM ( 2, X, NX/2 ) = 14.000000 Demonstrate with a matrix A: 11.000000 -12.000000 13.000000 -14.000000 -21.000000 22.000000 -23.000000 24.000000 31.000000 -32.000000 33.000000 -34.000000 -41.000000 42.000000 -43.000000 44.000000 51.000000 -52.000000 53.000000 -54.000000 DASUM(MA,A(1:MA,2),1) = 160.000000 DASUM(NA,A(2,1:NA),1) = 90.000000 DAXPY_TEST DAXPY adds a multiple of vector X to vector Y. X = 1 1.000000 2 2.000000 3 3.000000 4 4.000000 5 5.000000 6 6.000000 Y = 1 100.000000 2 200.000000 3 300.000000 4 400.000000 5 500.000000 6 600.000000 DAXPY ( N, 1.000000, X, 1, Y, 1 ) 1 101.000000 2 202.000000 3 303.000000 4 404.000000 5 505.000000 6 606.000000 DAXPY ( N, -2.000000, X, 1, Y, 1 ) 1 98.000000 2 196.000000 3 294.000000 4 392.000000 5 490.000000 6 588.000000 DAXPY ( 3, 3.000000, X, 2, Y, 1 ) 1 103.000000 2 209.000000 3 315.000000 4 400.000000 5 500.000000 6 600.000000 DAXPY ( 3, -4.000000, X, 1, Y, 2 ) 1 96.000000 2 200.000000 3 292.000000 4 400.000000 5 488.000000 6 600.000000 DCOPY_TEST DCOPY copies one vector into another. X = 1 1.000000 2 2.000000 3 3.000000 4 4.000000 5 5.000000 6 6.000000 7 7.000000 8 8.000000 9 9.000000 10 10.000000 Y = 1 10.000000 2 20.000000 3 30.000000 4 40.000000 5 50.000000 6 60.000000 7 70.000000 8 80.000000 9 90.000000 10 100.000000 A = 11.000000 12.000000 13.000000 14.000000 15.000000 21.000000 22.000000 23.000000 24.000000 25.000000 31.000000 32.000000 33.000000 34.000000 35.000000 41.000000 42.000000 43.000000 44.000000 45.000000 51.000000 52.000000 53.000000 54.000000 55.000000 DCOPY ( 5, X, 1, Y, 1 ) 1 1.000000 2 2.000000 3 3.000000 4 4.000000 5 5.000000 6 60.000000 7 70.000000 8 80.000000 9 90.000000 10 100.000000 DCOPY ( 3, X, 2, Y, 3 ) 1 1.000000 2 20.000000 3 30.000000 4 3.000000 5 50.000000 6 60.000000 7 5.000000 8 80.000000 9 90.000000 10 100.000000 DCOPY ( 5, X, 1, A, 1 ) A = 1.000000 12.000000 13.000000 14.000000 15.000000 2.000000 22.000000 23.000000 24.000000 25.000000 3.000000 32.000000 33.000000 34.000000 35.000000 4.000000 42.000000 43.000000 44.000000 45.000000 5.000000 52.000000 53.000000 54.000000 55.000000 DCOPY ( 5, X, 2, A, 5 ) A = 1.000000 3.000000 5.000000 7.000000 9.000000 21.000000 22.000000 23.000000 24.000000 25.000000 31.000000 32.000000 33.000000 34.000000 35.000000 41.000000 42.000000 43.000000 44.000000 45.000000 51.000000 52.000000 53.000000 54.000000 55.000000 DDOT_TEST DDOT computes the dot product of vectors. Dot product of X and Y is -55.000000 Product of row 2 of A and X is 85.000000 Product of column 2 of A and X is 85.000000 Matrix product computed with DDOT: 50.000000 30.000000 10.000000 -10.000000 -30.000000 60.000000 35.000000 10.000000 -15.000000 -40.000000 70.000000 40.000000 10.000000 -20.000000 -50.000000 80.000000 45.000000 10.000000 -25.000000 -60.000000 90.000000 50.000000 10.000000 -30.000000 -70.000000 DNRM2_TEST DNRM2 computes the Euclidean norm of a vector. The vector X: 1 1.000000 2 2.000000 3 3.000000 4 4.000000 5 5.000000 The 2-norm of X is 7.416198 The 2-norm of row 2 of A is 11.618950 The 2-norm of column 2 of A is 11.618950 DROT_TEST DROT carries out a Givens rotation. X and Y 1 1.000000 -11.000000 2 2.000000 -8.000000 3 3.000000 -3.000000 4 4.000000 4.000000 5 5.000000 13.000000 6 6.000000 24.000000 DROT ( N, X, 1, Y, 1, 0.500000, 0.866025 ) 1 -9.026279 -6.366025 2 -5.928203 -5.732051 3 -1.098076 -4.098076 4 5.464102 -1.464102 5 13.758330 2.169873 6 23.784610 6.803848 DROT ( N, X, 1, Y, 1, 0.090536, -0.995893 ) 1 11.045361 0.000000 2 8.148217 1.267500 3 3.259287 2.716072 4 -3.621430 4.345716 5 -12.493933 6.156431 6 -23.358222 8.148217 DROTG_TEST DROTG generates a real Givens rotation ( C S ) * ( A ) = ( R ) ( -S C ) ( B ) ( 0 ) A = 0.218418 B = 0.956318 C = 0.222661 S = 0.974896 R = 0.980943 Z = 4.491123 C*A+S*B = 0.980943 -S*A+C*B = 0.000000 A = 0.829509 B = 0.561695 C = 0.828025 S = 0.560691 R = 1.001792 Z = 0.560691 C*A+S*B = 1.001792 -S*A+C*B = 0.000000 A = 0.415307 B = 0.066119 C = 0.987563 S = 0.157224 R = 0.420537 Z = 0.157224 C*A+S*B = 0.420537 -S*A+C*B = 0.000000 A = 0.257578 B = 0.109957 C = 0.919705 S = 0.392611 R = 0.280066 Z = 0.392611 C*A+S*B = 0.280066 -S*A+C*B = 0.000000 A = 0.043829 B = 0.633966 C = 0.068970 S = 0.997619 R = 0.635479 Z = 14.499053 C*A+S*B = 0.635479 -S*A+C*B = 0.000000 DSCAL_TEST DSCAL multiplies a vector by a scalar. X = 1 1.000000 2 2.000000 3 3.000000 4 4.000000 5 5.000000 6 6.000000 DSCAL ( N, 5.000000, X, 1 ) 1 5.000000 2 10.000000 3 15.000000 4 20.000000 5 25.000000 6 30.000000 DSCAL ( 3, -2.000000, X, 2 ) 1 -2.000000 2 2.000000 3 -6.000000 4 4.000000 5 -10.000000 6 6.000000 DSWAP_TEST DSWAP swaps two vectors. X and Y 1 1.000000 100.000000 2 2.000000 200.000000 3 3.000000 300.000000 4 4.000000 400.000000 5 5.000000 500.000000 6 6.000000 600.000000 DSWAP ( N, X, 1, Y, 1 ) X and Y 1 100.000000 1.000000 2 200.000000 2.000000 3 300.000000 3.000000 4 400.000000 4.000000 5 500.000000 5.000000 6 600.000000 6.000000 DSWAP ( 3, X, 2, Y, 1 ) X and Y 1 100.000000 1.000000 2 2.000000 3.000000 3 200.000000 5.000000 4 4.000000 400.000000 5 300.000000 500.000000 6 6.000000 600.000000 IDAMAX_TEST IDAMAX returns the index of maximum magnitude; The vector X: 1 1.000000 2 -3.000000 3 4.000000 4 0.000000 5 -4.000000 6 3.000000 7 -1.000000 8 -5.000000 9 2.000000 10 -2.000000 11 -6.000000 The index of maximum magnitude = 11 BLAS1_D_TEST Normal end of execution. 02-Dec-2018 09:14:55