SphUtils.h

/*
 * This file contains inline utility functions and function declarations needed
 * for SPH
 * 
 * Based on the physics modules implemented by James Wadsley in GASOLINE.
 * Implemented by Tom Quinn in ChaNGa and refactored by Isaac Backus
 * Oct 14, 2016
 */
#ifndef SPHUTILS_H
#define SPHUTILS_H

// Include Sph.h for possible macro definitions in there
#include "Sph.h"

inline double presPdv(const double a, const double b){
    return a;
}
inline double presAcc(const double a, const double b){
    return a + b;
}

inline double switchCombine(GravityParticle *a, GravityParticle *b) {
#ifdef CULLENALPHA
    return 1.0;
#else
    return 0.5*(a->BalsaraSwitch()+b->BalsaraSwitch());
#endif
}

inline double rhoDivv(const double a, const double b) {
#ifdef GDFORCE
    return b;
#else
    return a;
#endif
}

inline double varAlpha(const double alpha, const GravityParticle *a, 
                       const GravityParticle *b) {
#if defined(VARALPHA) || defined(CULLENALPHA)
    return alpha*0.5*(a->CullenAlpha() + b->CullenAlpha());
#else
    return alpha;
#endif //VARALPHA
}

inline double varBeta(const double beta, GravityParticle *a, 
                      const GravityParticle *b) {
#if defined(VARALPHA) || defined(CULLENALPHA)
    return beta*0.5*(a->CullenAlpha() + b->CullenAlpha());
#else
    return beta;
#endif //VARALPHA
}

inline bool diffusionLimitTest(const double diffSum, const double dTime, 
                               GravityParticle *a, GravityParticle *b) {
#ifdef FEEDBACKDIFFLIMIT
    return (diffSum == 0 \
            || dTime < a->fTimeCoolIsOffUntil() \
            || dTime < b->fTimeCoolIsOffUntil());
#else
    return (diffSum == 0);
#endif
}

inline double massDiffFac(const GravityParticle *p) {
//#ifdef MASSDIFF /* compile-time flag */
//    return p->fMass;
//#else
    return 1.0;
//#endif //MASSDIFF
}

typedef struct PressSmoothUpdateStruct {
    double dvdotdr;
    double visc;
    double aFac;
    Vector3D<double> dx;
//    /* not implemented */
//    #ifdef DRHODT /* compile-time flag */
//        double fDivv_Corrector;
//    #endif
    #ifdef DIFFUSION
//        /* DIFFUSIONPRICE not implemented */
//        #ifdef DIFFUSIONPRICE /* compile-time flag */
//            double diffu;
//        #else
            #ifndef NODIFFUSIONTHERMAL /* compile-time flag */
                double diffu;
            #endif
//        #endif //DIFFUSIONPRICE
        /* not implemented */
//        #ifdef UNONCOOL
//            double diffuNc;
//        #endif
        double diffMetals;
        double diffMetalsOxygen;
        double diffMetalsIron;
//        /* not implemented */
//        #ifdef MASSDIFF /* compile-time flag */
//            double diffMass;
//            Vector3D<double> diffVelocity;
//        #endif
    #endif
} PressSmoothUpdate;

typedef struct PressSmoothParticleStruct {
    double rNorm;
    double PoverRho2;
    double PoverRho2f;
} PressSmoothParticle;

void updateParticle(GravityParticle *a, GravityParticle *b, 
                    PressSmoothUpdate *params, PressSmoothParticle *aParams, 
                    PressSmoothParticle *bParams, int sign);

#endif // SPHUTILS_H