inline std::ptrdiff_t sbev( const char jobz, const UpLo, const fortran_int_t n,
const fortran_int_t kd, float* ab, const fortran_int_t ldab, float* w,
float* z, const fortran_int_t ldz, float* work ) {
fortran_int_t info(0);
LAPACK_SSBEV( &jobz, &lapack_option< UpLo >::value, &n, &kd, ab, &ldab, w,
z, &ldz, work, &info );
return info;
}
inline
void hbev (char const jobz, char const uplo, int const n, int const kd,
float* ab, int const ldab, float* w, float* z, int const ldz,
float* work, int& info)
{
//for (int i=0; i<n*kd; ++i) std::cout << *(ab+i) << " " ;
//std::cout << "\n" ;
LAPACK_SSBEV (&jobz, &uplo, &n, &kd, ab, &ldab, w, z, &ldz,
work, &info);
}
inline void sbev( char const jobz, char const uplo, integer_t const n,
integer_t const kd, float* ab, integer_t const ldab, float* w,
float* z, integer_t const ldz, float* work, integer_t& info ) {
LAPACK_SSBEV( &jobz, &uplo, &n, &kd, ab, &ldab, w, z, &ldz, work,
&info );
}
