libflame  revision_anchor
Functions
FLA_LU_nopiv.h File Reference

(r)

Go to the source code of this file.

Functions

FLA_Error FLA_LU_nopiv_internal (FLA_Obj A, fla_lu_t *cntl)
 
FLA_Error FLA_LU_nopiv_solve (FLA_Obj A, FLA_Obj B, FLA_Obj X)
 
FLA_Error FLASH_LU_nopiv_solve (FLA_Obj A, FLA_Obj B, FLA_Obj X)
 

Function Documentation

◆ FLA_LU_nopiv_internal()

FLA_Error FLA_LU_nopiv_internal ( FLA_Obj  A,
fla_lu_t cntl 
)

References FLA_Check_error_level(), FLA_LU_nopiv_blk_var1(), FLA_LU_nopiv_blk_var2(), FLA_LU_nopiv_blk_var3(), FLA_LU_nopiv_blk_var4(), FLA_LU_nopiv_blk_var5(), fla_lu_nopiv_cntl_leaf, FLA_LU_nopiv_internal_check(), FLA_LU_nopiv_opt_var1(), FLA_LU_nopiv_opt_var2(), FLA_LU_nopiv_opt_var3(), FLA_LU_nopiv_opt_var4(), FLA_LU_nopiv_opt_var5(), FLA_LU_nopiv_unb_var1(), FLA_LU_nopiv_unb_var2(), FLA_LU_nopiv_unb_var3(), FLA_LU_nopiv_unb_var4(), FLA_LU_nopiv_unb_var5(), FLA_Obj_elemtype(), and FLASH_Queue_get_enabled().

Referenced by FLA_LU_nopiv(), FLA_LU_nopiv_blk_var1(), FLA_LU_nopiv_blk_var2(), FLA_LU_nopiv_blk_var3(), FLA_LU_nopiv_blk_var4(), FLA_LU_nopiv_blk_var5(), FLA_LU_nopiv_task(), and FLASH_LU_nopiv().

17 {
18  FLA_Error r_val = FLA_SUCCESS;
19 
20  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
21  FLA_LU_nopiv_internal_check( 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_LU_nopiv_internal( *FLASH_OBJ_PTR_AT( A ),
30  }
31  else if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
32  FLA_Obj_elemtype( A ) == FLA_SCALAR &&
34  {
35  // Enqueue
36  ENQUEUE_FLASH_LU_nopiv( A, cntl );
37  }
38  else
39  {
40  if ( FLA_Cntl_matrix_type( cntl ) == FLA_HIER &&
41  FLA_Obj_elemtype( A ) == FLA_SCALAR &&
43  {
44  // Execute leaf
46  }
47 
48  // Choose implementation.
49  if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT1 )
50  {
51  r_val = FLA_LU_nopiv_opt_var1( A );
52  }
53 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
54  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT2 )
55  {
56  r_val = FLA_LU_nopiv_opt_var2( A );
57  }
58  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT3 )
59  {
60  r_val = FLA_LU_nopiv_opt_var3( A );
61  }
62  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT4 )
63  {
64  r_val = FLA_LU_nopiv_opt_var4( A );
65  }
66 #endif
67  else if ( FLA_Cntl_variant( cntl ) == FLA_UNB_OPT_VARIANT5 )
68  {
69  r_val = FLA_LU_nopiv_opt_var5( A );
70  }
71 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
72  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT1 )
73  {
74  r_val = FLA_LU_nopiv_unb_var1( A );
75  }
76  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT2 )
77  {
78  r_val = FLA_LU_nopiv_unb_var2( A );
79  }
80  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT3 )
81  {
82  r_val = FLA_LU_nopiv_unb_var3( A );
83  }
84  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT4 )
85  {
86  r_val = FLA_LU_nopiv_unb_var4( A );
87  }
88 #endif
89  else if ( FLA_Cntl_variant( cntl ) == FLA_UNBLOCKED_VARIANT5 )
90  {
91  r_val = FLA_LU_nopiv_unb_var5( A );
92  }
93 #ifdef FLA_ENABLE_NON_CRITICAL_CODE
94  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT1 )
95  {
96  r_val = FLA_LU_nopiv_blk_var1( A, cntl );
97  }
98  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT2 )
99  {
100  r_val = FLA_LU_nopiv_blk_var2( A, cntl );
101  }
102  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT3 )
103  {
104  r_val = FLA_LU_nopiv_blk_var3( A, cntl );
105  }
106  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT4 )
107  {
108  r_val = FLA_LU_nopiv_blk_var4( A, cntl );
109  }
110 #endif
111  else if ( FLA_Cntl_variant( cntl ) == FLA_BLOCKED_VARIANT5 )
112  {
113  r_val = FLA_LU_nopiv_blk_var5( A, cntl );
114  }
115  else
116  {
117  FLA_Check_error_code( FLA_NOT_YET_IMPLEMENTED );
118  }
119  }
120 
121  return r_val;
122 }
FLA_Error FLA_LU_nopiv_opt_var3(FLA_Obj A)
Definition: FLA_LU_nopiv_opt_var3.c:15
FLA_Error FLA_LU_nopiv_blk_var3(FLA_Obj A, fla_lu_t *cntl)
Definition: FLA_LU_nopiv_blk_var3.c:15
fla_lu_t * fla_lu_nopiv_cntl_leaf
Definition: FLA_LU_nopiv_cntl_init.c:20
FLA_Error FLA_LU_nopiv_internal_check(FLA_Obj A, fla_lu_t *cntl)
Definition: FLA_LU_nopiv_internal_check.c:13
FLA_Error FLA_LU_nopiv_opt_var1(FLA_Obj A)
Definition: FLA_LU_nopiv_opt_var1.c:15
FLA_Error FLA_LU_nopiv_blk_var1(FLA_Obj A, fla_lu_t *cntl)
Definition: FLA_LU_nopiv_blk_var1.c:15
FLA_Error FLA_LU_nopiv_blk_var2(FLA_Obj A, fla_lu_t *cntl)
Definition: FLA_LU_nopiv_blk_var2.c:15
FLA_Error FLA_LU_nopiv_blk_var5(FLA_Obj A, fla_lu_t *cntl)
Definition: FLA_LU_nopiv_blk_var5.c:13
FLA_Error FLA_LU_nopiv_unb_var1(FLA_Obj A)
Definition: FLA_LU_nopiv_unb_var1.c:15
int FLA_Error
Definition: FLA_type_defs.h:47
FLA_Error FLA_LU_nopiv_unb_var2(FLA_Obj A)
Definition: FLA_LU_nopiv_unb_var2.c:15
FLA_Bool FLASH_Queue_get_enabled(void)
Definition: FLASH_Queue.c:171
FLA_Error FLA_LU_nopiv_unb_var5(FLA_Obj A)
Definition: FLA_LU_nopiv_unb_var5.c:15
FLA_Error FLA_LU_nopiv_opt_var2(FLA_Obj A)
Definition: FLA_LU_nopiv_opt_var2.c:15
FLA_Error FLA_LU_nopiv_opt_var4(FLA_Obj A)
Definition: FLA_LU_nopiv_opt_var4.c:15
FLA_Error FLA_LU_nopiv_opt_var5(FLA_Obj A)
Definition: FLA_LU_nopiv_opt_var5.c:13
FLA_Error FLA_LU_nopiv_blk_var4(FLA_Obj A, fla_lu_t *cntl)
Definition: FLA_LU_nopiv_blk_var4.c:15
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
FLA_Error FLA_LU_nopiv_unb_var3(FLA_Obj A)
Definition: FLA_LU_nopiv_unb_var3.c:15
fla_lu_t * flash_lu_nopiv_cntl
Definition: FLASH_LU_nopiv_cntl_init.c:17
FLA_Error FLA_LU_nopiv_internal(FLA_Obj A, fla_lu_t *cntl)
Definition: FLA_LU_nopiv_internal.c:16
FLA_Error FLA_LU_nopiv_unb_var4(FLA_Obj A)
Definition: FLA_LU_nopiv_unb_var4.c:15
FLA_Elemtype FLA_Obj_elemtype(FLA_Obj obj)
Definition: FLA_Query.c:51

◆ FLA_LU_nopiv_solve()

FLA_Error FLA_LU_nopiv_solve ( FLA_Obj  A,
FLA_Obj  B,
FLA_Obj  X 
)

References FLA_Check_error_level(), FLA_Copy_external(), FLA_LU_nopiv_solve_check(), 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_LU_nopiv_solve_check( A, B, X );
18 
19  if ( FLA_Obj_is_identical( B, X ) == FALSE )
20  FLA_Copy_external( B, X );
21 
22  FLA_Trsm_external( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE,
23  FLA_UNIT_DIAG, FLA_ONE, A, X );
24  FLA_Trsm_external( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE,
25  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
26 
27  return FLA_SUCCESS;
28 }
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_LU_nopiv_solve_check(FLA_Obj A, FLA_Obj B, FLA_Obj X)
Definition: FLA_LU_nopiv_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

◆ FLASH_LU_nopiv_solve()

FLA_Error FLASH_LU_nopiv_solve ( FLA_Obj  A,
FLA_Obj  B,
FLA_Obj  X 
)

References FLA_Check_error_level(), FLA_LU_nopiv_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_LU_nopiv_solve_check( A, B, X );
18 
19  FLASH_Copy( B, X );
20 
21  FLASH_Trsm( FLA_LEFT, FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE,
22  FLA_UNIT_DIAG, FLA_ONE, A, X );
23  FLASH_Trsm( FLA_LEFT, FLA_UPPER_TRIANGULAR, FLA_NO_TRANSPOSE,
24  FLA_NONUNIT_DIAG, FLA_ONE, A, X );
25 
26  return FLA_SUCCESS;
27 }
FLA_Obj FLA_ONE
Definition: FLA_Init.c:18
FLA_Error FLA_LU_nopiv_solve_check(FLA_Obj A, FLA_Obj B, FLA_Obj X)
Definition: FLA_LU_nopiv_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