BurgersAdvection.c File Reference

Contains the function to compute the hyperbolic flux for the Burgers equations over the domain. More...

#include <stdlib.h>
#include <basic.h>
#include <arrayfunctions.h>
#include <physicalmodels/burgers.h>
#include <hypar.h>

Go to the source code of this file.

Functions
int	BurgersAdvection (double *f, double *u, int dir, void *s, double t)

Detailed Description

Contains the function to compute the hyperbolic flux for the Burgers equations over the domain.

Author

John Loffeld

Definition in file BurgersAdvection.c.

Function Documentation

BurgersAdvection()

int BurgersAdvection	(	double *	f,
		double *	u,
		int	dir,
		void *	s,
		double	t
	)

Compute the hyperbolic flux over the local domain.

\begin{equation} {\bf F}\left({\bf u}\right) = 0.5 {\bf u}^2 \end{equation}

Parameters

f	Array to hold the computed flux (same size and layout as u)
u	Array containing the conserved solution
dir	Spatial dimension
s	Solver object of type HyPar
t	Current time

Definition at line 17 of file BurgersAdvection.c.

{
  HyPar     *solver = (HyPar*)   s;
  Burgers   *param  = (Burgers*) solver->physics;
  int        i, v;

  int *dim    = solver->dim_local;
  int ghosts  = solver->ghosts;
  int ndims   = solver->ndims;
  int nvars   = solver->nvars;
  int index[ndims], bounds[ndims], offset[ndims];

  /* set bounds for array index to include ghost points */
  _ArrayCopy1D_(dim,bounds,ndims);
  for (i=0; i<ndims; i++) bounds[i] += 2*ghosts;

  /* set offset such that index is compatible with ghost point arrangement */
  _ArraySetValue_(offset,ndims,-ghosts);

  int done = 0; _ArraySetValue_(index,ndims,0);
  while (!done) {
    int p; _ArrayIndex1DWO_(ndims,dim,index,offset,ghosts,p);
    for (v = 0; v < nvars; v++) f[nvars*p+v] = 0.5 * u[nvars*p+v] * u[nvars*p+v];
    _ArrayIncrementIndex_(ndims,bounds,index,done);
  }

  return(0);
}