origin.h

Go to the documentation of this file.

//
// Created by shuttle3468 on 8/2/17.
//

#ifndef COMPUTGEOALGOS_ORIGIN_H
#define COMPUTGEOALGOS_ORIGIN_H

#include <iostream>
#include <bits/stdc++.h>

#define ull unsigned long long 
#define pb push_back 
#define mp make_pair 
#define COLLINEAR 1 
#define CLOCKWISE 2 
#define ANTICLOCKWISE 3 
using namespace std;
vector<int> indices;

double calcEuclidDistance(pair<double, double> fpoint, pair<double, double> spoint) {
    double dist;
    dist = sqrt(pow((fpoint.first-spoint.first),2.0) + pow((fpoint.second-spoint.second),2.0));
    return dist;
}

int orientation(pair<double, double> a, pair<double, double> b, pair<double, double> c) {
    double m1, m2, dif;
    dif = (b.second - a.second)*(c.first-b.first) - (b.first-a.first)* (c.second - b.second);
    cout << "Difference in Slopes " << dif<<endl;
    if(dif == 0) {
        cout << "COLLINEAR\n";
        return COLLINEAR;
    } else if (dif > 0 ) {
        cout << "CLOCKWISE\n";
        return CLOCKWISE;
    } else {
        cout << "ANTICLOCKWISE\n";
        return ANTICLOCKWISE;
    }
}

bool orderedSort(pair<double,double> &f, pair<double,double> &s) {
    //cout << "Comparing (" << f.first << ", " << f.second <<")" << " (" << s.first << ", " << s.second <<")" << endl;
    if(f.first < s.first) {
        return true;
    }
    if (f.first > s.first) {
        return false;
    }
    return f.second<s.second;
}

bool orderedYSort(pair<double,double> &f, pair<double,double> &s) {
    if(f.second < s.second) {
        return true;
    }
    if (f.second > s.second) {
        return false;
    }
    return f.first<s.first;
}


pair<double ,double > P0;


bool orderByPolar(pair<double,double> &p1, pair<double,double> &p2){
    int oValue  = orientation(P0,p2,p1);
    cout << "comparing " << "(" << P0.first << ", " << P0.second <<")" <<  " "
         << "(" << p1.first << ", " << p1.second <<")" << " "
         << "(" << p2.first << ", " << p2.second <<")" <<endl;
    cout << " The orientation value is " << oValue << endl;
    if(oValue == 1) {
        cout << "Result :" << (calcEuclidDistance(P0,p2) >= calcEuclidDistance(P0,p1)) <<endl;
        return (calcEuclidDistance(P0,p2) < calcEuclidDistance(P0,p1));
    }
    return (oValue != 3);
}


void printVectorData(int len, vector<pair<double, double> > v, string s) {
    int i;
    cout << s << v.size() << endl;
    for(i=0;i<len;i++)
        cout << v[i].first <<" " << v[i].second <<endl;
}

vector<pair<double, double> > getData(char *filename) {
    vector<pair<double, double> > input;
    double a,b;
    ifstream in_file;
    in_file.open(filename);
    while (in_file.is_open()) {
        int n;
        in_file >> n;
        while(in_file >> a >> b) {
            input.pb(mp(a,b));
        }
        in_file.close();
    }
    return input;
}
#endif //COMPUTGEOALGOS_ORIGIN_H