320_LinearRegression/main.cpp at master · nilay036/320_LinearRegression · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
#include "Matrix.h"
#include <iostream>
#include <fstream>
#include <sstream>
#include <string>
#include <chrono>
using namespace std;

string time(string operation, Matrix<int>& a, Matrix<int>& b, int r, int c)
{
    // The "auto" type determines the correct type at compile-time
	Matrix<int> timeMe;
    	auto start = std::chrono::system_clock::now();
	if(operation == "+")
    		timeMe = a + b;

	if(operation == "-")
		timeMe = a - b;

	if(operation == "*")
		timeMe = a * b;

    auto end = std::chrono::system_clock::now();
    std::chrono::duration<double> elapsed_seconds = end-start;
    std::time_t end_time = std::chrono::system_clock::to_time_t(end);

    cout << "finished " << operation << " ing " << r << " #long matrixs" << endl
	 <<" at " << std::ctime(&end_time)
         << "elapsed time: " << elapsed_seconds.count() << " s" << endl;
	return timeMe.printMatrix();
}

int main1()
{
	//I was getting invalid free/delete from c++ trying to deconstruct x and y before they were pointers
	Matrix<int> x;// =  new Matrix<int>(size, size);
	Matrix<int> y;// = new Matrix<int>(size, size);

	string fills[] = {"fill", "diagonal", "upper-triangle", "lower-triangle", "random"};
	int numFills = 4;

	try{//multiplication timing

		for(int size = 2; size < 16; size*=2)
		{
			x = Matrix<int>(size, size);
			y = Matrix<int>(size, size);

			for(int f = 0; f< numFills; f++){
				x.fillMatrix(fills[f], 4);
				y.fillMatrix(fills[f], 4);
				cout << fills[f] << " matricies w/ size = " << size << endl;
				cout << time("*", x, y, size, size);
				cout << endl;
			}
		}
	}
	catch(string e)
	{
		cout << "multiplication error" << e << endl;
	}
	try{//addition timing

		for(int size = 2; size < 16; size*=2)
		{
			x = Matrix<int>(size, size);
			y = Matrix<int>(size, size);

			for(int f = 0; f< numFills; f++){
				x.fillMatrix(fills[f], 4);
				y.fillMatrix(fills[f], 4);
				cout << fills[f] << " matricies w/ size = " << size << endl;
				cout << time("+", x, y, size, size);
				cout << endl;
			}
		}
	}
	catch(string e)
	{
		cout << "addition error" << e << endl;
	}

	return 0;
}
Matrix<double> readFile(char* fileName)
{
	string line;
	int numLines = 0;
	Matrix<double> pointsMatrix;
	ifstream pointsFile (fileName);
	if(pointsFile.is_open()){
		while(getline(pointsFile, line))
			numLines++;
		pointsFile.clear();
		pointsFile.seekg(0, ios::beg);
		cout << numLines;
		pointsMatrix = Matrix<double>(numLines, 2);
		bool isX = true;
		double temp = 0;
		int i = 0;
		while(pointsFile >> temp && i < numLines)
		{
			if(isX)
			{
				pointsMatrix[i][0] = temp;
				isX = false;
			}
			else
			{
				pointsMatrix[i][1] = temp;
				isX = true;
			}
			i++;
		}
	pointsFile.close();
	}
	return pointsMatrix;

}
int main(int argc, char** argv)
{
	if(argc != 2)
	{
		cout << "this program only takes 1 parameter";
		return 0;
	}
	Matrix<double> pointsMatrix = readFile(argv[1]);
	cout << pointsMatrix;
}