#include <iostream>
#include <sstream>
#include <cstdlib>
#include <cmath>
#include <memory>

#include <TROOT.h>
#include <TGClient.h> 
#include <TGButton.h> 
#include <TGFrame.h> 
#include <TApplication.h>
#include <TCanvas.h>
#include <TH2F.h>
#include <TRandom.h>

#include "IDelaunay2D.hxx"

int main(int argc, char **argv) {
    int points = 2;
    if (argc>1) {
        std::istringstream arg(argv[1]);
        arg >> points;
    }
    
    int trials = 1;
    if (argc>2) {
        std::istringstream arg(argv[2]);
        arg >> trials;
    }
    
    TRandom ran;
    for (int i=0; i<trials; ++i) {
        std::cout << "trial: " << i << " " << points << std::endl;

        std::auto_ptr<COMET::IDelaunay2D> gDelaunay(new COMET::IDelaunay2D);

        gDelaunay->AddPoint(new COMET::IMeshXYPoint(10,10));
        
        gDelaunay->AddPoint(new COMET::IMeshXYPoint(4,10));
        
        gDelaunay->AddPoint(new COMET::IMeshXYPoint(10,4));
        
        gDelaunay->AddPoint(new COMET::IMeshXYPoint(12,12));
        
        // Add the requested number of points.
        for (int i=0; i<points; ++i) {
            double xmean = 5;
            double ymean = 5;
            if (ran.Integer(2)%2 == 1) xmean = -5; 
            if (ran.Integer(2)%2 == 1) ymean = -5; 
            xmean += 10;
            ymean += 10;
            double x = std::floor(10*ran.Gaus(xmean))/10.0;
            double y = std::floor(10*ran.Gaus(ymean))/10.0;
            gDelaunay->AddPoint(new COMET::IMeshXYPoint(x,y));
        }
        
        for (int i=0; i<points/10; ++i) {
            double x = ran.Uniform(1.0,19.0);
            double y = ran.Uniform(1.0,19.0);
            gDelaunay->AddPoint(new COMET::IMeshXYPoint(x,y));
        }
        
        gDelaunay->Close();
    }
}