eos.c File Reference

Implementation of the isothermal EOS. More...

#include "pluto.h"

Include dependency graph for eos.c:

Go to the source code of this file.

Functions
void	SoundSpeed2 (double **v, double *cs2, double *h, int beg, int end, int pos, Grid *grid)
void	Enthalpy (double **v, real *h, int beg, int end)
void	Entropy (double **v, double *s, int beg, int end)

Detailed Description

Implementation of the isothermal EOS.

Author

A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)

Date

April 14, 2014

Definition in file eos.c.

Function Documentation

void Enthalpy	(	double **	v,
		real *	h,
		int	beg,
		int	end
	)

Compute the enthalpy.

Parameters

[in]	v	1D array of primitive quantities
[in]	h	1D array of enthalpy values
[in]	beg	initial index of computation
[in]	end	final index of computation

Returns

This function has no return value.

Definition at line 63 of file eos.c.

{
  print ("! Enthalpy: enthalpy not defined in isothermal EOS\n");
  QUIT_PLUTO(1);
}

Here is the call graph for this function:

void Entropy	(	double **	v,
		double *	s,
		int	beg,
		int	end
	)

Compute the entropy.

Parameters

[in]	v	1D array of primitive quantities
[in]	s	1D array of entropy values
[in]	is	initial index of computation
[in]	ie	final index of computation

Returns

This function has no return value.

Definition at line 80 of file eos.c.

{
  print ("Entropy: entropy not defined for isothermal EOS\n");
  QUIT_PLUTO(1);
}

Here is the call graph for this function:

void SoundSpeed2	(	double **	v,
		double *	cs2,
		double *	h,
		int	beg,
		int	end,
		int	pos,
		Grid *	grid
	)

Define the square of the sound speed.

Parameters

[in]	v	1D array of primitive quantities
[out]	cs2	1D array containing the square of the sound speed
[in]	h	1D array of enthalpy values
[in]	beg	initial index of computation
[in]	end	final index of computation
[in]	pos	an integer specifying the spatial position inside the cell (only for spatially-dependent EOS)
[in]	grid	pointer to an array of Grid structures

Returns

This function has no return value.

Definition at line 13 of file eos.c.

{
  int  i, j,k;  /* -- used as multidimensional indices -- */
  double *x1, *x2, *x3;

  #if PHYSICS == HD || PHYSICS == MHD
   x1 = grid[IDIR].x;
   x2 = grid[JDIR].x;
   x3 = grid[KDIR].x;

   i = g_i; j = g_j; k = g_k;
    
   if (g_dir == IDIR) {

     x1 = (pos == FACE_CENTER ? grid[IDIR].xr : grid[IDIR].x);
     for (i = beg; i <= end; i++) cs2[i] = g_isoSoundSpeed*g_isoSoundSpeed;

   }else if (g_dir == JDIR){

     x2 = (pos == FACE_CENTER ? grid[JDIR].xr : grid[JDIR].x);
     for (j = beg; j <= end; j++) cs2[j] = g_isoSoundSpeed*g_isoSoundSpeed;

   }else if (g_dir == KDIR){

     x3 = (pos == FACE_CENTER ? grid[KDIR].xr : grid[KDIR].x);
     for (k = beg; k <= end; k++) cs2[k] = g_isoSoundSpeed*g_isoSoundSpeed;

  }
  #else
   print ("! SoundSpeed2: not defined for this EoS\n");
   QUIT_PLUTO(1);
  #endif
}

Here is the call graph for this function: