libflame  revision_anchor
Functions
FLA_Tridiag_unb_external.c File Reference

(r)

Functions

FLA_Error FLA_Tridiag_unb_external (FLA_Uplo uplo, FLA_Obj A, FLA_Obj t)
 
FLA_Error FLA_Tridiag_unb_ext (FLA_Uplo uplo, FLA_Obj A, FLA_Obj t)
 

Function Documentation

◆ FLA_Tridiag_unb_ext()

FLA_Error FLA_Tridiag_unb_ext ( FLA_Uplo  uplo,
FLA_Obj  A,
FLA_Obj  t 
)
124 {
125  return FLA_Tridiag_unb_external( uplo, A, t );
126 }
FLA_Error FLA_Tridiag_unb_external(FLA_Uplo uplo, FLA_Obj A, FLA_Obj t)
Definition: FLA_Tridiag_unb_external.c:13

References FLA_Tridiag_unb_external().

◆ FLA_Tridiag_unb_external()

FLA_Error FLA_Tridiag_unb_external ( FLA_Uplo  uplo,
FLA_Obj  A,
FLA_Obj  t 
)
14 {
15  int info = 0;
16 #ifdef FLA_ENABLE_EXTERNAL_LAPACK_INTERFACES
17  FLA_Datatype datatype;
18  int n_A, cs_A;
19  FLA_Obj d, e;
20  char blas_uplo;
21 
22  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
23  FLA_Tridiag_check( uplo, A, t );
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  FLA_Obj_create( FLA_Obj_datatype_proj_to_real( A ), n_A, 1, 0, 0, &d );
33  FLA_Obj_create( FLA_Obj_datatype_proj_to_real( A ), n_A - 1, 1, 0, 0, &e );
34 
35  FLA_Param_map_flame_to_netlib_uplo( uplo, &blas_uplo );
36 
37 
38  switch( datatype ){
39 
40  case FLA_FLOAT:
41  {
42  float* buff_A = ( float * ) FLA_FLOAT_PTR( A );
43  float* buff_d = ( float * ) FLA_FLOAT_PTR( d );
44  float* buff_e = ( float * ) FLA_FLOAT_PTR( e );
45  float* buff_t = ( float * ) FLA_FLOAT_PTR( t );
46 
47  F77_ssytd2( &blas_uplo,
48  &n_A,
49  buff_A, &cs_A,
50  buff_d,
51  buff_e,
52  buff_t,
53  &info );
54 
55  break;
56  }
57 
58  case FLA_DOUBLE:
59  {
60  double* buff_A = ( double * ) FLA_DOUBLE_PTR( A );
61  double* buff_d = ( double * ) FLA_DOUBLE_PTR( d );
62  double* buff_e = ( double * ) FLA_DOUBLE_PTR( e );
63  double* buff_t = ( double * ) FLA_DOUBLE_PTR( t );
64 
65  F77_dsytd2( &blas_uplo,
66  &n_A,
67  buff_A, &cs_A,
68  buff_d,
69  buff_e,
70  buff_t,
71  &info );
72 
73  break;
74  }
75 
76  case FLA_COMPLEX:
77  {
78  scomplex* buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
79  float* buff_d = ( float * ) FLA_FLOAT_PTR( d );
80  float* buff_e = ( float * ) FLA_FLOAT_PTR( e );
81  scomplex* buff_t = ( scomplex * ) FLA_COMPLEX_PTR( t );
82 
83  F77_chetd2( &blas_uplo,
84  &n_A,
85  buff_A, &cs_A,
86  buff_d,
87  buff_e,
88  buff_t,
89  &info );
90 
91  break;
92  }
93 
94  case FLA_DOUBLE_COMPLEX:
95  {
96  dcomplex* buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
97  double* buff_d = ( double * ) FLA_DOUBLE_PTR( d );
98  double* buff_e = ( double * ) FLA_DOUBLE_PTR( e );
99  dcomplex* buff_t = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( t );
100 
101  F77_zhetd2( &blas_uplo,
102  &n_A,
103  buff_A, &cs_A,
104  buff_d,
105  buff_e,
106  buff_t,
107  &info );
108 
109  break;
110  }
111 
112  }
113 
114  FLA_Obj_free( &d );
115  FLA_Obj_free( &e );
116 #else
117  FLA_Check_error_code( FLA_EXTERNAL_LAPACK_NOT_IMPLEMENTED );
118 #endif
119 
120  return info;
121 }
FLA_Error FLA_Tridiag_check(FLA_Uplo uplo, FLA_Obj A, FLA_Obj t)
Definition: FLA_Tridiag_check.c:13
int F77_chetd2(char *uplo, int *n, scomplex *a, int *lda, float *d, float *e, scomplex *tau, int *info)
int F77_dsytd2(char *uplo, int *n, double *a, int *lda, double *d, double *e, double *tau, int *info)
int F77_zhetd2(char *uplo, int *n, dcomplex *a, int *lda, double *d, double *e, dcomplex *tau, int *info)
int F77_ssytd2(char *uplo, int *n, float *a, int *lda, float *d, float *e, float *tau, 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
FLA_Datatype FLA_Obj_datatype_proj_to_real(FLA_Obj A)
Definition: FLA_Query.c:23
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
FLA_Error FLA_Obj_free(FLA_Obj *obj)
Definition: FLA_Obj.c:588
void FLA_Param_map_flame_to_netlib_uplo(FLA_Uplo uplo, void *blas_uplo)
Definition: FLA_Param.c:47
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_chetd2(), F77_dsytd2(), F77_ssytd2(), F77_zhetd2(), FLA_Check_error_level(), FLA_Obj_col_stride(), FLA_Obj_create(), FLA_Obj_datatype(), FLA_Obj_datatype_proj_to_real(), FLA_Obj_free(), FLA_Obj_has_zero_dim(), FLA_Obj_width(), FLA_Param_map_flame_to_netlib_uplo(), and FLA_Tridiag_check().

Referenced by FLA_Tridiag_unb_ext().