inline std::string aligned_write(const MatExpr &a) {
   //boost::io::ios_flags_saver coutstate(stream);
   std::ostringstream oss;
   oss << std::scientific;
   for (unsigned int i=0; i<a.size1(); ++i) {
     for (unsigned int j=0; j<a.size2(); ++j)
       oss << a(i,j) << "\t";
     if (i < a.size1()-1)
       oss << std::endl;
   }
   return oss.str();
 }