libflame  revision_anchor
Functions
FLASH_FS_incpiv_aux1.c File Reference

(r)

Functions

FLA_Error FLASH_FS_incpiv_aux1 (FLA_Obj A, FLA_Obj p, FLA_Obj L, FLA_Obj b, dim_t nb_alg)
 

Function Documentation

◆ FLASH_FS_incpiv_aux1()

FLA_Error FLASH_FS_incpiv_aux1 ( FLA_Obj  A,
FLA_Obj  p,
FLA_Obj  L,
FLA_Obj  b,
dim_t  nb_alg 
)

References FLA_Apply_pivots(), FLA_Cont_with_3x1_to_2x1(), FLA_Cont_with_3x3_to_2x2(), FLA_Obj_length(), FLA_Obj_width(), FLA_Part_2x1(), FLA_Part_2x2(), FLA_Repart_2x1_to_3x1(), FLA_Repart_2x2_to_3x3(), FLA_Trsv_external(), and FLASH_FS_incpiv_aux2().

Referenced by FLASH_FS_incpiv().

14 {
15  FLA_Obj ATL, ATR, A00, A01, A02,
16  ABL, ABR, A10, A11, A12,
17  A20, A21, A22;
18 
19  FLA_Obj pTL, pTR, p00, p01, p02,
20  pBL, pBR, p10, p11, p12,
21  p20, p21, p22;
22 
23  FLA_Obj LTL, LTR, L00, L01, L02,
24  LBL, LBR, L10, L11, L12,
25  L20, L21, L22;
26 
27  FLA_Obj bT, b0,
28  bB, b1,
29  b2;
30 
31  FLA_Obj p11_conf,
32  p11_rest;
33 
34  FLA_Part_2x2( A, &ATL, &ATR,
35  &ABL, &ABR, 0, 0, FLA_TL );
36 
37  FLA_Part_2x2( p, &pTL, &pTR,
38  &pBL, &pBR, 0, 0, FLA_TL );
39 
40  FLA_Part_2x2( L, &LTL, &LTR,
41  &LBL, &LBR, 0, 0, FLA_TL );
42 
43  FLA_Part_2x1( b, &bT,
44  &bB, 0, FLA_TOP );
45 
46  while ( FLA_Obj_length( ATL ) < FLA_Obj_length( A ) &&
47  FLA_Obj_width ( ATL ) < FLA_Obj_width ( A ) )
48  {
49  FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, /**/ &A01, &A02,
50  /* ************* */ /* ******************** */
51  &A10, /**/ &A11, &A12,
52  ABL, /**/ ABR, &A20, /**/ &A21, &A22,
53  1, 1, FLA_BR );
54 
55  FLA_Repart_2x2_to_3x3( pTL, /**/ pTR, &p00, /**/ &p01, &p02,
56  /* ************* */ /* ******************** */
57  &p10, /**/ &p11, &p12,
58  pBL, /**/ pBR, &p20, /**/ &p21, &p22,
59  1, 1, FLA_BR );
60 
61  FLA_Repart_2x2_to_3x3( LTL, /**/ LTR, &L00, /**/ &L01, &L02,
62  /* ************* */ /* ******************** */
63  &L10, /**/ &L11, &L12,
64  LBL, /**/ LBR, &L20, /**/ &L21, &L22,
65  1, 1, FLA_BR );
66 
67  FLA_Repart_2x1_to_3x1( bT, &b0,
68  /* ** */ /* ** */
69  &b1,
70  bB, &b2, 1, FLA_BOTTOM );
71 
72  /*------------------------------------------------------------*/
73 
74  FLA_Part_2x1( *FLASH_OBJ_PTR_AT( p11 ), &p11_conf,
75  &p11_rest,
76  FLA_Obj_length( *FLASH_OBJ_PTR_AT( b1 ) ), FLA_TOP );
77 
78 
79  FLA_Apply_pivots( FLA_LEFT, FLA_NO_TRANSPOSE,
80  p11_conf,
81  *FLASH_OBJ_PTR_AT( b1 ) );
82 
83  FLA_Trsv_external( FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_UNIT_DIAG,
84  *FLASH_OBJ_PTR_AT( A11 ),
85  *FLASH_OBJ_PTR_AT( b1 ) );
86 
88  A21, p21, b1,
89  b2, nb_alg );
90 
91  /*------------------------------------------------------------*/
92 
93  FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, A01, /**/ A02,
94  A10, A11, /**/ A12,
95  /* ************** */ /* ****************** */
96  &ABL, /**/ &ABR, A20, A21, /**/ A22,
97  FLA_TL );
98 
99  FLA_Cont_with_3x3_to_2x2( &pTL, /**/ &pTR, p00, p01, /**/ p02,
100  p10, p11, /**/ p12,
101  /* ************** */ /* ****************** */
102  &pBL, /**/ &pBR, p20, p21, /**/ p22,
103  FLA_TL );
104 
105  FLA_Cont_with_3x3_to_2x2( &LTL, /**/ &LTR, L00, L01, /**/ L02,
106  L10, L11, /**/ L12,
107  /* ************** */ /* ****************** */
108  &LBL, /**/ &LBR, L20, L21, /**/ L22,
109  FLA_TL );
110 
111  FLA_Cont_with_3x1_to_2x1( &bT, b0,
112  b1,
113  /* ** */ /* ** */
114  &bB, b2, FLA_TOP );
115  }
116 
117  return FLA_SUCCESS;
118 }
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_Repart_2x2_to_3x3(FLA_Obj ATL, FLA_Obj ATR, FLA_Obj *A00, FLA_Obj *A01, FLA_Obj *A02, FLA_Obj *A10, FLA_Obj *A11, FLA_Obj *A12, FLA_Obj ABL, FLA_Obj ABR, FLA_Obj *A20, FLA_Obj *A21, FLA_Obj *A22, dim_t mb, dim_t nb, FLA_Quadrant quadrant)
Definition: FLA_View.c:142
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_Part_2x2(FLA_Obj A, FLA_Obj *A11, FLA_Obj *A12, FLA_Obj *A21, FLA_Obj *A22, dim_t mb, dim_t nb, FLA_Quadrant quadrant)
Definition: FLA_View.c:17
FLA_Error FLASH_FS_incpiv_aux2(FLA_Obj L, FLA_Obj D, FLA_Obj p, FLA_Obj C, FLA_Obj E, dim_t nb_alg)
Definition: FLASH_FS_incpiv_aux2.c:13
Definition: FLA_type_defs.h:158
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
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
FLA_Error FLA_Cont_with_3x3_to_2x2(FLA_Obj *ATL, FLA_Obj *ATR, FLA_Obj A00, FLA_Obj A01, FLA_Obj A02, FLA_Obj A10, FLA_Obj A11, FLA_Obj A12, FLA_Obj *ABL, FLA_Obj *ABR, FLA_Obj A20, FLA_Obj A21, FLA_Obj A22, FLA_Quadrant quadrant)
Definition: FLA_View.c:304
FLA_Error FLA_Apply_pivots(FLA_Side side, FLA_Trans trans, FLA_Obj p, FLA_Obj A)
Definition: FLA_Apply_pivots.c:15
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