parameter_backup.h

Go to the documentation of this file.

 
#ifndef MILUPHPC_PARAMETER_H
#define MILUPHPC_PARAMETER_H
 
#include <limits>
#include <iostream>
 
#ifdef SINGLE_PRECISION
    typedef float real;
#else
    typedef double real;
#endif
typedef int integer;
typedef unsigned long keyType;
typedef int idInteger;
 
#define MAX_LEVEL 21
 
#define UNIT_TESTING 0
 
#define DEBUGGING 0
 
#define SAFETY_LEVEL 2
 
#define DIM 3
#define power_two(x) (1 << (x))
#define POW_DIM power_two(DIM)
 
#define SI_UNITS 1
 
#define CUBIC_DOMAINS 1
 
#define GRAVITY_SIM 1
 
#define SPH_SIM 1
 
#define INTEGRATE_ENERGY 0
 
#define INTEGRATE_DENSITY 1
 
#define INTEGRATE_SML 0
 
#define DECOUPLE_SML 0
 
#define VARIABLE_SML 1
 
#define SML_CORRECTION 0
 
#define SPH_EQU_VERSION 1
 
#define ARTIFICIAL_VISCOSITY 1
 
#define AVERAGE_KERNELS 0
#define DEAL_WITH_TOO_MANY_INTERACTIONS 0
#define SHEPARD_CORRECTION 0
#define SOLID 0
#define NAVIER_STOKES 0
#define ARTIFICIAL_STRESS 0
#define POROSITY 0
#define ZERO_CONSISTENCY 0
#define LINEAR_CONSISTENCY 0
#define FRAGMENTATION 0
#define PALPHA_POROSITY 0
#define PLASTICITY 0
#define KLEY_VISCOSITY 0
 
#define KEY_MAX ULONG_MAX
//#define DOMAIN_LIST_SIZE 512 // changed to be a runtime constant
#define MAX_DEPTH 128
#define MAX_NUM_INTERACTIONS 180
#define NUM_THREADS_LIMIT_TIME_STEP 256
#define NUM_THREADS_CALC_CENTER_OF_MASS 256
 
// (note that our sml is defined up to the zero of the kernel, not half of it)
#define COURANT_FACT 0.4
 
#define FORCES_FACT 0.2
 
constexpr real dbl_max = std::numeric_limits<real>::max();
#define DBL_MAX dbl_max;
 
namespace Constants {
    constexpr real G = 6.67430e-11;
}
 
typedef struct SimulationParameters {
    std::string directory;
    std::string logDirectory;
    int verbosity;
    bool timeKernels;
    int numOutputFiles;
    real timeStep;
    real maxTimeStep;
    real timeEnd;
    bool loadBalancing;
    int loadBalancingInterval;
    int loadBalancingBins;
    std::string inputFile;
    std::string materialConfigFile;
    int outputRank;
    bool performanceLog;
    bool particlesSent2H5;
    int sfcSelection;
    int integratorSelection;
//#if GRAVITY_SIM
    real theta;
    real smoothing;
    int gravityForceVersion;
//#endif
//#if SPH_SIM
    int smoothingKernelSelection;
    int sphFixedRadiusNNVersion;
//#endif
    bool removeParticles;
    int removeParticlesCriterion;
    real removeParticlesDimension;
    int bins;
    bool calculateAngularMomentum;
    bool calculateEnergy;
    bool calculateCenterOfMass;
    real particleMemoryContingent;
    int domainListSize;
} SimulationParameters;
 
struct To
{
    enum Target
    {
        host, device
    };
    Target t_;
    To(Target t) : t_(t) {}
    operator Target () const {return t_;}
private:
    template<typename T>
    operator T () const;
};
 
struct Smoothing
{
    enum Kernel
    {
        spiky, cubic_spline, wendlandc2, wendlandc4, wendlandc6
    };
    Kernel t_;
    Smoothing(Kernel t) : t_(t) {}
    operator Smoothing () const {return t_;}
private:
    template<typename T>
    operator T () const;
};
 
struct Execution
{
    enum Location
    {
        host, device
    };
    Location t_;
    Execution(Location t) : t_(t) {}
    operator Location () const {return t_;}
private:
    template<typename T>
    operator T () const;
};
 
struct Curve
{
    enum Type
    {
        lebesgue, hilbert
    };
    Type t_;
    Curve(Type t) : t_(t) {}
    operator Type () const {return t_;}
    //friend std::ostream& operator<<(std::ostream& out, const Curve::Type curveType);
private:
    template<typename T>
    operator T () const;
};
 
struct IntegratorSelection
{
    enum Type
    {
        explicit_euler, predictor_corrector_euler
    };
    Type t_;
    IntegratorSelection(Type t) : t_(t) {}
    operator Type () const {return t_;}
private:
    template<typename T>
    operator T () const;
};
 
 
enum EquationOfStates {
    //EOS_TYPE_ACCRETED = -2, // special flag for particles that got accreted by a gravitating point mass
    //EOS_TYPE_IGNORE = -1, // particle is ignored
    EOS_TYPE_POLYTROPIC_GAS = 0,
    //EOS_TYPE_MURNAGHAN = 1, // Murnaghan EOS for solid bodies, see Melosh "Impact Cratering", needs in material.cfg: rho_0, bulk_modulus, n
    //EOS_TYPE_TILLOTSON = 2, // Tillotson EOS for solid bodies, see Melosh "Impact Cratering", needs in material.cfg: till_rho_0, till_A, till_B, till_E_0, till_E_iv, till_E_cv, till_a, till_b, till_alpha, till_beta; bulk_modulus and shear_modulus are needed to calculate the sound speed and crack growth speed for FRAGMENTATION
    EOS_TYPE_ISOTHERMAL_GAS = 3,
    //EOS_TYPE_REGOLITH = 4, // The Bui et al. 2008 soil model
    //EOS_TYPE_JUTZI = 5, // Tillotson EOS with p-alpha model by Jutzi et al.
    //EOS_TYPE_JUTZI_MURNAGHAN = 6, // Murnaghan EOS with p-alpha model by Jutzi et al.
    //EOS_TYPE_ANEOS = 7, // ANEOS (or tabulated EOS in ANEOS format)
    //EOS_TYPE_VISCOUS_REGOLITH = 8, // describe regolith as a viscous material -> EXPERIMENTAL DO NOT USE
    EOS_TYPE_IDEAL_GAS = 9,
    //EOS_TYPE_SIRONO = 10, // Sirono EOS modifed by Geretshauser in 2009/10
    //EOS_TYPE_EPSILON = 11, // Tillotson EOS with epsilon-alpha model by Wuennemann, Collins et al.
    EOS_TYPE_LOCALLY_ISOTHERMAL_GAS = 12,
    //EOS_TYPE_JUTZI_ANEOS = 13// ANEOS EOS with p-alpha model by Jutzi et al.
};
 
struct Entry
{
    enum Name
    {
        x,
#if DIM > 1
        y,
#if DIM == 3
        z,
#endif
#endif
        mass
    };
    Name t_;
    Entry(Name t) : t_(t) {}
    operator Name () const {return t_;}
private:
    template<typename T>
    operator T () const;
};
 
#endif //MILUPHPC_PARAMETER_H