#include <stdlib.h>

#include <iostream>
#include <fstream>
using namespace std;

#include <fann.h>

int main(int narg, char *argv[])
{
	fann_type *calc_out;
	fann_type input[2];

	struct fann *ann = fann_create_from_file("mult.net");

	ofstream ofile("test_results.dat");
	ofile<<"# input1 input2 output"<<endl;

	srandom(1);
	for(int i=0; i<4000; i++){
		// Get 2 integer values between 0 and 3 inclusive
		double input1 = 3.0*(double)random()/(double)(0x7FFFFFFF);
		double input2 = 3.0*(double)random()/(double)(0x7FFFFFFF);

		// inputs need to be scaled to be between 0 and 1, but because the ANN
		// has harder time with near zero values, network inputs are shifted
		// bu 1.
		input[0] = (input1+1.0)/10.0;
		input[1] = (input2+1.0)/10.0;
		
		
		// Run through ANN
		calc_out = fann_run(ann, input);
		
		// Output needs to be scaled back up
		double output = calc_out[0]*10.0 - 1.0;
		
		char str[256];
		sprintf(str,"%3.2f", 100.0*(output-input1*input2)/(input1*input2));
	
		cout<<input1<<" * "<<input2<<" = "<<output<<" (actual="<<input1*input2<<")";
		//if(output!=(int)(input1*input2))cout<<" <====";
		cout<<"[err="<<str<<"%]";
		cout<<endl;
		
		ofile<<input1<<" "<<input2<<" "<<output<<endl;
	}
	
	ofile.close();

	fann_destroy(ann);

	return 0;
}