VlasovAdvectionCoeff.c File Reference

Returns the advection coefficient at a given grid index. More...

#include <float.h>
#include <math.h>
#include <basic.h>
#include <arrayfunctions.h>
#include <physicalmodels/vlasov.h>
#include <mpivars.h>
#include <hypar.h>

Go to the source code of this file.

Functions
double	VlasovAdvectionCoeff (int *idx, int dir, void *s)

Detailed Description

Returns the advection coefficient at a given grid index.

Author

John Loffeld

Definition in file VlasovAdvectionCoeff.c.

Function Documentation

VlasovAdvectionCoeff()

double VlasovAdvectionCoeff	(	int *	idx,
		int	dir,
		void *	s
	)

Returns the advection coefficient at a given grid index \(c\), where

\begin{equation} c = v_i, \end{equation}

if dir < Vlasov::ndims_x ( \(i = {\rm dir}\)), and

\begin{equation} c = E_i, \end{equation}

the electric field if dir >= Vlasov::ndims_x ( \(i = \) dir-Vlasov::ndims_x).

Note: this function assumes that the electric field has already been set.

Parameters

idx	grid index
dir	Spatial dimension
s	Solver object of type HyPar

Definition at line 28 of file VlasovAdvectionCoeff.c.

{

  HyPar  *solver = (HyPar*)  s;
  Vlasov *param  = (Vlasov*) solver->physics;

  int* dim    = solver->dim_local;
  int  ghosts = solver->ghosts;

  double retval = DBL_MAX;

  if (dir < param->ndims_x) {

    int veldim = dir + param->ndims_x;
    _GetCoordinate_(veldim,idx[veldim],dim,ghosts,solver->x,retval);

  }  else {

    int ndims_x = param->ndims_x;
    int dim_x[ndims_x]; _ArrayCopy1D_(dim, dim_x, ndims_x);

    int idx_x[ndims_x];
    _ArrayCopy1D_(idx, idx_x, ndims_x);
    int p; _ArrayIndex1D_(ndims_x, dim_x, idx_x, ghosts, p);
    retval = param->e_field[ndims_x*p+(dir-ndims_x)];
    
  }

  return retval;

}