main.cpp File Reference

#include <cassert>
#include <iostream>
#include <fstream>
#include <GL/glut.h>
#include <GL/gl.h>
#include <zpr.h>
#include <graphmisc.h>
#include <commandline.h>
#include <gobj.h>
#include <pointsdisplay.h>
#include <random.h>
#include <print.h>
#include <point.h>
#include <d3tess.h>
#include <d3tessdraw.h>
#include <d2func.h>
#include <d2circle.h>
#include <halfspaceD2.h>
#include <aclock.h>
#include <d3meshpointreader.h>
#include <d3minboundary.h>
#include <d3minrecursive.h>
#include <d3mincircle.h>
#include <d3mincentroid.h>
#include <d3mincentroid2.h>
#include <message.h>
#include <typedefs.h>
#include <trianglespace.h>

Include dependency graph for main.cpp:

Go to the source code of this file.

Defines
#define	SHOW(x)   #x << '=' << (x)
Typedefs
typedef point2< double >	pt2
typedef point2< double > const	pt2c
Functions
bool	enabletargetcurve (false)
void	init02 ()
string	getInputString ()
boolc	addpointinsidemesh ()
void	addrandompointtomeshptvector (pt2 &x)
void	addrandompointtomesh ()
void	timmingexperiment04 ()
void	timmingexperiment03 ()
void	filereadingmesh (d3meshpointreader &mr)
void	timmingexperiment02 ()
void	timmingexperiment01 ()
void	timmingexperiment05 ()
void	timmingexperiment06 ()
void	timmingexperiment00 ()
void	minimizermenu ()
void	searchStats (double &searchlongest, double &searchaverage, uintc n)
void	searchexp01 ()
void	searchexp02 ()
void	filewritingmesh ()
void	filereadingmesh ()
void	filemenu ()
void	searchmenu ()
void	timmingsmenu ()
void	displaymenu ()
void	mainmenu ()
void	keyboard (unsigned char key, int x, int y)
void	display ()
void	addpoint (pt2 const &x)
void	addpoint (doublec x, doublec y)
void	init00 ()
void	help ()
int	main (int argc, char **argv)
Variables
d3tess *	mesh
d3tessdraw *	meshdraw
pt2	targetcenter
doublec	targetradius = 0.6
d2func *	p2dfunc = new d2circle(targetcenter,targetradius)
gobjMyCircle	targetcircle
gobjMyCircleDraw	targetcurve (targetradius, pt3(targetcenter.x, targetcenter.y, 0.0), targetcircle)
gobjContainer	xGraphics
messagefile	error ("error.txt", true)
random11< double >	r

---

Define Documentation

#define SHOW

(

x

)

#x << '=' << (x)

Definition at line 41 of file main.cpp.

---

Typedef Documentation

typedef point2<double> pt2

Definition at line 42 of file main.cpp.

typedef point2<double> const pt2c

Definition at line 43 of file main.cpp.

---

Function Documentation

void addpoint	(	doublec	x,
		doublec	y
	)

Definition at line 999 of file main.cpp.

References d2func::make(), and d3tess::pt.

{
  mesh->pt.push_back(p2dfunc->make(x,y));
}

void addpoint

(

pt2 const &

x

)

Definition at line 994 of file main.cpp.

References d2func::make(), d3tess::pt, point2< T >::x, and point2< T >::y.

Referenced by partitionstest::addstar().

{
  mesh->pt.push_back(p2dfunc->make(x.x,x.y));
}

boolc addpointinsidemesh

(

)

Definition at line 101 of file main.cpp.

References d3tess::addpoint(), d2func::eval(), messageglobal(), d3tess::pt, gobjContainer::push(), r, gobjQuadric::radius, d3tess::searchminimizetowardspoint(), point2< T >::x, and point2< T >::y.

Referenced by keyboard().

{
  

  static gobjQuadric q4;
  q4.radius=0.005;

  xGraphics.push( new gobjglColor3f(1.0,0.0,0.0) );
                
  uint k= mesh->pt.size();

  pt2 x(-1.0+2.0*r(),-1.0+2.0*r());
  mesh->pt.push_back(pt3(x.x,x.y,p2dfunc->eval(x)));
  xGraphics.push( new gobjMySphereDraw( x, &q4 ) );

  if (mesh->searchminimizetowardspoint(k))
  {


messageglobal() << "addpointinsidemesh()" << "\n";
messageglobal() << *mesh << "\n\n";

messageglobal() << "k=" << k << "\n";

    mesh->addpoint(k);
messageglobal() << *mesh << "\n\n";
    return true;
  }

  return false;
}

void addrandompointtomesh

(

)

Definition at line 140 of file main.cpp.

References d3tess::addpoint(), addrandompointtomeshptvector(), d3tess::pt, gobjContainer::push(), and gobjQuadric::radius.

Referenced by keyboard(), searchStats(), timmingexperiment00(), and timmingexperiment01().

{
  static gobjQuadric q3;
  q3.radius=0.005;

  xGraphics.push( new gobjglColor3f(1.0,0.0,0.0) );
                
  uint k= mesh->pt.size();

  pt2 x;
  addrandompointtomeshptvector(x);
  //pt2 x(-1.0+2.0*r(),-1.0+2.0*r());
  //mesh->pt.push_back(pt3(x.x,x.y,p2dfunc->eval(x)));
  xGraphics.push( new gobjMySphereDraw( x, &q3 ) );

  mesh->addpoint(k);
}

void addrandompointtomeshptvector

(

pt2 &

x

)

Definition at line 133 of file main.cpp.

References d2func::eval(), d3tess::pt, r, point2< T >::x, and point2< T >::y.

Referenced by addrandompointtomesh(), and searchStats().

{
  x = pt2(-1.0+2.0*r(),-1.0+2.0*r());
  mesh->pt.push_back(pt3(x.x,x.y,p2dfunc->eval(x)));
}

void display

(

)

Definition at line 957 of file main.cpp.

References gobjMyCircleDraw::draw(), gobjContainer::draw(), d3tessdraw::draw(), enabletargetcurve(), glerrordisplay(), and targetcurve.

{ 
  //glMatrix temp;
  myglPushMatrix temp;

  glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);

  //axes(0.5);

  meshdraw->draw();

  xGraphics.draw();

  glerrordisplay();

  {
    glPushMatrix();

    glPushAttrib(GL_CURRENT_BIT);
    glPushAttrib(GL_LIGHTING_BIT);
 
    glColor3f(1.0,1.0,0.0);
    glEnable(GL_LIGHTING);
  
    if (enabletargetcurve)
      targetcurve.draw();

    glPopAttrib();
    glPopAttrib();

    glPopMatrix();
  }
  glerrordisplay();
  
  glutSwapBuffers();
}

void displaymenu

(

)

Definition at line 737 of file main.cpp.

References getInputString(), d3tessdraw::graphicsDeffered, d3tessdraw::graphicsImmediate, gobjContainerSwitch::gswitch, and d3tessdraw::meshupdate().

Referenced by mainmenu().

{
  cout << endl;
  cout << "Display Menu" << endl;
  cout << endl;
  cout << "  0:  points numbering toggle" << endl;
  cout << "  1:  triangles numbering toggle" << endl;
  cout << "  2:  winding toggle" << endl;
  cout << "  3:  grid toggle" << endl;
  cout << "  4:  cp toggle" << endl;
  cout << "  5:  surface toggle" << endl;
  cout << "  6:  cp circle toggle" << endl;
  cout << "  7:  circles through triangles" << endl;
  cout << "  8:  cp voronoi local polygons about base vertices" << endl;
  cout << "  9:  multicolored triangles" << endl;
  cout << " 10:  voronoi diagram" << endl;
    
  cout << "      Enter a choice  ";

  string const s = getInputString();
  if (s.empty())
  {
    cout << "No Choice Executed" << endl; 
    cout << "End of Display Menu" << endl;
    return;
  }

  if (s=="0")
  {
    meshdraw->graphicsDeffered.gswitch[0]->toggle();
    meshdraw->meshupdate();
  }

  if (s=="1")
  {
    meshdraw->graphicsDeffered.gswitch[1]->toggle();
    meshdraw->meshupdate();
  }

  if (s=="9")
  {
    meshdraw->graphicsDeffered.gswitch[2]->toggle();
    meshdraw->meshupdate();
  }

  if (s=="10")
  {
    meshdraw->graphicsDeffered.gswitch[3]->toggle();
    meshdraw->meshupdate();
  }

  if (s=="2")
    meshdraw->graphicsImmediate.gswitch[0]->toggle();

  if (s=="3")
    meshdraw->graphicsImmediate.gswitch[1]->toggle();
  
  if (s=="4")
    meshdraw->graphicsImmediate.gswitch[2]->toggle();

  if (s=="5")
    meshdraw->graphicsImmediate.gswitch[3]->toggle();

  if (s=="6")
    meshdraw->graphicsImmediate.gswitch[4]->toggle();

  if (s=="7")
    meshdraw->graphicsImmediate.gswitch[5]->toggle();

  if (s=="8")
    meshdraw->graphicsImmediate.gswitch[6]->toggle();



  cout << "End of Display Menu" << endl;
}

bool enabletargetcurve

(

false

)

Referenced by display(), and main().

void filemenu

(

)

Definition at line 632 of file main.cpp.

References filereadingmesh(), filewritingmesh(), getInputString(), timmingexperiment02(), timmingexperiment03(), and timmingexperiment04().

Referenced by mainmenu().

{
  cout << endl;
  cout << "File Menu" << endl;
  cout << endl;
  cout << "  0:  Read a tessellation in from a file." << endl;
  cout << "  1:  Write a tessellation out to a file." << endl;
  cout << "  2:  Read points from file." << endl;
  cout << "  3:  Create a new data set by writing 2D random data points to a file." << endl;
  cout << "  4:  Create a new data set by writing 3D random data points to a file." << endl;

  cout << "      Enter a choice  ";

  string const s = getInputString();
  if (s.empty())
  {
    cout << "No Choice Executed" << endl; 
    cout << "End of File Menu" << endl;
    return;
  }  

  if (s=="0")
    filereadingmesh();

  if (s=="1")
    filewritingmesh();

  if (s=="2")
    timmingexperiment02();

  if (s=="3")
    timmingexperiment03();

  if (s=="4")
    timmingexperiment04();

}

void filereadingmesh

(

)

Definition at line 608 of file main.cpp.

References getInputString().

{
  cout << endl;
  cout << "  Reading a mesh in." << endl;
  cout << endl;
  cout << "  Enter the filename: ";

  string fname(getInputString());

  if (fname.empty())
    return;

  ifstream targ(fname.c_str());
 
  if (targ.good()==false)
  {
    cout << "error:  could not read the file." << endl;
    return;
  }

  targ >> *mesh;
}

void filereadingmesh

(

d3meshpointreader &

mr

)

Definition at line 200 of file main.cpp.

References aclock::diff_s(), d3meshpointreader::eval(), and aclock::measure().

Referenced by filemenu(), main(), and timmingexperiment02().

{
    
  cout << "Starting insertion of points into the mesh" << endl;

  aclock c;
  c.measure();
  mr.eval();
  c.measure();

  cout << "Finished insertion of points into the mesh" << endl << endl;
  cout << "Time without graphics: ";
  cout << c.diff_s() << "s" << endl;
}

void filewritingmesh

(

)

Definition at line 585 of file main.cpp.

References getInputString().

Referenced by filemenu().

{
  cout << endl;
  cout << "  Writing the mesh out." << endl;
  cout << endl;
  cout << "  Enter the filename: ";

  string fname(getInputString());

  if (fname.empty())
    return;

  ofstream targ(fname.c_str());
 
  if (targ.good()==false)
  {
    cout << "error:  could not write the file out." << endl;
    return;
  }

  targ << *mesh;
}

string getInputString

(

)

Definition at line 75 of file main.cpp.

Referenced by displaymenu(), filemenu(), filereadingmesh(), filewritingmesh(), mainmenu(), minimizermenu(), searchmenu(), and timmingsmenu().

{
  string s;
  cin >> s;
  return s;
/*
  string s;
  char ch;
  for ( ; ; )
  {
    ch = cin.peek();
    cin.ignore();

    if (ch=='\n')
      return s;

    s += ch;
  }

  assert(false);
  return s;
*/
}

void help

(

)

Definition at line 1032 of file main.cpp.

Referenced by maze005::keyboard01(), simplexD2tessapp01::keyboard02(), delaunaysimpleD2test::keyboard02(), and main().

{
  cout << "Command Line Options" << endl;
  cout << "For boolean values true and 1, false and 0 are equivalent" << endl;
  cout << "file=filename  - read in points from a file" << endl;
  cout << "cp=true        - current pointer" << endl;
  cout << "winding=true   - display the triangle windings" << endl;
  cout << "surface=true   - draw the surface" << endl;
  cout << "grid=true      - draw the triangles as blue grid" << endl;
  cout << "points=true    - number each point with its integer index" << endl;
  cout << "debug=0        - only use when NDEBUG is undefined so checkdebug() is not excecuted" << endl;
  cout << "triangles=1    - number each triangle with its integer index" << endl;
  cout << "cpcircle=1     - draw a circle through the cp's simplex." << endl;
  cout << "circles=1      - draw circles through every simplex." << endl;
  cout << "cpvoronoi=1    - display voronoi polygons about the cp."  << endl;
  cout << "multcolor=1    - display the triangles with random colorings" << endl;
  cout << "voronoi=1      - display strict voronoi polygons" << endl;


  cout << endl;
  cout << "Keyboard Commands" << endl;
  cout << "'T' enters you into the main menu which is displayed at the" << endl;
  cout << "  command line rather than the application." << endl << endl;
  cout << "'a' 'A' - rotate anticlockwise or clockwise the cp triangle." << endl;
  cout << "'j' 'k' - move left, right and down on current triangle" << endl;
  cout << "  'm'" << endl;
  cout << "'b'     - orient the current triangle to the boundary if possible" << endl;
  cout << "'w' 'e' - First orient the cp to the boundary with 'b'.  Then move" << endl;
  cout << "         on the surface left and right respectively." << endl;
  cout << "'z'     - the flip operator" << endl;

  cout << endl;
  cout << "Minimizers" << endl;
  cout << "  Explore mesh balancing." << endl;
  cout << "  By default no mesh balancing is performed." << endl;
  cout << "minb=1        - boundary minimization on insertion" << endl;
  cout << "minbr=1       - boundary minimization with recursion" << endl;
  cout << "minc=1        - circle minimization on insertion" << endl;
  cout << "minDT=1       - Delauny Triangulation - recursive circle minimization" << endl;
  cout << "mincentroid=1 - circle by average points and radius, insertion" << endl;
  cout << "mincentroid2r=1 - circle by average points and min radius, recursive" << endl;
  //cout << "greedy=true    - very bad mesh generated." << endl;
  //cout << "greedy2=true   - very experimental, do I know what I am doing?" << endl;

}

void init00

(

)

Definition at line 1004 of file main.cpp.

References d2func::make(), d3tess::pt, and d3tess::vi.

{
  cout << endl;
  cout << "Test Mesh" << endl;

  mesh->pt.push_back(p2dfunc->make(0.0,0.0));
  mesh->pt.push_back(p2dfunc->make(1.0,0.0));
  mesh->pt.push_back(p2dfunc->make(1.0,1.0));

  mesh->vi.push_back( simplexD2linked(1,2,3, 0,0,0) );
}

void init02

(

)

Definition at line 1016 of file main.cpp.

References d4tess::debugcheck(), d3func::eval(), d3tess::initialize(), d2func::make(), d4tess::pt, d3tess::pt, r, d4tess::viadd(), point4< T >::x, point4< T >::y, and point4< T >::z.

Referenced by main(), mainmenu(), searchStats(), and timmingexperiment01().

{
  cout << endl;
  cout << "Test Mesh" << endl;
  cout << "Adding three random points to make the first" << endl;
  cout << "  simplex or triangle in the mesh." << endl;

  for (uint i=0; i<3; ++i)
    mesh->pt.push_back( 
      p2dfunc->make(-1.0+2.0*r(),-1.0+2.0*r()) );

  mesh->initialize();
}

void keyboard	(	unsigned char	key,
		int	x,
		int	y
	)

Definition at line 868 of file main.cpp.

References addpointinsidemesh(), addrandompointtomesh(), virtualtriangle::anticlockwise(), d3tess::boundaryorient(), virtualtriangle::clockwise(), d3tess::flip(), mainmenu(), d3tessdraw::meshupdate(), d3tess::movedown(), d3tess::moveleft(), d3tess::moveright(), d3tess::surfaceleft(), d3tess::surfaceright(), and d3tess::vs.

{
  switch (key)
  {
    case 27:
      delete mesh;
      mesh=0;
      delete meshdraw;
      meshdraw=0;
      delete p2dfunc;
      p2dfunc=0;
      
      exit(0);
      break;

    case 'a':
      mesh->vs.anticlockwise();
      break;

    case 'A':
      mesh->vs.clockwise();
      break;

    case 'j':
      mesh->moveleft();
      break;

    case 'k':
      mesh->moveright();
      break;

    case 'm':
      mesh->movedown();
      break;

    case 'b':
      mesh->boundaryorient();
      break;

    case 'w':
      mesh->surfaceleft();
      break;
    case 'e':
      mesh->surfaceright();
      break;

    case 'R':
      {
        for (; !addpointinsidemesh(); );

        meshdraw->meshupdate();
      }
      break;


    case 'r':
      addrandompointtomesh();
      meshdraw->meshupdate();
      break;

    case 't':
      for (uint i=0; i<100; ++i) 
        addrandompointtomesh();
      meshdraw->meshupdate();
      break;

    case 'T':
      mainmenu();
      break;

    case 'p':
      cout << endl;
      cout << *mesh;
      break;

    case 'z':
      mesh->flip();
      meshdraw->meshupdate();
      break;

    case 'l':
      cout << "l=" << mesh->cpbasemeasure() << endl;
      break;

  }

  glutPostRedisplay();
}

int main	(	int	argc,
		char **	argv
	)

Definition at line 1079 of file main.cpp.

References d3tessdraw::circle, d3tessdraw::circles, d3tessdraw::cp, d3tess::debugenable, display(), enabletargetcurve(), error, filereadingmesh(), message::globalset(), d3tessdraw::graphicsDeffered, d3tessdraw::graphicsImmediate, d3tessdraw::grid, gobjContainerSwitch::gswitch, help(), init02(), keyboard(), commandline::mapvar(), d3tessdraw::meshupdate(), d3tess::minimizerSet(), d3tessdraw::multicolor, d3tessdraw::points, random11< T, G >::rg, randomgenerator::seed(), gobjContainer::set(), d3tessdraw::simplexes, d3tessdraw::surface, d3tessdraw::voronoi, d3tessdraw::voronoip, and d3tessdraw::winding.

{
  glutInit(&argc,argv);
  glutInitDisplayMode(GLUT_DOUBLE | GLUT_RGB | GLUT_DEPTH);
  uintc wx=800;
  uintc wy=800;

  glutInitWindowSize(wx,wy);

  glutCreateWindow("");

  glutDisplayFunc(display);
  glutKeyboardFunc(keyboard);

  OpenGLinitialisation();

  glEnable(GL_DEPTH_TEST);

  //glDisable(GL_CULL_FACE);
  glEnable(GL_CULL_FACE);

  glEnable(GL_NORMALIZE);

  xGraphics.set();

  //error.setMessageGlobal();
  error.globalset();

  help();

  commandline cmd(argc,argv);

  uint N(1000);
  cmd.mapvar(N,"N");
  
  mesh = new d3tess(N);

  cmd.mapvar(mesh->debugenable,"debug");

  meshdraw = new d3tessdraw(*mesh);

  cmd.mapvar(meshdraw->graphicsDeffered.gswitch[d3tessdraw::points]->isdrawn,"points");
  cmd.mapvar(meshdraw->graphicsDeffered.gswitch[d3tessdraw::simplexes]->isdrawn,"triangles");
  cmd.mapvar(meshdraw->graphicsDeffered.gswitch[d3tessdraw::multicolor]->isdrawn,"multicolor");
  cmd.mapvar(meshdraw->graphicsDeffered.gswitch[d3tessdraw::voronoi]->isdrawn,"voronoi");


  cmd.mapvar(meshdraw->graphicsImmediate.gswitch[d3tessdraw::winding]->isdrawn,"winding");
  cmd.mapvar(meshdraw->graphicsImmediate.gswitch[d3tessdraw::grid]->isdrawn,"grid",true);
  cmd.mapvar(meshdraw->graphicsImmediate.gswitch[d3tessdraw::cp]->isdrawn,"cp",true);
  cmd.mapvar(meshdraw->graphicsImmediate.gswitch[d3tessdraw::surface]->isdrawn,"surface");
  cmd.mapvar(meshdraw->graphicsImmediate.gswitch[d3tessdraw::circle]->isdrawn,"cpcircle");
  cmd.mapvar(meshdraw->graphicsImmediate.gswitch[d3tessdraw::circles]->isdrawn,"circles");
  cmd.mapvar(meshdraw->graphicsImmediate.gswitch[d3tessdraw::voronoip]->isdrawn,"cpvoronoi");



  bool enableminb(false);
  cmd.mapvar(enableminb,"minb");
  if (enableminb)
    mesh->minimizerSet( new d3minboundary(*mesh) );

  bool enableminbr(false);
  cmd.mapvar(enableminbr,"minbr");
  if (enableminbr)
    mesh->minimizerSet
    ( 
      new d3minrecursiveoperator(*mesh,new d3minboundary(*mesh)) 
    );

  bool enableminc(false);
  cmd.mapvar(enableminc,"minc");
  if (enableminc)
    mesh->minimizerSet( new d3mincircle(*mesh) );

  bool enableminDT(false);
  cmd.mapvar(enableminDT,"minDT");
  if (enableminDT)
    mesh->minimizerSet
    ( 
      new d3minrecursiveoperator(*mesh,new d3mincircle(*mesh)) 
    );

  bool enablemincentroid(false);
  cmd.mapvar(enablemincentroid,"mincentroid");
  if (enablemincentroid)
    mesh->minimizerSet( new d3mincentroid(*mesh) );

  bool enablemincentroid2r(false);
  cmd.mapvar(enablemincentroid2r,"mincentroid2r");
  if (enablemincentroid2r)
    mesh->minimizerSet
    ( 
      new d3minrecursiveoperator(*mesh,new d3mincentroid2(*mesh)) 
    );

  cmd.mapvar(enabletargetcurve,"targetcurve");

  //test01();

  r.rg.seed(17);

  string fname;
  cmd.mapvar(fname,"file");
  bool fileread(false);

  if (fname.empty()==false)
  {
    d3meshpointreader mr(fileread,*mesh,fname);
 
    if (fileread)
      filereadingmesh(mr);
    else
      cout << "error:  could not read file " << fname << endl;
  }

  if (fileread==false)
    init02();

  meshdraw->meshupdate();

  zpr zz;

  glutPostRedisplay();

  glutMainLoop();

  return 0;
}

void mainmenu

(

)

Definition at line 815 of file main.cpp.

References displaymenu(), filemenu(), getInputString(), init02(), d3tessdraw::meshupdate(), minimizermenu(), d3tess::reset(), searchmenu(), and timmingsmenu().

Referenced by keyboard(), and menusystemtest::test02().

{
  for (; ; )
  {

    cout << endl;
    cout << "Main Menu" << endl;
    cout << endl;
    cout << "  0:  Timmings Menu" << endl;
    cout << "  1:  Display Menu" << endl;
    cout << "  2:  File Menu" << endl;

    cout << "  5:  Change the Minimizer." << endl;
    cout << "  6:  Reset the mesh." << endl;
    cout << "  9:  Search Menu" << endl;
    
    cout << "      Enter a choice  ";


    string const s = getInputString();
    if (s.empty())
    {
      cout << "No Choice Executed" << endl; 
      cout << "End of Main Menu" << endl;
      return;
    }

    if (s=="0")
      timmingsmenu();

    if (s=="1")
      displaymenu();

    if (s=="2")
      filemenu();

    if (s=="5")
      minimizermenu();

    if (s=="6")
    {
      mesh->reset();
      init02();
      meshdraw->meshupdate();
    }

    if (s=="9")
      searchmenu();

  }
}

void minimizermenu

(

)

Definition at line 418 of file main.cpp.

References getInputString(), and d3tess::minimizerSet().

Referenced by mainmenu().

{
  cout << endl;
  cout << "Minimizer Menu" << endl;
  cout << endl;
  cout << "The tessellator has a minimizer."  << endl;
  cout << "Setting the minimizer affects all future mesh insertions." << endl; 
  cout << endl;
  cout << "  0: No minimizer" << endl;
  cout << "  1: minb        - minimize the boundary on insertion" << endl;
  cout << "  2: minbr" << endl;
  cout << "  3: minc" << endl;
  cout << "  4: minDT   - Delaunay Triangulation" << endl;
  cout << "  5: mincentroid" << endl;
  cout << "  6: mincentroidr" << endl;

  cout << "      Enter a choice  ";

  string const s = getInputString();
  if (s.empty())
  {
    cout << "No Choice Executed" << endl; 
    cout << "End of Minimizer Menu" << endl;
    return;
  }

  if (s=="0")
    mesh->minimizerSet( new d3minoperator(*mesh) );

  if (s=="1")
    mesh->minimizerSet( new d3minboundary(*mesh) );
 
  if (s=="2")
    mesh->minimizerSet
    ( 
      new d3minrecursiveoperator(*mesh,new d3minboundary(*mesh)) 
    );

  if (s=="3")
    mesh->minimizerSet( new d3mincircle(*mesh) );

  if (s=="4")     
    mesh->minimizerSet
    ( 
      new d3minrecursiveoperator(*mesh,new d3mincircle(*mesh)) 
    );

  if (s=="5")
    mesh->minimizerSet( new d3mincentroid(*mesh) );

  if (s=="6")
    mesh->minimizerSet
    ( 
      new d3minrecursiveoperator(*mesh,new d3mincentroid2(*mesh)) 
    );

  cout << "End of Choice Menu" << endl;
}

void searchexp01

(

)

Definition at line 533 of file main.cpp.

References d3tessdraw::meshupdate(), searchStats(), and SHOW.

Referenced by searchmenu().

{
  cout << "Enter the number of points: ";
  uint n;
  cin >> n;

  double searchlongest;
  double searchaverage;

  searchStats(searchlongest,searchaverage,n);

  cout << SHOW(searchlongest) << endl;
  cout << SHOW(searchaverage) << endl;

  meshdraw->meshupdate();
}

void searchexp02

(

)

Definition at line 550 of file main.cpp.

References searchStats().

Referenced by searchmenu().

{
  uint nmax = 1000000;

  cout << "Writing search results to search.txt" << endl;
  cout << "The number of points is capped at " << nmax << endl;

  string fname("search.txt");
  ofstream f(fname.c_str());

  if (f.good()==false)
  {
    cout << "error:  can not open file: " << fname << endl;
    return;
  }

  uint n = 8;
  uint N = 3;

  cout << "Enter the initial number of points u and the number of times" << endl;
  cout << "  this is doubled." << endl;
  cin >> n;
  cin >> N;

  double searchlongest;
  double searchaverage;

  for (uint i=0; i<N; ++i)
  {
    searchStats(searchlongest,searchaverage,n);
    f << n << " " << searchaverage << endl;
    n *= 2;
  }
}

void searchmenu

(

)

Definition at line 671 of file main.cpp.

References getInputString(), searchexp01(), and searchexp02().

Referenced by mainmenu().

{
  cout << endl;
  cout << "Search Menu" << endl;
  cout << endl;
  cout << "  0: Calculate the shortes, longest and average lengths" << endl;
  cout << "       of a n point data set with n random inside searches." << endl;
  cout << "  1: Write the average search lengths to a file search.txt ." << endl;

  cout << "      Enter a choice  ";

  string const s = getInputString();
  if (s.empty())
  {
    cout << "No Choice Executed" << endl; 
    cout << "End of Search Menu" << endl;
    return;
  }

  if (s=="0")
    searchexp01();

  if (s=="1")
    searchexp02();

  cout << "End of Search Menu" << endl;
}

void searchStats	(	double &	searchlongest,
		double &	searchaverage,
		uintc	n
	)

Definition at line 479 of file main.cpp.

References addrandompointtomesh(), addrandompointtomeshptvector(), d3tess::cp, init02(), d3tess::movecounter, d3tess::pt, d3tess::reset(), and d3tess::searchminimizetowardspoint().

Referenced by searchexp01(), and searchexp02().

{
  // Create the mesh first.
  mesh->reset();
  init02();
  
  for (uint i=3; i<n; ++i) 
    addrandompointtomesh();

  // Create a vector of points inside the mesh
  vector<uint> v;

  pt2 x;
  uint k;
  for ( ; v.size()<n; )
  {
    k = mesh->pt.size();
    addrandompointtomeshptvector(x);
    if (mesh->searchminimizetowardspoint(k))
      v.push_back(k);
  }

  // Set the cp to some random point inside the mesh
  mesh->cp = (rand()% n)+1;

  mesh->searchminimizetowardspoint(v[0]);
  double c = mesh->movecounter; 
  searchlongest = c;
  searchaverage = c;

  for (uint i=1, sz=v.size(); i<sz; ++i)
  {
#ifndef NDEBUG
    assert( mesh->searchminimizetowardspoint(v[i]) );
#else
    mesh->searchminimizetowardspoint(v[i]);
#endif

    c = mesh->movecounter; 
    if (c>searchlongest)
      searchlongest = c;

    searchaverage += c; 
  }

  searchaverage *= 1.0/n;
}

void timmingexperiment00

(

)

Definition at line 399 of file main.cpp.

References addrandompointtomesh(), aclock::diff_s(), aclock::measure(), d3tessdraw::meshupdate(), and d3tess::vi.

Referenced by timmingsmenu().

{
  cout << "Timing experiment." << endl;
  cout << "Enter the number of points n: ";
  uint n;
  cin >> n;
  mesh->vi.reserve(n);

  aclock c;
  c.measure();
  for (uint i=0; i<n; ++i) 
    addrandompointtomesh();
  c.measure();
  cout << "Time without graphics: ";
  cout << c.diff_s() << "s" << endl;

  meshdraw->meshupdate();
}

void timmingexperiment01

(

)

Definition at line 240 of file main.cpp.

References addrandompointtomesh(), aclock::diff_s(), init02(), aclock::measure(), d3tess::reset(), SHOW, and d3tess::vi.

Referenced by timmingsmenu().

{
  cout << "Multiple timing experiments" << endl;
  cout << endl;
  cout << "Enter the filename: ";
  string fname;
  cin >> fname;
  
  ofstream fil(fname.c_str());

  cout << "Enter number of measurements M: ";
  uint M;
  cin >> M;

  uint n=100;

  aclock c;

  for (uint k=0; k<M; ++k)
  {
cout << SHOW(k) << endl;
cout << SHOW(n) << endl;
    mesh->vi.reserve(n);

    c.measure();
    for (uint i=0; i<n; ++i) 
      addrandompointtomesh();
    c.measure();
    fil << n << " " << c.diff_s() << endl;

    mesh->reset();

    init02();

    n *= 2;
  }
cout << "out of loop" << endl;
exit(0);

}

void timmingexperiment02

(

)

Definition at line 217 of file main.cpp.

References filereadingmesh().

Referenced by filemenu(), and keyboard().

{
  cout << "Reading Points from a File" << endl;
  cout << "Expecting the file to have 2 or 3 columns of numbers." << endl;
  cout << endl;
  cout << "  Enter the filename: ";
  string fname;

  //fname = "dataset2D08.txt";
  //cout << fname << endl;
  cin >> fname;

  bool res;
  d3meshpointreader mr(res,*mesh,fname);

  if (res==false)
    return;

  filereadingmesh(mr);
}

void timmingexperiment03

(

)

Definition at line 180 of file main.cpp.

References r.

Referenced by filemenu(), and keyboard().

{
  cout << "Creating a 2D Data Set" << endl;
  cout << "Writing Points to a File" << endl;
  cout << "  Enter the filename: ";
  string fname;
  cin >> fname;
  cout << "  Enter the number of points: ";
  uint n;
  cin >> n;

  ofstream targ(fname.c_str());

  for (uint i=0; i<n; ++i)
    targ << -1.0+2.0*r() << " " << -1.0+2.0*r() << endl;

  cout << "Finished writing file: " << fname << endl << endl;
}

void timmingexperiment04

(

)

Definition at line 158 of file main.cpp.

References d2func::eval(), r, point2< T >::x, and point2< T >::y.

Referenced by filemenu().

{
  cout << "Creating a 3D Data Set" << endl;
  cout << "Writing Points to a File" << endl;
  cout << "  Enter the filename: ";
  string fname;
  cin >> fname;
  cout << "  Enter the number of points: ";
  uint n;
  cin >> n;

  ofstream targ(fname.c_str());

  for (uint i=0; i<n; ++i)
  {
    pt2 x(-1.0+2.0*r(),-1.0+2.0*r());
    targ << x.x << " " << x.y << " " << p2dfunc->eval(x) << endl;
  }

  cout << "Finished writing file: " << fname << endl << endl;
}

void timmingexperiment05

(

)

Definition at line 281 of file main.cpp.

References d3meshpointreader::eval().

Referenced by timmingsmenu().

{
  cout << "Multiple Timing Experiments." << endl;
  cout << "The results are stored in timings.txt" << endl;

  cout << "Enter the filename with a list of 2D data files to be timed: ";
  string filenames;
  //cin >> filenames;
  filenames="dataset2Dlist.txt";
  cout << filenames << endl;

  ifstream filelist(filenames.c_str());

  ofstream filtim("timings.txt");

  if (filelist.good()==false)
  {
    cout << "error:  can not open file: " << filenames << endl;
    return;
  }

  while (filelist.eof()==false)
  {
    string f1;
    filelist >> f1;

    if (f1.empty())
      continue;

    cout << "Started reading file: " << f1 << endl;
    bool res;
    d3meshpointreader mr(res,*mesh,f1);
    if (res==false)
    {
      cout << "error:  could not read file:" << f1 << "." << endl;
      continue;
    }
    cout << "Finished reading file: " << f1 << endl;
    cout << "Timing file:" << f1 << "." << endl;

    cout << "Starting insertion of points into the mesh" << endl;
    mr.eval(filtim);
    cout << "Finished insertion of points into the mesh" << endl << endl;
  }
}

void timmingexperiment06

(

)

Definition at line 330 of file main.cpp.

References d3meshpointreader::eval().

Referenced by timmingsmenu().

{
  cout << "Simulating Multiple Timing Experiments." << endl;
  cout << "The results are stored in timings.txt" << endl;
  cout << "This is testing an algorithm not yet build." << endl;
  cout << "  It has the same complexity as this one." << endl;
  cout << "  A O(n) operation is applied after this step." << endl;
  cout << endl;
  cout << "A bounding box is first inserted.  All the points " << endl;
  cout << "  should lie well within the bounding box." << endl;
  cout << endl;
  cout << "For a preordered data set I am expecting linear results" << endl;


  cout << "Enter the filename with a list of 2D data files to be timed: ";
  string filenames;
  //cin >> filenames;
  filenames="datasetlist.txt";
  cout << filenames << endl;

  ifstream filelist(filenames.c_str());

  ofstream filtim("timings.txt");

  if (filelist.good()==false)
  {
    cout << "error:  can not open file: " << filenames << endl;
    return;
  }

  // Points defining the bounding box.
  vector< pt3 > vbox;
  vbox.push_back( pt3(-2.0,-2.0,0.0) );
  vbox.push_back( pt3(2.0,-2.0,0.0) );
  vbox.push_back( pt3(2.0,2.0,0.0) );
  vbox.push_back( pt3(-2.0,2.0,0.0) );


  while (filelist.eof()==false)
  {
    string f1;
    filelist >> f1;

    if (f1.empty())
      continue;

    cout << "Started reading file: " << f1 << endl;
    bool res;
    d3meshpointreader mr(res,*mesh,vbox,f1);
    if (res==false)
    {
      cout << "error:  could not read file:" << f1 << "." << endl;
      continue;
    }
    cout << "Finished reading file: " << f1 << endl;
    cout << "Timing file:" << f1 << "." << endl;

    cout << "Starting insertion of points into the mesh" << endl;
    mr.eval(filtim);
    cout << "Finished insertion of points into the mesh" << endl << endl;
  }
}

void timmingsmenu

(

)

Definition at line 699 of file main.cpp.

References getInputString(), timmingexperiment00(), timmingexperiment01(), timmingexperiment05(), and timmingexperiment06().

Referenced by mainmenu().

{
  cout << endl;
  cout << "Timmings Menu" << endl;
  cout << endl;
  cout << "  0:  One tessellation with n points measured." << endl;
  cout << "  1:  Multiple timing measurements to file." << endl;
  cout << "  2:  Multiple timmings on 2D data files" << endl;
  cout << "  3:  Simulating Multiple timmings on my Bounded Algorithm" << endl;
  cout << "      Enter a choice  ";


  string const s = getInputString();
  if (s.empty())
  {
    cout << "No Choice Executed" << endl; 
    cout << "End of Timmings Menu" << endl;
    return;
  }

  if (s=="0")
    timmingexperiment00();

  if (s=="1")
    timmingexperiment01();

  if (s=="2")
    timmingexperiment05();

  if (s=="3")
    timmingexperiment06();

  cout << "End of Timmings Menu" << endl;
}

---

Variable Documentation

messagefile error("error.txt", true)

Referenced by main(), minexpdimtest::test01(), minexpdimtest::test03(), desystest::test03(), minexpdimtest::test04(), desystest::test04(), and minexpdimtest::test05().

d3tess* mesh

Definition at line 48 of file main.cpp.

d3tessdraw* meshdraw

Definition at line 49 of file main.cpp.

d2func* p2dfunc = new d2circle(targetcenter,targetradius)

Definition at line 53 of file main.cpp.

random11<double> r

Definition at line 69 of file main.cpp.

Referenced by protractor::addAngleRuler(), addpointinsidemesh(), protractor::addRadiusRuler(), addrandompointtomeshptvector(), pointgrid3D::createIndexedTriangles(), d3writefile(), tessD1draw01linesmulticolor< TESS, PT >::draw(), triangledisplaygergonnepoint< TR >::draw(), triangledisplayincenter< TR >::draw(), triangledisplaycircumcenter< TR >::draw(), tessD2draw02circles< TESS, PT, INDX >::draw(), regionD2tessdisplaymesh< TS, Indx >::draw(), writemulticolorobj::draw(), d4mingreedy2::eval(), trianglevisualize< T >::eval(), tessD2disp01< PT >::eval(), polytopeD2tessdisp01< VPTS, VPOLY >::eval(), generateRandomPointsInCircle< T, D, RG >::generateRandomPointsInCircle(), generateRandomPointsInSphere< T, D, RG >::generateRandomPointsInSphere(), visconvex::handlecommand(), mazematrixD2createmaze< T >::impropermaze01(), init02(), particle::isIntersecting(), keyboard(), integration< D, F >::LJcalc(), d2homogeneous::matrixMultiply(), zprmouse::mouseratio(), d2homogeneous::operator*=(), halfspaceD2draw::rotate(), d2homogeneous::setRotateAboutPoint(), gausselim< T >::substitute(), randomtest::test01(), histogramtest::test01(), delaunaysimpleD2test::test01(), test01(), spiralindextest::test02(), randomtest::test02(), delaunaysimpleD2test::test02(), test02(), randomtest::test03(), halfspaceD3test::test03(), triangles3Tdisplaytest::test03(), test03(), spiralindextest::test04(), randomtest::test04(), triangles3Tdisplaytest::test04(), partitionstest::test06(), test13(), timmingexperiment03(), timmingexperiment04(), triangles3Tdisplaymulticolored< T, W >::triangles3Tdisplaymulticolored(), pointgrid3Ddraw::update(), pointgrid3Dbilineardraw::update(), and rpnp_to_c::visit().

pt2 targetcenter

Definition at line 51 of file main.cpp.

gobjMyCircle targetcircle

Definition at line 57 of file main.cpp.

gobjMyCircleDraw targetcurve(targetradius, pt3(targetcenter.x, targetcenter.y, 0.0), targetcircle)

Referenced by display().

doublec targetradius = 0.6

Definition at line 52 of file main.cpp.

gobjContainer xGraphics

Definition at line 65 of file main.cpp.