void save_matrix( const char * filename, const vnl_matrix<double>& x) { if (strlen(filename)>0) { std::ofstream outfile(filename,std::ios_base::out); x.print(outfile); } }