/*
 *
 */

#define ENABLE_TIMER

#include "framework/DynamicExtension.h"
#include "shard/VPTree.h"
#include "query/knn.h"
#include "framework/interface/Record.h"
#include "file_util.h"
#include "standard_benchmarks.h"

#include <gsl/gsl_rng.h>

#include "psu-util/timer.h"


typedef Word2VecRec Rec;

typedef de::VPTree<Rec, 100, true> Shard;
typedef de::knn::Query<Rec, Shard> Q;
typedef de::DynamicExtension<Rec, Shard, Q, de::LayoutPolicy::TEIRING, de::DeletePolicy::TAGGING, de::SerialScheduler> Ext;
typedef de::DynamicExtension<Rec, Shard, Q, de::LayoutPolicy::LEVELING, de::DeletePolicy::TAGGING, de::SerialScheduler> Ext2;
typedef de::knn::Parms<Rec> QP;

void usage(char *progname) {
    fprintf(stderr, "%s reccnt datafile queryfile\n", progname);
}

int main(int argc, char **argv) {

    if (argc < 4) {
        usage(argv[0]);
        exit(EXIT_FAILURE);
    }

    size_t n = atol(argv[1]);
    std::string d_fname = std::string(argv[2]);
    std::string q_fname = std::string(argv[3]);

    gsl_rng * rng = gsl_rng_alloc(gsl_rng_mt19937);
    
    auto data = read_vector_file<Rec, 300>(d_fname, n);

    std::vector<size_t> to_delete(n * delete_proportion);
    size_t j=0;
    for (size_t i=0; i<data.size() && j<to_delete.size(); i++) {
        if (gsl_rng_uniform(rng) <= delete_proportion) {
            to_delete[j++] = i;
        } 
    }
    auto queries = read_knn_queries<QP>(q_fname, 10);


    const std::vector<de::LayoutPolicy> policies = {de::LayoutPolicy::LEVELING, de::LayoutPolicy::TEIRING};
    const std::vector<size_t> buffer_sizes = {100, 400, 800, 1200, 1500, 2000};
    const std::vector<size_t> scale_factors = {2, 4, 6, 8, 10, 12};

    for (const auto &bs : buffer_sizes) {
        for (const auto &sf : scale_factors) {
            auto extension = new Ext(bs, bs, sf, 0, 64);

            /* warmup structure w/ 10% of records */
            size_t warmup = .1 * n;
            size_t delete_idx = 0;
            insert_records<Ext, Rec>(extension, 0, warmup, data, to_delete, delete_idx, false, rng);

            extension->await_next_epoch();

            TIMER_INIT();

            TIMER_START();
            insert_records<Ext, Rec>(extension, warmup, data.size(), data, to_delete, delete_idx, true, rng);
            TIMER_STOP();

            auto insert_latency = TIMER_RESULT();
            size_t insert_throughput = (size_t) ((double) (n - warmup) / (double) insert_latency * 1e9);

            TIMER_START();
            run_queries<Ext, QP>(extension, queries);
            TIMER_STOP();

            auto query_latency = TIMER_RESULT() / queries.size();

            auto ext_size = extension->get_memory_usage() + extension->get_aux_memory_usage();

            fprintf(stdout, "TIERING\t%ld\t%ld\t%ld\t%ld\t%ld\n", bs, sf, insert_throughput, query_latency, ext_size);
            delete extension;
        }
    }

    for (const auto &bs : buffer_sizes) {
        for (const auto &sf : scale_factors) {
            auto extension = new Ext2(bs, bs, sf, 0, 64);

            /* warmup structure w/ 10% of records */
            size_t warmup = .1 * n;
            size_t delete_idx = 0;
            insert_records<Ext2, Rec>(extension, 0, warmup, data, to_delete, delete_idx, false, rng);

            extension->await_next_epoch();

            TIMER_INIT();

            TIMER_START();
            insert_records<Ext2, Rec>(extension, warmup, data.size(), data, to_delete, delete_idx, true, rng);
            TIMER_STOP();

            auto insert_latency = TIMER_RESULT();
            size_t insert_throughput = (size_t) ((double) (n - warmup) / (double) insert_latency * 1e9);

            TIMER_START();
            run_queries<Ext2, QP>(extension, queries);
            TIMER_STOP();

            auto query_latency = TIMER_RESULT() / queries.size();

            auto ext_size = extension->get_memory_usage() + extension->get_aux_memory_usage();

            fprintf(stdout, "LEVELING\t%ld\t%ld\t%ld\t%ld\t%ld\n", bs, sf, insert_throughput, query_latency, ext_size);
            delete extension;
        }
    }

    gsl_rng_free(rng);
    fflush(stderr);
}