libflame revision_anchor
|
Functions | |
void | bli_sconjmr (char uplo, int m, int n, float *a, int a_rs, int a_cs) |
void | bli_dconjmr (char uplo, int m, int n, double *a, int a_rs, int a_cs) |
void | bli_cconjmr (char uplo, int m, int n, scomplex *a, int a_rs, int a_cs) |
void | bli_zconjmr (char uplo, int m, int n, dcomplex *a, int a_rs, int a_cs) |
void bli_cconjmr | ( | char | uplo, |
int | m, | ||
int | n, | ||
scomplex * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
References bli_is_row_storage(), bli_is_upper(), bli_sm1(), bli_sscal(), and bli_zero_dim2().
Referenced by bli_chemm(), bli_ctrmm(), bli_ctrsm(), and FLA_Conjugate_r().
{ float m1 = bli_sm1(); float* a_conj; int lda, inca; int n_iter; int n_elem_max; int n_elem; int j; // Return early if possible. if ( bli_zero_dim2( m, n ) ) return; // We initialize for column-major. n_iter = n; n_elem_max = m; lda = a_cs; inca = a_rs; // An optimization: if A is row-major, then let's access the matrix // by rows instead of by columns to increase spatial locality. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( n_iter, n_elem_max ); bli_swap_ints( lda, inca ); bli_toggle_uplo( uplo ); } if ( bli_is_upper( uplo ) ) { for ( j = 0; j < n_iter; ++j ) { n_elem = bli_min( j + 1, n_elem_max ); a_conj = ( float* )( a + j*lda ) + 1; bli_sscal( n_elem, &m1, a_conj, 2*inca ); } } else // if ( bli_is_lower( uplo ) ) { for ( j = 0; j < n_iter; ++j ) { n_elem = bli_max( 0, n_elem_max - j ); a_conj = ( float* )( a + j*lda + j*inca ) + 1; if ( n_elem <= 0 ) break; bli_sscal( n_elem, &m1, a_conj, 2*inca ); } } }
void bli_dconjmr | ( | char | uplo, |
int | m, | ||
int | n, | ||
double * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
{
return;
}
void bli_sconjmr | ( | char | uplo, |
int | m, | ||
int | n, | ||
float * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
{
return;
}
void bli_zconjmr | ( | char | uplo, |
int | m, | ||
int | n, | ||
dcomplex * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
References bli_dm1(), bli_dscal(), bli_is_row_storage(), bli_is_upper(), and bli_zero_dim2().
Referenced by bli_zhemm(), bli_ztrmm(), bli_ztrsm(), and FLA_Conjugate_r().
{ double m1 = bli_dm1(); double* a_conj; int lda, inca; int n_iter; int n_elem_max; int n_elem; int j; // Return early if possible. if ( bli_zero_dim2( m, n ) ) return; // We initialize for column-major. n_iter = n; n_elem_max = m; lda = a_cs; inca = a_rs; // An optimization: if A is row-major, then let's access the matrix // by rows instead of by columns to increase spatial locality. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( n_iter, n_elem_max ); bli_swap_ints( lda, inca ); bli_toggle_uplo( uplo ); } if ( bli_is_upper( uplo ) ) { for ( j = 0; j < n_iter; ++j ) { n_elem = bli_min( j + 1, n_elem_max ); a_conj = ( double* )( a + j*lda ) + 1; bli_dscal( n_elem, &m1, a_conj, 2*inca ); } } else // if ( bli_is_lower( uplo ) ) { for ( j = 0; j < n_iter; ++j ) { n_elem = bli_max( 0, n_elem_max - j ); a_conj = ( double* )( a + j*lda + j*inca ) + 1; if ( n_elem <= 0 ) break; bli_dscal( n_elem, &m1, a_conj, 2*inca ); } } }