plm_coeffs.c File Reference

Compute linear reconstruction coefficients. More...

#include "pluto.h"

Include dependency graph for plm_coeffs.c:

Go to the source code of this file.

Functions
void	PLM_CoefficientsSet (Grid *grid)
void	PLM_CoefficientsGet (PLM_Coeffs *plm_coeffs, int dir)

Variables
static double **	cp3D
static double **	cm3D
static double **	wp3D
static double **	wm3D
static double **	dp3D
static double **	dm3D

Detailed Description

Compute linear reconstruction coefficients.

Reconstruction coefficients needed for linear interpolation are computed in the PLM_CoefficientsSet() each time a new grid is created. There's a different set of 1D coefficients in each direction.

The function PLM_CoefficientsGet() can be used to obtain a set of coefficients along a desired direction.

References

• "High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates" A. Mignone, JCP (2014), 270, 784.

Authors

A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)

Date

May 19, 2014

Definition in file plm_coeffs.c.

Function Documentation

void PLM_CoefficientsGet	(	PLM_Coeffs *	plm_coeffs,
		int	dir
	)

Retrieve reconstruction coefficients in the PLM_Coeffs structure. This function can be called only if the previous one has been completed already.

Parameters

[out]	plm_coeffs	a pointer to a PLM_Coeffs structure containing the 1D coefficients needed for reconstruction
[in]	dir	the desired direction

Definition at line 79 of file plm_coeffs.c.

{
  if (cp3D == NULL) {
    print1 ("! PLM_Coefficients: coefficients not set.\n");
    QUIT_PLUTO(1);
  }

  plm_coeffs->cp = cp3D[dir];
  plm_coeffs->cm = cm3D[dir];

  plm_coeffs->wp = wp3D[dir];
  plm_coeffs->wm = wm3D[dir];

  plm_coeffs->dp = dp3D[dir];
  plm_coeffs->dm = dm3D[dir];

}

Here is the call graph for this function:

Here is the caller graph for this function:

void PLM_CoefficientsSet

(

Grid *

grid

)

Compute interpolation coefficients for linear reconstruction.

Parameters

[in]

grid

pointer to array of grid structure

Returns

This function has no return value

Definition at line 30 of file plm_coeffs.c.

{
  int    i, d, beg, end;
  double *dx, *xr, *xgc;

  if (cp3D == NULL) {
    cp3D = ARRAY_2D(DIMENSIONS, NMAX_POINT, double);
    cm3D = ARRAY_2D(DIMENSIONS, NMAX_POINT, double);
    dp3D = ARRAY_2D(DIMENSIONS, NMAX_POINT, double);
    dm3D = ARRAY_2D(DIMENSIONS, NMAX_POINT, double);
    wp3D = ARRAY_2D(DIMENSIONS, NMAX_POINT, double);
    wm3D = ARRAY_2D(DIMENSIONS, NMAX_POINT, double);
  }

/* -----------------------------------------------------
    Compute interpolation coefficients.
    This must be the first call to this function in
    order to store coefficients into memory.
   ----------------------------------------------------- */

  for (d = 0; d < DIMENSIONS; d++){
    dx  = grid[d].dx;
    xgc = grid[d].xgc;
    xr  = grid[d].xr;

  /* -- first and last zone are excluded -- */
    
    beg = 1;
    end = grid[d].np_tot - 2;
    for (i = beg; i <= end; i++){
      wp3D[d][i] = dx[i]/(xgc[i+1] - xgc[i]); /* coeff. for dQF in Eq. [29] */
      wm3D[d][i] = dx[i]/(xgc[i] - xgc[i-1]); /* coeff. for dQB in Eq. [29] */

      cp3D[d][i] = (xgc[i+1] - xgc[i])/(xr[i] - xgc[i]);   /* Eq. [33], cF */
      cm3D[d][i] = (xgc[i] - xgc[i-1])/(xgc[i] - xr[i-1]); /* Eq. [33], cB */

      dp3D[d][i] = (xr[i] - xgc[i])/dx[i];     /* Eq. [30], plus sign */
      dm3D[d][i] = (xgc[i] - xr[i-1])/dx[i];   /* Eq. [30], minus sign */
    }
  }
}

Here is the caller graph for this function:

Variable Documentation

double ** cm3D

static

Definition at line 25 of file plm_coeffs.c.

double** cp3D

static

Definition at line 25 of file plm_coeffs.c.

double ** dm3D

static

Definition at line 27 of file plm_coeffs.c.

double** dp3D

static

Definition at line 27 of file plm_coeffs.c.

double ** wm3D

static

Definition at line 26 of file plm_coeffs.c.

double** wp3D

static

Definition at line 26 of file plm_coeffs.c.