|
libflame
revision_anchor
|
Go to the source code of this file.
| FLA_Error FLA_Apply_Q_UT_create_workspace | ( | FLA_Obj | T, |
| FLA_Obj | B, | ||
| FLA_Obj * | W | ||
| ) |
References FLA_Obj_create(), FLA_Obj_datatype(), FLA_Obj_length(), and FLA_Obj_width().
Referenced by FLA_LQ_UT_form_Q(), FLA_LQ_UT_solve(), FLA_QR_UT_form_Q(), and FLA_QR_UT_solve().
{
FLA_Datatype datatype;
dim_t m_W, n_W;
datatype = FLA_Obj_datatype( T );
m_W = FLA_Obj_length( T );
n_W = FLA_Obj_width( B );
FLA_Obj_create( datatype, m_W, n_W, 0, 0, W );
return FLA_SUCCESS;
}
| FLA_Error FLA_Apply_Q_UT_internal | ( | FLA_Side | side, |
| FLA_Trans | trans, | ||
| FLA_Direct | direct, | ||
| FLA_Store | storev, | ||
| FLA_Obj | A, | ||
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_internal(), FLA_Apply_Q_UT_internal_check(), FLA_Apply_Q_UT_lhbc(), FLA_Apply_Q_UT_lhbr(), FLA_Apply_Q_UT_lhfc(), FLA_Apply_Q_UT_lhfr(), FLA_Apply_Q_UT_lnbc(), FLA_Apply_Q_UT_lnbr(), FLA_Apply_Q_UT_lnfc(), FLA_Apply_Q_UT_lnfr(), FLA_Apply_Q_UT_rhbc(), FLA_Apply_Q_UT_rhbr(), FLA_Apply_Q_UT_rhfc(), FLA_Apply_Q_UT_rhfr(), FLA_Apply_Q_UT_rnbc(), FLA_Apply_Q_UT_rnbr(), FLA_Apply_Q_UT_rnfc(), FLA_Apply_Q_UT_rnfr(), fla_apqut_cntl_leaf, FLA_Check_error_level(), FLA_Obj_elemtype(), and FLASH_Queue_get_enabled().
Referenced by FLA_Apply_Q_UT(), FLA_Apply_Q_UT_inc_lhfc_blk_var1(), FLA_Apply_Q_UT_inc_lnfc_blk_var1(), FLA_Apply_Q_UT_internal(), FLA_Apply_Q_UT_lhbc_blk_var2(), FLA_Apply_Q_UT_lhbc_task(), FLA_Apply_Q_UT_lhbr_blk_var2(), FLA_Apply_Q_UT_lhbr_task(), FLA_Apply_Q_UT_lhfc_blk_var2(), FLA_Apply_Q_UT_lhfc_task(), FLA_Apply_Q_UT_lhfr_blk_var2(), FLA_Apply_Q_UT_lhfr_task(), FLA_Apply_Q_UT_lnbc_blk_var2(), FLA_Apply_Q_UT_lnbc_task(), FLA_Apply_Q_UT_lnbr_blk_var2(), FLA_Apply_Q_UT_lnbr_task(), FLA_Apply_Q_UT_lnfc_blk_var2(), FLA_Apply_Q_UT_lnfc_task(), FLA_Apply_Q_UT_lnfr_blk_var2(), FLA_Apply_Q_UT_lnfr_task(), FLA_Apply_Q_UT_rhbc_blk_var2(), FLA_Apply_Q_UT_rhbc_task(), FLA_Apply_Q_UT_rhbr_blk_var2(), FLA_Apply_Q_UT_rhbr_task(), FLA_Apply_Q_UT_rhfc_blk_var2(), FLA_Apply_Q_UT_rhfc_task(), FLA_Apply_Q_UT_rhfr_blk_var2(), FLA_Apply_Q_UT_rhfr_task(), FLA_Apply_Q_UT_rnbc_blk_var2(), FLA_Apply_Q_UT_rnbc_task(), FLA_Apply_Q_UT_rnbr_blk_var2(), FLA_Apply_Q_UT_rnbr_task(), FLA_Apply_Q_UT_rnfc_blk_var2(), FLA_Apply_Q_UT_rnfc_task(), FLA_Apply_Q_UT_rnfr_blk_var2(), FLA_Apply_Q_UT_rnfr_task(), FLA_Apply_Q_UT_task(), FLA_LQ_UT_blk_var1(), FLA_LQ_UT_blk_var2(), FLA_LQ_UT_blk_var3(), FLA_QR_UT_blk_var1(), FLA_QR_UT_blk_var2(), FLA_QR_UT_blk_var3(), FLA_QR_UT_inc_blk_var1(), FLA_QR_UT_inc_blk_var2(), and FLASH_Apply_Q_UT().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
FLA_Apply_Q_UT_internal_check( side, trans, direct, storev, A, T, W, B, cntl );
if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
FLA_Obj_elemtype( A ) == FLA_MATRIX &&
FLA_Cntl_variant( cntl ) == FLA_SUBPROBLEM )
{
// Recurse
r_val = FLA_Apply_Q_UT_internal( side,
trans,
direct,
storev,
*FLASH_OBJ_PTR_AT( A ),
*FLASH_OBJ_PTR_AT( T ),
*FLASH_OBJ_PTR_AT( W ),
*FLASH_OBJ_PTR_AT( B ),
flash_apqut_cntl );
}
else if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
FLA_Obj_elemtype( A ) == FLA_SCALAR &&
FLASH_Queue_get_enabled( ) )
{
// Enqueue
ENQUEUE_FLASH_Apply_Q_UT( side, trans, direct, storev, A, T, W, B, cntl );
}
else
{
if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
FLA_Obj_elemtype( A ) == FLA_SCALAR &&
!FLASH_Queue_get_enabled( ) )
{
// Execute leaf.
cntl = fla_apqut_cntl_leaf;
}
if ( side == FLA_LEFT )
{
if ( trans == FLA_NO_TRANSPOSE )
{
if ( direct == FLA_FORWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_lnfc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_lnfr( A, T, W, B, cntl );
}
else if ( direct == FLA_BACKWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_lnbc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_lnbr( A, T, W, B, cntl );
}
}
else if ( trans == FLA_TRANSPOSE || trans == FLA_CONJ_TRANSPOSE )
{
if ( direct == FLA_FORWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_lhfc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_lhfr( A, T, W, B, cntl );
}
else if ( direct == FLA_BACKWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_lhbc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_lhbr( A, T, W, B, cntl );
}
}
}
else if ( side == FLA_RIGHT )
{
if ( trans == FLA_NO_TRANSPOSE )
{
if ( direct == FLA_FORWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_rnfc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_rnfr( A, T, W, B, cntl );
}
else if ( direct == FLA_BACKWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_rnbc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_rnbr( A, T, W, B, cntl );
}
}
else if ( trans == FLA_TRANSPOSE || trans == FLA_CONJ_TRANSPOSE )
{
if ( direct == FLA_FORWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_rhfc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_rhfr( A, T, W, B, cntl );
}
else if ( direct == FLA_BACKWARD )
{
if ( storev == FLA_COLUMNWISE )
r_val = FLA_Apply_Q_UT_rhbc( A, T, W, B, cntl );
else if ( storev == FLA_ROWWISE )
r_val = FLA_Apply_Q_UT_rhbr( A, T, W, B, cntl );
}
}
}
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lhbc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lhbc_blk_var1(), FLA_Apply_Q_UT_lhbc_blk_var2(), and FLA_Apply_Q_UT_lhbc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lhbc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lhbc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lhbc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lhbr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lhbr_blk_var1(), FLA_Apply_Q_UT_lhbr_blk_var2(), and FLA_Apply_Q_UT_lhbr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lhbr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lhbr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lhbr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lhfc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lhfc_blk_var1(), FLA_Apply_Q_UT_lhfc_blk_var2(), and FLA_Apply_Q_UT_lhfc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lhfc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lhfc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lhfc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lhfr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lhfr_blk_var1(), FLA_Apply_Q_UT_lhfr_blk_var2(), and FLA_Apply_Q_UT_lhfr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lhfr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lhfr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lhfr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lnbc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lnbc_blk_var1(), FLA_Apply_Q_UT_lnbc_blk_var2(), and FLA_Apply_Q_UT_lnbc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lnbc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lnbc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lnbc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lnbr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lnbr_blk_var1(), FLA_Apply_Q_UT_lnbr_blk_var2(), and FLA_Apply_Q_UT_lnbr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lnbr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lnbr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lnbr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lnfc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lnfc_blk_var1(), FLA_Apply_Q_UT_lnfc_blk_var2(), and FLA_Apply_Q_UT_lnfc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lnfc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lnfc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lnfc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_lnfr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_lnfr_blk_var1(), FLA_Apply_Q_UT_lnfr_blk_var2(), and FLA_Apply_Q_UT_lnfr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_lnfr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_lnfr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_lnfr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rhbc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rhbc_blk_var1(), FLA_Apply_Q_UT_rhbc_blk_var2(), and FLA_Apply_Q_UT_rhbc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rhbc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rhbc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rhbc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rhbr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rhbr_blk_var1(), FLA_Apply_Q_UT_rhbr_blk_var2(), and FLA_Apply_Q_UT_rhbr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rhbr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rhbr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rhbr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rhfc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rhfc_blk_var1(), FLA_Apply_Q_UT_rhfc_blk_var2(), and FLA_Apply_Q_UT_rhfc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rhfc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rhfc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rhfc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rhfr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rhfr_blk_var1(), FLA_Apply_Q_UT_rhfr_blk_var2(), and FLA_Apply_Q_UT_rhfr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rhfr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rhfr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rhfr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rnbc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rnbc_blk_var1(), FLA_Apply_Q_UT_rnbc_blk_var2(), and FLA_Apply_Q_UT_rnbc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rnbc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rnbc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rnbc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rnbr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rnbr_blk_var1(), FLA_Apply_Q_UT_rnbr_blk_var2(), and FLA_Apply_Q_UT_rnbr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rnbr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rnbr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rnbr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rnfc | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rnfc_blk_var1(), FLA_Apply_Q_UT_rnfc_blk_var2(), and FLA_Apply_Q_UT_rnfc_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rnfc_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rnfc_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rnfc_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLA_Apply_Q_UT_rnfr | ( | FLA_Obj | A, |
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B, | ||
| fla_apqut_t * | cntl | ||
| ) |
References FLA_Apply_Q_UT_rnfr_blk_var1(), FLA_Apply_Q_UT_rnfr_blk_var2(), and FLA_Apply_Q_UT_rnfr_blk_var3().
Referenced by FLA_Apply_Q_UT_internal().
{
FLA_Error r_val = FLA_SUCCESS;
if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
{
r_val = FLA_Apply_Q_UT_rnfr_blk_var1( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
{
r_val = FLA_Apply_Q_UT_rnfr_blk_var2( A, T, W, B, cntl );
}
else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
{
r_val = FLA_Apply_Q_UT_rnfr_blk_var3( A, T, W, B, cntl );
}
else
{
FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
}
return r_val;
}
| FLA_Error FLASH_Apply_Q_UT | ( | FLA_Side | side, |
| FLA_Trans | trans, | ||
| FLA_Direct | direct, | ||
| FLA_Store | storev, | ||
| FLA_Obj | A, | ||
| FLA_Obj | T, | ||
| FLA_Obj | W, | ||
| FLA_Obj | B | ||
| ) |
References FLA_Abort(), FLA_Apply_Q_UT_check(), FLA_Apply_Q_UT_internal(), FLA_Blocksize_set(), FLA_Check_error_level(), FLA_Print_message(), FLASH_Obj_scalar_length_tl(), FLASH_Obj_scalar_width_tl(), FLASH_Queue_begin(), and FLASH_Queue_end().
{
FLA_Error r_val;
dim_t b_alg;
// Check parameters.
if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
FLA_Apply_Q_UT_check( side, trans, direct, storev, A, T, W, B );
// Inspect the length of TTL to get the blocksize used by the QR/LQ
// factorization, which will be our inner blocksize for Apply_Q_UT.
b_alg = FLASH_Obj_scalar_length_tl( T );
// The traditional (non-incremental) Apply_Q_UT algorithm-by-blocks
// requires that the algorithmic blocksize be equal to the storage
// blocksize.
if ( b_alg != FLASH_Obj_scalar_width_tl( T ) )
{
FLA_Print_message( "FLASH_Apply_Q_UT() requires that b_alg == b_store",
__FILE__, __LINE__ );
FLA_Abort();
}
// Adjust the blocksize of the control tree node for the flat subproblem.
if ( FLA_Cntl_blocksize( fla_apqut_cntl_leaf ) != NULL )
FLA_Blocksize_set( FLA_Cntl_blocksize( fla_apqut_cntl_leaf ),
b_alg, b_alg, b_alg, b_alg );
// Begin a parallel region.
FLASH_Queue_begin();
// Invoke FLA_Apply_Q_UT_internal() with the standard control tree.
r_val = FLA_Apply_Q_UT_internal( side, trans, direct, storev, A, T, W, B,
flash_apqut_cntl_blas );
// End the parallel region.
FLASH_Queue_end();
return r_val;
}
| FLA_Error FLASH_Apply_Q_UT_create_workspace | ( | FLA_Obj | TW, |
| FLA_Obj | B, | ||
| FLA_Obj * | W | ||
| ) |
References FLA_Abort(), FLA_Obj_datatype(), FLA_Print_message(), FLASH_Obj_create_ext(), FLASH_Obj_depth(), FLASH_Obj_scalar_length_tl(), FLASH_Obj_scalar_width(), and FLASH_Obj_scalar_width_tl().
Referenced by FLASH_LQ_UT_solve(), and FLASH_QR_UT_solve().
{
FLA_Datatype datatype;
dim_t depth;
dim_t b_alg;
dim_t b_flash;
dim_t m, n;
// Query the depth.
depth = FLASH_Obj_depth( TW );
// *** The current Apply_Q_UT algorithm implemented assumes that
// the matrix has a hierarchical depth of 1. We check for that here
// because we anticipate that we'll use a more general algorithm in the
// future, and we don't want to forget to remove the constraint. ***
if ( depth != 1 )
{
FLA_Print_message( "FLASH_Apply_Q_UT() currently only supports matrices of depth 1",
__FILE__, __LINE__ );
FLA_Abort();
}
// Query the datatype of matrix TW.
datatype = FLA_Obj_datatype( TW );
// Inspect the dimensions of a the top-left element of TW to get the
// algorithmic/storage blocksize we'll use throughout the Apply_Q_UT
// algorithm.
b_alg = FLASH_Obj_scalar_length_tl( TW );
b_flash = FLASH_Obj_scalar_width_tl( TW );
// The traditional (non-incremental) Apply_Q_UT algorithm-by-blocks
// requires that the algorithmic blocksize be equal to the storage
// blocksize.
if ( b_alg != b_flash )
{
FLA_Print_message( "FLASH_Apply_Q_UT() requires that b_alg == b_store",
__FILE__, __LINE__ );
FLA_Abort();
}
// The scalar length of W should be the algorithmc/storage blocksize
// encoded in TW.
m = b_alg;
// Query the scalar (not element) width of the right-hand side
// matrix B.
n = FLASH_Obj_scalar_width( B );
// Create hierarchical matrix W.
FLASH_Obj_create_ext( datatype, m, n,
depth, &b_alg, &b_flash,
W );
return FLA_SUCCESS;
}
1.7.6.1