DcpseInterpolation.hpp

//
// Created by Abhinav Singh on 03.11.21.
//
 
#ifndef OPENFPM_PDATA_DCPSEINTERPOLATION_HPP
#define OPENFPM_PDATA_DCPSEINTERPOLATION_HPP
#include "DCPSE/Dcpse.hpp"
 
template<typename particlesFrom_type, typename particlesTo_type>
class PPInterpolation 
{
 
    void *dcpse;
 
    particlesFrom_type & particlesFrom;
    particlesTo_type & particlesTo;
 
public:
    PPInterpolation(particlesFrom_type &particlesFrom,particlesTo_type &particlesTo, unsigned int ord, typename particlesFrom_type::stype rCut,
                      double oversampling_factor = dcpse_oversampling_factor,
                      support_options opt = support_options::RADIUS)
    :particlesFrom(particlesFrom),particlesTo(particlesTo)
    {
        Point<particlesFrom_type::dims, unsigned int> p;
        p.zero();
        dcpse = new Dcpse<particlesFrom_type::dims, particlesFrom_type,particlesTo_type>(particlesFrom,particlesTo, p, ord, rCut, oversampling_factor, opt);
    }
 
    void deallocate() {
        delete (Dcpse<particlesFrom_type::dims, particlesFrom_type, particlesTo_type> *) dcpse;
    }
 
   /* template<typename operand_type>
    vector_dist_expression_op<operand_type, Dcpse_type<operand_type::vtype::dims, typename operand_type::vtype>, VECT_DCPSE>
    operator()(operand_type arg) {
        typedef Dcpse_type<operand_type::vtype::dims, typename operand_type::vtype> dcpse_type;
        return vector_dist_expression_op<operand_type, dcpse_type, VECT_DCPSE>(arg, *(dcpse_type *) dcpse);
    }*/
 
   template<unsigned int prp1,unsigned int prp2>
   void p2p() {
       auto dcpse_temp = (Dcpse<particlesFrom_type::dims, particlesFrom_type, particlesTo_type>*) dcpse;
       dcpse_temp->template p2p<prp1,prp2>();
 
   }
 
    // template<unsigned int prp, typename particles_type>
    // void DrawKernel(particles_type &particles, int k) {
    //     auto dcpse_temp = (Dcpse_type<particlesFrom_type::dims, particlesFrom_type, particlesTo_type> *) dcpse;
    //     dcpse_temp->template DrawKernel<prp>(particles, k);
 
    // }
 
    // template<unsigned int prp, typename particles_type>
    // void DrawKernelNN(particles_type &particles, int k) {
    //     auto dcpse_temp = (Dcpse_type<particlesFrom_type::dims, particlesFrom_type,particlesTo_type> *) dcpse;
    //     dcpse_temp->template DrawKernelNN<prp>(particles, k);
 
    // }
 
    // template<typename particles_type>
    // void checkMomenta(particles_type &particles) {
    //     auto dcpse_temp = (Dcpse_type<particles_type::dims, particlesFrom_type, particlesTo_type> *) dcpse;
    //     dcpse_temp->checkMomenta(particles);
 
    // }
 
    void update() {
        auto dcpse_temp = (Dcpse<particlesFrom_type::dims, particlesFrom_type, particlesTo_type> *) dcpse;
        dcpse_temp->initializeUpdate(particlesFrom,particlesTo);
 
    }
 
};
 
 
 
#endif //OPENFPM_PDATA_DCPSEINTERPOLATION_HPP