Hi,
I would like to put in a feature request for MKL:
Currently, the c interface of mkl does not provide for const pointers. It only does for pointers to const ( i.e. double const * and not double const * const ) - the latter meaning that the pointer itself is not changing as well and not only what it points to (through it, to be exact).
This results in a lot of const (un)casting when one tries to write wrapper classes, since any sensible c++ programmer would provide with 2 versions of the (example, schematically) container data:
double * data() { return ptr_ ; }
double const * const data() const { return ptr_ ;}
This creates a logical havoc when writing lapack/blas wrappers. Would you please consider correcting it ?