Vector_util.hpp

/*
 * Vector_util.hpp
 *
 *  Created on: Dec 7, 2015
 *      Author: i-bird
 */

#ifndef OPENFPM_NUMERICS_SRC_VECTOR_VECTOR_UTIL_HPP_
#define OPENFPM_NUMERICS_SRC_VECTOR_VECTOR_UTIL_HPP_

#include "Grid/grid_dist_key.hpp"
#include "Space/Shape/HyperCube.hpp"
#include "util/mul_array_extents.hpp"

template<typename copy_type, typename T, typename Ev, typename Eqs_sys, int sa>
struct copy_ele_sca_array
{
    template<typename Grid> inline static void copy(Grid & grid_dst, const grid_dist_key_dx<Eqs_sys::dims> & key, const Ev & x,size_t lin_id, size_t base_id, size_t gs_size)
    {
        grid_dst.template get<T::value>(key) = x(lin_id * Eqs_sys::nvar + base_id);
    }
};

template<typename copy_type, typename T, typename Ev, typename Eqs_sys>
struct copy_ele_sca_array<copy_type,T,Ev,Eqs_sys,1>
{
    template<typename Grid> inline static void copy(Grid & grid_dst, const grid_dist_key_dx<Eqs_sys::dims> & key, const Ev & x,size_t lin_id, size_t base_id, size_t gs_size)
    {
        for (size_t i = 0 ; i < std::extent<copy_type>::value ; i++)
        {
            grid_dst.template get<T::value>(key)[i] = x(lin_id * Eqs_sys::nvar + base_id + i);
        }
    }
};


template<typename Eqs_sys, typename S, typename Ev>
struct copy_ele
{
    const grid_dist_key_dx<Eqs_sys::dims> key;

    S & grid_dst;

    size_t lin_id;

    size_t prp_id;

    size_t gs_size;

    const Ev & x;

    inline copy_ele(const grid_dist_key_dx<Eqs_sys::dims> & key, S & grid_dst, const Ev & x, size_t lin_id, size_t gs_size)
    :key(key),grid_dst(grid_dst),lin_id(lin_id),prp_id(0),gs_size(gs_size),x(x){};


#ifdef SE_CLASS1

    inline copy_ele(grid_dist_key_dx<Eqs_sys::dims> & key, S & grid_dst, Ev && x)
    :key(key),grid_dst(grid_dst),x(x)
    {std::cerr << "Error: " <<__FILE__ << ":" << __LINE__ << " Passing a temporal object";};
#endif

    template<typename T>
    inline void operator()(T& t)
    {
        // This is the type of the object we have to copy
        typedef typename boost::mpl::at_c<typename S::value_type::type,T::value>::type copy_type;

        copy_ele_sca_array<copy_type,T,Ev,Eqs_sys,std::is_array<copy_type>::value>::copy(grid_dst,key,x,lin_id,prp_id,gs_size);
        prp_id += array_extents<copy_type>::mul();
    }
};


#endif /* OPENFPM_NUMERICS_SRC_VECTOR_VECTOR_UTIL_HPP_ */