proteus  1.8.1
C/C++/Fortran libraries
subgridError.h
Go to the documentation of this file.
1 #ifndef SUBGRIDERROR_H
2 #define SUBGRIDERROR_H
3 #include <math.h>
4 #include <stdio.h> /*mwf for debugging*/
5 #include <assert.h>
18  int nElements_global,
19  int nQuadraturePoints_element,
20  int nDOF_trial_element,
21  double *tau,
22  double *pdeResidual,
23  double *dpdeResidual,
24  double *subgridError,
25  double *dsubgridError
26 );
27 extern void calculateSubgridError_ADR_generic_tau(int nElements_global,
28  int nQuadraturePoints_element,
29  int nSpace,
30  double* inverseJ,
31  double* dmt,
32  double* df,
33  double* a,
34  double* da,
35  double* grad_phi,
36  double* dphi,
37  double* dr,
38  double* pe,
39  double* cfl,
40  double* tau);
41 extern void calculateSubgridErrorNavierStokes2D_generic_tau(int nElements_global,
42  int nQuadraturePoints_element,
43  int nSpace,
44  double* inverseJ,
45  double* dmt,
46  double* dm,
47  double* f,
48  double* a,
49  double* tau0,
50  double* tau1,
51  double* cfl);
53  int nElements_global,
54  int nQuadraturePoints_element,
55  int nSpace,
56  double *elementDiameter,
57  double *dmt,
58  double *df,
59  double *a,
60  double *da,
61  double *grad_phi,
62  double *dphi,
63  double *dr,
64  double *pe,
65  double *cfl,
66  double *tau
67 );
69  int nElements_global,
70  int nQuadraturePoints_element,
71  int nSpace,
72  double *elementDiameter,
73  double *dmt,
74  double *df,
75  double *a,
76  double *da,
77  double *grad_phi,
78  double *dphi,
79  double *dr,
80  double *pe,
81  double *cfl,
82  double *tau
83 );
85  int nElements_global,
86  int nQuadraturePoints_element,
87  int nSpace,
88  double *elementDiameter,
89  double *dmt,
90  double *df,
91  double *a,
92  double *da,
93  double *grad_phi,
94  double *dphi,
95  double *dr,
96  double *pe,
97  double *cfl,
98  double *tau
99 );
101  int nElements_global,
102  int nQuadraturePoints_element,
103  int nSpace,
104  char stabilization,
105  double *elementDiameter,
106  double *dmt,
107  double *df,
108  double *a,
109  double *da,
110  double *grad_phi,
111  double *dphi,
112  double *dr,
113  double *pe,
114  double *cfl,
115  double *tau
116 );
118  int nElements_global,
119  int nQuadraturePoints_element,
120  int nSpace,
121  double *elementDiameter,
122  double *dmt,
123  double *df,
124  double *cfl,
125  double *tau
126 );
128  int nElements_global,
129  int nQuadraturePoints_element,
130  int nSpace,
131  double *elementDiameter,
132  double *dmt,
133  double *df,
134  double *cfl,
135  double *tau
136 );
137 extern void calculateSubgridError_A_tau(
138  int nElements_global,
139  int nQuadraturePoints_element,
140  int nSpace,
141  char stabilization,
142  double *elementDiameter,
143  double *dmt,
144  double *df,
145  double *cfl,
146  double *tau
147 );
149  int nElements_global,
150  int nQuadraturePoints_element,
151  int nDOF_trial_element,
152  int nSpace,
153  double *elementDiameter,
154  double *a,
155  double *pdeResidualU,
156  double *dpdeResidualU_dp,
157  double *dpdeResidualU_du,
158  double *pdeResidualV,
159  double *dpdeResidualV_dp,
160  double *dpdeResidualV_dv,
161  double *subgridErrorU,
162  double *dsubgridErrorU_dp,
163  double *dsubgridErrorU_du,
164  double *subgridErrorV,
165  double *dsubgridErrorV_dp,
166  double *dsubgridErrorV_dv
167 );
169  int nElements_global,
170  int nQuadraturePoints_element,
171  int nDOF_trial_element,
172  int nSpace,
173  double *elementDiameter,
174  double *a,
175  double *pdeResidualU,
176  double *dpdeResidualU_dp,
177  double *dpdeResidualU_du,
178  double *pdeResidualV,
179  double *dpdeResidualV_dp,
180  double *dpdeResidualV_dv,
181  double *pdeResidualW,
182  double *dpdeResidualW_dp,
183  double *dpdeResidualW_dw,
184  double *subgridErrorU,
185  double *dsubgridErrorU_dp,
186  double *dsubgridErrorU_du,
187  double *subgridErrorV,
188  double *dsubgridErrorV_dp,
189  double *dsubgridErrorV_dv,
190  double *subgridErrorW,
191  double *dsubgridErrorW_dp,
192  double *dsubgridErrorW_dw
193 );
195  int nElements_global,
196  int nQuadraturePoints_element,
197  int nDOF_trial_element,
198  int nSpace,
199  double *elementDiameter,
200  double *a,
201  double *pdeResidualP,
202  double *dpdeResidualP_du,
203  double *dpdeResidualP_dv,
204  double *pdeResidualU,
205  double *dpdeResidualU_dp,
206  double *dpdeResidualU_du,
207  double *pdeResidualV,
208  double *dpdeResidualV_dp,
209  double *dpdeResidualV_dv,
210  double *subgridErrorP,
211  double *dsubgridErrorP_du,
212  double *dsubgridErrorP_dv,
213  double *subgridErrorU,
214  double *dsubgridErrorU_dp,
215  double *dsubgridErrorU_du,
216  double *subgridErrorV,
217  double *dsubgridErrorV_dp,
218  double *dsubgridErrorV_dv
219 );
221  int nElements_global,
222  int nQuadraturePoints_element,
223  int nDOF_trial_element,
224  int nSpace,
225  double *elementDiameter,
226  double *a,
227  double *pdeResidualP,
228  double *dpdeResidualP_du,
229  double *dpdeResidualP_dv,
230  double *dpdeResidualP_dw,
231  double *pdeResidualU,
232  double *dpdeResidualU_dp,
233  double *dpdeResidualU_du,
234  double *pdeResidualV,
235  double *dpdeResidualV_dp,
236  double *dpdeResidualV_dv,
237  double *pdeResidualW,
238  double *dpdeResidualW_dp,
239  double *dpdeResidualW_dw,
240  double *subgridErrorP,
241  double *dsubgridErrorP_du,
242  double *dsubgridErrorP_dv,
243  double *dsubgridErrorP_dw,
244  double *subgridErrorU,
245  double *dsubgridErrorU_dp,
246  double *dsubgridErrorU_du,
247  double *subgridErrorV,
248  double *dsubgridErrorV_dp,
249  double *dsubgridErrorV_dv,
250  double *subgridErrorW,
251  double *dsubgridErrorW_dp,
252  double *dsubgridErrorW_dw
253 );
255  int nElements_global,
256  int nQuadraturePoints_element,
257  int nSpace,
258  double hFactor,
259  double *elementDiameter,
260  double *dmt,
261  double *dm,
262  double *f,
263  double *a,
264  double *tau0,
265  double *tau1,
266  double *cfl
267 );
269  int nElements_global,
270  int nQuadraturePoints_element,
271  int nSpace,
272  double *elementDiameter,
273  double *dH,
274  double *a,
275  double *tau0,
276  double *tau1
277 );
279  int nElements_global,
280  int nQuadraturePoints_element,
281  int nSpace,
282  double *elementDiameter,
283  double *dH,
284  double *a,
285  double *tau0,
286  double *tau1
287 );
289  int nElements_global,
290  int nQuadraturePoints_element,
291  int nDOF_trial_element,
292  int nSpace,
293  double *tau0,
294  double *tau1,
295  double *pdeResidualP,
296  double *dpdeResidualP_du,
297  double *dpdeResidualP_dv,
298  double *pdeResidualU,
299  double *dpdeResidualU_dp,
300  double *dpdeResidualU_du,
301  double *pdeResidualV,
302  double *dpdeResidualV_dp,
303  double *dpdeResidualV_dv,
304  double *subgridErrorP,
305  double *dsubgridErrorP_du,
306  double *dsubgridErrorP_dv,
307  double *subgridErrorU,
308  double *dsubgridErrorU_dp,
309  double *dsubgridErrorU_du,
310  double *subgridErrorV,
311  double *dsubgridErrorV_dp,
312  double *dsubgridErrorV_dv
313 );
314 
315 extern void calculateSubgridErrorStokes3D_GLS_tauRes(int nElements_global,
316  int nQuadraturePoints_element,
317  int nDOF_trial_element,
318  int nSpace,
319  double* tau0,
320  double* tau1,
321  double* pdeResidualP,
322  double* dpdeResidualP_du,
323  double* dpdeResidualP_dv,
324  double* dpdeResidualP_dw,
325  double* pdeResidualU,
326  double* dpdeResidualU_dp,
327  double* dpdeResidualU_du,
328  double* pdeResidualV,
329  double* dpdeResidualV_dp,
330  double* dpdeResidualV_dv,
331  double* pdeResidualW,
332  double* dpdeResidualW_dp,
333  double* dpdeResidualW_dw,
334  double* subgridErrorP,
335  double* dsubgridErrorP_du,
336  double* dsubgridErrorP_dv,
337  double* dsubgridErrorP_dw,
338  double* subgridErrorU,
339  double* dsubgridErrorU_dp,
340  double* dsubgridErrorU_du,
341  double* subgridErrorV,
342  double* dsubgridErrorV_dp,
343  double* dsubgridErrorV_dv,
344  double* subgridErrorW,
345  double* dsubgridErrorW_dp,
346  double* dsubgridErrorW_dw);
347 
348 
350  int nElements_global,
351  int nQuadraturePoints_element,
352  int nDOF_trial_element,
353  int nSpace,
354  double *elementDiameter,
355  double *dm,
356  double *f,
357  double *a,
358  double *pdeResidualP,
359  double *dpdeResidualP_du,
360  double *dpdeResidualP_dv,
361  double *dpdeResidualP_dw,
362  double *pdeResidualU,
363  double *dpdeResidualU_dp,
364  double *dpdeResidualU_du,
365  double *pdeResidualV,
366  double *dpdeResidualV_dp,
367  double *dpdeResidualV_dv,
368  double *pdeResidualW,
369  double *dpdeResidualW_dp,
370  double *dpdeResidualW_dw,
371  double *subgridErrorP,
372  double *dsubgridErrorP_du,
373  double *dsubgridErrorP_dv,
374  double *dsubgridErrorP_dw,
375  double *subgridErrorU,
376  double *dsubgridErrorU_dp,
377  double *dsubgridErrorU_du,
378  double *subgridErrorV,
379  double *dsubgridErrorV_dp,
380  double *dsubgridErrorV_dv,
381  double *subgridErrorW,
382  double *dsubgridErrorW_dp,
383  double *dsubgridErrorW_dw
384 );
386  int nElements_global,
387  int nQuadraturePoints_element,
388  int nDOF_trial_element,
389  int nSpace,
390  double *elementDiameter,
391  double *u,
392  double *v,
393  double *a,
394  double *pdeResidualP,
395  double *dpdeResidualP_du,
396  double *dpdeResidualP_dv,
397  double *pdeResidualU,
398  double *dpdeResidualU_dp,
399  double *dpdeResidualU_du,
400  double *pdeResidualV,
401  double *dpdeResidualV_dp,
402  double *dpdeResidualV_dv,
403  double *subgridErrorP,
404  double *dsubgridErrorP_dp,
405  double *dsubgridErrorP_du,
406  double *dsubgridErrorP_dv,
407  double *subgridErrorU,
408  double *dsubgridErrorU_dp,
409  double *dsubgridErrorU_du,
410  double *dsubgridErrorU_dv,
411  double *subgridErrorV,
412  double *dsubgridErrorV_dp,
413  double *dsubgridErrorV_du,
414  double *dsubgridErrorV_dv
415 );
417 #endif
calculateSubgridError_A_tau
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.
Definition: subgridError.c:971
calculateSubgridError_ADR_tau_1
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 ...
Definition: subgridError.c:224
calculateSubgridErrorStokes2D_GLS_tauRes
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 ...
Definition: subgridError.c:2685
calculateSubgridError_tauRes
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.
Definition: subgridError.c:12
f
Double f
Definition: Headers.h:64
calculateSubgridErrorStokes_GLS_tau
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 ...
Definition: subgridError.c:2632
calculateSubgridErrorNavierStokes3D_GLS_velocity_pressure
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 ...
Definition: subgridError.c:3163
df
double df(double C, double b, double a, int q, int r)
Definition: analyticalSolutions.c:2209
calculateSubgridErrorStokes3D_GLS_velocity_pressure
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 ...
Definition: subgridError.c:1617
calculateSubgridErrorStokes3D_GLS_tauRes
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
Definition: subgridError.c:2797
calculateSubgridErrorStokes_GLS_tau_sd
void calculateSubgridErrorStokes_GLS_tau_sd(int nElements_global, int nQuadraturePoints_element, int nSpace, double *elementDiameter, double *dH, double *a, double *tau0, double *tau1)
Definition: subgridError.c:2657
calculateSubgridErrorStokes3D_GLS_velocity
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.
Definition: subgridError.c:1157
v
Double v
Definition: Headers.h:95
calculateSubgridError_ADR_tau_2
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 ...
Definition: subgridError.c:394
pe
Double pe
Definition: Headers.h:75
u
Double u
Definition: Headers.h:89
calculateSubgridErrorStokes2D_1
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 ...
Definition: subgridError.c:3464
calculateSubgridError_ADR_tau_p
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 ...
Definition: subgridError.c:49
calculateSubgridErrorStokes2D_GLS_velocity
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.
Definition: subgridError.c:990
calculateSubgridError_ADR_generic_tau
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)
Definition: subgridError.c:580
calculateSubgridErrorNavierStokes2D_GLS_tau
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 ...
Definition: subgridError.c:1914
calculateSubgridError_A_tau_2
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 ...
Definition: subgridError.c:923
calculateSubgridErrorStokes2D_GLS_velocity_pressure
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 ...
Definition: subgridError.c:1382
calculateSubgridError_ADR_tau
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.
Definition: subgridError.c:826
calculateSubgridErrorNavierStokes2D_generic_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)
Definition: subgridError.c:2125
calculateSubgridError_A_tau_1
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 ...
Definition: subgridError.c:880