libflame  revision_anchor
Functions
FLA_Chol_solve.c File Reference

(r)

Functions

FLA_Error FLA_Chol_solve (FLA_Uplo uplo, FLA_Obj A, FLA_Obj B, FLA_Obj X)
 

Function Documentation

◆ FLA_Chol_solve()

FLA_Error FLA_Chol_solve ( FLA_Uplo  uplo,
FLA_Obj  A,
FLA_Obj  B,
FLA_Obj  X 
)

References FLA_Check_error_level(), FLA_Chol_solve_check(), FLA_Copy_external(), FLA_Obj_is_identical(), FLA_ONE, and FLA_Trsm_external().

14 {
15  // Check parameters.
16  if ( FLA_Check_error_level() >= FLA_MIN_ERROR_CHECKING )
17  FLA_Chol_solve_check( uplo, A, B, X );
18 
19  if ( FLA_Obj_is_identical( B, X ) == FALSE )
20  FLA_Copy_external( B, X );
21 
22  if ( uplo == FLA_LOWER_TRIANGULAR )
23  {
24  FLA_Trsm_external( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE,
25  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
26  FLA_Trsm_external( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_CONJ_TRANSPOSE,
27  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
28  }
29  else // if ( uplo == FLA_UPPER_TRIANGULAR )
30  {
31  FLA_Trsm_external( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_CONJ_TRANSPOSE,
32  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
33  FLA_Trsm_external( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE,
34  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
35  }
36 
37  return FLA_SUCCESS;
38 }
FLA_Obj FLA_ONE
Definition: FLA_Init.c:18
FLA_Error FLA_Copy_external(FLA_Obj A, FLA_Obj B)
Definition: FLA_Copy_external.c:13
FLA_Error FLA_Chol_solve_check(FLA_Uplo uplo, FLA_Obj A, FLA_Obj B, FLA_Obj X)
Definition: FLA_Chol_solve_check.c:13
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
FLA_Bool FLA_Obj_is_identical(FLA_Obj A, FLA_Obj B)
Definition: FLA_Query.c:470
FLA_Error FLA_Trsm_external(FLA_Side side, FLA_Uplo uplo, FLA_Trans trans, FLA_Diag diag, FLA_Obj alpha, FLA_Obj A, FLA_Obj B)
Definition: FLA_Trsm_external.c:13