SIDM.h

class SIDMSmoothParams : public SmoothParams
{
  protected:
    double a, H; // Cosmological parameters
        double dTime;
        double dDelta;
        double dSIDMSigma;
        double dSIDMVariable;
        double iSIDMSelect;
    virtual void fcnSmooth(GravityParticle *p, int nSmooth, pqSmoothNode *nList);
        virtual int isSmoothActive(GravityParticle *p);
        virtual void initTreeParticle(GravityParticle *p);
        virtual void postTreeParticle(GravityParticle *p) {}
        virtual void initSmoothParticle(GravityParticle *p);
        virtual void initSmoothCache(GravityParticle *p);
        virtual void combSmoothCache(GravityParticle *p1,ExternalSmoothParticle *p2);

  public:
    SIDMSmoothParams() {} //empty constructor
    SIDMSmoothParams(int _iType, int am, CSM csm, double _dTime, double _dSIDMSigma, double _dSIDMVariable, int _iSIDMSelect, double _dDelta) {
    iType = _iType;
    activeRung = am;
    dTime = _dTime;
    dDelta = _dDelta;
    dSIDMSigma= _dSIDMSigma;
    dSIDMVariable= _dSIDMVariable;
    iSIDMSelect= _iSIDMSelect;

        if(csm->bComove) {
            H = csmTime2Hub(csm,dTime);
            a = csmTime2Exp(csm,dTime);
            }
        else {
            H = 0.0;
            a = 1.0;
            }

    }
    PUPable_decl(SIDMSmoothParams);
    SIDMSmoothParams(CkMigrateMessage *m) : SmoothParams(m) {}
    virtual void pup(PUP::er &p) {
        SmoothParams::pup(p);//Call  base class
        p|dTime;
        p|dDelta;
        p|dSIDMSigma;
        p|dSIDMVariable;
        p|iSIDMSelect;
        p|a;
        p|H;
        }
    };