//

// Pointer class implementation

//

#include "point.h"

#include <iostream>

using namespace std;

// Constructor

Point::Point(int dimen) {

// Could not accept dimension less then 2, don't return an error, but set default

if (dimen < 2)

dimen= 2;

coord.resize(dimen);

for (int i = 0; i < dimen; ++i)

coord[i] = 0.0;

} // Constructor

Point::Point(vector<double> p_coord)

{

coord = p_coord;

size_t dimen = coord.size();

if (dimen < 1)

{

coord.resize(2);

if (dimen == 0)

coord[0] = 0.0;

coord[1] = 0.0;

}

}

// Function to calculate distance to other point

double Point::dist(Point other) {

double sum = 0;

for (int i = 0; i < coord.size(); ++i)

sum += pow(coord[i] - other.coord[i], 2);

return sqrt(sum);

}

Point operator+ (Point& p1, Point& p2) {

int majorDim = (p1.coord.size() > p2.coord.size() ? p1.coord.size() : p2.coord.size());

Point sum(majorDim); // Create a variable coordinates 0s with greater dimension

for (int i=0; i < majorDim; i++) {

if ((i+1) <= p1.coord.size())

sum.coord[i] += p1.coord[i];

if ((i+1) <= p2.coord.size())

sum.coord[i] += p2.coord[i];

}

return sum;

}

Point operator* (double m1, Point& p1) {

Point prod(p1.coord.size()); // Create a variable coordinates 0s

for (int i=0; i < p1.coord.size(); i++)

prod.coord[i] = m1 * p1.coord[i];

return prod;

}

ostream& operator<< (ostream&out, Point& p) {

out << "(";

for (int i = 0; i < p.coord.size() - 1; ++i)

out << p.coord[i] << ", ";

out << p.coord[p.coord.size() - 1] << ")";

return out;

}