source: GTP/trunk/App/Demos/Illum/pathmap/DelaunayMesh.cpp @ 2197

// DelaunayMesh.cpp: implementation of the DelaunayMesh class.
//
//////////////////////////////////////////////////////////////////////

#include "dxstdafx.h"
#include "DelaunayMesh.h"

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

DelaunayMesh::DelaunayMesh(Vector* c,double encloseRad)
{
        encloseRadius = encloseRad;
        encloseCentre = *c;
//to see how we crashed for testing
        crashed = 0;
//WARNING should contain everything, now still hardwired for testing
        tetrahedra = new Woctreebranch(c->x,c->y,c->z,
                encloseRad*5,encloseRad*5,encloseRad*5);

//enclosing dodecaeder
        double radius = encloseRadius*4;
//WARNING relative coords to c 
        Vector* top = new Vector(c->x, c->y, c->z + radius);
        nodeArray.push_back( top );
        Vector* bottom = new Vector(c->x, c->y, c->z - radius);
        nodeArray.push_back( bottom );
        
        Vector* highcircle[5];
        Vector* lowcircle[5];
        for(int i=0;i<5;i++)
        {
                highcircle[i] = new Vector(c->x + radius*sin(i*1.257),c->y + radius*cos(i*1.257),c->z + radius/3.0);
                nodeArray.push_back( highcircle[i] );
                lowcircle[i] = new Vector(c->x + radius*sin(i*1.257+0.628),c->y + radius*cos(i*1.257+0.628),c->z - radius/3.0);
                nodeArray.push_back( lowcircle[i] );
        }
        nodeArray.push_back( c );

        cavvie = new WsubMesh();
        for(int s=0;s<5;s++)
        {
                cavvie->add(new Wsimplex(top, highcircle[s], highcircle[(s+1)%5],1));
                cavvie->add(new Wsimplex(bottom, lowcircle[s], lowcircle[(s+1)%5],1));
                cavvie->add(new Wsimplex(highcircle[s], lowcircle[(s+4)%5], lowcircle[s],1));
                cavvie->add(new Wsimplex(lowcircle[s], highcircle[s], highcircle[(s+1)%5],1));
        }
        ::std::vector<Wsimplex*> a = cavvie->complete;
        ::std::vector<Wsimplex*>::const_iterator b = a.begin();

        while(b != a.end())
        {
                Wsimplex* pyramid = new Wsimplex(c,*b);
                (*b)->inFrontOf = pyramid;
                tetrahedra->add(pyramid);
                simplexArray.push_back(pyramid);
                b++;
        }
}


DelaunayMesh::~DelaunayMesh()
{
        delete tetrahedra;
        delete cavvie;
}

char DelaunayMesh::add(Vector * proxy,Wsimplex* starter)
{
        if(crashed) return 'C';
        Wsimplex* first;
        if(starter)
                first = starter;
        else
                first = tetrahedra->search(proxy);
        if(!first)
                return 'I';

        WsubMesh* suff = first->cavity(proxy);

        if(!suff->incomplete.empty())
        {
                crashed = 1;
                return 'C';
        }

        ::std::vector<Wsimplex*>::iterator a = simplexArray.begin();
        while(a != simplexArray.end())
        {
                if((*a)->deleted)
                {
                        delete(*a);
                        simplexArray.erase(a);
                }
                else
                        a++;
        }

        a = suff->complete.begin();
        while(a != suff->complete.end())
        {
                Wsimplex* pyramid = new Wsimplex(proxy,*a);
                (*a)->inFrontOf = pyramid;
                tetrahedra->add(pyramid);
                simplexArray.push_back(pyramid);
                a++;
        }

        nodeArray.push_back(proxy);

        delete suff;
        return 'O';
}

int DelaunayMesh::getSimplexCount()
{
        return simplexArray.size();
}

Wsimplex* DelaunayMesh::search(Vector* pixie)
{
        return tetrahedra->search(pixie);
}