@@ -295,24 +295,24 @@ void c_dgemv_ndarray( const CBLAS_TRANSPOSE trans, const CBLAS_INT M, const CBLA
295
295
296
296
int main ( void ) {
297
297
// 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 };
301
301
302
302
// Specify the number of elements along each dimension of `A`:
303
303
const int M = 3;
304
304
const int N = 3;
305
305
306
306
// 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 );
308
308
309
309
// Print the result:
310
310
for ( int i = 0; i < N; i++ ) {
311
311
printf( "y[ %i ] = %lf\n", i, y[ i ] );
312
312
}
313
313
314
314
// 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 );
316
316
317
317
// Print the result:
318
318
for ( int i = 0; i < N; i++ ) {
0 commit comments