fluxes.c File Reference

Compute the flux for the relativistic hydro (RHD) equations. More...

#include "pluto.h"

Include dependency graph for fluxes.c:

Go to the source code of this file.

Functions
void	Flux (double **u, double **v, double *a2, double **fx, double *p, int beg, int end)

Detailed Description

Compute the flux for the relativistic hydro (RHD) equations.

Compute the flux of the conservative RHD equations in the direction given by g_dir. This function defines the component of the hyperbolic flux tensor of the standard RHD equations.
In what follows:

• VXn, MXn are the velocity, momentum components in the direction given by g_dir (normal, "n")
• VXt, MXt and VXb, MXb are the transverse components (tangent "t" and bi-tangent "b").

Author

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

Date

July 5, 2015

Definition in file fluxes.c.

Function Documentation

void Flux	(	double **	u,
		double **	v,
		double *	a2,
		double **	fx,
		double *	p,
		int	beg,
		int	end
	)

Parameters

[in]	u	1D array of conserved quantities
[in]	w	1D array of primitive quantities
[in]	a2	1D array of sound speeds
[out]	fx	1D array of fluxes (total pressure excluded)
[out]	p	1D array of pressure values
[in]	beg	initial index of computation
[in]	end	final index of computation

Returns

This function has no return value.

Definition at line 23 of file fluxes.c.

{
  int    nv, i;
  double vn;

  for (i = beg ; i <= end; i++) {
    vn = v[i][VXn];

    fx[i][RHO]  = u[i][RHO]*vn;
    EXPAND(fx[i][MX1] = u[i][MX1]*vn;  ,
           fx[i][MX2] = u[i][MX2]*vn;  ,
           fx[i][MX3] = u[i][MX3]*vn;)
    fx[i][ENG] = u[i][MXn];
    p[i] = v[i][PRS];
  }
}