libflame revision_anchor
|
Go to the source code of this file.
References FLA_Obj_create(), FLA_Obj_datatype(), FLA_Obj_length(), and FLA_Obj_width().
Referenced by FLA_LQ_UT_solve(), 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_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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }
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; }