libflame  revision_anchor
Functions
FLA_Hess_blk_external.c File Reference

(r)

Functions

FLA_Error FLA_Hess_blk_external (FLA_Obj A, FLA_Obj t, int ilo, int ihi)
 
FLA_Error FLA_Hess_blk_ext (FLA_Obj A, FLA_Obj t, int ilo, int ihi)
 

Function Documentation

◆ FLA_Hess_blk_ext()

FLA_Error FLA_Hess_blk_ext ( FLA_Obj  A,
FLA_Obj  t,
int  ilo,
int  ihi 
)
125 {
126  return FLA_Hess_blk_external( A, t, ilo, ihi );
127 }
FLA_Error FLA_Hess_blk_external(FLA_Obj A, FLA_Obj t, int ilo, int ihi)
Definition: FLA_Hess_blk_external.c:13

References FLA_Hess_blk_external().

◆ FLA_Hess_blk_external()

FLA_Error FLA_Hess_blk_external ( FLA_Obj  A,
FLA_Obj  t,
int  ilo,
int  ihi 
)
14 {
15  int info = 0;
16 #ifdef FLA_ENABLE_EXTERNAL_LAPACK_INTERFACES
17  FLA_Datatype datatype;
18  int 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_Hess_check( A, t, ilo, ihi );
24 
25  if ( FLA_Obj_has_zero_dim( A ) ) return FLA_SUCCESS;
26 
27  datatype = FLA_Obj_datatype( A );
28 
29  n_A = FLA_Obj_width( A );
30  cs_A = FLA_Obj_col_stride( A );
31 
32  lwork = n_A * FLA_Query_blocksize( datatype, FLA_DIMENSION_MIN );
33 
34  FLA_Obj_create( datatype, lwork, 1, 0, 0, &work_obj );
35 
36  // Shift ilo and ihi from zero-based indexing to one-based indexing.
37  ilo += 1;
38  ihi += 1;
39 
40  switch( datatype ){
41 
42  case FLA_FLOAT:
43  {
44  float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
45  float *buff_t = ( float * ) FLA_FLOAT_PTR( t );
46  float *buff_work = ( float * ) FLA_FLOAT_PTR( work_obj );
47 
48  F77_sgehrd( &n_A,
49  &ilo,
50  &ihi,
51  buff_A, &cs_A,
52  buff_t,
53  buff_work,
54  &lwork,
55  &info );
56 
57  break;
58  }
59 
60  case FLA_DOUBLE:
61  {
62  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
63  double *buff_t = ( double * ) FLA_DOUBLE_PTR( t );
64  double *buff_work = ( double * ) FLA_DOUBLE_PTR( work_obj );
65 
66  F77_dgehrd( &n_A,
67  &ilo,
68  &ihi,
69  buff_A, &cs_A,
70  buff_t,
71  buff_work,
72  &lwork,
73  &info );
74 
75  break;
76  }
77 
78  case FLA_COMPLEX:
79  {
80  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
81  scomplex *buff_t = ( scomplex * ) FLA_COMPLEX_PTR( t );
82  scomplex *buff_work = ( scomplex * ) FLA_COMPLEX_PTR( work_obj );
83 
84  F77_cgehrd( &n_A,
85  &ilo,
86  &ihi,
87  buff_A, &cs_A,
88  buff_t,
89  buff_work,
90  &lwork,
91  &info );
92 
93  break;
94  }
95 
96  case FLA_DOUBLE_COMPLEX:
97  {
98  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
99  dcomplex *buff_t = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( t );
100  dcomplex *buff_work = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( work_obj );
101 
102  F77_zgehrd( &n_A,
103  &ilo,
104  &ihi,
105  buff_A, &cs_A,
106  buff_t,
107  buff_work,
108  &lwork,
109  &info );
110 
111  break;
112  }
113 
114  }
115 
116  FLA_Obj_free( &work_obj );
117 #else
118  FLA_Check_error_code( FLA_EXTERNAL_LAPACK_NOT_IMPLEMENTED );
119 #endif
120 
121  return info;
122 }
FLA_Error FLA_Hess_check(FLA_Obj A, FLA_Obj t, int ilo, int ihi)
Definition: FLA_Hess_check.c:13
int F77_zgehrd(int *n, int *ilo, int *ihi, dcomplex *a, int *lda, dcomplex *tau, dcomplex *work, int *lwork, int *info)
int F77_dgehrd(int *n, int *ilo, int *ihi, double *a, int *lda, double *tau, double *work, int *lwork, int *info)
int F77_cgehrd(int *n, int *ilo, int *ihi, scomplex *a, int *lda, scomplex *tau, scomplex *work, int *lwork, int *info)
int F77_sgehrd(int *n, int *ilo, int *ihi, 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_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_cgehrd(), F77_dgehrd(), F77_sgehrd(), F77_zgehrd(), FLA_Check_error_level(), FLA_Hess_check(), FLA_Obj_col_stride(), FLA_Obj_create(), FLA_Obj_datatype(), FLA_Obj_free(), FLA_Obj_has_zero_dim(), FLA_Obj_width(), and FLA_Query_blocksize().

Referenced by FLA_Hess_blk_ext().