#include <../BLAS/c_interface_base.h>
#include <complex>

extern "C" {
#include "../BLAS/blas.h"

//void sgesv_(int*, int*, float *, int*, int*, float *, int*, int*);
//void dgesv_(int*, int*, double*, int*, int*, double*, int*, int*);

void sgels_(char*, int*, int*, int*, float *, int*, float *, int*, float *, int*, int*);
void dgels_(char*, int*, int*, int*, double*, int*, double*, int*, double*, int*, int*);

//void sgetrf_(int*, int*, float *, int*, int*, int*);
//void dgetrf_(int*, int*, double*, int*, int*, int*);

//void spotrf_(char*, int*, float *, int*, int*);
//void dpotrf_(char*, int*, double*, int*, int*);

void ssyev_(char*, char*, int*, float *, int*, float *, float *, int*, int*);
void dsyev_(char*, char*, int*, double*, int*, double*, double*, int*, int*);
}


#define MAKE_STRING2(S) #S
#define MAKE_STRING(S) MAKE_STRING2(S)

#define CAT2(A,B) A##B
#define CAT(A,B) CAT2(A,B)

template <typename real> class lapack_interface;


static char notrans = 'N';
static char trans = 'T';
static char nonunit = 'N';
static char lower = 'L';
static char right = 'R';
static char left = 'L';
static int intone = 1;
static int zeroint = 0;


#define SCALAR        float
#define SCALAR_PREFIX s
#include "lapack_interface_impl.hh"
#undef SCALAR
#undef SCALAR_PREFIX

#define SCALAR        double
#define SCALAR_PREFIX d
#include "lapack_interface_impl.hh"
#undef SCALAR
#undef SCALAR_PREFIX