flat_states.c File Reference

#include "pluto.h"

Include dependency graph for flat_states.c:

Go to the source code of this file.

Functions
void	States (const State_1D *state, int beg, int end, Grid *grid)

Function Documentation

void States	(	const State_1D *	state,
		int	beg,
		int	end,
		Grid *	grid
	)

Compute 1D left and right interface states using piecewise linear reconstruction and the characteristic decomposition of the quasi-linear form of the equations.

This is done by first extrapolating the cell center value to the interface using piecewise limited linear reconstruction on the characteristic variables.

Left and right states are then evolved for the half time step using characteristic tracing if necessary.

Definition at line 4 of file flat_states.c.

{
  int nv, i;

  #if TIME_STEPPING != EULER
   #error FLAT Reconstruction must be used with EULER integration only
  #endif
  
  for (i = beg; i <= end; i++) {
  for (nv = 0; nv < NVAR; nv++) {
    state->vm[i][nv] = state->vp[i][nv] = state->v[i][nv];
  }}
  PrimToCons(state->vm, state->um, beg, end);
  PrimToCons(state->vp, state->up, beg, end);
  
/*  -------------------------------------------
      Assign face-centered magnetic field
    -------------------------------------------  */

  #ifdef STAGGERED_MHD
   for (i = beg; i <= end-1; i++) {
     state->vR[i][BXn] = state->vL[i][BXn] = state->bn[i];
   }
  #endif

/* -------------------------------------------
    compute states in conservative variables
   ------------------------------------------- */

  PrimToCons (state->vp, state->up, beg, end);
  PrimToCons (state->vm, state->um, beg, end);
}

Here is the call graph for this function: