libflame  revision_anchor
Functions
FLA_Trsm_rlh_unb_var3.c File Reference

(r)

Functions

FLA_Error FLA_Trsm_rlh_unb_var3 (FLA_Diag diagA, FLA_Obj alpha, FLA_Obj A, FLA_Obj B)
 

Function Documentation

◆ FLA_Trsm_rlh_unb_var3()

FLA_Error FLA_Trsm_rlh_unb_var3 ( FLA_Diag  diagA,
FLA_Obj  alpha,
FLA_Obj  A,
FLA_Obj  B 
)

References FLA_Cont_with_3x1_to_2x1(), FLA_Obj_length(), FLA_Part_2x1(), FLA_Repart_2x1_to_3x1(), FLA_Scal_external(), and FLA_Trsv_external().

Referenced by FLA_Trsm_rlh().

16 {
17  FLA_Obj BT, B0,
18  BB, b1t,
19  B2;
20 
21  FLA_Scal_external( alpha, B );
22 
23  FLA_Part_2x1( B, &BT,
24  &BB, 0, FLA_TOP );
25 
26  while ( FLA_Obj_length( BT ) < FLA_Obj_length( B ) ){
27 
28  FLA_Repart_2x1_to_3x1( BT, &B0,
29  /* ** */ /* *** */
30  &b1t,
31  BB, &B2, 1, FLA_BOTTOM );
32 
33  /*------------------------------------------------------------*/
34 
35  /* b1t = b1t / tril( A' ); */
36  FLA_Trsv_external( FLA_LOWER_TRIANGULAR, FLA_CONJ_NO_TRANSPOSE, diagA, A, b1t );
37 
38  /*------------------------------------------------------------*/
39 
40  FLA_Cont_with_3x1_to_2x1( &BT, B0,
41  b1t,
42  /* ** */ /* *** */
43  &BB, B2, FLA_TOP );
44 
45  }
46 
47  return FLA_SUCCESS;
48 }
FLA_Error FLA_Repart_2x1_to_3x1(FLA_Obj AT, FLA_Obj *A0, FLA_Obj *A1, FLA_Obj AB, FLA_Obj *A2, dim_t mb, FLA_Side side)
Definition: FLA_View.c:226
FLA_Error FLA_Scal_external(FLA_Obj alpha, FLA_Obj A)
Definition: FLA_Scal_external.c:13
FLA_Error FLA_Cont_with_3x1_to_2x1(FLA_Obj *AT, FLA_Obj A0, FLA_Obj A1, FLA_Obj *AB, FLA_Obj A2, FLA_Side side)
Definition: FLA_View.c:428
Definition: FLA_type_defs.h:158
FLA_Error FLA_Part_2x1(FLA_Obj A, FLA_Obj *A1, FLA_Obj *A2, dim_t mb, FLA_Side side)
Definition: FLA_View.c:76
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
FLA_Error FLA_Trsv_external(FLA_Uplo uplo, FLA_Trans transa, FLA_Diag diag, FLA_Obj A, FLA_Obj x)
Definition: FLA_Trsv_external.c:13