libflame  revision_anchor
Functions
FLA_LQ_blk_external.c File Reference

(r)

Functions

FLA_Error FLA_LQ_blk_external (FLA_Obj A, FLA_Obj t)
 
FLA_Error FLA_LQ_blk_ext (FLA_Obj A, FLA_Obj t)
 

Function Documentation

◆ FLA_LQ_blk_ext()

FLA_Error FLA_LQ_blk_ext ( FLA_Obj  A,
FLA_Obj  t 
)
113 {
114  return FLA_LQ_blk_external( A, t );
115 }
FLA_Error FLA_LQ_blk_external(FLA_Obj A, FLA_Obj t)
Definition: FLA_LQ_blk_external.c:13

References FLA_LQ_blk_external().

◆ FLA_LQ_blk_external()

FLA_Error FLA_LQ_blk_external ( FLA_Obj  A,
FLA_Obj  t 
)
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  int lwork;
20  FLA_Obj work_obj;
21 
22  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
23  FLA_LQ_check( A, t );
24 
25  if ( FLA_Obj_has_zero_dim( A ) ) return FLA_SUCCESS;
26 
27  datatype = FLA_Obj_datatype( A );
28 
29  m_A = FLA_Obj_length( A );
30  n_A = FLA_Obj_width( A );
31  cs_A = FLA_Obj_col_stride( A );
32  lwork = m_A * FLA_Query_blocksize( datatype, FLA_DIMENSION_MIN );
33 
34  FLA_Obj_create( datatype, lwork, 1, 0, 0, &work_obj );
35 
36  switch( datatype ){
37 
38  case FLA_FLOAT:
39  {
40  float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
41  float *buff_t = ( float * ) FLA_FLOAT_PTR( t );
42  float *buff_work = ( float * ) FLA_FLOAT_PTR( work_obj );
43 
44  F77_sgelqf( &m_A,
45  &n_A,
46  buff_A, &cs_A,
47  buff_t,
48  buff_work, &lwork,
49  &info );
50 
51  break;
52  }
53 
54  case FLA_DOUBLE:
55  {
56  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
57  double *buff_t = ( double * ) FLA_DOUBLE_PTR( t );
58  double *buff_work = ( double * ) FLA_DOUBLE_PTR( work_obj );
59 
60  F77_dgelqf( &m_A,
61  &n_A,
62  buff_A, &cs_A,
63  buff_t,
64  buff_work, &lwork,
65  &info );
66 
67  break;
68  }
69 
70  case FLA_COMPLEX:
71  {
72  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
73  scomplex *buff_t = ( scomplex * ) FLA_COMPLEX_PTR( t );
74  scomplex *buff_work = ( scomplex * ) FLA_COMPLEX_PTR( work_obj );
75 
76  F77_cgelqf( &m_A,
77  &n_A,
78  buff_A, &cs_A,
79  buff_t,
80  buff_work, &lwork,
81  &info );
82 
83  break;
84  }
85 
86  case FLA_DOUBLE_COMPLEX:
87  {
88  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
89  dcomplex *buff_t = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( t );
90  dcomplex *buff_work = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( work_obj );
91 
92  F77_zgelqf( &m_A,
93  &n_A,
94  buff_A, &cs_A,
95  buff_t,
96  buff_work, &lwork,
97  &info );
98 
99  break;
100  }
101 
102  }
103 
104  FLA_Obj_free( &work_obj );
105 #else
106  FLA_Check_error_code( FLA_EXTERNAL_LAPACK_NOT_IMPLEMENTED );
107 #endif
108 
109  return info;
110 }
FLA_Error FLA_LQ_check(FLA_Obj A, FLA_Obj t)
Definition: FLA_LQ_check.c:13
int F77_cgelqf(int *m, int *n, scomplex *a, int *lda, scomplex *tau, scomplex *work, int *lwork, int *info)
int F77_zgelqf(int *m, int *n, dcomplex *a, int *lda, dcomplex *tau, dcomplex *work, int *lwork, int *info)
int F77_dgelqf(int *m, int *n, double *a, int *lda, double *tau, double *work, int *lwork, int *info)
int F77_sgelqf(int *m, int *n, float *a, int *lda, float *tau, float *work, int *lwork, int *info)
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
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
FLA_Bool FLA_Obj_has_zero_dim(FLA_Obj A)
Definition: FLA_Query.c:400
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
dim_t FLA_Query_blocksize(FLA_Datatype dt, FLA_Dimension dim)
Definition: FLA_Blocksize.c:161
FLA_Error FLA_Obj_free(FLA_Obj *obj)
Definition: FLA_Obj.c:588
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
int FLA_Datatype
Definition: FLA_type_defs.h:49
Definition: FLA_type_defs.h:159
Definition: blis_type_defs.h:138
Definition: blis_type_defs.h:133

References F77_cgelqf(), F77_dgelqf(), F77_sgelqf(), F77_zgelqf(), FLA_Check_error_level(), FLA_LQ_check(), 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_Query_blocksize().

Referenced by FLA_LQ_blk_ext().