libflame revision_anchor
Functions
FLA_Apply_QUD_UT.h File Reference

(r)

Go to the source code of this file.

Functions

FLA_Error FLA_Apply_QUD_UT (FLA_Side side, FLA_Trans trans, FLA_Direct direct, FLA_Store storev, FLA_Obj T, FLA_Obj W, FLA_Obj R, FLA_Obj U, FLA_Obj C, FLA_Obj V, FLA_Obj D)
FLA_Error FLA_Apply_QUD_UT_internal (FLA_Side side, FLA_Trans trans, FLA_Direct direct, FLA_Store storev, FLA_Obj T, FLA_Obj W, FLA_Obj R, FLA_Obj U, FLA_Obj C, FLA_Obj V, FLA_Obj D, fla_apqudut_t *cntl)
FLA_Error FLA_Apply_QUD_UT_lhfc (FLA_Obj T, FLA_Obj W, FLA_Obj R, FLA_Obj U, FLA_Obj C, FLA_Obj V, FLA_Obj D, fla_apqudut_t *cntl)
FLA_Error FLA_Apply_QUD_UT_create_workspace (FLA_Obj T, FLA_Obj R, FLA_Obj *W)

Function Documentation

FLA_Error FLA_Apply_QUD_UT ( FLA_Side  side,
FLA_Trans  trans,
FLA_Direct  direct,
FLA_Store  storev,
FLA_Obj  T,
FLA_Obj  W,
FLA_Obj  R,
FLA_Obj  U,
FLA_Obj  C,
FLA_Obj  V,
FLA_Obj  D 
)

References FLA_Apply_QUD_UT_check(), FLA_Apply_QUD_UT_internal(), and FLA_Check_error_level().

Referenced by FLA_UDdate_UT_update_rhs().

{
  FLA_Error r_val;

  // Check parameters.
  if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
    FLA_Apply_QUD_UT_check( side, trans, direct, storev, T, W, R, U, C, V, D );

  // Invoke _internal() back-end with the standard control tree.
  r_val = FLA_Apply_QUD_UT_internal( side, trans, direct, storev,
                                     T, W, R, U, C, V, D, fla_apqudut_cntl_leaf );

  return r_val;
}
FLA_Error FLA_Apply_QUD_UT_create_workspace ( FLA_Obj  T,
FLA_Obj  R,
FLA_Obj W 
)

References FLA_Obj_create(), FLA_Obj_datatype(), FLA_Obj_length(), and FLA_Obj_width().

Referenced by FLA_UDdate_UT_update_rhs().

{
    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( R );

    FLA_Obj_create( datatype, m_W, n_W, 0, 0, W );

    return FLA_SUCCESS;
}
FLA_Error FLA_Apply_QUD_UT_internal ( FLA_Side  side,
FLA_Trans  trans,
FLA_Direct  direct,
FLA_Store  storev,
FLA_Obj  T,
FLA_Obj  W,
FLA_Obj  R,
FLA_Obj  U,
FLA_Obj  C,
FLA_Obj  V,
FLA_Obj  D,
fla_apqudut_t cntl 
)

References FLA_Apply_QUD_UT_internal(), FLA_Apply_QUD_UT_internal_check(), FLA_Apply_QUD_UT_lhfc(), fla_apqudut_cntl_leaf, FLA_Check_error_level(), FLA_Obj_elemtype(), and FLASH_Queue_get_enabled().

Referenced by FLA_Apply_QUD_UT(), FLA_Apply_QUD_UT_inc_lhfc_blk_var1(), FLA_Apply_QUD_UT_internal(), FLA_Apply_QUD_UT_lhfc_blk_var2(), FLA_Apply_QUD_UT_lhfc_blk_var3(), FLA_Apply_QUD_UT_lhfc_task(), FLA_Apply_QUD_UT_task(), FLA_UDdate_UT_blk_var1(), and FLA_UDdate_UT_inc_blk_var1().

{
    FLA_Error r_val = FLA_SUCCESS;

    if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
        FLA_Apply_QUD_UT_internal_check( side, trans, direct, storev,
                                         T, W, R, U, C, V, D, cntl );

    if      ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
              FLA_Obj_elemtype( R ) == FLA_MATRIX &&
              FLA_Cntl_variant( cntl ) == FLA_SUBPROBLEM )
    {
        // Recurse
        r_val = FLA_Apply_QUD_UT_internal( side,
                                           trans,
                                           direct,
                                           storev,
                                           *FLASH_OBJ_PTR_AT( T ),
                                           *FLASH_OBJ_PTR_AT( W ),
                                           *FLASH_OBJ_PTR_AT( R ),
                                           *FLASH_OBJ_PTR_AT( U ),
                                           *FLASH_OBJ_PTR_AT( C ),
                                           *FLASH_OBJ_PTR_AT( V ),
                                           *FLASH_OBJ_PTR_AT( D ),
                                           flash_apqudut_cntl );
    }
    else if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
              FLA_Obj_elemtype( R ) == FLA_SCALAR &&
              FLASH_Queue_get_enabled( ) )
    {
        // Enqueue
        ENQUEUE_FLASH_Apply_QUD_UT( side, trans, direct, storev,
                                    T, W, R, U, C, V, D, cntl );
    }
    else
    {
        if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
             FLA_Obj_elemtype( R ) == FLA_SCALAR &&
             !FLASH_Queue_get_enabled( ) )
        {
            // Execute leaf.
            cntl = fla_apqudut_cntl_leaf;
        }

        if      ( side == FLA_LEFT )
        {
            if      ( trans == FLA_NO_TRANSPOSE )
            {
                if      ( direct == FLA_FORWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
                else if ( direct == FLA_BACKWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
            }
            else if ( trans == FLA_TRANSPOSE || trans == FLA_CONJ_TRANSPOSE )
            {
                if      ( direct == FLA_FORWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        r_val = FLA_Apply_QUD_UT_lhfc( T, W, R, U, C, V, D, cntl );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
                else if ( direct == FLA_BACKWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
            }
        }
        else if ( side == FLA_RIGHT )
        {
            if      ( trans == FLA_NO_TRANSPOSE )
            {
                if      ( direct == FLA_FORWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
                else if ( direct == FLA_BACKWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
            }
            else if ( trans == FLA_TRANSPOSE || trans == FLA_CONJ_TRANSPOSE )
            {
                if      ( direct == FLA_FORWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
                else if ( direct == FLA_BACKWARD )
                {
                    if      ( storev == FLA_COLUMNWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                    else if ( storev == FLA_ROWWISE )
                        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
                }
            }
        }
    }

    return r_val;
}
FLA_Error FLA_Apply_QUD_UT_lhfc ( FLA_Obj  T,
FLA_Obj  W,
FLA_Obj  R,
FLA_Obj  U,
FLA_Obj  C,
FLA_Obj  V,
FLA_Obj  D,
fla_apqudut_t cntl 
)

References FLA_Apply_QUD_UT_lhfc_blk_var1(), FLA_Apply_QUD_UT_lhfc_blk_var2(), and FLA_Apply_QUD_UT_lhfc_blk_var3().

Referenced by FLA_Apply_QUD_UT_internal().

{
    FLA_Error r_val = FLA_SUCCESS;
    
    if      ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
    {
        r_val = FLA_Apply_QUD_UT_lhfc_blk_var1( T, W, R, U, C, V, D, cntl );
    }
    else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
    {
        r_val = FLA_Apply_QUD_UT_lhfc_blk_var2( T, W, R, U, C, V, D, cntl );
    }
    else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
    {
        r_val = FLA_Apply_QUD_UT_lhfc_blk_var3( T, W, R, U, C, V, D, cntl );
    }
    else
    {
        FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
    }

    return r_val;
}