DcpseInterpolation.hpp

//
// Created by Abhinav Singh on 03.11.21.
// Surface interpolation: Lennt Schulze and Ale Foggia on 06.03.24
//
 
#ifndef OPENFPM_PDATA_DCPSEINTERPOLATION_HPP
#define OPENFPM_PDATA_DCPSEINTERPOLATION_HPP
#include "DCPSE/Dcpse.hpp"
 
template<typename particlesSupport_type, typename particlesDomain_type, typename VerletList_type, size_t NORMAL_ID = INT_MAX>
class PPInterpolation 
{
 
    void *dcpse;
 
    particlesSupport_type & particlesSupport;
    particlesDomain_type & particlesDomain;
 
 
public:
    PPInterpolation(
        particlesSupport_type &particlesSupport,
        particlesDomain_type &particlesDomain,
        VerletList_type& verletList,
        unsigned int ord,
        typename particlesSupport_type::stype rCut,
        support_options opt = support_options::RADIUS
    ):
        particlesSupport(particlesSupport),
        particlesDomain(particlesDomain)
    {
        Point<particlesSupport_type::dims, unsigned int> p; p.zero();
 
        dcpse = new Dcpse<
            particlesSupport_type::dims,
            VerletList_type,
            particlesSupport_type,
            particlesDomain_type>
        (
            particlesSupport,
            particlesDomain,
            verletList,
            p,
            ord,
            rCut,
            opt
        );
    }
 
    PPInterpolation(
        particlesSupport_type &particlesSupport,
        particlesDomain_type &particlesDomain,
        VerletList_type& verletList,
        unsigned int ord,
        typename particlesSupport_type::stype rCut,
        typename particlesSupport_type::stype nSpacing,
        bool isSurfaceInterpolation,
        support_options opt = support_options::RADIUS
    ):
        particlesSupport(particlesSupport),
        particlesDomain(particlesDomain)
    {
        Point<particlesSupport_type::dims, unsigned int> p;
        p.zero();
        dcpse = new SurfaceDcpse<particlesSupport_type::dims,VerletList_type,particlesSupport_type,particlesDomain_type>(particlesSupport, particlesDomain, verletList, p, ord, rCut, nSpacing, static_cast<unsigned int>(rCut/nSpacing), value_t<NORMAL_ID>(), opt);
    }
 
    void deallocate() {
        delete (Dcpse<particlesSupport_type::dims, VerletList_type, particlesSupport_type, particlesDomain_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 propSupport,unsigned int propDomain>
   void p2p() {
       auto dcpse_temp = (Dcpse<particlesSupport_type::dims, VerletList_type, particlesSupport_type, particlesDomain_type>*) dcpse;
       dcpse_temp->template p2p<propSupport,propDomain>();
   }
 
  // foggia 16.09.24
  template<unsigned int propSupport,unsigned int propDomain, size_t N1>
  void p2p() {
    auto dcpse_temp = (Dcpse<particlesSupport_type::dims, VerletList_type, particlesSupport_type, particlesDomain_type>*) dcpse;
    dcpse_temp->template p2p<propSupport,propDomain,N1>();
  }
 
    // template<unsigned int prp, typename particles_type>
    // void DrawKernel(particles_type &particles, int k) {
    //     auto dcpse_temp = (Dcpse_type<particlesSupport_type::dims, VerletList_type, particlesSupport_type, particlesDomain_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<particlesSupport_type::dims, particlesSupport_type,particlesDomain_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, particlesSupport_type, particlesDomain_type> *) dcpse;
    //     dcpse_temp->checkMomenta(particles);
 
    // }
 
    void update() {
        auto dcpse_temp = (Dcpse<particlesSupport_type::dims, VerletList_type, particlesSupport_type, particlesDomain_type> *) dcpse;
        dcpse_temp->initializeUpdate(particlesSupport,particlesDomain);
    }
};
 
 
 
#endif //OPENFPM_PDATA_DCPSEINTERPOLATION_HPP