proteus
1.8.1
C/C++/Fortran libraries
|
Go to the documentation of this file.
19 int nQuadraturePoints_element,
20 int nDOF_trial_element,
28 int nQuadraturePoints_element,
42 int nQuadraturePoints_element,
54 int nQuadraturePoints_element,
56 double *elementDiameter,
70 int nQuadraturePoints_element,
72 double *elementDiameter,
86 int nQuadraturePoints_element,
88 double *elementDiameter,
101 int nElements_global,
102 int nQuadraturePoints_element,
105 double *elementDiameter,
118 int nElements_global,
119 int nQuadraturePoints_element,
121 double *elementDiameter,
128 int nElements_global,
129 int nQuadraturePoints_element,
131 double *elementDiameter,
138 int nElements_global,
139 int nQuadraturePoints_element,
142 double *elementDiameter,
149 int nElements_global,
150 int nQuadraturePoints_element,
151 int nDOF_trial_element,
153 double *elementDiameter,
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
169 int nElements_global,
170 int nQuadraturePoints_element,
171 int nDOF_trial_element,
173 double *elementDiameter,
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
195 int nElements_global,
196 int nQuadraturePoints_element,
197 int nDOF_trial_element,
199 double *elementDiameter,
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
221 int nElements_global,
222 int nQuadraturePoints_element,
223 int nDOF_trial_element,
225 double *elementDiameter,
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
255 int nElements_global,
256 int nQuadraturePoints_element,
259 double *elementDiameter,
269 int nElements_global,
270 int nQuadraturePoints_element,
272 double *elementDiameter,
279 int nElements_global,
280 int nQuadraturePoints_element,
282 double *elementDiameter,
289 int nElements_global,
290 int nQuadraturePoints_element,
291 int nDOF_trial_element,
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
316 int nQuadraturePoints_element,
317 int nDOF_trial_element,
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);
350 int nElements_global,
351 int nQuadraturePoints_element,
352 int nDOF_trial_element,
354 double *elementDiameter,
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
386 int nElements_global,
387 int nQuadraturePoints_element,
388 int nDOF_trial_element,
390 double *elementDiameter,
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
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.
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 ...
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 ...
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.
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 ...
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 ...
double df(double C, double b, double a, int q, int r)
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 ...
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
void calculateSubgridErrorStokes_GLS_tau_sd(int nElements_global, int nQuadraturePoints_element, int nSpace, double *elementDiameter, double *dH, double *a, double *tau0, double *tau1)
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.
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 ...
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 ...
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 ...
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.
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_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 ...
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 ...
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 ...
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.
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_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 ...