C implementations of subgrid error calculations. More...

#include <math.h>
#include <stdio.h>
#include <assert.h>

Include dependency graph for subgridError.h:

This graph shows which files directly or indirectly include this file:

Functions
void	calculateSubgridError_tauRes (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, double tau, double pdeResidual, double dpdeResidual, double subgridError, double *dsubgridError)
	Calculate the ASGS subgrid error given tau and the strong residual. More...

void	calculateSubgridError_ADR_generic_tau (int nElements_global, int nQuadraturePoints_element, int nSpace, double inverseJ, double dmt, double df, double a, double da, double grad_phi, double dphi, double dr, double pe, double cfl, double *tau)

void	calculateSubgridErrorNavierStokes2D_generic_tau (int nElements_global, int nQuadraturePoints_element, int nSpace, double inverseJ, double dmt, double dm, double f, double a, double tau0, double tau1, double cfl)

void	calculateSubgridError_ADR_tau_p (int nElements_global, int nQuadraturePoints_element, int nSpace, double elementDiameter, double dmt, double df, double a, double da, double grad_phi, double dphi, double dr, double pe, double cfl, double *tau)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the peclet number formula. More...

void	calculateSubgridError_ADR_tau_1 (int nElements_global, int nQuadraturePoints_element, int nSpace, double elementDiameter, double dmt, double df, double a, double da, double grad_phi, double dphi, double dr, double pe, double cfl, double *tau)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

void	calculateSubgridError_ADR_tau_2 (int nElements_global, int nQuadraturePoints_element, int nSpace, double elementDiameter, double dmt, double df, double a, double da, double grad_phi, double dphi, double dr, double pe, double cfl, double *tau)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_2 norm" formula. More...

void	calculateSubgridError_ADR_tau (int nElements_global, int nQuadraturePoints_element, int nSpace, char stabilization, double elementDiameter, double dmt, double df, double a, double da, double grad_phi, double dphi, double dr, double pe, double cfl, double *tau)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation. More...

void	calculateSubgridError_A_tau_1 (int nElements_global, int nQuadraturePoints_element, int nSpace, double elementDiameter, double dmt, double df, double cfl, double *tau)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_2 norm" formula. More...

void	calculateSubgridError_A_tau_2 (int nElements_global, int nQuadraturePoints_element, int nSpace, double elementDiameter, double dmt, double df, double cfl, double *tau)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_2 norm" formula. More...

void	calculateSubgridError_A_tau (int nElements_global, int nQuadraturePoints_element, int nSpace, char stabilization, double elementDiameter, double dmt, double df, double cfl, double *tau)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation. More...

void	calculateSubgridErrorStokes2D_GLS_velocity (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double elementDiameter, double a, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_dv)
	Calculate the subgrid error for velocity in 2D Stokes equation with a GLS-like formula. More...

void	calculateSubgridErrorStokes3D_GLS_velocity (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double elementDiameter, double a, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double pdeResidualW, double dpdeResidualW_dp, double dpdeResidualW_dw, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_dv, double subgridErrorW, double dsubgridErrorW_dp, double dsubgridErrorW_dw)
	Calculate the subgrid error for velocity in 3D Stokes equation with a GLS-like formula. More...

void	calculateSubgridErrorStokes2D_GLS_velocity_pressure (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double elementDiameter, double a, double pdeResidualP, double dpdeResidualP_du, double dpdeResidualP_dv, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double subgridErrorP, double dsubgridErrorP_du, double dsubgridErrorP_dv, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_dv)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

void	calculateSubgridErrorStokes3D_GLS_velocity_pressure (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double elementDiameter, double a, double pdeResidualP, double dpdeResidualP_du, double dpdeResidualP_dv, double dpdeResidualP_dw, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double pdeResidualW, double dpdeResidualW_dp, double dpdeResidualW_dw, double subgridErrorP, double dsubgridErrorP_du, double dsubgridErrorP_dv, double dsubgridErrorP_dw, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_dv, double subgridErrorW, double dsubgridErrorW_dp, double dsubgridErrorW_dw)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

void	calculateSubgridErrorNavierStokes2D_GLS_tau (int nElements_global, int nQuadraturePoints_element, int nSpace, double hFactor, double elementDiameter, double dmt, double dm, double f, double a, double tau0, double tau1, double cfl)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

void	calculateSubgridErrorStokes_GLS_tau (int nElements_global, int nQuadraturePoints_element, int nSpace, double elementDiameter, double dH, double a, double tau0, double *tau1)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

void	calculateSubgridErrorStokes_GLS_tau_sd (int nElements_global, int nQuadraturePoints_element, int nSpace, double elementDiameter, double dH, double a, double tau0, double *tau1)

void	calculateSubgridErrorStokes2D_GLS_tauRes (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double tau0, double tau1, double pdeResidualP, double dpdeResidualP_du, double dpdeResidualP_dv, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double subgridErrorP, double dsubgridErrorP_du, double dsubgridErrorP_dv, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_dv)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

void	calculateSubgridErrorStokes3D_GLS_tauRes (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double tau0, double tau1, double pdeResidualP, double dpdeResidualP_du, double dpdeResidualP_dv, double dpdeResidualP_dw, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double pdeResidualW, double dpdeResidualW_dp, double dpdeResidualW_dw, double subgridErrorP, double dsubgridErrorP_du, double dsubgridErrorP_dv, double dsubgridErrorP_dw, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_dv, double subgridErrorW, double dsubgridErrorW_dp, double dsubgridErrorW_dw)
	3D version of Stokes GLS tau More...

void	calculateSubgridErrorNavierStokes3D_GLS_velocity_pressure (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double elementDiameter, double dm, double f, double a, double pdeResidualP, double dpdeResidualP_du, double dpdeResidualP_dv, double dpdeResidualP_dw, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double pdeResidualW, double dpdeResidualW_dp, double dpdeResidualW_dw, double subgridErrorP, double dsubgridErrorP_du, double dsubgridErrorP_dv, double dsubgridErrorP_dw, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_dv, double subgridErrorW, double dsubgridErrorW_dp, double dsubgridErrorW_dw)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

void	calculateSubgridErrorStokes2D_1 (int nElements_global, int nQuadraturePoints_element, int nDOF_trial_element, int nSpace, double elementDiameter, double u, double v, double a, double pdeResidualP, double dpdeResidualP_du, double dpdeResidualP_dv, double pdeResidualU, double dpdeResidualU_dp, double dpdeResidualU_du, double pdeResidualV, double dpdeResidualV_dp, double dpdeResidualV_dv, double subgridErrorP, double dsubgridErrorP_dp, double dsubgridErrorP_du, double dsubgridErrorP_dv, double subgridErrorU, double dsubgridErrorU_dp, double dsubgridErrorU_du, double dsubgridErrorU_dv, double subgridErrorV, double dsubgridErrorV_dp, double dsubgridErrorV_du, double *dsubgridErrorV_dv)
	Calculate the stabilization parameter for the scalar advection-diffusion-reaction equation using the "l_1 norm formula". More...

Functions