main.cpp

Go to the documentation of this file.

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

// TODO clean this up! 

#include <GL/glut.h>
#include <GL/gl.h>

#include <arc.h>
#include <arcdraw.h>
#include <arcsconnected.h>
#include <arcsdef.h>
#include <arcprob.h>
#include <commandline.h>
#include <farcmin01.h>
#include <gobj.h>
#include <graphmisc.h>
#include <openwindowresource.h>
#include <pathlinesegtest.h>
#include <zpr.h>

// Instanciated in main as pathlineseg.cpp needs it but cannot instantiate it.
template<>
double zero<double>::val = 1E-10;


void print(vector<arc> const & v)
{
  cout << endl << endl;
  cout << "Printing the arcs" << endl;
  for (uint i=0; i<v.size(); ++i)
    v[i].print();
}


void writearcgeometry(int argc, char** & argv)
{
  commandline cmd(argc,argv);
  string filearcs = "pa03.txt";
  cmd.mapvar(filearcs,"filearcs");

  if (filearcs.empty())
  {
    cout << "error:  no filearcs=<filename> command";
    cout << " line argument" << endl;
    return;
  }

  vector<arc> vi;
  arc::arcreader(vi,filearcs.c_str());

  for (uint i=0; i<vi.size(); ++i)
  {
    arcdraw *ad = new arcdraw(vi[i]);
    gobj::global->push(ad);
  }

  print(vi);

}





void prog01(int argc, char** & argv)
{

  openwindowresource window(argc,argv);

  writearcgeometry(argc,argv);


}

void test01()
{
  cout << "Looking at the atan2 function" << endl;
  cout << "Really asking the question that given a point" << endl;
  cout << "on a unit circle what is its angle." << endl;
  cout << endl;
  double phi;
  cout << "Enter phi(degrees): ";
  cin >> phi;
  phi *= PI/180.0;

  pt2 p0(cos(phi),sin(phi));
  p0.normalize();
  cout << p0 << endl;
  double theta;
  theta = atan2(p0.y,p0.x);
  cout << theta << endl;
  cout << theta*180.0/PI << endl;
}


void test02()
{
  arc a(1.0,pt2(0.0,1.0),pt2(1.0,0.0));
  cout << "center=" << a.center << endl;
}

void test03()
{
  vector<arc> v;

  arc::arcreader(v,"pa03.txt");

  for (uint i=0; i<v.size(); ++i)
  {
    cout << v[i].p0 << " " << v[i].p1 << " ";
    cout << v[i].radius << endl;  
  }

}

void test04()
{
  vector<pt2> pts;
  pts.push_back( pt2(0.0,0.0) );
  pts.push_back( pt2(1.0,1.1) );
  pts.push_back( pt2(1.5,2.8) );


  arc a1(2.5,pts[0],pts[1]);

  cout << "a1" << endl;
  a1.print();

  double theta0 = a1.phi1 + PI;
  if (theta0 > 2.0*PI)
    theta0 -= 2.0*PI;

cout << SHOW(a1.phi1*radtodeg) << endl;
cout << SHOW(theta0*radtodeg) << endl;

  arc a2;
  a2.constructPhi0TwoPoints(theta0,pts[1],pts[2]);

  //cout << "a2" << endl;
  //a2.print();

}

void test05(int argc, char** & argv)
{
  cout << "Constructing connected arcs" << endl;
  openwindowresource window(argc,argv);

  vector<pt2> pts;
  pts.push_back( pt2(0.0,0.0) );
  pts.push_back( pt2(1.0,1.1) );
  pts.push_back( pt2(1.5,2.8) );
  pts.push_back( pt2(2.8,4.2) );

  pts.push_back( pt2(4.8,4.0) );
  pts.push_back( pt2(6.8,4.2) );
  pts.push_back( pt2(8.1,5.5) );


  arcsconnected ac(3,pts.size()-1);

  ac.constructR0(2.5,pts);

  uintc N = ac.vi.size();
  vector<arc> & vi(ac.vi);


  for (uint i=0; i<N; ++i)
  {
    arcdraw *ad = new arcdraw(vi[i]);
    gobj::global->push(ad);
  }

  //print(vi);
}


void test06()
{
  cout << "Testing the distance from any point to an arc" << endl;

  pt2 p0(1.0,0.0);
  pt2 p1(0.0,1.0);

  arc a1(1.0,p0,p1);
  a1.print();

  cout << "Enter a point" << endl;
  pt2 p2(0.0,1.0);
  cin >> p2;
  double dist(0.0);
  a1.distance(dist,p2);
  cout << SHOW(dist) << endl;

}

void test07(int argc, char** & argv)
{
  string fileprob="pb01.txt";
  arcprob a1(fileprob);

  string filearcs="pa04.txt";

//cout << "arcprob constructed???" << endl;

  assert(a1.ac!=0);

  a1.ac->arcwriter("pa04.txt");

cout << endl << endl;
cout << "Look at the arc produced" << endl;

  vector<arc> & vi(a1.ac->vi);
  for (uint i=0; i<vi.size(); ++i)
    vi[i].print();
}

void displayarcsfileformat()
{
  cout << "Displaying the Arcs File Format" << endl;
  cout << endl;
  cout << "  Format for filearcs:  This is a text file containing a list" << endl;
  cout << "  of points and arcs. An arc is described by three parameters" << endl; 
  cout << "  i j radius where i and j are indexes into the list of points." << endl; 
  cout << "  The arc is on the right of line ij if the radius is +ve, " << endl;
  cout << "  left if the radius is -ve." << endl;
  cout << endl;
  cout << "  Here is an example of the filearc file." << endl;
  cout << endl;

  cout << "3" << endl;
  cout << "0.0 0.0" << endl;
  cout << "1.0 1.1" << endl;
  cout << "1.5 2.8" << endl;
  cout << "2" << endl;
  cout << "0 1 2.5" << endl;
  cout << "1 2 -2.5" << endl;
  cout << endl;
    
  cout << "  The first number is the number of points to be read in." << endl;
  cout << "  After the 2D points comes the number of arcs to be read" << endl;
  cout << "  Each arc indexes to two points and has a signed radius" << endl;
  cout << endl;
}

void test08()
{
  pt2 p0(0.0,0.0);
  //pt2 p1(1.5,2.8);
  pt2 p1(1.5,0.8);
  double theta=315.0*degtorad;

  arc a1;
  a1.constructPhi0TwoPoints(theta,p0,p1);
  a1.print();


}

int main(int argc, char** argv)
{
  commandline cmd(argc,argv);
  uint prog(0);
  cmd.mapvar(prog,"prog");

  bool help(false);
  cmd.mapvar(help,"help");
  if (help==true)
  {
    displayarcsfileformat();
  }

  switch (prog)
  {
    case 0:
      cout << "Test programs for arc workspace" << endl << endl;

      cout << "$./main prog=1 filearcs=<pa03.txt> " << endl;
      cout << "  Display 2D arcs in file filearcs." << endl;

      cout << "$./main prog=5" << endl;
      cout << "  Constructing connected arcs" << endl;
      cout << "$./main prog=6"  << endl;
      cout << "  Testing the distance from a point to an arc." << endl;
      cout << "$/.main prog=7 fileprob=<pb01.txt> filearcs=<pa04.txt>" << endl;
      cout << "  Testing initial condition solution" << endl;
      cout << "$./main prog=8" << endl;
      cout << "  Looking at initial angle and 2 points construction." << endl;
      cout << "$./main prog=30 " << pathlinesegtest::doc[1] << endl;
      cout << "$./main prog=31 " << pathlinesegtest::doc[2] << endl;
      cout << "$./main prog=32 " << pathlinesegtest::doc[3] << endl;
      cout << "$./main prog=33 " << pathlinesegtest::doc[4] << endl;


      cout << "./main help=true" << endl;
      cout << "  Display info such as the arc file formats." << endl;

      break;
  
    case 1: prog01(argc,argv); break;

    case 4: test04(); break; 
    case 5: test05(argc,argv); break;
    case 6: test06(); break;
    case 7: test07(argc,argv); break; 
    case 8: test08(); break;

    case 30: pathlinesegtest::test01(argc,argv); break;
    case 31: pathlinesegtest::test02(); break; 
    case 32: return pathlinesegtest::unittest01(argc,argv); 
    case 33: pathlinesegtest::test03(argc,argv); break;

    default: cout << "error:  No case handled." << endl; return 1;
  }


  return 0;
}