FLA_LU_piv.h File Reference

(r)

Go to the source code of this file.

Functions
FLA_Error	FLA_LU_piv_internal (FLA_Obj A, FLA_Obj p, fla_lu_t *cntl)
FLA_Error	FLA_LU_piv_solve (FLA_Obj A, FLA_Obj p, FLA_Obj B, FLA_Obj X)
FLA_Error	FLASH_LU_piv_solve (FLA_Obj A, FLA_Obj p, FLA_Obj B, FLA_Obj X)

---

Function Documentation

FLA_Error FLA_LU_piv_internal	(	FLA_Obj	A,
		FLA_Obj	p,
		fla_lu_t *	cntl
	)

References FLA_Check_error_level(), FLA_Check_null_pointer(), FLA_LU_piv_blk_ext(), FLA_LU_piv_blk_var3(), FLA_LU_piv_blk_var4(), FLA_LU_piv_blk_var5(), FLA_LU_piv_macro_task(), FLA_LU_piv_opt_var3(), FLA_LU_piv_opt_var4(), FLA_LU_piv_opt_var5(), FLA_LU_piv_unb_ext(), FLA_LU_piv_unb_var3(), FLA_LU_piv_unb_var4(), FLA_LU_piv_unb_var5(), and FLASH_Queue_get_enabled().

Referenced by FLA_LU_piv(), FLA_LU_piv_blk_var3(), FLA_LU_piv_blk_var4(), FLA_LU_piv_blk_var5(), FLA_LU_piv_task(), and FLASH_LU_piv().

{
    FLA_Error r_val = FLA_SUCCESS;

    if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
    {
        FLA_Error e_val = FLA_Check_null_pointer( ( void* ) cntl );
        FLA_Check_error_code( e_val );
    }

    if      ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
              FLA_Cntl_variant( cntl ) == FLA_SUBPROBLEM &&
              FLASH_Queue_get_enabled( ) )
    {
        // Enqueue
        ENQUEUE_FLASH_LU_piv_macro( A, *FLASH_OBJ_PTR_AT( p ), cntl );
    }
    else
    {
/*
        if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
             FLA_Obj_elemtype( A ) == FLA_SCALAR &&
             !FLASH_Queue_get_enabled( ) )
        {
            // Execute leaf
            cntl = fla_lu_piv_cntl_leaf;
        }
*/
        // Choose implementation.
        if      ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_EXTERN ||
                  FLA_Cntl_variant( cntl ) == FLA_SUBPROBLEM )
        {
            if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER )
            {
                // When performing an unblocked operation on a hierarhical
                // object, we assume it's because we need to perform the
                // the operation on a macro block.
                r_val = FLA_LU_piv_macro_task( A, *FLASH_OBJ_PTR_AT( p ), cntl );
            }
            else
            {
                r_val = FLA_LU_piv_blk_ext( A, p );
            }
        }
        else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_EXTERN )
        {
            if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER )
            {
                // When performing an unblocked operation on a hierarhical
                // object, we assume it's because we need to perform the
                // the operation on a macro block.
                r_val = FLA_LU_piv_macro_task( A, *FLASH_OBJ_PTR_AT( p ), cntl );
            }
            else
            {
                r_val = FLA_LU_piv_unb_ext( A, p );
            }
        }
#ifdef FLA_ENABLE_NON_CRITICAL_CODE
        else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT3 )
        {
            r_val = FLA_LU_piv_unb_var3( A, p );
        }
        else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT4 )
        {
            r_val = FLA_LU_piv_unb_var4( A, p );
        }
#endif
        else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT5 )
        {
            r_val = FLA_LU_piv_unb_var5( A, p );
        }
#ifdef FLA_ENABLE_NON_CRITICAL_CODE
        else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT3 )
        {
            r_val = FLA_LU_piv_opt_var3( A, p );
        }
        else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT4 )
        {
            r_val = FLA_LU_piv_opt_var4( A, p );
        }
#endif
        else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT5 )
        {
            r_val = FLA_LU_piv_opt_var5( A, p );
        }
#ifdef FLA_ENABLE_NON_CRITICAL_CODE
        else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
        {
            r_val = FLA_LU_piv_blk_var3( A, p, cntl );
        }
        else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT4 )
        {
            r_val = FLA_LU_piv_blk_var4( A, p, cntl );
        }
#endif
        else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT5 )
        {
            r_val = FLA_LU_piv_blk_var5( A, p, cntl );
        }
        else
        {
            FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
        }
    }

    return r_val;
}

FLA_Error FLA_LU_piv_solve	(	FLA_Obj	A,
		FLA_Obj	p,
		FLA_Obj	B,
		FLA_Obj	X
	)

References FLA_Apply_pivots(), FLA_Check_error_level(), FLA_Copy_external(), FLA_LU_piv_solve_check(), FLA_ONE, and FLA_Trsm_external().

{
  // Check parameters.
  if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
    FLA_LU_piv_solve_check( A, p, B, X );

  FLA_Copy_external( B, X );

  FLA_Apply_pivots( FLA_LEFT, FLA_NO_TRANSPOSE, p, X );

  FLA_Trsm_external( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE,
                     FLA_UNIT_DIAG, FLA_ONE, A, X );
  FLA_Trsm_external( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE,
                     FLA_NONUNIT_DIAG, FLA_ONE, A, X );

  return FLA_SUCCESS;
}

FLA_Error FLASH_LU_piv_solve	(	FLA_Obj	A,
		FLA_Obj	p,
		FLA_Obj	B,
		FLA_Obj	X
	)

References FLA_Check_error_level(), FLA_LU_piv_solve_check(), FLA_ONE, FLASH_Apply_pivots(), FLASH_Copy(), and FLASH_Trsm().

{
  // Check parameters.
  if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
    FLA_LU_piv_solve_check( A, p, B, X );

  FLASH_Copy( B, X );

  FLASH_Apply_pivots( FLA_LEFT, FLA_NO_TRANSPOSE, p, X );

  FLASH_Trsm( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE,
              FLA_UNIT_DIAG, FLA_ONE, A, X );
  FLASH_Trsm( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE,
              FLA_NONUNIT_DIAG, FLA_ONE, A, X );

  return FLA_SUCCESS;
}