file2mat.cc

#include "cuttraj.h"
#include <iostream>
#include <armadillo>

using namespace arma;
using namespace std;


mat file2mat(const string& filename, mat& y) {
	mat yTemp;
	yTemp.load(filename, raw_ascii);

	
	y = zeros<mat>(yTemp.n_rows, yTemp.n_cols);

	
	double dTemp;
	for (int i = 0; i < y.n_rows; ++i) {
		for (int j = 0; j < y.n_cols; ++j) {
			dTemp = yTemp(i,j);
			y(i,j) = dTemp;
		}
	}
	cout << y.n_rows << endl;
	cout << y.n_cols << endl;
	
	y = ones<mat>(yTemp.n_rows, yTemp.n_cols);
	
}