libflame  revision_anchor
Functions
FLA_Her2k_lh_unb_var9.c File Reference

(r)

Functions

FLA_Error FLA_Her2k_lh_unb_var9 (FLA_Obj alpha, FLA_Obj A, FLA_Obj B, FLA_Obj beta, FLA_Obj C)
 

Function Documentation

◆ FLA_Her2k_lh_unb_var9()

FLA_Error FLA_Her2k_lh_unb_var9 ( FLA_Obj  alpha,
FLA_Obj  A,
FLA_Obj  B,
FLA_Obj  beta,
FLA_Obj  C 
)

References FLA_Cont_with_3x1_to_2x1(), FLA_Her2c_external(), FLA_Obj_length(), FLA_Part_2x1(), FLA_Repart_2x1_to_3x1(), and FLA_Scalr_external().

Referenced by FLA_Her2k_lh().

16 {
17  FLA_Obj AT, A0,
18  AB, a1t,
19  A2;
20 
21  FLA_Obj BT, B0,
22  BB, b1t,
23  B2;
24 
25  FLA_Scalr_external( FLA_LOWER_TRIANGULAR, beta, C );
26 
27  FLA_Part_2x1( A, &AT,
28  &AB, 0, FLA_TOP );
29 
30  FLA_Part_2x1( B, &BT,
31  &BB, 0, FLA_TOP );
32 
33  while ( FLA_Obj_length( AT ) < FLA_Obj_length( A ) ){
34 
35 
36  FLA_Repart_2x1_to_3x1( AT, &A0,
37  /* ** */ /* ** */
38  &a1t,
39  AB, &A2, 1, FLA_BOTTOM );
40 
41  FLA_Repart_2x1_to_3x1( BT, &B0,
42  /* ** */ /* ** */
43  &b1t,
44  BB, &B2, 1, FLA_BOTTOM );
45 
46  /*------------------------------------------------------------*/
47 
48  /* C = C + a1t' * b1t + b1t' * a1t */
49  FLA_Her2c_external( FLA_LOWER_TRIANGULAR, FLA_CONJUGATE, alpha, a1t, b1t, C );
50 
51  /*------------------------------------------------------------*/
52 
53  FLA_Cont_with_3x1_to_2x1( &AT, A0,
54  a1t,
55  /* ** */ /* ** */
56  &AB, A2, FLA_TOP );
57 
58  FLA_Cont_with_3x1_to_2x1( &BT, B0,
59  b1t,
60  /* ** */ /* ** */
61  &BB, B2, FLA_TOP );
62 
63  }
64 
65  return FLA_SUCCESS;
66 }
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_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
FLA_Error FLA_Scalr_external(FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A)
Definition: FLA_Scalr_external.c:13
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_Her2c_external(FLA_Uplo uplo, FLA_Conj conj, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj A)
Definition: FLA_Her2c_external.c:13