arcprob.cpp

Go to the documentation of this file.

#include <cmath>
#include <fstream>
using namespace std;

#include <arcprob.h>


arcprob::arcprob( stringc & filename )
  : W(0), angle0(0.0), ac(0)
{
  reader(filename);
}

arcprob::~arcprob()
{
  delete ac;
}

void arcprob::reader( stringc & filename )
{
  ifstream targ(filename.c_str());
  assert(targ.good()==true);
  if (targ.good()==false)
    return;

  targ >> angle0;
  angle0 *= degtorad;
  uint N;
  targ >> N;
  targ >> W;

cout << SHOW(angle0) << endl;
cout << SHOW(N) << endl;
cout << SHOW(W) << endl;

  assert(W>N);
  assert(N>1);
  
  qi.resize(W);
  for (uint i=0; i<W; ++i)
    targ >> qi[i];

  if (ac!=0)
    delete ac;

  ac = new arcsconnected(5,N);

for (uint i=0; i<qi.size(); ++i)
  cout << SHOW(qi[i]) << endl;
    
  solveInitialCondition();
}

void arcprob::solveInitialCondition()
{
  assert(ac!=0);
  if (ac==0)
    return;

  uintc N(ac->N);

  // Let the arcs initially pass through the
  // control points.
  vector< pt2 > u(N+1);
  doublec dx = (double)(W-1) / (double)N;
  uint temp;
  for (uint i=1; i<N; ++i)
  {
    temp = (uint)(floor(dx*i));
    u[i] = qi[temp];
  }
  u[0] = qi[0];
  u[N] = qi[W-1];

  ac->constructPhi0(angle0,u);
}