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

//REAL f(const realmat& b) {b(1,1)=1.0; return b(1,1);}  //This is an error.
REAL f(const realmat& b) {realmat a=b; a(1,1)=1.0; return a(1,1);} //This is 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;
}