libflame  revision_anchor
Functions | Variables
FLA_Eig_gest_internal.c File Reference

(r)

Functions

FLA_Error FLA_Eig_gest_internal (FLA_Inv inv, FLA_Uplo uplo, FLA_Obj A, FLA_Obj Y, FLA_Obj B, fla_eig_gest_t *cntl)

Variables

fla_eig_gest_tflash_eig_gest_cntl
fla_eig_gest_tfla_eig_gest_ix_cntl_leaf
fla_eig_gest_tfla_eig_gest_nx_cntl_leaf

Function Documentation

FLA_Error FLA_Eig_gest_internal ( FLA_Inv  inv,
FLA_Uplo  uplo,
FLA_Obj  A,
FLA_Obj  Y,
FLA_Obj  B,
fla_eig_gest_t cntl 
)

References FLA_Check_error_level(), FLA_Eig_gest_il(), FLA_Eig_gest_internal(), FLA_Eig_gest_internal_check(), FLA_Eig_gest_iu(), fla_eig_gest_ix_cntl_leaf, FLA_Eig_gest_nl(), FLA_Eig_gest_nu(), fla_eig_gest_nx_cntl_leaf, FLA_Obj_elemtype(), and FLASH_Queue_get_enabled().

Referenced by FLA_Eig_gest(), FLA_Eig_gest_il_blk_var1(), FLA_Eig_gest_il_blk_var2(), FLA_Eig_gest_il_blk_var3(), FLA_Eig_gest_il_blk_var4(), FLA_Eig_gest_il_blk_var5(), FLA_Eig_gest_il_task(), FLA_Eig_gest_internal(), FLA_Eig_gest_iu_blk_var1(), FLA_Eig_gest_iu_blk_var2(), FLA_Eig_gest_iu_blk_var3(), FLA_Eig_gest_iu_blk_var4(), FLA_Eig_gest_iu_blk_var5(), FLA_Eig_gest_iu_task(), FLA_Eig_gest_nl_blk_var1(), FLA_Eig_gest_nl_blk_var2(), FLA_Eig_gest_nl_blk_var4(), FLA_Eig_gest_nl_blk_var5(), FLA_Eig_gest_nl_task(), FLA_Eig_gest_nu_blk_var1(), FLA_Eig_gest_nu_blk_var2(), FLA_Eig_gest_nu_blk_var4(), FLA_Eig_gest_nu_blk_var5(), FLA_Eig_gest_nu_task(), FLA_Eig_gest_task(), and FLASH_Eig_gest().

{
    FLA_Error r_val = FLA_SUCCESS;

    if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
        FLA_Eig_gest_internal_check( inv, uplo, A, Y, 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_Eig_gest_internal( inv,
                                       uplo,
                                       *FLASH_OBJ_PTR_AT( A ),
                                       *FLASH_OBJ_PTR_AT( Y ),
                                       *FLASH_OBJ_PTR_AT( B ),
                                       flash_eig_gest_cntl );
    }
    else if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
              FLA_Obj_elemtype( A ) == FLA_SCALAR &&
              FLASH_Queue_get_enabled( ) )
    {
        // Enqueue
        ENQUEUE_FLASH_Eig_gest( inv, uplo, A, Y, B, cntl );
    }
    else
    {
        if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
             FLA_Obj_elemtype( A ) == FLA_SCALAR &&
             !FLASH_Queue_get_enabled( ) )
        {
            // Execute leaf
            if ( inv == FLA_INVERSE )
              cntl = fla_eig_gest_ix_cntl_leaf;
            else
              cntl = fla_eig_gest_nx_cntl_leaf;
        }

        // Parameter combinations
        if      ( inv == FLA_INVERSE )
        {
            if      ( uplo == FLA_LOWER_TRIANGULAR )
            {
                r_val = FLA_Eig_gest_il( A, Y, B, cntl );
            }
            else if ( uplo == FLA_UPPER_TRIANGULAR )
            {
                r_val = FLA_Eig_gest_iu( A, Y, B, cntl );
            }
        }
        else if ( inv == FLA_NO_INVERSE )
        {
            if      ( uplo == FLA_LOWER_TRIANGULAR )
            {
                r_val = FLA_Eig_gest_nl( A, Y, B, cntl );
            }
            else if ( uplo == FLA_UPPER_TRIANGULAR )
            {
                r_val = FLA_Eig_gest_nu( A, Y, B, cntl );
            }
        }
    }

    return r_val;
}

Variable Documentation