testOnCurve.cc

#include "modelGen2d.h"
#include "curveReader.h"

int main(int argc, char **argv) {
  const int numExpectedArgs = 5;
  if (argc != numExpectedArgs) {
    std::cerr << "Usage: <input csv file> <ptIdx> <isOnCurve> <onCurveAngleTol>\n";
    std::cerr << "input csv file with the following columns: "
                 "x,y,z,isOnCurve,angle,isMdlVtx\n";
    std::cerr << "ptIdx: the point to check\n";
    std::cerr << "isOnCurve: the expected result for the point being checked\n";
    std::cerr << "onCurveAngleTol: the angle being passed to onCurve(...)\n";
    return 1;
  }
  assert(argc == numExpectedArgs);

  std::string filename = argv[1];
  const auto ptIdx = std::atoi(argv[2]);
  const auto isOnCurve = std::atoi(argv[3]);
  const auto onCurveAngleTol = std::atof(argv[4]);
  std::cout << "input csv file: " << argv[1] << " "
            << "onCurveAngleTol : " << onCurveAngleTol << "\n";

  const auto debug = true;

  auto curveInfo = CurveReader::readCurveInfo(filename);
  CurveReader::printCurveInfo(curveInfo);

  OnCurve onCurve(onCurveAngleTol);
  const double tc_m1 = curveInfo.tcAngle.at(ptIdx-1);
  const double tc = curveInfo.tcAngle.at(ptIdx);
  const double tc_p1 = curveInfo.tcAngle.at(ptIdx+1);
  const auto on = onCurve(tc_m1, tc, tc_p1);
  std::cerr << "pt " << ptIdx 
            << " (" << curveInfo.x.at(ptIdx) << ", " << curveInfo.y.at(ptIdx) << ")"
            << " isOnCurve = " << on << "\n";
  if( on != isOnCurve ) {
    return 1;
  } else {
    return 0;
  }
}