libflame  revision_anchor
Functions
FLA_QR_unb_external.c File Reference

(r)

Functions

FLA_Error FLA_QR_unb_external (FLA_Obj A, FLA_Obj t)
 
FLA_Error FLA_QR_unb_ext (FLA_Obj A, FLA_Obj t)
 

Function Documentation

◆ FLA_QR_unb_ext()

FLA_Error FLA_QR_unb_ext ( FLA_Obj  A,
FLA_Obj  t 
)

References FLA_QR_unb_external().

111 {
112  return FLA_QR_unb_external( A, t );
113 }
FLA_Error FLA_QR_unb_external(FLA_Obj A, FLA_Obj t)
Definition: FLA_QR_unb_external.c:13

◆ FLA_QR_unb_external()

FLA_Error FLA_QR_unb_external ( FLA_Obj  A,
FLA_Obj  t 
)

References F77_cgeqr2(), F77_dgeqr2(), F77_sgeqr2(), F77_zgeqr2(), FLA_Check_error_level(), FLA_Obj_col_stride(), FLA_Obj_create(), FLA_Obj_datatype(), FLA_Obj_free(), FLA_Obj_has_zero_dim(), FLA_Obj_length(), FLA_Obj_width(), and FLA_QR_check().

Referenced by FLA_QR_unb_ext().

14 {
15  int info = 0;
16 #ifdef FLA_ENABLE_EXTERNAL_LAPACK_INTERFACES
17  FLA_Datatype datatype;
18  int m_A, n_A, cs_A;
19  FLA_Obj work_obj;
20 
21  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
22  FLA_QR_check( A, t );
23 
24  if ( FLA_Obj_has_zero_dim( A ) ) return FLA_SUCCESS;
25 
26  datatype = FLA_Obj_datatype( A );
27 
28  m_A = FLA_Obj_length( A );
29  n_A = FLA_Obj_width( A );
30  cs_A = FLA_Obj_col_stride( A );
31 
32  FLA_Obj_create( datatype, n_A, 1, 0, 0, &work_obj );
33 
34  switch( datatype ){
35 
36  case FLA_FLOAT:
37  {
38  float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
39  float *buff_t = ( float * ) FLA_FLOAT_PTR( t );
40  float *buff_work = ( float * ) FLA_FLOAT_PTR( work_obj );
41 
42  F77_sgeqr2( &m_A,
43  &n_A,
44  buff_A, &cs_A,
45  buff_t,
46  buff_work,
47  &info );
48 
49  break;
50  }
51 
52  case FLA_DOUBLE:
53  {
54  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
55  double *buff_t = ( double * ) FLA_DOUBLE_PTR( t );
56  double *buff_work = ( double * ) FLA_DOUBLE_PTR( work_obj );
57 
58  F77_dgeqr2( &m_A,
59  &n_A,
60  buff_A, &cs_A,
61  buff_t,
62  buff_work,
63  &info );
64 
65  break;
66  }
67 
68  case FLA_COMPLEX:
69  {
70  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
71  scomplex *buff_t = ( scomplex * ) FLA_COMPLEX_PTR( t );
72  scomplex *buff_work = ( scomplex * ) FLA_COMPLEX_PTR( work_obj );
73 
74  F77_cgeqr2( &m_A,
75  &n_A,
76  buff_A, &cs_A,
77  buff_t,
78  buff_work,
79  &info );
80 
81  break;
82  }
83 
84  case FLA_DOUBLE_COMPLEX:
85  {
86  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
87  dcomplex *buff_t = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( t );
88  dcomplex *buff_work = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( work_obj );
89 
90  F77_zgeqr2( &m_A,
91  &n_A,
92  buff_A, &cs_A,
93  buff_t,
94  buff_work,
95  &info );
96 
97  break;
98  }
99 
100  }
101 
102  FLA_Obj_free( &work_obj );
103 #else
104  FLA_Check_error_code( FLA_EXTERNAL_LAPACK_NOT_IMPLEMENTED );
105 #endif
106 
107  return info;
108 }
FLA_Error FLA_Obj_create(FLA_Datatype datatype, dim_t m, dim_t n, dim_t rs, dim_t cs, FLA_Obj *obj)
Definition: FLA_Obj.c:55
int F77_zgeqr2(int *m, int *n, dcomplex *a, int *lda, dcomplex *tau, dcomplex *work, int *info)
FLA_Error FLA_Obj_free(FLA_Obj *obj)
Definition: FLA_Obj.c:588
int F77_cgeqr2(int *m, int *n, scomplex *a, int *lda, scomplex *tau, scomplex *work, int *info)
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
FLA_Error FLA_QR_check(FLA_Obj A, FLA_Obj t)
Definition: FLA_QR_check.c:13
Definition: FLA_type_defs.h:158
FLA_Bool FLA_Obj_has_zero_dim(FLA_Obj A)
Definition: FLA_Query.c:400
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
int F77_dgeqr2(int *m, int *n, double *a, int *lda, double *tau, double *work, int *info)
Definition: blis_type_defs.h:132
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
int FLA_Datatype
Definition: FLA_type_defs.h:49
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
int F77_sgeqr2(int *m, int *n, float *a, int *lda, float *tau, float *work, int *info)
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
Definition: blis_type_defs.h:137