#include <iostream>
#include <fstream>
#include <string>
#include <sstream>
#include <math.h> 
#include <iostream>
#include <vector>
using namespace std;
#include "function.h"


	//part 1
	void read(string fileName)
	{
		ifstream myfile;
		myfile.open(fileName);
		string line;
		while(getline(myfile,line))
		{
			cout <<line <<endl;
		}
		myfile.close();
	}

	vector<double> observation;
	vector<double> modelOutput;
	//part 2
	void readAndStore(string fileName)
	{
		

		ifstream myfile;
		myfile.open(fileName);
		string line;

		while(getline(myfile,line))
		{
			stringstream totalSString(line);
			string token;
			getline( totalSString, token, '\t' );
			observation.push_back(atof(token.c_str()));
			getline( totalSString, token, '\t' );
			modelOutput.push_back(atof(token.c_str()));
		
		
		}
		myfile.close();
	}
	//part 3
	void print()
	{
		for(int j=0;j<observation.size();j=j+1)
		{
			cout<<observation[j]<<endl;
		}
		for(int j=0;j<modelOutput.size();j=j+1)
		{
			cout<<modelOutput[j]<<endl;
		}
	}
	void substractAndWrie(string output)
	{
		for(int j=0;j<observation.size();j=j+1)
		{
			observation[j]=observation[j]-0.2;
		}
		ofstream outputfile;
		outputfile.open	(output );
		for(int j=0;j<observation.size();j=j+1)
		{
			outputfile <<observation[j]<<'\t'<<modelOutput[j]<<endl;

		}
		outputfile.close();
	}
	//5
	void readNewInfo(string fileName)
	{
	ifstream myfile;
	myfile.open(fileName);

		int i=0;
		string line;
		while(getline(myfile,line))
		{

			stringstream totalSString(line);
			string token;
			getline( totalSString, token, '\t' );
			observation[i]=atof(token.c_str());
			getline( totalSString, token, '\t' );
			modelOutput[i]=atof(token.c_str());
			i++;
		
		}
		myfile.close();
		
	}
	//MSE
	double MSE()
	{
		double result=0;

		for(int j=0;j<modelOutput.size();j=j+1)
		{
			result=result+pow(modelOutput[j]-observation[j],2)/modelOutput.size();


		}
			double finalResult=sqrt(result);
			return finalResult;

	}


	
	//7
	void avarageTemp(string input,string ouput)
	{
	ifstream myfile;
	myfile.open(input);
	int oldMonth=1;
    int oldDay = 0;
	double tmp=0;
	ofstream outputfile;
	outputfile.open(ouput);
	string line;
	while(getline(myfile,line))
	{
		cout<<line<<endl;
		stringstream totalSString(line);
		int value;
		totalSString>> value;
		totalSString>> value;
		int month=value;
		totalSString>> value;
		int days=value;
		double mesurement;
		totalSString>> mesurement;
	

		if (oldMonth == month)
                    {
                        tmp += mesurement;
                        oldDay = days;
                        oldMonth = month;

                    }
                    else
                    {
                        double avarege = tmp / oldDay;
 
                        outputfile<<avarege<<endl;;
                        tmp = mesurement;
                        oldDay = days;
                        oldMonth = month;
                    }
		outputfile.close();

	
	}
	}

	double calculateR(string intput)
	{
		ifstream myfile;
		myfile.open(intput);

		//8

		vector<double> observations;
		vector<double> modelOutputs;
		string line;
		while(getline(myfile,line))
		{

			stringstream totalSString(line);
			double value;
			totalSString>> value;
			observations.push_back(value);
			totalSString>> value;
			modelOutputs.push_back(value);
	
		
		}
		myfile.close();
		double firstSum=0;
		for(int j =0;j<observations.size();j=j+1)
		{
			firstSum=firstSum+observations[j]-modelOutputs[j];
		}
		double SumOfObs=0;
		double SumOfMessured=0;
		for(int j =0;j<observations.size();j=j+1)
		{
			SumOfObs=SumOfObs+observations[j];
			SumOfMessured=SumOfMessured+modelOutputs[j];
		}
		double SumOfObsPower=0;
		double SumOfMessuredPower=0;
		for(int j =0;j<observations.size();j=j+1)
		{
			SumOfObsPower=SumOfObsPower+pow(observations[j],2);
			SumOfMessuredPower=SumOfMessuredPower+pow(modelOutputs[j],2);
		}
		double SumOfObsPowerLast=0;
		double SumOfMessuredPowerLast=0;
		for(int j =0;j<observations.size();j=j+1)
		{
			SumOfObsPower=SumOfObsPower+pow(observations[j],2)/observations.size();
			SumOfMessuredPower=SumOfMessuredPower+pow(modelOutputs[j],2)/observations.size();
		}
		double result=(firstSum*((SumOfObs*SumOfMessured)/observations.size()))/sqrt((SumOfObsPower-SumOfObsPowerLast)*(SumOfMessuredPower-SumOfMessuredPower));
		return result;
	}