#include <iostream>
#include <octave/oct.h>
#include<octave/octave.h>
#include<octave/parse.h>
#include <octave/toplev.h>
#include <octave/ov-struct.h>
#include <octave/Cell.h>

using namespace std;

int
main(void)
{
	vector<string> file_names;
	
	string_vector argv(2);
	argv(0)="embedded";
	argv(1)="-q";
	octave_main(2,argv.c_str_vec(),1);

	octave_value_list dir_stuff;
	dir_stuff(0) = "*.saeC";
	octave_value_list dir_return = feval("dir", dir_stuff, 1); //calls dir function returns a octave_map
	Octave_map my_map = dir_return(0).map_value(); //returns the map_value to a map instance
	Octave_map::iterator p1 = my_map.seek("name"); //looks for the the key "name"
	
	for(int i = 0; i < my_map.contents(p1).numel(); i++) //loop over the times where "name" is found
	{
		//my_map.contents(p1) returns a cell array with the value
		//my_map.contents(p1)(i) returns the ith cell array
		//my_map.contents(p1)(i).char_matrix_value() returns an octave::charMatrix
		//my_map.contents(p1)(i).char_matrix_value().row_as_string(0)  returns a std::string
		cout<<my_map.contents(p1)(i).char_matrix_value().row_as_string(0)<<endl; 
		file_names.push_back(my_map.contents(p1)(i).char_matrix_value().row_as_string(0)); 
		cout<<"--------------------"<<endl;
	}
	
	for(int i = 0; i < file_names.size(); i++)
	{
		cout<<file_names[i]<<endl;
		cout<<"XXXXXXXXXXXXXXXXXXXXXXXXXX"<<endl;
	}
		
	cout<<"!!!!!!!!!!!!!!!!!!!!!!!!!!!!"<<endl;
	cout<<endl;
	//read data file and store in matrix
	
	octave_value_list dlm_in;
	octave_value_list grid_in;
	octave_value_list datafile;
	Matrix read_data;
	ColumnVector Mdot;
	ColumnVector PR;
	ColumnVector Eff;
	octave_value_list out;
	
	for(int i = 0; i < file_names.size(); i++)
	{
		cout<<file_names[i]<<endl;
		
	
		dlm_in(0) = file_names[i]; 
		dlm_in(1) = ""; 
		dlm_in(2) = 7;
		dlm_in(3) = 0;
		
		datafile = feval("dlmread",dlm_in,1);
		
		read_data = datafile(0).matrix_value();
			
		//get the columns needed for griddata
		Mdot = read_data.column(2);
		PR = read_data.column(3);
		Eff = read_data.column(4);
		
		grid_in(0) = octave_value(Mdot);
		grid_in(1) = octave_value(PR);
		grid_in(2) = octave_value(Eff);
		grid_in(3) = octave_value(.5);
		grid_in(4) = octave_value(3);
		
		out=feval("griddata",grid_in,1);
		
		cout<<out(0).double_value()<<endl;
		cout<<"XXXXXXXXXXXXXXXXXXXXXXXXXX"<<endl;
	}
	cout<<"Here we go"<<endl;
    do_octave_atexit(); //quit the interpreter
    cout<<"there we went"<<endl;
    return 0;
}