amd/blis
1
0
Fork 0
mirror of https://github.com/amd/blis.git synced 2026-05-11 17:50:00 +00:00
Code Issues Packages Projects Releases Wiki Activity

Merge branch 'master' of https://github.com/flame/blis

Browse Source
This commit is contained in:
Tyler Smith
2014-07-23 15:07:30 -05:00
parent cdb9413e14 47b243ef08
commit acff74041b
2 changed files with 18 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
frame/3/her2k/bli_her2k_front.c
View File

@@ -52,10 +52,12 @@ void bli_her2k_front( obj_t* alpha,
if ( bli_error_checking_is_enabled() )
bli_her2k_check( alpha, a, b, beta, c );
// If alpha is zero, scale by beta and return.
// If alpha is zero, scale by beta, zero the imaginary components of
// the diagonal elements, and return.
if ( bli_obj_equals( alpha, &BLIS_ZERO ) )
{
bli_scalm( beta, c );
bli_setid( &BLIS_ZERO, c );
return;
}
@@ -139,7 +141,12 @@ void bli_her2k_front( obj_t* alpha,
#endif
// Explicitly set diagonal elements' imaginary components to zero.
// The Hermitian rank-2k product was computed as A*B'+B*A', even for
// the diagonal elements. Mathematically, the imaginary components of
// diagonal elements of a Hermitian rank-2k product should always be
// zero. However, in practice, they sometimes accumulate meaningless
// non-zero values. To prevent this, we explicitly set those values
// to zero before returning.
bli_setid( &BLIS_ZERO, &c_local );
}

11
frame/3/herk/bli_herk_front.c
View File

@@ -48,10 +48,12 @@ void bli_herk_front( obj_t* alpha,
if ( bli_error_checking_is_enabled() )
bli_herk_check( alpha, a, beta, c );
// If alpha is zero, scale by beta and return.
// If alpha is zero, scale by beta, zero the imaginary components of
// the diagonal elements, and return.
if ( bli_obj_equals( alpha, &BLIS_ZERO ) )
{
bli_scalm( beta, c );
bli_setid( &BLIS_ZERO, c );
return;
}
@@ -93,7 +95,12 @@ void bli_herk_front( obj_t* alpha,
bli_herk_thrinfo_free_paths( infos, n_threads );
// Explicitly set diagonal elements' imaginary components to zero.
// The Hermitian rank-k product was computed as A*A', even for the
// diagonal elements. Mathematically, the imaginary components of
// diagonal elements of a Hermitian rank-k product should always be
// zero. However, in practice, they sometimes accumulate meaningless
// non-zero values. To prevent this, we explicitly set those values
// to zero before returning.
bli_setid( &BLIS_ZERO, &c_local );
}