#include "numericalFlux.h"

Include dependency graph for numericalFlux.c:

Macros
#define	TR_ALPHA 0.5

#define	TR_ALPHA_EXT 1.0

Functions
void	calculateInteriorNumericalAdvectiveFluxConvexOneSonicPoint (double sonicPoint, double sonicFlux, int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double f, double df, double flux, double dflux_left, double dflux_right)
	Calculate the advective flux at at interior element boundaries for simple scalar nonlinear hyperbolic pdes. More...

void	calculateInteriorNumericalAdvectiveFluxRusanov (double safetyFactor, int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_element, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double f, double df, double df_element, double flux, double dflux_left, double *dflux_right)

void	calculateExteriorNumericalAdvectiveFluxRusanov (double safetyFactor, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double df_element, double flux, double *dflux)

void	calculateGlobalExteriorNumericalAdvectiveFluxRusanov (double safetyFactor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double df_element, double flux, double *dflux)

void	calculateInteriorNumericalAdvectiveFluxRusanovWithEigenvalueBound (double safetyFactor, int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_element, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double f, double lambda_bar_element, double flux)

void	calculateGlobalExteriorNumericalAdvectiveFluxRusanovWithEigenvalueBound (double safetyFactor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double u, double f, double lambda_bar, double flux)

void	calculateExteriorNumericalAdvectiveFluxRusanovWithEigenvalueBound (double safetyFactor, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double u, double f, double lambda_bar, double flux)

void	calculateInteriorNumericalDiffusiveFlux (int scale_penalty, double penalty_floor, int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double grad_phi, double u, double penalty, double *flux)
	Calculate the diffusive flux at interior element boundary quadrature points. More...

void	calculateInteriorNumericalDiffusiveFlux_sd (int scale_penalty, double penalty_floor, int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double grad_phi, double u, double penalty, double *flux)

void	updateInteriorNumericalDiffusiveFluxJacobian (int scale_penalty, double penalty_floor, int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int l2g, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double *fluxJacobian)
	Calculate the diffusive flux Jacobian at interior element boundary quadrature points. More...

void	updateInteriorNumericalDiffusiveFluxJacobian_sd (int scale_penalty, double penalty_floor, int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int l2g, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double *fluxJacobian)

void	calculateExteriorNumericalDiffusiveFlux (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)
	Calculate the diffusive flux at exterior element boundary quadrature points. More...

void	calculateExteriorNumericalDiffusiveFlux_sd (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)

void	calculateGlobalExteriorNumericalDiffusiveFlux (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)
	Calculate the diffusive flux at exterior element boundary quadrature points. More...

void	calculateGlobalExteriorNumericalDiffusiveFlux_sd (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)

void	calculateExteriorNumericalDiffusiveFlux_free (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)

void	calculateExteriorNumericalDiffusiveFlux_free_sd (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)

void	calculateGlobalExteriorNumericalDiffusiveFlux_free (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)

void	calculateGlobalExteriorNumericalDiffusiveFlux_free_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double *flux)

void	updateExteriorNumericalDiffusiveFluxJacobian (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)
	Update the diffusive flux Jacobian at exterior element boundary quadrature points. More...

void	updateExteriorNumericalDiffusiveFluxJacobian_sd (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)

void	updateGlobalExteriorNumericalDiffusiveFluxJacobian (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)
	Update the diffusive flux Jacobian at exterior element boundary quadrature points. More...

void	updateGlobalExteriorNumericalDiffusiveFluxJacobian_sd (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)

void	updateExteriorNumericalDiffusiveFluxJacobian_free (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)

void	updateExteriorNumericalDiffusiveFluxJacobian_free_sd (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)

void	updateGlobalExteriorNumericalDiffusiveFluxJacobian_free (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)

void	updateGlobalExteriorNumericalDiffusiveFluxJacobian_free_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double fluxJacobian)

void	calculateInteriorNumericalAdvectiveFlux (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double f, double df, double flux, double dflux_left, double dflux_right)
	Calculate the advective flux at at interior element boundaries. More...

void	updateInteriorNumericalAdvectiveFluxJacobian (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double dflux_left, double dflux_right, double v, double *fluxJacobian)
	Calculate the advective flux at at interior element boundaries. More...

void	updateInteriorTwoSidedNumericalFluxJacobian (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double dflux_left, double dflux_right, double v, double *fluxJacobian_2sided)
	Calculate the two-sided flux jacobian at at interior element boundaries. More...

void	calculateInteriorNumericalAdvectiveFlux_average (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double f, double df, double flux, double dflux_left, double dflux_right)
	Calculate the advective flux at at interior element boundaries. More...

void	calculateExteriorNumericalAdvectiveFlux_NoBC (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int inflowFlag, double n, double f, double df, double flux, double *dflux_left)
	Update the advective flux at exterior element boundaries. More...

void	calculateGlobalExteriorNumericalAdvectiveFlux_NoBC (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int inflowFlag, double n, double f, double df, double flux, double *dflux_left)
	Update the advective flux at exterior element boundaries. More...

void	calculateExteriorNumericalAdvectiveFlux (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double flux, double dflux)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateGlobalExteriorNumericalAdvectiveFlux (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double flux, double dflux)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateExteriorNumericalAdvectiveFlux_free (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double flux, double dflux)

void	calculateGlobalExteriorNumericalAdvectiveFlux_free (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double flux, double dflux)

void	calculateExteriorNumericalAdvectiveFluxStokesP2D (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, double n, double bc_f, double bc_fpu, double bc_fpv, double f, double fpu, double fpv, double df_du, double df_dv, double dfpu_dp, double dfpv_dp, double flux, double fluxpu, double fluxpv, double dflux_du, double dflux_dv, double dfluxpu_dp, double dfluxpv_dp)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateExteriorNumericalAdvectiveFluxNavierStokes2D (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, double n, double bc_p, double bc_f_mass, double bc_f_umom, double bc_f_vmom, double p, double f_mass, double f_umom, double f_vmom, double df_mass_du, double df_mass_dv, double df_umom_du, double df_umom_dv, double df_vmom_du, double df_vmom_dv, double flux_mass, double flux_umom, double flux_vmom, double dflux_mass_du, double dflux_mass_dv, double dflux_umom_dp, double dflux_umom_du, double dflux_umom_dv, double dflux_vmom_dp, double dflux_vmom_du, double dflux_vmom_dv)

void	calculateGlobalExteriorNumericalAdvectiveFluxNavierStokes2D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, double n, double bc_p, double bc_f_mass, double bc_f_umom, double bc_f_vmom, double p, double oneByRho, double f_mass, double f_umom, double f_vmom, double df_mass_du, double df_mass_dv, double df_umom_dp, double df_umom_du, double df_umom_dv, double df_vmom_dp, double df_vmom_du, double df_vmom_dv, double flux_mass, double flux_umom, double flux_vmom, double dflux_mass_dp, double dflux_mass_du, double dflux_mass_dv, double dflux_umom_dp, double dflux_umom_du, double dflux_umom_dv, double dflux_vmom_dp, double dflux_vmom_du, double dflux_vmom_dv, double *velocity)

void	calculateGlobalExteriorNumericalAdvectiveFluxNavierStokes3D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, int isDOFBoundary_w, double n, double bc_p, double bc_f_mass, double bc_f_umom, double bc_f_vmom, double bc_f_wmom, double p, double f_mass, double f_umom, double f_vmom, double f_wmom, double df_mass_du, double df_mass_dv, double df_mass_dw, double df_umom_dp, double df_umom_du, double df_umom_dv, double df_umom_dw, double df_vmom_dp, double df_vmom_du, double df_vmom_dv, double df_vmom_dw, double df_wmom_dp, double df_wmom_du, double df_wmom_dv, double df_wmom_dw, double flux_mass, double flux_umom, double flux_vmom, double flux_wmom, double dflux_mass_du, double dflux_mass_dv, double dflux_mass_dw, double dflux_umom_dp, double dflux_umom_du, double dflux_umom_dv, double dflux_umom_dw, double dflux_vmom_dp, double dflux_vmom_du, double dflux_vmom_dv, double dflux_vmom_dw, double dflux_wmom_dp, double dflux_wmom_du, double dflux_wmom_dv, double dflux_wmom_dw, double *velocity)

void	calculateGlobalExteriorNumericalAdvectiveFluxStokesP2D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, double n, double bc_f, double bc_fpu, double bc_fpv, double f, double fpu, double fpv, double df_du, double df_dv, double dfpu_dp, double dfpv_dp, double flux, double fluxpu, double fluxpv, double dflux_du, double dflux_dv, double dfluxpu_dp, double dfluxpv_dp)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateExteriorNumericalAdvectiveFluxStokesP3D (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, int isDOFBoundary_w, double n, double bc_f, double bc_fpu, double bc_fpv, double bc_fpw, double f, double fpu, double fpv, double fpw, double df_du, double df_dv, double df_dw, double dfpu_dp, double dfpv_dp, double dfpw_dp, double flux, double fluxpu, double fluxpv, double fluxpw, double dflux_du, double dflux_dv, double dflux_dw, double dfluxpu_dp, double dfluxpv_dp, double dfluxpw_dp)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateGlobalExteriorNumericalAdvectiveFluxStokesP3D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, int isDOFBoundary_w, double n, double bc_f, double bc_fpu, double bc_fpv, double bc_fpw, double f, double fpu, double fpv, double fpw, double df_du, double df_dv, double df_dw, double dfpu_dp, double dfpv_dp, double dfpw_dp, double flux, double fluxpu, double fluxpv, double fluxpw, double dflux_du, double dflux_dv, double dflux_dw, double dfluxpu_dp, double dfluxpv_dp, double dfluxpw_dp)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateGlobalExteriorNumericalAdvectiveFluxStokes2D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, double n, double bc_p, double bc_f_mass, double p, double f_mass, double df_mass_du, double df_mass_dv, double flux_mass, double flux_umom, double flux_vmom, double dflux_mass_du, double dflux_mass_dv, double dflux_umom_dp, double dflux_vmom_dp, double *velocity)
	Apply basic pressure boundary penalty term for Stokes. More...

void	calculateGlobalExteriorNumericalAdvectiveFluxStokes3D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_p, int isDOFBoundary_u, int isDOFBoundary_v, int isDOFBoundary_w, double n, double bc_p, double bc_f_mass, double p, double f_mass, double df_mass_du, double df_mass_dv, double df_mass_dw, double flux_mass, double flux_umom, double flux_vmom, double flux_wmom, double dflux_mass_du, double dflux_mass_dv, double dflux_mass_dw, double dflux_umom_dp, double dflux_vmom_dp, double dflux_wmom_dp, double velocity)

void	calculateExteriorNumericalAdvectiveFlux_average (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double flux, double dflux)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateGlobalExteriorNumericalAdvectiveFlux_average (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_f, double bc_df, double u, double f, double df, double flux, double dflux)
	Calculate the advective flux at at exterior element boundaries. More...

void	calculateGlobalExteriorInflowNumericalAdvectiveFlux (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int inflowFlag, double inflowFlux, double n, double f, double df, double flux, double dflux_left)
	Update the advective flux at exterior inflow element boundaries. More...

void	updateExteriorNumericalAdvectiveFluxJacobian (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int inflowFlag, double dflux_left, double v, double *fluxJacobian)
	Calculate the advective flux at exterior element boundaries. More...

void	updateGlobalExteriorNumericalAdvectiveFluxJacobian (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int inflowFlag, double dflux_left, double v, double *fluxJacobian)
	Calculate the advective flux at global exterior element boundaries. More...

void	updateExteriorNumericalAdvectiveFluxJacobian_free (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int inflowFlag, double dflux_left, double v, double *fluxJacobian)

void	updateGlobalExteriorNumericalAdvectiveFluxJacobian_free (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int inflowFlag, double dflux_left, double v, double *fluxJacobian)

void	setInflowFlux (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int exteriorElementBoundaries, double inflowFlux, double *flux)
	Set the advective flux boundary condition at exterior element boundaries from the current exterior flux. More...

void	calculateInteriorNumericalDiffusiveFlux_LDG_upwind (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double a, double phi, double V, double penalty, double flux)
	Calculate the advective flux at at interior element boundaries. More...

void	calculateInteriorNumericalDiffusiveFlux_LDG_upwind_sd (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double a, double phi, double V, double penalty, double flux)

void	updateInteriorNumericalDiffusiveFluxJacobian_LDG_upwind (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double dphi, double V, double DV, double DV_eb, double v, double penalty, double fluxJacobian, double fluxJacobian_eb)
	Update the advective flux at at interior element boundaries. More...

void	updateInteriorNumericalDiffusiveFluxJacobian_LDG_upwind_sd (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double dphi, double V, double DV, double DV_eb, double v, double penalty, double fluxJacobian, double fluxJacobian_eb)

void	calculateExteriorNumericalDiffusiveFlux_LDG_upwind (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double a, double phi_bc, double phi, double V, double penalty, double *flux)
	Calculate the advective flux at at interior element boundaries. More...

void	calculateExteriorNumericalDiffusiveFlux_LDG_upwind_sd (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double a, double phi_bc, double phi, double V, double penalty, double *flux)

void	calculateGlobalExteriorNumericalDiffusiveFlux_LDG_upwind (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double a, double phi_bc, double phi, double V, double penalty, double *flux)
	Calculate the advective flux at at global exterior element boundaries. More...

void	calculateGlobalExteriorNumericalDiffusiveFlux_LDG_upwind_sd (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double a, double phi_bc, double phi, double V, double penalty, double *flux)

void	updateExteriorNumericalDiffusiveFluxJacobian_LDG_upwind (int isDiffusiveFluxBoundary, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double dphi, double V, double DV, double DV_eb, double v, double penalty, double fluxJacobian, double *fluxJacobian_eb)
	update the flux Jacobian with the advective flux contribution at at exterior element boundaries More...

void	updateExteriorNumericalDiffusiveFluxJacobian_LDG_upwind_sd (int isDiffusiveFluxBoundary, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double dphi, double V, double DV, double DV_eb, double v, double penalty, double fluxJacobian, double *fluxJacobian_eb)

void	updateGlobalExteriorNumericalDiffusiveFluxJacobian_LDG_upwind (int isDiffusiveFluxBoundary, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double dphi, double V, double DV, double DV_eb, double v, double penalty, double fluxJacobian_exterior, double *fluxJacobian_eb)
	update the flux Jacobian with the advective flux contribution at at exterior element boundaries More...

void	updateGlobalExteriorNumericalDiffusiveFluxJacobian_LDG_upwind_sd (int isDiffusiveFluxBoundary, int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double a, double da, double dphi, double V, double DV, double DV_eb, double v, double penalty, double fluxJacobian_exterior, double *fluxJacobian_eb)

void	calculateDiffusionMatrixSplittings_LDG_sd (int aSplit, int nElements_global, int nElementBoundaries_element, int nQuadraturePoints_element, int nQuadraturePoints_elementBoundary, int nSpace, const int rowptr, const int colind, const double ebq_a, const double q_a, double eb_aHat, double eb_aTilde, double aHat, double aTilde)

void	calculateInteriorLesaintRaviartNumericalFlux (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int speedEvalFlag, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double H, double dH, double flux, double dflux_left, double dflux_right)

void	calculateExteriorLesaintRaviartNumericalFlux (int nExteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nSpace, int speedEvalFlag, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_H, double bc_dH, double u, double H, double dH, double flux, double dflux)

void	calculateGlobalExteriorLesaintRaviartNumericalFlux (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int speedEvalFlag, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int inflowFlag, double n, double bc_u, double bc_H, double bc_dH, double u, double H, double dH, double flux, double dflux)

void	calculateGlobalExteriorNumericalFluxDarcyFCFF (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_um, const double n, const double bc_f_m, const double bc_a_wm, const double bc_a_mw, const double bc_a_mm, const double bc_grad_phi_w, const double bc_grad_phi_m, const double bc_u_w, const double bc_u_m, const double f_m, const double df_m_dw, const double a_wm, const double a_mw, const double a_mm, const double grad_phi_w, const double grad_phi_m, const double u_w, const double u_m, const double penalty_w, const double penalty_m, double advectiveFlux_m, double dadvectiveFlux_m_dw, double diffusiveFlux_wm, double diffusiveFlux_mw, double diffusiveFlux_mm)

void	calculateGlobalExteriorNumericalFluxDarcyFCFF_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr_wm, int colind_wm, int rowptr_mw, int colind_mw, int rowptr_mm, int colind_mm, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_um, const double n, const double bc_f_m, const double bc_a_wm, const double bc_a_mw, const double bc_a_mm, const double bc_grad_phi_w, const double bc_grad_phi_m, const double bc_u_w, const double bc_u_m, const double f_m, const double df_m_dw, const double a_wm, const double a_mw, const double a_mm, const double grad_phi_w, const double grad_phi_m, const double u_w, const double u_m, const double penalty_w, const double penalty_m, double advectiveFlux_m, double dadvectiveFlux_m_dw, double diffusiveFlux_wm, double diffusiveFlux_mw, double diffusiveFlux_mm)

void	calculateGlobalExteriorNumericalFluxDarcyFCFF_diffusiveFluxJacobian (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, const int l2g, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_um, const double n, const double f_m, const double df_m_dw, const double a_wm, const double da_wm_dw, const double da_wm_dm, const double a_mw, const double da_mw_dw, const double da_mw_dm, const double a_mm, const double da_mm_dw, const double da_mm_dm, const double grad_phi_w, const double grad_phi_m, const double dphi_w_w, const double dphi_w_m, const double dphi_m_w, const double dphi_m_m, const double u_w, const double u_m, const double v, const double grad_v, const double penalty_w, const double penalty_m, double fluxJacobian_ww, double fluxJacobian_wm, double fluxJacobian_mw, double fluxJacobian_mm)

void	calculateGlobalExteriorNumericalFluxDarcyFCFF_diffusiveFluxJacobian_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, int rowptr_wm, int colind_wm, int rowptr_mw, int colind_mw, int rowptr_mm, int colind_mm, const int l2g, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_um, const double n, const double f_m, const double df_m_dw, const double a_wm, const double da_wm_dw, const double da_wm_dm, const double a_mw, const double da_mw_dw, const double da_mw_dm, const double a_mm, const double da_mm_dw, const double da_mm_dm, const double grad_phi_w, const double grad_phi_m, const double dphi_w_w, const double dphi_w_m, const double dphi_m_w, const double dphi_m_m, const double u_w, const double u_m, const double v, const double grad_v, const double penalty_w, const double penalty_m, double fluxJacobian_ww, double fluxJacobian_wm, double fluxJacobian_mw, double fluxJacobian_mm)

void	calculateGlobalExteriorNumericalFluxDarcyFC (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double bc_a_ww, const double bc_a_nn, const double bc_grad_phi_w, const double bc_grad_phi_n, const double bc_s_w, const double bc_psi_w, const double bc_psi_n, const double a_ww, const double a_nn, const double grad_phi_w, const double grad_phi_n, const double s_w, const double psi_w, const double psi_n, const double penalty_w, const double penalty_n, double diffusiveFlux_ww, double diffusiveFlux_nn)

void	calculateGlobalExteriorNumericalFluxDarcyFC_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr_ww, int colind_ww, int rowptr_nn, int colind_nn, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double bc_a_ww, const double bc_a_nn, const double bc_grad_phi_w, const double bc_grad_phi_n, const double bc_s_w, const double bc_psi_w, const double bc_psi_n, const double a_ww, const double a_nn, const double grad_phi_w, const double grad_phi_n, const double s_w, const double psi_w, const double psi_n, const double penalty_w, const double penalty_n, double diffusiveFlux_ww, double diffusiveFlux_nn)

void	calculateGlobalExteriorNumericalFluxDarcyFC_diffusiveFluxJacobian (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, const int l2g, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double a_ww, const double da_ww_dw, const double da_ww_dn, const double a_nn, const double da_nn_dw, const double da_nn_dn, const double grad_phi_w, const double grad_phi_n, const double dphi_w_w, const double dphi_w_n, const double dphi_n_w, const double dphi_n_n, const double s_w, const double psi_w, const double psi_n, const double dpsi_n_dsw, const double dpsi_n_dpsiw, const double v, const double grad_v, const double penalty_w, const double penalty_n, double fluxJacobian_ww, double fluxJacobian_wn, double fluxJacobian_nw, double fluxJacobian_nn)

void	calculateGlobalExteriorNumericalFluxDarcyFC_diffusiveFluxJacobian_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, int rowptr_ww, int colind_ww, int rowptr_nn, int colind_nn, const int l2g, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double a_ww, const double da_ww_dw, const double da_ww_dn, const double a_nn, const double da_nn_dw, const double da_nn_dn, const double grad_phi_w, const double grad_phi_n, const double dphi_w_w, const double dphi_w_n, const double dphi_n_w, const double dphi_n_n, const double s_w, const double psi_w, const double psi_n, const double dpsi_n_dsw, const double dpsi_n_dpsiw, const double v, const double grad_v, const double penalty_w, const double penalty_n, double fluxJacobian_ww, double fluxJacobian_wn, double fluxJacobian_nw, double fluxJacobian_nn)

void	calculateGlobalExteriorNumericalAdvectiveFlux_DarcyFC (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_sw, int isDOFBoundary_psiw, double n, double bc_sw, double bc_psiw, double bc_fw, double bc_dfw_dsw, double bc_dfw_dpsiw, double bc_fn, double bc_dfn_dsw, double bc_dfn_dpsiw, double sw, double psiw, double fw, double dfw_dsw, double dfw_dpsiw, double fn, double dfn_dsw, double dfn_dpsiw, double fluxw, double dfluxw_dsw, double dfluxw_dpsiw, double fluxn, double dfluxn_dsw, double dfluxn_dpsiw)
	Calculate the advective (gravity) flux at at exterior element boundaries for Darcy FC does not upwind right now. More...

void	calculateGlobalExteriorNumericalFluxDarcyFCPP (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double bc_a_ww, const double bc_a_nn, const double bc_grad_phi_w, const double bc_grad_phi_n, const double bc_psi_w, const double bc_psi_c, const double bc_psi_n, const double a_ww, const double a_nn, const double grad_phi_w, const double grad_phi_n, const double psi_w, const double psi_c, const double psi_n, const double penalty_w, const double penalty_n, double diffusiveFlux_ww, double diffusiveFlux_nn)

void	calculateGlobalExteriorNumericalFluxDarcyFCPP_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr_ww, int colind_ww, int rowptr_nn, int colind_nn, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double bc_a_ww, const double bc_a_nn, const double bc_grad_phi_w, const double bc_grad_phi_n, const double bc_psi_w, const double bc_psi_c, const double bc_psi_n, const double a_ww, const double a_nn, const double grad_phi_w, const double grad_phi_n, const double psi_w, const double psi_c, const double psi_n, const double penalty_w, const double penalty_n, double diffusiveFlux_ww, double diffusiveFlux_nn)

void	calculateGlobalExteriorNumericalFluxDarcyFCPP_diffusiveFluxJacobian (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, const int l2g, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double a_ww, const double da_ww_dw, const double da_ww_dn, const double a_nn, const double da_nn_dw, const double da_nn_dn, const double grad_phi_w, const double grad_phi_n, const double dphi_w_w, const double dphi_w_n, const double dphi_n_w, const double dphi_n_n, const double psi_w, const double psi_c, const double psi_n, const double dpsi_n_dpsiw, const double dpsi_n_dpsic, const double v, const double grad_v, const double penalty_w, const double penalty_n, double fluxJacobian_ww, double fluxJacobian_wn, double fluxJacobian_nw, double fluxJacobian_nn)

void	calculateGlobalExteriorNumericalFluxDarcyFCPP_diffusiveFluxJacobian_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int nDOF_trial_element, int rowptr_ww, int colind_ww, int rowptr_nn, int colind_nn, const int l2g, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_uw, const int isDOFBoundary_un, int fluxBoundaryFlag_uw, int fluxBoundaryFlag_un, const double n, const double a_ww, const double da_ww_dw, const double da_ww_dn, const double a_nn, const double da_nn_dw, const double da_nn_dn, const double grad_phi_w, const double grad_phi_n, const double dphi_w_w, const double dphi_w_n, const double dphi_n_w, const double dphi_n_n, const double psi_w, const double psi_c, const double psi_n, const double dpsi_n_dpsiw, const double dpsi_n_dpsic, const double v, const double grad_v, const double penalty_w, const double penalty_n, double fluxJacobian_ww, double fluxJacobian_wn, double fluxJacobian_nw, double fluxJacobian_nn)

void	calculateGlobalExteriorNumericalFluxDarcySplitPressure (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_u, const double n, const double bc_a, const double bc_grad_phi, const double bc_psi_w, const double bc_psi_n, const double a, const double grad_phi, const double psi_w, const double psi_n, const double penalty, double *diffusiveFlux)

void	calculateGlobalExteriorNumericalFluxDarcySplitPressure_sd (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, const int rowptr, const int colind, const int exteriorElementBoundaries, const int elementBoundaryElements, const int elementBoundaryLocalElementBoundaries, const int isDOFBoundary_u, const double n, const double bc_a, const double bc_grad_phi, const double bc_psi_w, const double bc_psi_n, const double a, const double grad_phi, const double psi_w, const double psi_n, const double penalty, double *diffusiveFlux)

void	shallowWater_phi (double g, double h_l, double h_r, double u_l, double u_r, double c_l, double c_r, double h, double phi, double dphi, double u_f, int w_1, int *w_2)

void	shallowWater_Riemann (int verbose, double h_eps, double tol_u, double g, double h_l, double h_r, double hu_l, double hu_r, double h_G, double u_G)

void	calculateInteriorNumericalFluxShallowWater_1D (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, double h_eps, double tol_u, double g, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double h, double hu, double flux_h, double flux_hu)

void	calculateExteriorNumericalFluxShallowWater_1D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, double h_eps, double tol_u, double g, double n, double h_lv, double hu_lv, double h_rv, double hu_rv, double flux_h, double *flux_hu)

void	calculateInteriorNumericalFluxShallowWater_2D (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, double h_eps, double tol_u, double g, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double h, double hu, double hv, double flux_h, double flux_hu, double flux_hv)

void	calculateExteriorNumericalFluxShallowWater_2D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, double h_eps, double tol_u, double g, double n, double h_lq, double hu_lq, double hv_lq, double h_rq, double hu_rq, double hv_rq, double flux_h, double flux_hu, double flux_hv)

void	calculateInteriorNumericalFluxShallowWaterHLL_1D (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, double h_eps, double tol_u, double g, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double h, double hu, double flux_h, double flux_hu)

void	calculateExteriorNumericalFluxShallowWaterHLL_1D (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, double h_eps, double tol_u, double g, double n, double h_lv, double hu_lv, double h_rv, double hu_rv, double flux_h, double *flux_hu)

void	calculateInteriorChengShuNumericalFlux (int nInteriorElementBoundaries_global, int nElementBoundaries_element, int nQuadraturePoints_elementBoundary, int nQuadraturePoints_element, int nSpace, int speedEvalFlag, int interiorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, double n, double u, double H, double dH, double H_element, double dH_element, double flux, double dflux_left, double dflux_right)

double	smoothedHeaviside (double eps, double phi)

double	smoothedHeaviside_integral (double eps, double phi)

double	smoothedDirac (double eps, double phi)

void	applySeepageFace (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isSeepageFace, int isDOFBoundary, double epsFact, double elementDiameters, double g, double n, double grad_u, double u, double advectiveFlux, double diffusiveFlux)

void	calculateExteriorNumericalFluxRichards_sd (int rowptr, int colind, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int isSeepageFace, int isDOFBoundary, double n, double bc_u, double K, double grad_psi, double u, double K_rho_g, double penalty, double diffusiveFlux)

void	calculateExteriorNumericalFluxJacobianRichards_sd (int rowptr, int colind, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int isDOFBoundary, double n, double bc_u, double K, double dK, double grad_psi, double grad_v, double u, double dK_rho_g, double v, double penalty, double fluxJacobian)

void	applySeepageFaceJacobian (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isSeepageFace, double epsFact, double elementDiameters, double g, double n, double grad_u, double u, double advectiveFlux, double diffusiveFlux, double v, double *fluxJacobian)

void	calculateGlobalExteriorNumericalStressFlux (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_u, int isDOFBoundary_v, int isDOFBoundary_w, double n, double bc_u, double bc_v, double bc_w, double sigma, double u, double v, double w, double penalty, double stressFlux_u, double stressFlux_v, double stressFlux_w)

void	updateExteriorNumericalStressFluxJacobian (int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary_u, int isDOFBoundary_v, int isDOFBoundary_w, int isStressBoundary_u, int isStressBoundary_v, int isStressBoundary_w, double n, double dstress_u_u, double dstress_u_v, double dstress_u_w, double dstress_v_u, double dstress_v_v, double dstress_v_w, double dstress_w_u, double dstress_w_v, double dstress_w_w, double v, double grad_v, double penalty, double fluxJacobian_u_u, double fluxJacobian_u_v, double fluxJacobian_u_w, double fluxJacobian_v_u, double fluxJacobian_v_v, double fluxJacobian_v_w, double fluxJacobian_w_u, double fluxJacobian_w_v, double *fluxJacobian_w_w)

void	calculateGlobalExteriorNumericalDiffusiveFluxWithUpwinding_sd (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nSpace, int rowptr, int colind, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int fluxBoundaryFlag, double n, double bc_a, double bc_grad_phi, double bc_u, double a, double grad_phi, double u, double penalty, double flux)

void	updateGlobalExteriorNumericalDiffusiveFluxWithUpwindingJacobian_sd (int scale_penalty, double penalty_floor, int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, int nSpace, int rowptr, int colind, int l2g, int exteriorElementBoundaries, int elementBoundaryElements, int elementBoundaryLocalElementBoundaries, int isDOFBoundary, int fluxBoundaryFlag, double n, double a, double da, double grad_phi, double dphi, double v, double grad_v, double penalty, double *fluxJacobian)

Macro Definition Documentation

◆ TR_ALPHA

#define TR_ALPHA 0.5

Definition at line 3 of file numericalFlux.c.

◆ TR_ALPHA_EXT

#define TR_ALPHA_EXT 1.0

Definition at line 4 of file numericalFlux.c.

Function Documentation

◆ applySeepageFace()

void applySeepageFace	(	int	nExteriorElementBoundaries_global,
		int	nQuadraturePoints_elementBoundary,
		int	nSpace,
		int *	exteriorElementBoundaries,
		int *	elementBoundaryElements,
		int *	elementBoundaryLocalElementBoundaries,
		int *	isSeepageFace,
		int *	isDOFBoundary,
		double	epsFact,
		double *	elementDiameters,
		double *	g,
		double *	n,
		double *	grad_u,
		double *	u,
		double *	advectiveFlux,
		double *	diffusiveFlux
	)

Definition at line 14211 of file numericalFlux.c.

◆ applySeepageFaceJacobian()

void applySeepageFaceJacobian	(	int	nExteriorElementBoundaries_global,
		int	nQuadraturePoints_elementBoundary,
		int	nDOF_trial_element,
		int	nSpace,
		int *	exteriorElementBoundaries,
		int *	elementBoundaryElements,
		int *	elementBoundaryLocalElementBoundaries,
		int *	isSeepageFace,
		double	epsFact,
		double *	elementDiameters,
		double *	g,
		double *	n,
		double *	grad_u,
		double *	u,
		double *	advectiveFlux,
		double *	diffusiveFlux,
		double *	v,
		double *	fluxJacobian
	)

Definition at line 14454 of file numericalFlux.c.

◆ calculateGlobalExteriorNumericalDiffusiveFluxWithUpwinding_sd()

void calculateGlobalExteriorNumericalDiffusiveFluxWithUpwinding_sd	(	int	scale_penalty,
		double	penalty_floor,
		int	nExteriorElementBoundaries_global,
		int	nQuadraturePoints_elementBoundary,
		int	nSpace,
		int *	rowptr,
		int *	colind,
		int *	exteriorElementBoundaries,
		int *	elementBoundaryElements,
		int *	elementBoundaryLocalElementBoundaries,
		int *	isDOFBoundary,
		int *	fluxBoundaryFlag,
		double *	n,
		double *	bc_a,
		double *	bc_grad_phi,
		double *	bc_u,
		double *	a,
		double *	grad_phi,
		double *	u,
		double *	penalty,
		double *	flux
	)

Allow upwinding on the diffusion term if this actually is part of an advective process. Comes up in fully coupled Darcy flow and Transport formulations

Definition at line 14855 of file numericalFlux.c.

◆ shallowWater_phi()

void shallowWater_phi	(	double	g,
		double	h_l,
		double	h_r,
		double	u_l,
		double	u_r,
		double	c_l,
		double	c_r,
		double	h,
		double *	phi,
		double *	dphi,
		double *	u_f,
		int *	w_1,
		int *	w_2
	)

Definition at line 13163 of file numericalFlux.c.

◆ shallowWater_Riemann()

void shallowWater_Riemann	(	int	verbose,
		double	h_eps,
		double	tol_u,
		double	g,
		double	h_l,
		double	h_r,
		double	hu_l,
		double	hu_r,
		double *	h_G,
		double *	u_G
	)

Definition at line 13202 of file numericalFlux.c.

Macros

Functions