Reference/internal-api-ref/grid3D_8h_source.html

#ifndef GRID3D_H

#define GRID3D_H


#ifdef MPI_CHOLLA

  #include "../mpi/mpi_routines.h"

#endif /*MPI_CHOLLA*/


#include <stdio.h>


#include <functional>


#include "../global/global.h"  // declares Parameter and forward-declares ParameterMap

#include "../global/global_cuda.h"

#include "../grid/field_info.h"

#include "../io/FnameTemplate.h"

#include "../model/model_collection.h"


#ifdef HDF5

  #include <hdf5.h>

#endif


#ifdef GRAVITY

  #include "../gravity/grav3D.h"

#endif


#ifdef PARTICLES

  #include "../particles/particles_3D.h"

#endif


#ifdef COSMOLOGY

  #include "../cosmology/cosmology.h"

#endif


#ifdef COOLING_GRACKLE

  #include "../cooling_grackle/cool_grackle.h"

#endif


#ifdef CPU_TIME

  #include "../utils/timing_functions.h"

#endif


#ifdef CHEMISTRY_GPU

  #include "chemistry_gpu/chemistry_gpu.h"

#endif


#ifdef ANALYSIS

  #include "../analysis/analysis.h"

#endif


class AttrRecorderInterface;


struct Header {

  int n_cells;


  int n_ghost;


  int nx;


  int ny;


  int nz;


  int nx_real;


  int ny_real;


  int nz_real;


  /*  \brief Global domain x-direction minimum */

  Real xbound;


  /*  \brief Global domain y-direction minimum */

  Real ybound;


  /*  \brief Global domain z-direction minimum */

  Real zbound;


  /*  \brief Local domain x-direction minimum */

  Real xblocal;


  /*  \brief Local domain y-direction minimum */

  Real yblocal;


  /*  \brief Local domain z-direction minimum */

  Real zblocal;


  /*  \brief Local domain x-direction maximum */

  Real xblocal_max;


  /*  \brief Local domain y-direction maximum */

  Real yblocal_max;


  /*  \brief Local domain z-direction maximum */

  Real zblocal_max;


  /*  \brief Global domain length in x-direction */

  Real xdglobal;


  /*  \brief Global domain length in y-direction */

  Real ydglobal;


  /*  \brief Global domain length in z-direction */

  Real zdglobal;


  Real dx;


  Real dy;


  Real dz;


  Real t;


  Real dt;


  Real min_dt_slow;


  Real t_wall;


  int n_step;


  int n_fields;


  int custom_grav;


  // Values for lower limit for density and temperature

  Real temperature_floor;

  Real density_floor;

  Real scalar_floor;


  Real Ekin_avrg;


  // Flag to indicate when to transfer the Conserved boundaries

  bool TRANSFER_HYDRO_BOUNDARIES;


  // Parameters For Spherical Colapse Problem

  Real sphere_density;

  Real sphere_radius;

  Real sphere_background_density;

  Real sphere_center_x;

  Real sphere_center_y;

  Real sphere_center_z;


  // only meaningful when GRAVITY and GRAVITY_ANALYTIC_COMP are defined

  bool gas_only_use_static_grav;


#ifdef GRAVITY

  int n_ghost_potential_offset;

#endif


#ifdef COSMOLOGY

  bool OUTPUT_SCALE_FACOR;

#endif


  bool Output_Now;

  bool Output_Initial;


  bool Output_Complete_Data;


#ifdef SCALAR

  #ifdef DUST

  Real grain_radius;

  #endif

#endif

};


class Grid3D

{

 public:

  int flag_init;


  struct Header H;


  FieldInfo field_info;


  ModelCollection model_collection;


#ifdef GRAVITY

  // Object that contains data for gravity

  Grav3D Grav;

#endif


#ifdef PARTICLES

  // Object that contains data for particles

  Particles3D Particles;

#endif


#ifdef COSMOLOGY

  // Object that contains data for cosmology

  Cosmology Cosmo;

#endif


#ifdef COOLING_GRACKLE

  // Object that contains data for Grackle cooling

  Cool_GK Cool;

#endif


#ifdef CPU_TIME

  Time Timer;

#endif


#ifdef CHEMISTRY_GPU

  // Object that contains data for the GPU chemistry solver

  Chem_GPU Chem;

#endif


#ifdef ANALYSIS

  AnalysisModule Analysis;

#endif


  struct Conserved {

    Real *host;


    Real *density;


    Real *momentum_x;


    Real *momentum_y;


    Real *momentum_z;


    Real *Energy;


#ifdef SCALAR

    Real *scalar;

  #ifdef BASIC_SCALAR

    Real *basic_scalar;

  #endif

  #ifdef DUST

    Real *dust_density;

  #endif

#endif  // SCALAR


#ifdef MHD

    Real *magnetic_x;


    Real *magnetic_y;


    Real *magnetic_z;

#endif  // MHD


#ifdef DE

    Real *GasEnergy;

#endif  // DE


    Real *Grav_potential;


#ifdef CHEMISTRY_GPU

    Real *HI_density;

    Real *HII_density;

    Real *HeI_density;

    Real *HeII_density;

    Real *HeIII_density;

    Real *e_density;

#endif


    Real *device;

    Real *d_density, *d_momentum_x, *d_momentum_y, *d_momentum_z, *d_Energy, *d_scalar, *d_basic_scalar,

        *d_dust_density, *d_magnetic_x, *d_magnetic_y, *d_magnetic_z, *d_GasEnergy;


    Real *d_Grav_potential;

  } C;


  Grid3D(void);


  void Initialize(struct Parameters *P);


  void AllocateMemory(void);


  void Set_Initial_Conditions(Parameters P, const ParameterMap &pmap);


  void Get_Position(long i, long j, long k, Real *xpos, Real *ypos, Real *zpos) const;


  Real Calc_Inverse_Timestep();


  void Set_Domain_Properties(struct Parameters P);


  void set_dt(Real dti);


#ifdef GRAVITY

  void set_dt_Gravity();

#endif


  void Execute_Hydro_Integrator(void);


  Real Update_Hydro_Grid(std::function<void(Grid3D &)> &feedback_callback,

                         std::function<void(Grid3D &)> &chemistry_callback);


  void Update_Time();


  void Write_Header(AttrRecorderInterface &attr_recorder) const;


  void Read_Grid(struct Parameters P);


  void Read_Grid_Cat(struct Parameters P);


#ifdef HDF5

  void Read_Grid_HDF5(hid_t file_id, struct Parameters P);

#endif


#if defined(PRINT_INITIAL_STATS) && defined(COSMOLOGY)

  void Print_Grid_Stats(void);

#endif


  void Reset(void);


  void FreeMemory(void);


  void Constant(Parameters const &P);


  void Sound_Wave(Parameters const &P);


  void Linear_Wave(Parameters const &P);


  void Square_Wave(Parameters const &P);


  void Riemann(Parameters const &P);


  void Shu_Osher();


  void Blast_1D();


  void KH();


  void KH_res_ind();


  void Rayleigh_Taylor();


  void Gresho();


  void Implosion_2D();


  void Explosion_2D();


  void Noh_2D();


  void Noh_3D();


  void Disk_2D();


  void Disk_3D(Parameters P);


  void Set_Boundary_Conditions(Parameters P);


  void Set_Boundary_Conditions_Grid(Parameters P);


  int Check_Custom_Boundary(int *flags, struct Parameters P);


  void Set_Boundaries(int dir, int flags[]);


  void Set_Boundary_Extents(int dir, int *imin, int *imax);


  void Custom_Boundary(char bcnd[MAXLEN]);


  void Wind_Boundary();


  void Noh_Boundary();


  void Spherical_Overpressure_3D();


  void Spherical_Overdensity_3D();


  void Clouds();


  void Uniform_Grid();


  void Zeldovich_Pancake(struct Parameters P);


  void Adiabatic_Expansion(struct Parameters P);


  void Chemistry_Test(struct Parameters P);


#ifdef MHD

  void Circularly_Polarized_Alfven_Wave(struct Parameters const P);


  void Advecting_Field_Loop(struct Parameters const P);


  void MHD_Spherical_Blast(struct Parameters const P);


  void Orszag_Tang_Vortex();

#endif  // MHD


#ifdef MPI_CHOLLA

  void Set_Boundaries_MPI(struct Parameters P);

  void Set_Boundaries_MPI_BLOCK(int *flags, struct Parameters P);

  void Load_and_Send_MPI_Comm_Buffers(int dir, int *flags);

  void Wait_and_Unload_MPI_Comm_Buffers(int dir, int *flags);

  void Unload_MPI_Comm_Buffers(int index);


  int Load_Hydro_DeviceBuffer_X0(Real *buffer);

  int Load_Hydro_DeviceBuffer_X1(Real *buffer);

  int Load_Hydro_DeviceBuffer_Y0(Real *buffer);

  int Load_Hydro_DeviceBuffer_Y1(Real *buffer);

  int Load_Hydro_DeviceBuffer_Z0(Real *buffer);

  int Load_Hydro_DeviceBuffer_Z1(Real *buffer);


  void Unload_Hydro_DeviceBuffer_X0(Real *buffer);

  void Unload_Hydro_DeviceBuffer_X1(Real *buffer);

  void Unload_Hydro_DeviceBuffer_Y0(Real *buffer);

  void Unload_Hydro_DeviceBuffer_Y1(Real *buffer);

  void Unload_Hydro_DeviceBuffer_Z0(Real *buffer);

  void Unload_Hydro_DeviceBuffer_Z1(Real *buffer);

#endif /*MPI_CHOLLA*/


  void Initialize_Models(ParameterMap &pmap) { model_collection = ModelCollection(pmap); }


  const ModelCollection &models() const { return model_collection; }


#ifdef GRAVITY

  void Initialize_Gravity(struct Parameters *P);

  void Compute_Gravitational_Potential(struct Parameters *P);

  void Copy_Hydro_Density_to_Gravity_Function(int g_start, int g_end);

  void Copy_Hydro_Density_to_Gravity();

  void Extrapolate_Grav_Potential_Function(int g_start, int g_end);

  void Extrapolate_Grav_Potential();

  void Set_Potential_Boundaries_Periodic(int direction, int side, int *flags);

  int Load_Gravity_Potential_To_Buffer(int direction, int side, Real *buffer, int buffer_start);

  void Unload_Gravity_Potential_from_Buffer(int direction, int side, Real *buffer, int buffer_start);

  void Set_Potential_Boundaries_Isolated(int direction, int side, int *flags);

  void Compute_Potential_Boundaries_Isolated(int dir, struct Parameters *P);

  void Compute_Potential_Isolated_Boundary(int direction, int side, int bc_potential_type);

  #ifdef SOR

  void Get_Potential_SOR(Real Grav_Constant, Real dens_avrg, Real current_a, struct Parameters *P);

  int Load_Poisson_Boundary_To_Buffer(int direction, int side, Real *buffer);

  void Unload_Poisson_Boundary_From_Buffer(int direction, int side, Real *buffer_host);

  #endif

  #ifdef GRAVITY_GPU

  void Copy_Hydro_Density_to_Gravity_GPU();

  void Extrapolate_Grav_Potential_GPU();

  int Load_Gravity_Potential_To_Buffer_GPU(int direction, int side, Real *buffer, int buffer_start);

  void Unload_Gravity_Potential_from_Buffer_GPU(int direction, int side, Real *buffer, int buffer_start);

  void Set_Potential_Boundaries_Isolated_GPU(int direction, int side, int *flags);

  void Set_Potential_Boundaries_Periodic_GPU(int direction, int side, int *flags);

  #endif


#endif  // GRAVITY


#ifdef GRAVITY_ANALYTIC_COMP

  void Add_Analytic_Potential();

  void Add_Analytic_Potential(int g_start, int g_end);

  void Setup_Analytic_Potential(struct Parameters *P);

  void Setup_Analytic_Galaxy_Potential(int g_start, int g_end, const DiskGalaxy &gal);

  #ifdef GRAVITY_GPU

  void Add_Analytic_Potential_GPU();

  #endif

#endif  // GRAVITY_ANALYTIC_COMP


#ifdef PARTICLES

  void Initialize_Particles(struct Parameters *P);

  void Initialize_Uniform_Particles();

  void Copy_Particles_Density_function(int g_start, int g_end);

  void Copy_Particles_Density();

  void Copy_Particles_Density_to_Gravity(struct Parameters P);

  void Set_Particles_Density_Boundaries_Periodic(int direction, int side);

  void Transfer_Particles_Boundaries(struct Parameters P);

  Real Update_Grid_and_Particles_KDK(struct Parameters P);

  void Set_Particles_Boundary(int dir, int side);

  #ifdef PARTICLES_CPU

  void Set_Particles_Open_Boundary_CPU(int dir, int side);

  #endif

  #ifdef MPI_CHOLLA

  int Load_Particles_Density_Boundary_to_Buffer(int direction, int side, Real *buffer);

  void Unload_Particles_Density_Boundary_From_Buffer(int direction, int side, Real *buffer);

  void Load_and_Send_Particles_X0(int ireq_n_particles, int ireq_particles_transfer);

  void Load_and_Send_Particles_X1(int ireq_n_particles, int ireq_particles_transfer);

  void Load_and_Send_Particles_Y0(int ireq_n_particles, int ireq_particles_transfer);

  void Load_and_Send_Particles_Y1(int ireq_n_particles, int ireq_particles_transfer);

  void Load_and_Send_Particles_Z0(int ireq_n_particles, int ireq_particles_transfer);

  void Load_and_Send_Particles_Z1(int ireq_n_particles, int ireq_particles_transfer);

  void Unload_Particles_from_Buffer_X0(int *flags);

  void Unload_Particles_from_Buffer_X1(int *flags);

  void Unload_Particles_from_Buffer_Y0(int *flags);

  void Unload_Particles_from_Buffer_Y1(int *flags);

  void Unload_Particles_from_Buffer_Z0(int *flags);

  void Unload_Particles_from_Buffer_Z1(int *flags);

  void Wait_NTransfer_and_Request_Recv_Particles_Transfer_BLOCK(int dir, int *flags);

  void Load_NTtransfer_and_Request_Receive_Particles_Transfer(int index, int *ireq_particles_transfer);

  void Wait_and_Unload_MPI_Comm_Particles_Buffers_BLOCK(int dir, int *flags);

  void Unload_Particles_From_Buffers_BLOCK(int index, int *flags);

  void Finish_Particles_Transfer();

  #endif  // MPI_CHOLLA

  void Transfer_Particles_Density_Boundaries(struct Parameters P);

  void Copy_Particles_Density_Buffer_Device_to_Host(int direction, int side, Real *buffer_d, Real *buffer_h);

  // void Transfer_Particles_Boundaries( struct Parameters P );

  void WriteData_Particles(struct Parameters P, int nfile, const FnameTemplate &fname_template);

  void OutputData_Particles(struct Parameters P, int nfile, const FnameTemplate &fname_template);

  void Load_Particles_Data(struct Parameters P);

  #ifdef HDF5


  void Write_Particles_Header(AttrRecorderInterface &attr_recorder);

  void Write_Particles_Data_HDF5(hid_t file_id);

  void Load_Particles_Data_HDF5(hid_t file_id, int nfile);

  #endif  // HDF5

  void Get_Gravity_Field_Particles_function(int g_start, int g_end);

  void Get_Gravity_Field_Particles();

  void Get_Gravity_CIC_function(part_int_t p_start, part_int_t p_end);

  void Get_Gravity_CIC();

  void Advance_Particles_KDK_Step1();

  void Advance_Particles_KDK_Step2();

  void Advance_Particles_KDK_Step1_function(part_int_t p_start, part_int_t p_end);

  void Advance_Particles_KDK_Step2_function(part_int_t p_start, part_int_t p_end);

  void Get_Particles_Acceleration();

  void Advance_Particles(int N_KDK_step);

  Real Calc_Particles_dt_function(part_int_t p_start, part_int_t p_end);

  Real Calc_Particles_dt();

  #ifdef PARTICLES_GPU

  Real Calc_Particles_dt_GPU();

  void Advance_Particles_KDK_Step1_GPU();

  void Advance_Particles_KDK_Step2_GPU();

  void Set_Particles_Boundary_GPU(int dir, int side);

  void Set_Particles_Density_Boundaries_Periodic_GPU(int direction, int side);

  int Load_Particles_Density_Boundary_to_Buffer_GPU(int direction, int side, Real *buffer);

  void Unload_Particles_Density_Boundary_From_Buffer_GPU(int direction, int side, Real *buffer);

  #endif  // PARTICLES_GPU

  #ifdef GRAVITY_GPU

  void Copy_Potential_From_GPU();

  void Copy_Particles_Density_to_GPU();

  void Copy_Particles_Density_GPU();

  #endif  // GRAVITY_GPU

#endif    // PARTICLES


#ifdef COSMOLOGY

  void Initialize_Cosmology(struct Parameters *P);

  void Change_DM_Frame_System(bool forward);

  void Change_GAS_Frame_System(bool forward);

  void Change_GAS_Frame_System_GPU(bool forward);

  void Change_Cosmological_Frame_System(bool forward);

  void Advance_Particles_KDK_Cosmo_Step1_function(part_int_t p_start, part_int_t p_end);

  void Advance_Particles_KDK_Cosmo_Step2_function(part_int_t p_start, part_int_t p_end);

  Real Calc_Particles_dt_Cosmo_function(part_int_t p_start, part_int_t p_end);

  Real Calc_Particles_dt_Cosmo();

  #ifdef PARTICLES_GPU

  void Advance_Particles_KDK_Cosmo_Step1_GPU();

  void Advance_Particles_KDK_Cosmo_Step2_GPU();

  #endif  // PARTICLES_GPU

#endif    // COSMOLOGY


#ifdef COOLING_GRACKLE

  void Initialize_Grackle(struct Parameters *P);

  void Allocate_Memory_Grackle();

  void Initialize_Fields_Grackle();

  void Copy_Fields_To_Grackle_function(int g_start, int g_end);

  void Copy_Fields_To_Grackle();

  void Update_Internal_Energy_function(int g_start, int g_end);

  void Update_Internal_Energy();

  void Do_Cooling_Step_Grackle();

#endif


#ifdef CHEMISTRY_GPU

  void Initialize_Chemistry(struct Parameters *P);

  void Compute_Gas_Temperature(Real *temperature, bool convert_cosmo_units);

  void Update_Chemistry();

#endif


#ifdef ANALYSIS

  void Initialize_AnalysisModule(struct Parameters *P);

  void Compute_and_Output_Analysis(struct Parameters *P);

  void Output_Analysis(struct Parameters *P);

  void Write_Analysis_Header_HDF5(hid_t file_id);

  void Write_Analysis_Data_HDF5(hid_t file_id);


  #ifdef PHASE_DIAGRAM

  void Compute_Phase_Diagram();

  #endif


  #ifdef LYA_STATISTICS

  void Populate_Lya_Skewers_Local(int axis);

  void Compute_Transmitted_Flux_Skewer(int skewer_id, int axis);

  void Compute_Lya_Statistics();

  void Compute_Flux_Power_Spectrum_Skewer(int skewer_id, int axis);

  void Initialize_Power_Spectrum_Measurements(int axis);

    #ifdef OUTPUT_SKEWERS

  void Output_Skewers_File(struct Parameters *P);

  void Write_Skewers_Header_HDF5(hid_t file_id);

  void Write_Skewers_Data_HDF5(hid_t file_id);

    #endif

  #endif  // LYA_STATISTICS

#endif    // ANALYSIS

};


// typedef for Grid3D_PointerMemberFunction

typedef void (Grid3D::*Grid3D_PMF_UnloadHydroBuffer)(Real *);

typedef void (Grid3D::*Grid3D_PMF_UnloadGravityPotential)(int, int, Real *, int);

typedef void (Grid3D::*Grid3D_PMF_UnloadParticleDensity)(int, int, Real *);


#endif  // GRID3D_H

AttrRecorderInterface
Definition AttrRecorderInterface.h:36

DiskGalaxy
Definition disk_galaxy.h:39

FieldInfo
Definition field_info.h:49

FnameTemplate
Definition FnameTemplate.h:43

Grav3D
Definition grav3D.h:38

Grid3D
Class to create a the gravity object.
Definition grid3D.h:229

Grid3D::Grid3D
Grid3D(void)
Constructor for the grid.
Definition grid3D.cpp:43

Grid3D::Explosion_2D
void Explosion_2D()
Explosion test described in Liska, 2003.

Grid3D::Wind_Boundary
void Wind_Boundary()
Apply a constant wind to the -x boundary.
Definition boundary_conditions.cpp:525

Grid3D::Initialize_Models
void Initialize_Models(ParameterMap &pmap)
initialize the model collection
Definition grid3D.h:672

Grid3D::Noh_2D
void Noh_2D()
Noh test described in Liska, 2003.
Definition initial_conditions.cpp:1076

Grid3D::Write_Header
void Write_Header(AttrRecorderInterface &attr_recorder) const
Definition io.cpp:150

Grid3D::H
struct Header H
Header for the grid.
Definition grid3D.h:237

Grid3D::Blast_1D
void Blast_1D()
Initialize the grid with two interacting blast waves. See Stone 2008, Section 8.1.
Definition initial_conditions.cpp:635

Grid3D::Update_Hydro_Grid
Real Update_Hydro_Grid(std::function< void(Grid3D &)> &feedback_callback, std::function< void(Grid3D &)> &chemistry_callback)
Do all steps to update the hydro.
Definition grid3D.cpp:437

Grid3D::Gresho
void Gresho()
Initialize the grid with the 2D Gresho problem described in LW03.
Definition initial_conditions.cpp:932

Grid3D::Riemann
void Riemann(Parameters const &P)
Initialize the grid with a Riemann problem.
Definition initial_conditions.cpp:509

Grid3D::Noh_3D
void Noh_3D()
Noh test described in Stone, 2008.
Definition initial_conditions.cpp:1107

Grid3D::models
const ModelCollection & models() const
get the model collection
Definition grid3D.h:675

Grid3D::Initialize
void Initialize(struct Parameters *P)
Initialize the grid.
Definition grid3D.cpp:114

Grid3D::Spherical_Overpressure_3D
void Spherical_Overpressure_3D()
Initialize the grid with a 3D spherical overdensity and overpressue.
Definition initial_conditions.cpp:1203

Grid3D::Shu_Osher
void Shu_Osher()
Initialize the grid with the Shu-Osher shock tube problem. See Stone 2008, Section 8....
Definition initial_conditions.cpp:597

Grid3D::Custom_Boundary
void Custom_Boundary(char bcnd[MAXLEN])
Select appropriate custom boundary function.
Definition boundary_conditions.cpp:509

Grid3D::Rayleigh_Taylor
void Rayleigh_Taylor()
Initialize the grid with a 2D Rayleigh-Taylor instability.
Definition initial_conditions.cpp:883

Grid3D::Set_Domain_Properties
void Set_Domain_Properties(struct Parameters P)
Set local domain properties.
Definition initial_conditions.cpp:109

Grid3D::Read_Grid
void Read_Grid(struct Parameters P)
Read in grid data from 1-per-process output files.
Definition io.cpp:522

Grid3D::Set_Boundary_Conditions_Grid
void Set_Boundary_Conditions_Grid(Parameters P)
Set the boundary conditions for all components based on info in the parameters structure.
Definition boundary_conditions.cpp:19

Grid3D::field_info
FieldInfo field_info
Definition grid3D.h:240

Grid3D::Set_Boundaries
void Set_Boundaries(int dir, int flags[])
Apply boundary conditions to the grid.
Definition boundary_conditions.cpp:165

Grid3D::Check_Custom_Boundary
int Check_Custom_Boundary(int *flags, struct Parameters P)
Check for custom boundary conditions.
Definition boundary_conditions.cpp:133

Grid3D::KH_res_ind
void KH_res_ind()
Initialize the grid with a Kelvin-Helmholtz instability whose modes are resolution independent.
Definition initial_conditions.cpp:749

Grid3D::FreeMemory
void FreeMemory(void)
Free the memory for the density array.
Definition grid3D.cpp:624

Grid3D::Set_Boundary_Extents
void Set_Boundary_Extents(int dir, int *imin, int *imax)
Set the extents of the ghost region we are initializing.
Definition boundary_conditions.cpp:419

Grid3D::Square_Wave
void Square_Wave(Parameters const &P)
Square wave density perturbation with amplitude A*rho in pressure equilibrium.
Definition initial_conditions.cpp:437

Grid3D::Read_Grid_Cat
void Read_Grid_Cat(struct Parameters P)
Read in grid data from a single concatenated output file.
Definition io_parallel.cpp:133

Grid3D::Noh_Boundary
void Noh_Boundary()
Apply analytic boundary conditions to +x, +y (and +z) faces, as per the Noh problem in Liska,...
Definition boundary_conditions.cpp:544

Grid3D::Set_Boundary_Conditions
void Set_Boundary_Conditions(Parameters P)
Set the boundary conditions based on info in the parameters structure.
Definition boundary_conditions.cpp:54

Grid3D::Constant
void Constant(Parameters const &P)
Constant gas properties.
Definition initial_conditions.cpp:180

Grid3D::Disk_3D
void Disk_3D(Parameters P)
Initialize the grid with a 3D disk following a Miyamoto-Nagai profile.
Definition disk_ICs.cu:829

Grid3D::flag_init
int flag_init
Initialization flag.
Definition grid3D.h:233

Grid3D::Reset
void Reset(void)
Reset the Grid3D class.
Definition grid3D.cpp:613

Grid3D::Set_Initial_Conditions
void Set_Initial_Conditions(Parameters P, const ParameterMap &pmap)
Definition initial_conditions.cpp:27

Grid3D::Disk_2D
void Disk_2D()
Initialize the grid with a 2D disk following a Kuzmin profile.
Definition initial_conditions.cpp:1139

Grid3D::KH
void KH()
Initialize the grid with a Kelvin-Helmholtz instability with a discontinuous interface.
Definition initial_conditions.cpp:677

Grid3D::Sound_Wave
void Sound_Wave(Parameters const &P)
Sine wave perturbation.
Definition initial_conditions.cpp:244

Grid3D::Implosion_2D
void Implosion_2D()
Implosion test described in Liska, 2003.
Definition initial_conditions.cpp:1033

Grid3D::AllocateMemory
void AllocateMemory(void)
Allocate memory for the d, m, E arrays.
Definition grid3D.cpp:230

Grid3D::set_dt
void set_dt(Real dti)
Calculate the timestep.
Definition grid3D.cpp:311

Grid3D::Linear_Wave
void Linear_Wave(Parameters const &P)
Initialize the grid with a simple linear wave.
Definition initial_conditions.cpp:302

Grid3D::model_collection
ModelCollection model_collection
Definition grid3D.h:243

ModelCollection
Holds a collection of models.
Definition model_collection.h:215

ParameterMap
A class that provides map-like access to parameter files.
Definition ParameterMap.h:74

Time
Definition timing_functions.h:42

Grid3D::Conserved
Definition grid3D.h:278

Grid3D::Conserved::device
Real * device
Definition grid3D.h:365

Grid3D::Conserved::d_Grav_potential
Real * d_Grav_potential
Definition grid3D.h:372

Grid3D::Conserved::momentum_x
Real * momentum_x
Array containing the momentum in the x direction of each cell in the grid.
Definition grid3D.h:289

Grid3D::Conserved::momentum_y
Real * momentum_y
Array containing the momentum in the y direction of each cell in the grid.
Definition grid3D.h:294

Grid3D::Conserved::Energy
Real * Energy
Array containing the total Energy of each cell in the grid.
Definition grid3D.h:303

Grid3D::Conserved::host
Real * host
Definition grid3D.h:280

Grid3D::Conserved::momentum_z
Real * momentum_z
Array containing the momentum in the z direction of each cell in the grid.
Definition grid3D.h:299

Grid3D::Conserved::density
Real * density
Array containing the density of each cell in the grid.
Definition grid3D.h:284

Header
Definition grid3D.h:55

Header::min_dt_slow
Real min_dt_slow
Cells that introduce timesteps shorter than will be averaged with neighboring cells....
Definition grid3D.h:159

Header::nz
int nz
Total number of cells in the z-dimension.
Definition grid3D.h:74

Header::ny
int ny
Total number of cells in the y-dimension.
Definition grid3D.h:70

Header::nx_real
int nx_real
Number of real cells in the x-dimension.
Definition grid3D.h:78

Header::dy
Real dy
y-width of cells
Definition grid3D.h:142

Header::dx
Real dx
x-width of cells
Definition grid3D.h:138

Header::n_cells
int n_cells
Total number of cells in the grid (including ghost cells)
Definition grid3D.h:58

Header::n_fields
int n_fields
Number of fields (conserved variables, scalars, etc.)
Definition grid3D.h:171

Header::dz
Real dz
z-width of cells
Definition grid3D.h:146

Header::n_step
int n_step
Number of timesteps taken.
Definition grid3D.h:167

Header::t
Real t
Simulation time.
Definition grid3D.h:150

Header::t_wall
Real t_wall
Wall time.
Definition grid3D.h:163

Header::n_ghost
int n_ghost
Number of ghost cells on each side of the grid.
Definition grid3D.h:62

Header::nx
int nx
Total number of cells in the x-dimension.
Definition grid3D.h:66

Header::Output_Now
bool Output_Now
Flag set to true when data has to be written to file.
Definition grid3D.h:211

Header::custom_grav
int custom_grav
Flag to set specific static gravity field.
Definition grid3D.h:175

Header::dt
Real dt
Length of the current timestep.
Definition grid3D.h:154

Header::Output_Complete_Data
bool Output_Complete_Data
Flag set to true when all the data will be written to file (Restart File )
Definition grid3D.h:217

Parameters
Definition global.h:217