VlasovFunctions.c File Reference

Misc. functions for the Vlasov equations. More...

#include <physicalmodels/vlasov.h>
#include <hypar.h>

Go to the source code of this file.

Functions
int	VlasovEField (double *, void *, double)
int	VlasovPreStep (double *u, void *s, void *m, double waqt)
int	VlasovPostStage (double *u, void *s, void *m, double waqt)

Detailed Description

Misc. functions for the Vlasov equations.

Author

Debojyoti Ghosh

Definition in file VlasovFunctions.c.

Function Documentation

int VlasovEField	(	double *	u,
		void *	s,
		double	t
	)

Compute electric field.

Parameters

u	Conserved solution
s	Solver object of type HyPar
t	Current time

Definition at line 258 of file VlasovEField.c.

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

  int ierr;

  if (param->self_consistent_electric_field) {
    ierr = SetEFieldSelfConsistent(u, s, t);
  } else {
    ierr = SetEFieldPrescribed(u, s, t);
  }

  return ierr;
}

int VlasovPreStep	(	double *	u,
		void *	s,
		void *	m,
		double	waqt
	)

Vlasov-specific function called at the beginning of each time-step: Calls the function to set the electric field

Parameters

u	Solution (conserved variables)
s	Solver object of type HyPar
m	MPI object of type MPIVariables
waqt	Current solution time

Definition at line 14 of file VlasovFunctions.c.

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

  int ierr = VlasovEField(u, solver, waqt);
  if (ierr) return ierr;

  return 0;
}

int VlasovPostStage	(	double *	u,
		void *	s,
		void *	m,
		double	waqt
	)

Vlasov-specific function called at the end of each stage in a multistage time integrator: Calls the function to set the electric field

Parameters

u	Solution (conserved variables)
s	Solver object of type HyPar
m	MPI object of type MPIVariables
waqt	Current solution time

Definition at line 33 of file VlasovFunctions.c.

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

  int ierr = VlasovEField(u, solver, waqt);
  if (ierr) return ierr;

  return 0;
}