libflame  revision_anchor
Functions
FLA_Chol.h File Reference

(r)

Go to the source code of this file.

Functions

FLA_Error FLA_Chol_internal (FLA_Uplo uplo, FLA_Obj A, fla_chol_t *cntl)
 
FLA_Error FLA_Chol_l (FLA_Obj A, fla_chol_t *cntl)
 
FLA_Error FLA_Chol_u (FLA_Obj A, fla_chol_t *cntl)
 
FLA_Error FLA_Chol_solve (FLA_Uplo uplo, FLA_Obj A, FLA_Obj B, FLA_Obj X)
 
FLA_Error FLASH_Chol_solve (FLA_Uplo uplo, FLA_Obj A, FLA_Obj B, FLA_Obj X)
 

Function Documentation

◆ FLA_Chol_internal()

FLA_Error FLA_Chol_internal ( FLA_Uplo  uplo,
FLA_Obj  A,
fla_chol_t cntl 
)

References FLA_Check_error_level(), fla_chol_cntl_leaf, FLA_Chol_internal_check(), FLA_Chol_l(), FLA_Chol_u(), FLA_Obj_elemtype(), and FLASH_Queue_get_enabled().

Referenced by FLA_Chol(), FLA_Chol_l_blk_var1(), FLA_Chol_l_blk_var2(), FLA_Chol_l_blk_var3(), FLA_Chol_l_task(), FLA_Chol_task(), FLA_Chol_u_blk_var1(), FLA_Chol_u_blk_var2(), FLA_Chol_u_blk_var3(), FLA_Chol_u_task(), FLA_SPDinv_internal(), and FLASH_Chol().

17 {
18  FLA_Error r_val = FLA_SUCCESS;
19 
20  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
21  FLA_Chol_internal_check( uplo, A, cntl );
22 
23  if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
24  FLA_Obj_elemtype( A ) == FLA_MATRIX &&
25  FLA_Cntl_variant( cntl ) == FLA_SUBPROBLEM )
26  {
27  // Recurse
28  r_val = FLA_Chol_internal( uplo,
29  *FLASH_OBJ_PTR_AT( A ),
31  }
32  else if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
33  FLA_Obj_elemtype( A ) == FLA_SCALAR &&
35  {
36  // Enqueue
37  ENQUEUE_FLASH_Chol( uplo, A, cntl );
38  }
39  else
40  {
41  if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
42  FLA_Obj_elemtype( A ) == FLA_SCALAR &&
44  {
45  // Execute leaf
46  cntl = fla_chol_cntl_leaf;
47  }
48 
49  // Parameter combinations
50  if ( uplo == FLA_LOWER_TRIANGULAR )
51  {
52  r_val = FLA_Chol_l( A, cntl );
53  }
54  else if ( uplo == FLA_UPPER_TRIANGULAR )
55  {
56  r_val = FLA_Chol_u( A, cntl );
57  }
58  }
59 
60  return r_val;
61 }
fla_chol_t * flash_chol_cntl
Definition: FLASH_Chol_cntl_init.c:17
FLA_Error FLA_Chol_internal_check(FLA_Uplo uplo, FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_internal_check.c:13
int FLA_Error
Definition: FLA_type_defs.h:47
FLA_Bool FLASH_Queue_get_enabled(void)
Definition: FLASH_Queue.c:171
FLA_Error FLA_Chol_l(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_l.c:13
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
fla_chol_t * fla_chol_cntl_leaf
Definition: FLA_Chol_cntl_init.c:20
FLA_Error FLA_Chol_u(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_u.c:13
FLA_Error FLA_Chol_internal(FLA_Uplo uplo, FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_internal.c:16
FLA_Elemtype FLA_Obj_elemtype(FLA_Obj obj)
Definition: FLA_Query.c:51

◆ FLA_Chol_l()

FLA_Error FLA_Chol_l ( FLA_Obj  A,
fla_chol_t cntl 
)

References FLA_Chol_l_blk_ext(), FLA_Chol_l_blk_var1(), FLA_Chol_l_blk_var2(), FLA_Chol_l_blk_var3(), FLA_Chol_l_opt_var1(), FLA_Chol_l_opt_var2(), FLA_Chol_l_opt_var3(), FLA_Chol_l_unb_ext(), FLA_Chol_l_unb_var1(), FLA_Chol_l_unb_var2(), and FLA_Chol_l_unb_var3().

Referenced by FLA_Chol_internal().

14 {
15  FLA_Error r_val = FLA_SUCCESS;
16 
17  if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_EXTERN )
18  {
19  r_val = FLA_Chol_l_blk_ext( A );
20  }
21  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_EXTERN )
22  {
23  r_val = FLA_Chol_l_unb_ext( A );
24  }
25 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
26  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT1 )
27  {
28  r_val = FLA_Chol_l_opt_var1( A );
29  }
30 #endif
31  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT2 )
32  {
33  r_val = FLA_Chol_l_opt_var2( A );
34  }
35  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT3 )
36  {
37  r_val = FLA_Chol_l_opt_var3( A );
38  }
39 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
40  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
41  {
42  r_val = FLA_Chol_l_blk_var1( A, cntl );
43  }
44  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
45  {
46  r_val = FLA_Chol_l_blk_var2( A, cntl );
47  }
48 #endif
49  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
50  {
51  r_val = FLA_Chol_l_blk_var3( A, cntl );
52  }
53 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
54  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT1 )
55  {
56  r_val = FLA_Chol_l_unb_var1( A );
57  }
58  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT2 )
59  {
60  r_val = FLA_Chol_l_unb_var2( A );
61  }
62  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT3 )
63  {
64  r_val = FLA_Chol_l_unb_var3( A );
65  }
66 #endif
67  else
68  {
69  FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
70  }
71 
72  return r_val;
73 }
FLA_Error FLA_Chol_l_blk_var1(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_l_blk_var1.c:15
FLA_Error FLA_Chol_l_blk_var3(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_l_blk_var3.c:13
FLA_Error FLA_Chol_l_unb_var1(FLA_Obj A)
Definition: FLA_Chol_l_unb_var1.c:15
FLA_Error FLA_Chol_l_blk_var2(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_l_blk_var2.c:15
int FLA_Error
Definition: FLA_type_defs.h:47
FLA_Error FLA_Chol_l_blk_ext(FLA_Obj A)
Definition: FLA_Chol_blk_external.c:98
FLA_Error FLA_Chol_l_unb_var2(FLA_Obj A)
Definition: FLA_Chol_l_unb_var2.c:15
FLA_Error FLA_Chol_l_opt_var3(FLA_Obj A)
Definition: FLA_Chol_l_opt_var3.c:13
FLA_Error FLA_Chol_l_opt_var1(FLA_Obj A)
Definition: FLA_Chol_l_opt_var1.c:15
FLA_Error FLA_Chol_l_unb_var3(FLA_Obj A)
Definition: FLA_Chol_l_unb_var3.c:15
FLA_Error FLA_Chol_l_opt_var2(FLA_Obj A)
Definition: FLA_Chol_l_opt_var2.c:13
FLA_Error FLA_Chol_l_unb_ext(FLA_Obj A)
Definition: FLA_Chol_unb_external.c:98

◆ 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

◆ FLA_Chol_u()

FLA_Error FLA_Chol_u ( FLA_Obj  A,
fla_chol_t cntl 
)

References FLA_Chol_u_blk_ext(), FLA_Chol_u_blk_var1(), FLA_Chol_u_blk_var2(), FLA_Chol_u_blk_var3(), FLA_Chol_u_opt_var1(), FLA_Chol_u_opt_var2(), FLA_Chol_u_opt_var3(), FLA_Chol_u_unb_ext(), FLA_Chol_u_unb_var1(), FLA_Chol_u_unb_var2(), and FLA_Chol_u_unb_var3().

Referenced by FLA_Chol_internal().

14 {
15  FLA_Error r_val = FLA_SUCCESS;
16 
17  if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_EXTERN )
18  {
19  r_val = FLA_Chol_u_blk_ext( A );
20  }
21  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_EXTERN )
22  {
23  r_val = FLA_Chol_u_unb_ext( A );
24  }
25 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
26  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT1 )
27  {
28  r_val = FLA_Chol_u_opt_var1( A );
29  }
30 #endif
31  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT2 )
32  {
33  r_val = FLA_Chol_u_opt_var2( A );
34  }
35  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT3 )
36  {
37  r_val = FLA_Chol_u_opt_var3( A );
38  }
39 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
40  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
41  {
42  r_val = FLA_Chol_u_blk_var1( A, cntl );
43  }
44  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
45  {
46  r_val = FLA_Chol_u_blk_var2( A, cntl );
47  }
48 #endif
49  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
50  {
51  r_val = FLA_Chol_u_blk_var3( A, cntl );
52  }
53 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
54  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT1 )
55  {
56  r_val = FLA_Chol_u_unb_var1( A );
57  }
58  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT2 )
59  {
60  r_val = FLA_Chol_u_unb_var2( A );
61  }
62  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT3 )
63  {
64  r_val = FLA_Chol_u_unb_var3( A );
65  }
66 #endif
67  else
68  {
69  FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
70  }
71 
72  return r_val;
73 }
FLA_Error FLA_Chol_u_opt_var3(FLA_Obj A)
Definition: FLA_Chol_u_opt_var3.c:13
FLA_Error FLA_Chol_u_unb_var1(FLA_Obj A)
Definition: FLA_Chol_u_unb_var1.c:15
FLA_Error FLA_Chol_u_unb_ext(FLA_Obj A)
Definition: FLA_Chol_unb_external.c:103
FLA_Error FLA_Chol_u_unb_var3(FLA_Obj A)
Definition: FLA_Chol_u_unb_var3.c:15
int FLA_Error
Definition: FLA_type_defs.h:47
FLA_Error FLA_Chol_u_blk_var1(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_u_blk_var1.c:15
FLA_Error FLA_Chol_u_blk_var2(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_u_blk_var2.c:15
FLA_Error FLA_Chol_u_opt_var2(FLA_Obj A)
Definition: FLA_Chol_u_opt_var2.c:13
FLA_Error FLA_Chol_u_unb_var2(FLA_Obj A)
Definition: FLA_Chol_u_unb_var2.c:15
FLA_Error FLA_Chol_u_opt_var1(FLA_Obj A)
Definition: FLA_Chol_u_opt_var1.c:15
FLA_Error FLA_Chol_u_blk_var3(FLA_Obj A, fla_chol_t *cntl)
Definition: FLA_Chol_u_blk_var3.c:13
FLA_Error FLA_Chol_u_blk_ext(FLA_Obj A)
Definition: FLA_Chol_blk_external.c:103

◆ FLASH_Chol_solve()

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

References FLA_Check_error_level(), FLA_Chol_solve_check(), FLA_ONE, FLASH_Copy(), and FLASH_Trsm().

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  FLASH_Copy( B, X );
20 
21  if ( uplo == FLA_LOWER_TRIANGULAR )
22  {
23  FLASH_Trsm( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE,
24  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
25  FLASH_Trsm( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_CONJ_TRANSPOSE,
26  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
27  }
28  else // if ( uplo == FLA_UPPER_TRIANGULAR )
29  {
30  FLASH_Trsm( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_CONJ_TRANSPOSE,
31  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
32  FLASH_Trsm( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE,
33  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
34  }
35 
36  return FLA_SUCCESS;
37 }
FLA_Obj FLA_ONE
Definition: FLA_Init.c:18
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_Error FLASH_Copy(FLA_Obj A, FLA_Obj B)
Definition: FLASH_Copy.c:15
FLA_Error FLASH_Trsm(FLA_Side side, FLA_Uplo uplo, FLA_Trans trans, FLA_Diag diag, FLA_Obj alpha, FLA_Obj A, FLA_Obj B)
Definition: FLASH_Trsm.c:15