libflame
revision_anchor
|
Functions | |
void | bli_smaxabsm (int m, int n, float *a, int a_rs, int a_cs, float *maxabs) |
void | bli_dmaxabsm (int m, int n, double *a, int a_rs, int a_cs, double *maxabs) |
void | bli_cmaxabsm (int m, int n, scomplex *a, int a_rs, int a_cs, float *maxabs) |
void | bli_zmaxabsm (int m, int n, dcomplex *a, int a_rs, int a_cs, double *maxabs) |
void bli_cmaxabsm | ( | int | m, |
int | n, | ||
scomplex * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
float * | maxabs | ||
) |
References bli_cmaxabsv(), bli_is_row_storage(), bli_s0(), and bli_zero_dim2().
Referenced by FLA_Max_abs_value().
{ float zero = bli_s0(); scomplex* a_begin; float maxabs_cand; float maxabs_temp; int inca, lda; int n_iter; int n_elem; int j; // Return early if possible. if ( bli_zero_dim2( m, n ) ) { *maxabs = zero; return; } // Initialize with optimal values for column-major storage. inca = a_rs; lda = a_cs; n_iter = n; n_elem = m; // An optimization: if A is row-major, then let's access the matrix by // rows instead of by columns for increased spatial locality. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( n_iter, n_elem ); bli_swap_ints( lda, inca ); } // Initialize the maximum absolute value candidate to the first element. bli_csabsval2( a, &maxabs_cand ); for ( j = 0; j < n_iter; j++ ) { a_begin = a + j*lda; bli_cmaxabsv( n_elem, a_begin, inca, &maxabs_temp ); if ( maxabs_temp > maxabs_cand ) maxabs_cand = maxabs_temp; } *maxabs = maxabs_cand; }
void bli_dmaxabsm | ( | int | m, |
int | n, | ||
double * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
double * | maxabs | ||
) |
References bli_d0(), bli_dmaxabsv(), bli_is_row_storage(), and bli_zero_dim2().
Referenced by FLA_Max_abs_value().
{ double zero = bli_d0(); double* a_begin; double maxabs_cand; double maxabs_temp; int inca, lda; int n_iter; int n_elem; int j; // Return early if possible. if ( bli_zero_dim2( m, n ) ) { *maxabs = zero; return; } // Initialize with optimal values for column-major storage. inca = a_rs; lda = a_cs; n_iter = n; n_elem = m; // An optimization: if A is row-major, then let's access the matrix by // rows instead of by columns for increased spatial locality. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( n_iter, n_elem ); bli_swap_ints( lda, inca ); } // Initialize the maximum absolute value candidate to the first element. bli_dabsval2( a, &maxabs_cand ); for ( j = 0; j < n_iter; j++ ) { a_begin = a + j*lda; bli_dmaxabsv( n_elem, a_begin, inca, &maxabs_temp ); if ( maxabs_temp > maxabs_cand ) maxabs_cand = maxabs_temp; } *maxabs = maxabs_cand; }
void bli_smaxabsm | ( | int | m, |
int | n, | ||
float * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
float * | maxabs | ||
) |
References bli_is_row_storage(), bli_s0(), bli_smaxabsv(), and bli_zero_dim2().
Referenced by FLA_Max_abs_value().
{ float zero = bli_s0(); float* a_begin; float maxabs_cand; float maxabs_temp; int inca, lda; int n_iter; int n_elem; int j; // Return early if possible. if ( bli_zero_dim2( m, n ) ) { *maxabs = zero; return; } // Initialize with optimal values for column-major storage. inca = a_rs; lda = a_cs; n_iter = n; n_elem = m; // An optimization: if A is row-major, then let's access the matrix by // rows instead of by columns for increased spatial locality. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( n_iter, n_elem ); bli_swap_ints( lda, inca ); } // Initialize the maximum absolute value candidate to the first element. bli_sabsval2( a, &maxabs_cand ); for ( j = 0; j < n_iter; j++ ) { a_begin = a + j*lda; bli_smaxabsv( n_elem, a_begin, inca, &maxabs_temp ); if ( maxabs_temp > maxabs_cand ) maxabs_cand = maxabs_temp; } *maxabs = maxabs_cand; }
void bli_zmaxabsm | ( | int | m, |
int | n, | ||
dcomplex * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
double * | maxabs | ||
) |
References bli_d0(), bli_is_row_storage(), bli_zero_dim2(), and bli_zmaxabsv().
Referenced by FLA_Max_abs_value().
{ double zero = bli_d0(); dcomplex* a_begin; double maxabs_cand; double maxabs_temp; int inca, lda; int n_iter; int n_elem; int j; // Return early if possible. if ( bli_zero_dim2( m, n ) ) { *maxabs = zero; return; } // Initialize with optimal values for column-major storage. inca = a_rs; lda = a_cs; n_iter = n; n_elem = m; // An optimization: if A is row-major, then let's access the matrix by // rows instead of by columns for increased spatial locality. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( n_iter, n_elem ); bli_swap_ints( lda, inca ); } // Initialize the maximum absolute value candidate to the first element. bli_zdabsval2( a, &maxabs_cand ); for ( j = 0; j < n_iter; j++ ) { a_begin = a + j*lda; bli_zmaxabsv( n_elem, a_begin, inca, &maxabs_temp ); if ( maxabs_temp > maxabs_cand ) maxabs_cand = maxabs_temp; } *maxabs = maxabs_cand; }