interppoint.para.c File Reference

Parabolic/PPM Spatial Interpolation for fluxes based upon providing each point / parabolic interpolation subroutin / ref. Colella && Woodward's paper / Colella, P., & Woodward, P. R. 1984, J. Comput. Phys., 54, 174-201 / / using zone-centered value of 5 continuous zones / to get left and right value of the middle zone. / / / Latest JCM version is para4()/parapl(): 02/25/08. More...

#include "decs.h"
#include "para_and_paraenohybrid.h"

Include dependency graph for interppoint.para.c:

Go to the source code of this file.

Macros
#define	DO4MONO   0
#define	OVERRIDEWITHMINM   1
#define	PARAGENMINMOD(a, b)   MINMOD(a,b)

Functions
void	para (FTYPE *y, FTYPE *lout, FTYPE *rout)
	lout/rout is left and right sides of cell note how used in step_ch.c to get appropriate interface value given by Xiaoyue Guan to Scott Noble on Nov 9, 2004, given to me Jan 7, 2005 More...
void	para4 (int realisinterp, int pl, FTYPE *y, FTYPE *lout, FTYPE *rout)
	doesn't use dq More...
void	paracontsmooth (int pl, FTYPE *y, FTYPE *facecont, int *smooth)
	PARA starts with continuous 3rd order polynomial (4th order through averages if used) using monotonized slopes for second derivative term facecont[i] = 0.5*(y[i]+y[i-1]) - ddq[i] where ddq[0] = dq[0]-dq[-1]. More...
void	parapl (int i, int j, int k, int loc, int realisinterp, int dir, FTYPE **yrealpl, FTYPE **ypl, FTYPE *loutpl, FTYPE *routpl)
	used when lim=PARAFLAT More...
FTYPE	divftilde (int dir, int shift, FTYPE *Vabs, FTYPE *Pabs)
void	getPressure (int whicheom, int i, int j, int k, int loc, FTYPE **yrealpl, FTYPE *P)
	Get pressure Note this is quite inefficient since operating per-point get same pressure for entire line multiple times (SUPERGODMARK: Also no accounting of magnetic field) More...

Detailed Description

Parabolic/PPM Spatial Interpolation for fluxes based upon providing each point / parabolic interpolation subroutin / ref. Colella && Woodward's paper / Colella, P., & Woodward, P. R. 1984, J. Comput. Phys., 54, 174-201 / / using zone-centered value of 5 continuous zones / to get left and right value of the middle zone. / / / Latest JCM version is para4()/parapl(): 02/25/08.

Definition in file interppoint.para.c.

Macro Definition Documentation

#define DO4MONO   0

Definition at line 405 of file interppoint.para.c.

#define OVERRIDEWITHMINM   1

Definition at line 407 of file interppoint.para.c.

#define PARAGENMINMOD	(		a,
			b
	)		MINMOD(a,b)

Definition at line 604 of file interppoint.para.c.

Function Documentation

FTYPE divftilde	(	int	dir,
		int	shift,
		FTYPE *	Vabs,
		FTYPE *	Pabs
	)

Definition at line 1371 of file interppoint.para.c.

void getPressure	(	int	whicheom,
		int	i,
		int	j,
		int	k,
		int	loc,
		FTYPE **	yrealpl,
		FTYPE *	P
	)

Get pressure Note this is quite inefficient since operating per-point get same pressure for entire line multiple times (SUPERGODMARK: Also no accounting of magnetic field)

Definition at line 1433 of file interppoint.para.c.

Here is the call graph for this function:

Here is the caller graph for this function:

void para	(	FTYPE *	y,
		FTYPE *	lout,
		FTYPE *	rout
	)

lout/rout is left and right sides of cell note how used in step_ch.c to get appropriate interface value given by Xiaoyue Guan to Scott Noble on Nov 9, 2004, given to me Jan 7, 2005

Definition at line 56 of file interppoint.para.c.

Here is the caller graph for this function:

void para4	(	int	realisinterp,
		int	pl,
		FTYPE *	y,
		FTYPE *	lout,
		FTYPE *	rout
	)

doesn't use dq

Definition at line 571 of file interppoint.para.c.

Here is the caller graph for this function:

void paracontsmooth	(	int	pl,
		FTYPE *	y,
		FTYPE *	facecont,
		int *	smooth
	)

PARA starts with continuous 3rd order polynomial (4th order through averages if used) using monotonized slopes for second derivative term facecont[i] = 0.5*(y[i]+y[i-1]) - ddq[i] where ddq[0] = dq[0]-dq[-1].

Definition at line 985 of file interppoint.para.c.

void parapl	(	int	i,
		int	j,
		int	k,
		int	loc,
		int	realisinterp,
		int	dir,
		FTYPE **	yrealpl,
		FTYPE **	ypl,
		FTYPE *	loutpl,
		FTYPE *	routpl
	)

used when lim=PARAFLAT

Definition at line 1205 of file interppoint.para.c.

Here is the call graph for this function:

Here is the caller graph for this function: