#include "libscl.h"
using namespace scl;  
using namespace std;

//REAL f(const realmat* bptr) {realmat* a=bptr; return (*a)(1,1);}   //Error.
//REAL f(const realmat* bptr) {const realmat* a=bptr; return (*a)(1,1);} //Ok.
REAL f(const realmat* bptr) {realmat a=*bptr; a(1,1)=1; return a(1,1);} //Ok.

int main(int argc, char** argp, char** envp)
{
  realmat a(5,5,0.0);
  REAL r = f(&a);
  cout << "a(1,1) = " << a(1,1) << '\n';
  cout << "r = " << r  << '\n';
  return 0;
}