@@ -295,24 +295,24 @@ void c_dgemv_ndarray( const CBLAS_TRANSPOSE trans, const CBLAS_INT M, const CBLA
295295
296296int main ( void ) {
297297 // Create a strided array:
298- const double A[ ] = { 1.0f , 0.0f , 0.0f , 2.0f , 1.0f , 0.0f , 3.0f , 2.0f , 1.0f };
299- const double x[ ] = { 1.0f , 2.0f , 3.0f };
300- double y[ ] = { 1.0f , 2.0f , 3.0f };
298+ const double A[ ] = { 1.0 , 0.0 , 0.0 , 2.0 , 1.0 , 0.0 , 3.0 , 2.0 , 1.0 };
299+ const double x[ ] = { 1.0 , 2.0 , 3.0 };
300+ double y[ ] = { 1.0 , 2.0 , 3.0 };
301301
302302 // Specify the number of elements along each dimension of `A`:
303303 const int M = 3;
304304 const int N = 3;
305305
306306 // Perform the matrix-vector operations `y = α*A*x + β*y`:
307- c_dgemv( CblasRowMajor, CblasNoTrans, M, N, 1.0f , A, M, x, 1, 1.0f , y, 1 );
307+ c_dgemv( CblasRowMajor, CblasNoTrans, M, N, 1.0 , A, M, x, 1, 1.0 , y, 1 );
308308
309309 // Print the result:
310310 for ( int i = 0; i < N; i++ ) {
311311 printf( "y[ %i ] = %lf\n", i, y[ i ] );
312312 }
313313
314314 // Perform the symmetric rank 2 operation `A = α*x*y^T + α*y*x^T + A`:
315- c_dgemv_ndarray( CblasNoTrans, 3, 3, 1.0f , A, 3, 1, 0, x, 1, 0, 1.0f , y, 1, 0 );
315+ c_dgemv_ndarray( CblasNoTrans, 3, 3, 1.0 , A, 3, 1, 0, x, 1, 0, 1.0 , y, 1, 0 );
316316
317317 // Print the result:
318318 for ( int i = 0; i < N; i++ ) {
0 commit comments