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_Tridiag_unb_external
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_Tridiag_check
FLA_Error FLA_Tridiag_check(FLA_Uplo uplo, FLA_Obj A, FLA_Obj t)
Definition FLA_Tridiag_check.c:13
F77_chetd2
int F77_chetd2(char *uplo, int *n, scomplex *a, int *lda, float *d, float *e, scomplex *tau, int *info)
F77_dsytd2
int F77_dsytd2(char *uplo, int *n, double *a, int *lda, double *d, double *e, double *tau, int *info)
F77_zhetd2
int F77_zhetd2(char *uplo, int *n, dcomplex *a, int *lda, double *d, double *e, dcomplex *tau, int *info)
F77_ssytd2
int F77_ssytd2(char *uplo, int *n, float *a, int *lda, float *d, float *e, float *tau, int *info)
FLA_Obj_width
dim_t FLA_Obj_width(FLA_Obj obj)
Definition FLA_Query.c:123
FLA_Obj_create
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_Obj_has_zero_dim
FLA_Bool FLA_Obj_has_zero_dim(FLA_Obj A)
Definition FLA_Query.c:400
FLA_Obj_datatype_proj_to_real
FLA_Datatype FLA_Obj_datatype_proj_to_real(FLA_Obj A)
Definition FLA_Query.c:23
FLA_Obj_col_stride
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition FLA_Query.c:174
FLA_Check_error_level
unsigned int FLA_Check_error_level(void)
Definition FLA_Check.c:18
FLA_Obj_free
FLA_Error FLA_Obj_free(FLA_Obj *obj)
Definition FLA_Obj.c:588
FLA_Param_map_flame_to_netlib_uplo
void FLA_Param_map_flame_to_netlib_uplo(FLA_Uplo uplo, void *blas_uplo)
Definition FLA_Param.c:47
FLA_Obj_datatype
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition FLA_Query.c:13
FLA_Datatype
int FLA_Datatype
Definition FLA_type_defs.h:49
FLA_Obj_view
Definition FLA_type_defs.h:159
dcomplex
Definition blis_type_defs.h:138
scomplex
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().

Generated on Thu Sep 23 2021 17:46:17 for libflame by doxygen 1.12.0