allvars.h

Go to the documentation of this file.

/*******************************************************************************
 * \copyright   This file is part of the GADGET4 N-body/SPH code developed
 * \copyright   by Volker Springel. Copyright (C) 2014-2020 by Volker Springel
 * \copyright   (vspringel@mpa-garching.mpg.de) and all contributing authors.
 *******************************************************************************/
#ifndef ALLVARS_H
#define ALLVARS_H
 
#include "gadgetconfig.h"
 
#include <math.h>
 
#include "../data/constants.h"
#include "../data/dtypes.h"
#include "../data/macros.h"
#include "../io/parameters.h"
 
struct global_data_all_processes : public parameters
{
#if defined(COOLING)
  char TreecoolFile[255];
#endif
 
#ifdef INDIVIDUAL_GRAVITY_SOFTENING
  double AvgType1Mass;
#endif
 
  double TopNodeFactor;
 
  int ICFormat; 
  int SnapFormat; 
  int NumFilesPerSnapshot;      
  int MaxFilesWithConcurrentIO; 
  double TreeAllocFactor; 
  double TopNodeAllocFactor; 
  double NgbTreeAllocFactor; 
  double ForeignNodeAllocFactor;
 
  double ForeignPointAllocFactor;
 
  enum restart_options RestartFlag; 
  int MaxMemSize; 
  /* some SPH parameters */
 
  int DesNumNgb; 
#ifdef LEAN
  MyDouble PartMass;
#endif
 
#ifdef TIMEDEP_ART_VISC
  double AlphaMin; 
#endif
 
#ifdef SUBFIND
  int DesLinkNgb;
#endif
 
  long long GlobalNSynchronizedHydro;
  long long GlobalNSynchronizedGravity;
 
  double MaxNumNgbDeviation; 
  double ArtBulkViscConst; 
  double InitGasTemp;      
  double InitGasU;         
  double MinEgySpec;       
  /* some force counters  */
 
  long long TotNumOfForces;     
  long long TotNumDirectForces; 
  long long TotNumDensity;      
  long long TotNumHydro;        
  /* various cosmological factors that are only a function of the current scale factor, and in non-comoving runs are set to 1 */
  double cf_atime, cf_atime2, cf_ainv, cf_a2inv, cf_a3inv, cf_afac1, cf_afac2, cf_afac3, cf_hubble_a, cf_atime_hubble_a,
      cf_atime2_hubble_a, cf_redshift;
  /* Hubble rate at the current time, valid both for comoving and non-comoving intergation */
  double cf_H;
 
  double accel_normalize_fac; /* used in I/O to normalize accelerations if reduced precision storage is used */
 
  /* system of units  */
  double UnitTime_in_s;            
  double UnitMass_in_g;            
  double UnitVelocity_in_cm_per_s; 
  double UnitLength_in_cm;         
  double UnitPressure_in_cgs;      
  double UnitDensity_in_cgs;       
  double UnitCoolingRate_in_cgs;   
  double UnitEnergy_in_cgs;        
  double UnitTime_in_Megayears;    
  double UnitTime_in_years;        
  double GravityConstantInternal;  
  double G;                        
  /* Cosmology */
 
  double Hubble;      
  double Omega0;      
  double OmegaLambda; 
  double OmegaBaryon; 
  double HubbleParam; 
  double BoxSize; 
  /* Code options */
 
  int ComovingIntegrationOn;  
  int TypeOfOpeningCriterion; 
  int OutputListOn;           
  int LowestActiveTimeBin;
  int HighestActiveTimeBin;
  int LowestOccupiedTimeBin;
  int HighestOccupiedTimeBin;
  int LowestOccupiedGravTimeBin;
  int HighestOccupiedGravTimeBin;
  int HighestSynchronizedTimeBin;
  int SmallestTimeBinWithDomainDecomposition;
  double ActivePartFracForNewDomainDecomp;
  /* parameters determining output frequency */
 
#if defined(TREEPM_NOTIMESPLIT) || defined(PLACEHIGHRESREGION)
  double ActivePartFracForPMinsteadOfEwald;
#endif
 
  int SnapshotFileCount;        
  double TimeBetSnapshot;       
  double TimeOfFirstSnapshot;   
  double CpuTimeBetRestartFile; 
  double TimeLastRestartFile;   
  double TimeBetStatistics;     
  double TimeLastStatistics;    
  int NumCurrentTiStep;         
  /* Current time of the simulation, global step, and end of simulation */
 
  double Time;      
  double TimeBegin; 
  double TimeStep;  
  double TimeMax;   
  double TimeOld;   
  /* variables for organizing discrete timeline */
 
  double Timebase_interval;  
  integertime Ti_Current;    
  integertime Ti_nextoutput; 
  integertime Ti_lastoutput;
 
#if defined(PMGRID) && !defined(TREEPM_NOTIMESPLIT)
  integertime PM_Ti_endstep, PM_Ti_begstep;
#endif
 
#if defined(EVALPOTENTIAL) && defined(PMGRID) && defined(PERIODIC)
  double TotalMass;
#endif
 
  inline double get_absolutetime_from_integertime(integertime ti)
  {
    if(ComovingIntegrationOn)
      return TimeBegin * exp(ti * Timebase_interval);
    else
      return TimeBegin + ti * Timebase_interval;
  }
 
  char DumpFlag_nextoutput;
 
  integertime Ti_begstep[TIMEBINS]; 
#ifdef FORCE_EQUAL_TIMESTEPS
  integertime GlobalTimeStep;
#endif
 
  /* variables that keep track of CPU consumption */
 
  double TimeLimitCPU;
 
  double CPUForLastPMExecution;
 
  /* tree code opening criterion */
 
  double ErrTolTheta;    
  double ErrTolThetaMax; 
  double ErrTolForceAcc; 
  char RelOpeningCriterionInUse; 
  /* adjusts accuracy of time-integration */
 
  double ErrTolIntAccuracy; 
  double MinSizeTimestep; 
  double MaxSizeTimestep; 
  double CourantFac; 
  int CPU_TimeBinCountMeasurements[TIMEBINS];
  double CPU_TimeBinMeasurements[TIMEBINS][NUMBER_OF_MEASUREMENTS_TO_RECORD];
 
  /* gravitational and hydrodynamical softening lengths (given in terms of an `equivalent' Plummer softening
   * length)
   *
   * five groups of particles are supported 0=gas,1=halo,2=disk,3=bulge,4=stars
   */
 
  int SofteningClassOfPartType[NTYPES];
 
  double SofteningComoving[NSOFTCLASSES]; 
  double SofteningMaxPhys[NSOFTCLASSES];  
  double SofteningTable[NSOFTCLASSES +
                        NSOFTCLASSES_HYDRO]; 
  double ForceSoftening[NSOFTCLASSES + NSOFTCLASSES_HYDRO +
                        2]; 
#ifdef ADAPTIVE_HYDRO_SOFTENING
  double MinimumComovingHydroSoftening;
  double AdaptiveHydroSofteningSpacing;
  double GasSoftFactor;
#endif
 
  double MassTable[NTYPES];
 
  /* some filenames */
  char InitCondFile[MAXLEN_PATH], OutputDir[MAXLEN_PATH], SnapshotFileBase[MAXLEN_PATH], OutputListFilename[MAXLEN_PATH];
 
  double OutputListTimes[MAXLEN_OUTPUTLIST];
  char OutputListFlag[MAXLEN_OUTPUTLIST];
  int OutputListLength; 
#ifdef SECOND_ORDER_LPT_ICS
  double LptScalingfactor;
#endif
 
#if defined(PMGRID) && !defined(TREEPM_NOTIMESPLIT)
  double DtDisplacement;
#endif
 
#ifdef LIGHTCONE
 
#ifdef LIGHTCONE_PARTICLES
  char LightConeDefinitionFile[MAXLEN_PATH];
  int LightconeFileCount;
#endif
 
#ifdef LIGHTCONE_MASSMAPS
  int LightConeMassMapsNside;
  double LightConeMassMapThickness;
  double LightConeMassMapMaxRedshift;
  int CurrentMassMapBoundary;
#endif
 
  int LightConeImageConeNr;
  double LightConeImageLengthX;
  double LightConeImageLengthY;
  double LightConeImageCornerX;
  double LightConeImageCornerY;
  int LightConeImagePixelsX;
  int LightConeImagePixelsY;
  char LightConeImagePicName[MAXLEN_PATH];
  int LightConeImageFirstConeDir;
  int LightConeImageLastConeDir;
#endif
 
#ifdef STARFORMATION /* star formation and feedback sector */
  double CritOverDensity;
  double CritPhysDensity;
  double OverDensThresh;
  double PhysDensThresh;
  double EgySpecSN;
  double EgySpecCold;
  double FactorEVP;
  double TempSupernova;
  double TempClouds;
  double MaxSfrTimescale;
  double FactorSN;
  MyIDType MaxID;
#endif
 
#ifdef REDUCE_FLUSH
  double FlushCpuTimeDiff;
  double FlushLast;
#endif
 
#ifdef NGENIC
  int NSample;
  int SphereMode;
  int PowerSpectrumType;
  int ReNormalizeInputSpectrum;
  double PrimordialIndex;
  double ShapeGamma;
  double Sigma8;
  char PowerSpectrumFile[MAXLEN_PATH];
  double InputSpectrum_UnitLength_in_cm;
  int NgenicSeed;
#endif
 
#ifdef CREATE_GRID
  int GridSize;
#endif
 
#ifdef EXTERNALGRAVITY_STATICHQ
  double A_StaticHQHalo;
  double Mass_StaticHQHalo;
#endif
 
  void set_cosmo_factors_for_current_time(void);
  void register_parameters(void);
  void read_outputlist(char *fname);
  void some_parameter_checks(void);
 
  inline char *get_data_ptr(void) { return (char *)this + sizeof(parameters); }
 
  inline size_t get_data_size(void) { return sizeof(global_data_all_processes) - sizeof(parameters); }
};
 
extern global_data_all_processes All;
 
#endif