incremental.cxx

/***************************************************************
  *  Copyright (c) 2013, Tsinghua University.
  *  This is a source file of C-Coupler.
  *  This file was initially finished by Haoyu Yang,
  *  supervised by Dr. Li Liu. If you have any problem,
  *  please contact Dr. Li Liu via liuli-cess@tsinghua.edu.cn
  ***************************************************************/

#include <cstdio>
#include <sys/time.h>
#include <cstdlib>
#include <cstring>
#include "common_utils.h"
#include "delaunay.h"

char cube_fine[][64] = {
    "CUBE_grid_0.1_0.flat",
    "CUBE_grid_0.1_1.flat",
    "CUBE_grid_0.1_2.flat",
    "CUBE_grid_0.1_3.flat",
    "CUBE_grid_0.1_4.flat",
    "CUBE_grid_0.1_5.flat",
    "CUBE_grid_0.1_6.flat",
    "CUBE_grid_0.1_7.flat",
    "CUBE_grid_0.1_8.flat",
    "CUBE_grid_0.1_9.flat",
};
char cube_midd[][64] = {
    "CUBE_grid_0.3_0.flat",
    "CUBE_grid_0.3_1.flat",
    "CUBE_grid_0.3_2.flat",
    "CUBE_grid_0.3_3.flat",
    "CUBE_grid_0.3_4.flat",
    "CUBE_grid_0.3_5.flat",
    "CUBE_grid_0.3_6.flat",
    "CUBE_grid_0.3_7.flat",
    "CUBE_grid_0.3_8.flat",
    "CUBE_grid_0.3_9.flat",
};
char cube_cors[][64] = {
    "CUBE_grid_1_0.flat",
    "CUBE_grid_1_1.flat",
    "CUBE_grid_1_2.flat",
    "CUBE_grid_1_3.flat",
    "CUBE_grid_1_4.flat",
    "CUBE_grid_1_5.flat",
    "CUBE_grid_1_6.flat",
    "CUBE_grid_1_7.flat",
    "CUBE_grid_1_8.flat",
    "CUBE_grid_1_9.flat",
};
char rect_fine[][64] = {
    "lonlat_grid_0.1_0.flat",
    "lonlat_grid_0.1_1.flat",
    "lonlat_grid_0.1_2.flat",
    "lonlat_grid_0.1_3.flat",
    "lonlat_grid_0.1_4.flat",
    "lonlat_grid_0.1_5.flat",
    "lonlat_grid_0.1_6.flat",
    "lonlat_grid_0.1_7.flat",
    "lonlat_grid_0.1_8.flat",
    "lonlat_grid_0.1_9.flat",
};
char rect_midd[][64] = {
    "lonlat_grid_0.3_0.flat",
    "lonlat_grid_0.3_1.flat",
    "lonlat_grid_0.3_2.flat",
    "lonlat_grid_0.3_3.flat",
    "lonlat_grid_0.3_4.flat",
    "lonlat_grid_0.3_5.flat",
    "lonlat_grid_0.3_6.flat",
    "lonlat_grid_0.3_7.flat",
    "lonlat_grid_0.3_8.flat",
    "lonlat_grid_0.3_9.flat",
};
char rect_cors[][64] = {
    "lonlat_grid_1_0.flat",
    "lonlat_grid_1_1.flat",
    "lonlat_grid_1_2.flat",
    "lonlat_grid_1_3.flat",
    "lonlat_grid_1_4.flat",
    "lonlat_grid_1_5.flat",
    "lonlat_grid_1_6.flat",
    "lonlat_grid_1_7.flat",
    "lonlat_grid_1_8.flat",
    "lonlat_grid_1_9.flat",
};
char rand_fine[][64] = {
    "lonlat_random_10000000_0.flat",
    "lonlat_random_10000000_1.flat",
    "lonlat_random_10000000_2.flat",
    "lonlat_random_10000000_3.flat",
    "lonlat_random_10000000_4.flat",
    "lonlat_random_10000000_5.flat",
    "lonlat_random_10000000_6.flat",
    "lonlat_random_10000000_7.flat",
    "lonlat_random_10000000_8.flat",
    "lonlat_random_10000000_9.flat",
};
char rand_midd[][64] = {
    "lonlat_random_1000000_0.flat",
    "lonlat_random_1000000_1.flat",
    "lonlat_random_1000000_2.flat",
    "lonlat_random_1000000_3.flat",
    "lonlat_random_1000000_4.flat",
    "lonlat_random_1000000_5.flat",
    "lonlat_random_1000000_6.flat",
    "lonlat_random_1000000_7.flat",
    "lonlat_random_1000000_8.flat",
    "lonlat_random_1000000_9.flat",
};
char rand_cors[][64] = {
    "lonlat_random_100000_0.flat",
    "lonlat_random_100000_1.flat",
    "lonlat_random_100000_2.flat",
    "lonlat_random_100000_3.flat",
    "lonlat_random_100000_4.flat",
    "lonlat_random_100000_5.flat",
    "lonlat_random_100000_6.flat",
    "lonlat_random_100000_7.flat",
    "lonlat_random_100000_8.flat",
    "lonlat_random_100000_9.flat",
};

char path[] = "grid/expand/%s";

int main(int argc, char* argv[])
{
    bool non_incremental = false;

    Delaunay_Voronoi* triangulation = new Delaunay_Voronoi();

    for (int i = 0; i < 10; i++)
    {
        char filepath[64];
        snprintf(filepath, 64, path, cube_fine[i]);

        if (i > 0 && non_incremental)
            strncat(filepath, ".png.flat", 64);

        FILE* fp = fopen(filepath, "r");
        if (!fp) {
            perror("No such file\n");
            return -1;
        }

        int num_points = 0;
        fscanf(fp, "%d", &num_points);
        if (num_points < 1) {
            perror("Wrong number of points\n");
            return -1;
        }

        double* x = new double[num_points];
        double* y = new double[num_points];

        for (int i = 0; i < num_points; i++)
            fscanf(fp, "%lf, %lf", &x[i], &y[i]);
        fclose(fp);

        if (non_incremental)
        {
            delete triangulation;
            triangulation = new Delaunay_Voronoi();
        }

        timeval start, end;
        gettimeofday(&start, NULL);

        triangulation->add_points(x, y, num_points);
        triangulation->triangulate();

        gettimeofday(&end, NULL);
        fprintf(stderr, "time: %ld us, points: %d\n", (end.tv_sec - start.tv_sec) * 1000000 + (end.tv_usec - start.tv_usec), num_points);

        delete[] x;
        delete[] y;
    }

    return 0;
}