#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "variables_global.h"
#include "prototypes.h"

Functions
double	max_states (double a, double b)
int	states (int nmin, int nmax, int flag, double pre_1D, double rho_1D, double dx, double Cdtdx, double fCdtdx, double plft, double pl, double dp, double p6, double ulft, double ul, double du, double u6, double rlft, double rl, double dr, double r6, double prgh, double urgh, double rrgh)

Function Documentation

double max_states	(	double	a,
		double	b
	)

inline

Inline function to calculate the maximum of two doubles

Definition at line 17 of file states.c.

Referenced by states().

                                             {
    double tmp;
    if (a < b) tmp = b;
    if (a == b)tmp = b;
    if (a > b) tmp = a;
    return tmp;
}

int states	(	int	nmin,
		int	nmax,
		int	flag,
		double *	pre_1D,
		double *	rho_1D,
		double *	dx,
		double *	Cdtdx,
		double *	fCdtdx,
		double *	plft,
		double *	pl,
		double *	dp,
		double *	p6,
		double *	ulft,
		double *	ul,
		double *	du,
		double *	u6,
		double *	rlft,
		double *	rl,
		double *	dr,
		double *	r6,
		double *	prgh,
		double *	urgh,
		double *	rrgh
	)

Get hydroquantities at cell centres and faces.

Definition at line 30 of file states.c.

References dt, Gamma1, grav_acc, gravity_on_off, max_states(), smallp, and smallr.

Referenced by ppm_step().

                                                              {
    int i, ii;
    double hdt, svel;
    double fourthd;
    fourthd = 4.0 / 3.0;
    svel = 0.0;
    hdt = 0.5 * dt;
    for (i = nmin - 4; i <= nmax + 4; i++) {
        Cdtdx[i] = sqrt(Gamma1 * pre_1D[i] / rho_1D[i]) / (dx[i]);
        svel = max_states(svel, Cdtdx[i]);
        Cdtdx[i] = Cdtdx[i] * hdt;
        fCdtdx[i] = 1.0 - fourthd * Cdtdx[i];
    }
    for (i = nmin - 4; i <= nmax + 4; i++) {
        ii = i + 1;
        plft[ii] = pl[i] + dp[i] - Cdtdx[i]*(dp[i] - fCdtdx[i] * p6[i]);
        ulft[ii] = ul[i] + du[i] - Cdtdx[i]*(du[i] - fCdtdx[i] * u6[i]);
        rlft[ii] = rl[i] + dr[i] - Cdtdx[i]*(dr[i] - fCdtdx[i] * r6[i]);
        plft[ii] = max_states(smallp, plft[ii]);
        rlft[ii] = max_states(smallr, rlft[ii]);
        if (flag == 0) {
            ulft[ii] = ulft[ii] + hdt * grav_acc;
        } else {
            ulft[ii] = ulft[ii];
        }
        prgh[i] = pl[i] + Cdtdx[i]*(dp[i] + fCdtdx[i] * p6[i]);
        urgh[i] = ul[i] + Cdtdx[i]*(du[i] + fCdtdx[i] * u6[i]);
        rrgh[i] = rl[i] + Cdtdx[i]*(dr[i] + fCdtdx[i] * r6[i]);
        prgh[i] = max_states(smallp, prgh[i]);
        rrgh[i] = max_states(smallr, rrgh[i]);
        if (flag == 1) {
            if (gravity_on_off == 1) urgh[i] = urgh[i] + hdt * grav_acc;
            if (gravity_on_off == 0) urgh[i] = urgh[i];
        } else {
            urgh[i] = urgh[i];
        }
    }
    
  
    return 0;
}

Functions

Function Documentation