VlasovComputeCFL.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*, int, void*);

double VlasovComputeCFL( void    *s, 
                         void    *m, 
                         double  dt, 
                         double  t   
                       )
{
  HyPar   *solver = (HyPar*)  s;
  Vlasov  *params = (Vlasov*) solver->physics;

  int     ndims  = solver->ndims;
  int     ghosts = solver->ghosts;
  int     *dim   = solver->dim_local;
  double  *u     = solver->u;

  double max_cfl = 0;
  int done = 0; 
  int index[ndims];
  _ArraySetValue_(index,ndims,0);

  while (!done) {

    for (int dir=0; dir<ndims; dir++) {
      double dxinv;
      _GetCoordinate_(dir,index[dir],dim,ghosts,solver->dxinv,dxinv);
      double eig = VlasovAdvectionCoeff(index, dir, solver);
      double local_cfl = eig*dt*dxinv; 
      if (local_cfl > max_cfl) max_cfl = local_cfl;
    }

    _ArrayIncrementIndex_(ndims,dim,index,done);
  }

  return(max_cfl);
}