VlasovAdvectionCoeff.c

Go to the documentation of this file.

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

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

  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;

}