libflame  revision_anchor
Functions
FLA_Eig_gest_unb_external.c File Reference

(r)

Functions

FLA_Error FLA_Eig_gest_unb_external (FLA_Inv inv, FLA_Uplo uplo, FLA_Obj A, FLA_Obj B)
 
FLA_Error FLA_Eig_gest_il_unb_ext (FLA_Obj A, FLA_Obj B)
 
FLA_Error FLA_Eig_gest_iu_unb_ext (FLA_Obj A, FLA_Obj B)
 
FLA_Error FLA_Eig_gest_nl_unb_ext (FLA_Obj A, FLA_Obj B)
 
FLA_Error FLA_Eig_gest_nu_unb_ext (FLA_Obj A, FLA_Obj B)
 

Function Documentation

◆ FLA_Eig_gest_il_unb_ext()

FLA_Error FLA_Eig_gest_il_unb_ext ( FLA_Obj  A,
FLA_Obj  B 
)

References FLA_Eig_gest_unb_external().

Referenced by FLA_Eig_gest_il().

115 {
116  return FLA_Eig_gest_unb_external( FLA_INVERSE, FLA_LOWER_TRIANGULAR, A, B );
117 }
FLA_Error FLA_Eig_gest_unb_external(FLA_Inv inv, FLA_Uplo uplo, FLA_Obj A, FLA_Obj B)
Definition: FLA_Eig_gest_unb_external.c:13

◆ FLA_Eig_gest_iu_unb_ext()

FLA_Error FLA_Eig_gest_iu_unb_ext ( FLA_Obj  A,
FLA_Obj  B 
)

References FLA_Eig_gest_unb_external().

Referenced by FLA_Eig_gest_iu().

120 {
121  return FLA_Eig_gest_unb_external( FLA_INVERSE, FLA_UPPER_TRIANGULAR, A, B );
122 }
FLA_Error FLA_Eig_gest_unb_external(FLA_Inv inv, FLA_Uplo uplo, FLA_Obj A, FLA_Obj B)
Definition: FLA_Eig_gest_unb_external.c:13

◆ FLA_Eig_gest_nl_unb_ext()

FLA_Error FLA_Eig_gest_nl_unb_ext ( FLA_Obj  A,
FLA_Obj  B 
)

References FLA_Eig_gest_unb_external().

Referenced by FLA_Eig_gest_nl().

125 {
126  return FLA_Eig_gest_unb_external( FLA_NO_INVERSE, FLA_LOWER_TRIANGULAR, A, B );
127 }
FLA_Error FLA_Eig_gest_unb_external(FLA_Inv inv, FLA_Uplo uplo, FLA_Obj A, FLA_Obj B)
Definition: FLA_Eig_gest_unb_external.c:13

◆ FLA_Eig_gest_nu_unb_ext()

FLA_Error FLA_Eig_gest_nu_unb_ext ( FLA_Obj  A,
FLA_Obj  B 
)

References FLA_Eig_gest_unb_external().

Referenced by FLA_Eig_gest_nu().

130 {
131  return FLA_Eig_gest_unb_external( FLA_NO_INVERSE, FLA_UPPER_TRIANGULAR, A, B );
132 }
FLA_Error FLA_Eig_gest_unb_external(FLA_Inv inv, FLA_Uplo uplo, FLA_Obj A, FLA_Obj B)
Definition: FLA_Eig_gest_unb_external.c:13

◆ FLA_Eig_gest_unb_external()

FLA_Error FLA_Eig_gest_unb_external ( FLA_Inv  inv,
FLA_Uplo  uplo,
FLA_Obj  A,
FLA_Obj  B 
)

References F77_chegs2(), F77_dsygs2(), F77_ssygs2(), F77_zhegs2(), FLA_Check_error_level(), FLA_Eig_gest_check(), FLA_Obj_col_stride(), FLA_Obj_datatype(), FLA_Obj_length(), and FLA_Param_map_flame_to_netlib_uplo().

Referenced by FLA_Eig_gest_il_unb_ext(), FLA_Eig_gest_iu_unb_ext(), FLA_Eig_gest_nl_unb_ext(), and FLA_Eig_gest_nu_unb_ext().

14 {
15  FLA_Error r_val = FLA_SUCCESS;
16 #ifdef FLA_ENABLE_EXTERNAL_LAPACK_INTERFACES
17  int itype;
18  int info;
19  FLA_Datatype datatype;
20  int m_A, cs_A;
21  int cs_B;
22  char blas_uplo;
23 
24  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
25  FLA_Eig_gest_check( inv, uplo, A, B );
26 
27 // if ( FLA_Obj_has_zero_dim( A ) ) return FLA_SUCCESS;
28 
29  if ( inv == FLA_INVERSE )
30  itype = 1;
31  else
32  itype = 2;
33 
34  datatype = FLA_Obj_datatype( A );
35 
36  m_A = FLA_Obj_length( A );
37  cs_A = FLA_Obj_col_stride( A );
38  cs_B = FLA_Obj_col_stride( B );
39 
40  FLA_Param_map_flame_to_netlib_uplo( uplo, &blas_uplo );
41 
42 
43  switch( datatype ){
44 
45  case FLA_FLOAT:
46  {
47  float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
48  float *buff_B = ( float * ) FLA_FLOAT_PTR( B );
49 
50  F77_ssygs2( &itype,
51  &blas_uplo,
52  &m_A,
53  buff_A, &cs_A,
54  buff_B, &cs_B,
55  &info );
56 
57  break;
58  }
59 
60  case FLA_DOUBLE:
61  {
62  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
63  double *buff_B = ( double * ) FLA_DOUBLE_PTR( B );
64 
65  F77_dsygs2( &itype,
66  &blas_uplo,
67  &m_A,
68  buff_A, &cs_A,
69  buff_B, &cs_B,
70  &info );
71 
72  break;
73  }
74 
75  case FLA_COMPLEX:
76  {
77  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
78  scomplex *buff_B = ( scomplex * ) FLA_COMPLEX_PTR( B );
79 
80  F77_chegs2( &itype,
81  &blas_uplo,
82  &m_A,
83  buff_A, &cs_A,
84  buff_B, &cs_B,
85  &info );
86 
87  break;
88  }
89 
90  case FLA_DOUBLE_COMPLEX:
91  {
92  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
93  dcomplex *buff_B = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( B );
94 
95  F77_zhegs2( &itype,
96  &blas_uplo,
97  &m_A,
98  buff_A, &cs_A,
99  buff_B, &cs_B,
100  &info );
101 
102  break;
103  }
104 
105  }
106 
107 #else
108  FLA_Check_error_code( FLA_EXTERNAL_LAPACK_NOT_IMPLEMENTED );
109 #endif
110 
111  return r_val;
112 }
void FLA_Param_map_flame_to_netlib_uplo(FLA_Uplo uplo, void *blas_uplo)
Definition: FLA_Param.c:47
int FLA_Error
Definition: FLA_type_defs.h:47
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
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
int F77_zhegs2(int *itype, char *uplo, int *n, dcomplex *a, int *lda, dcomplex *b, int *ldb, int *info)
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
Definition: blis_type_defs.h:137
FLA_Error FLA_Eig_gest_check(FLA_Inv inv, FLA_Uplo uplo, FLA_Obj A, FLA_Obj B)
Definition: FLA_Eig_gest_check.c:13
int F77_ssygs2(int *itype, char *uplo, int *n, float *a, int *lda, float *b, int *ldb, int *info)
int F77_dsygs2(int *itype, char *uplo, int *n, double *a, int *lda, double *b, int *ldb, int *info)
int F77_chegs2(int *itype, char *uplo, int *n, scomplex *a, int *lda, scomplex *b, int *ldb, int *info)