vector_dist_performance_util.hpp

/*
 * vector_dist_performance_util.hpp
 *
 *  Created on: Jun 17, 2016
 *      Author: Yaroslav Zaluzhnyi
 */

#ifndef SRC_VECTOR_VECTOR_DIST_PERFORMANCE_UTIL_HPP_
#define SRC_VECTOR_VECTOR_DIST_PERFORMANCE_UTIL_HPP_

#include <boost/property_tree/ptree.hpp>
#include "Plot/GoogleChart.hpp"
#include "vector_dist_performance_common.hpp"

// Number of tests
#define N_VERLET_TEST 3
#define N_STAT_TEST 30


template<unsigned int dim, size_t prp = 0, typename T, typename V> double benchmark_calc_forces(T & NN, V & vd, float r_cut)
{
    //Timer
    timer t;
    t.start();

    calc_forces<dim,prp>(NN,vd,r_cut);

    t.stop();

    return t.getwct();
}

template<typename V> double benchmark_reorder(V & vd, size_t m)
{
    //Timer
    timer t;
    t.start();

    //Reorder
    vd.reorder(m);

    t.stop();

    return t.getwct();
}

template<typename T, typename V>
double benchmark_get_celllist(T & NN, V & vd, float r_cut)
{
    // Cell list timer
    timer t;
    t.start();

    //get cell list
    NN = vd.getCellList(r_cut);

    t.stop();

    return t.getwct();
}

template<typename V> double benchmark_get_verlet(V & vd, float r_cut)
{
    openfpm::vector<openfpm::vector<size_t>> verlet;
    //Timer
    timer t;
    t.start();

    //get verlet
    auto vr = vd.getVerlet(r_cut);

    t.stop();

    return t.getwct();
}



template<unsigned int dim, size_t prp = 0, typename T, typename V> double benchmark_calc_forces_hilb(T & NN, V & vd, float r_cut)
{
    //Forces timer
    timer t;
    t.start();

    calc_forces_hilb<dim,prp>(NN,vd,r_cut);

    t.stop();

    return t.getwct();
}

template<typename T, typename V> double benchmark_get_celllist_hilb(T & NN, V & vd, float r_cut)
{
    // Cell list timer
    timer t;
    t.start();

    //get cell list
    NN = vd.getCellList_hilb(r_cut);

    t.stop();

    return t.getwct();
}

template<unsigned int dim, typename v_dist> void move_particles(v_dist & vd, double dist)
{
    double phi;
    double theta;
    const double pi = 3.14159;

    auto it = vd.getDomainIterator();

    while (it.isNext())
    {
        phi = rand()/double(RAND_MAX);
        theta = rand()/double(RAND_MAX);
        phi *= 2.0*pi;
        theta *= pi;

        auto key = it.get();

        if(dim == 1)
            vd.getPos(key)[0] += dist*sin(theta)*cos(phi);
        else if(dim == 2)
        {
            vd.getPos(key)[0] += dist*sin(theta)*cos(phi);
            vd.getPos(key)[1] += dist*sin(theta)*sin(phi);
        }
        else if(dim == 3)
        {
            vd.getPos(key)[0] += dist*sin(theta)*cos(phi);
            vd.getPos(key)[1] += dist*sin(theta)*sin(phi);
            vd.getPos(key)[2] += dist*cos(phi);
        }

        ++it;
    }
}


extern void StandardPerformanceGraph(std::string file_mean,
                              std::string file_var,
                              std::string file_mean_save,
                              std::string file_var_save,
                              GoogleChart & cg,
                              openfpm::vector<size_t> & xp,
                              openfpm::vector<openfpm::vector<openfpm::vector<double>>> & yp_mean,
                              openfpm::vector<openfpm::vector<openfpm::vector<double>>> & yp_dev,
                              openfpm::vector<std::string> & names,
                              openfpm::vector<std::string> & gnames,
                              std::string x_string,
                              std::string y_string,
                              bool use_log);

#endif /* SRC_VECTOR_VECTOR_DIST_PERFORMANCE_UTIL_HPP_ */