| /libCEED/examples/fluids/qfunctions/ |
| H A D | eulervortex.h | 87 CeedScalar rho, P, T, E, u[3] = {0.}; in Exact_Euler() local
96 rho = pow(T / S_vortex, 1 / (gamma - 1.)); in Exact_Euler()
97 P = rho * T; in Exact_Euler()
102 q[0] = rho; in Exact_Euler()
103 q[1] = rho * u[0]; in Exact_Euler()
104 q[2] = rho * u[1]; in Exact_Euler()
105 q[3] = rho * u[2]; in Exact_Euler()
106 q[4] = P / (gamma - 1.) + rho * (u[0] * u[0] + u[1] * u[1]) / 2.; in Exact_Euler()
109 rho = 1.; in Exact_Euler()
113 q[0] = rho; in Exact_Euler()
[all …]
|
| H A D | shocktube.h | 90 CeedScalar rho, P, u[3] = {0.}; in Exact_ShockTube() local
94 rho = rho_high; in Exact_ShockTube()
97 rho = rho_low; in Exact_ShockTube()
102 q[0] = rho; in Exact_ShockTube()
103 q[1] = rho * u[0]; in Exact_ShockTube()
104 q[2] = rho * u[1]; in Exact_ShockTube()
105 q[3] = rho * u[2]; in Exact_ShockTube()
106 q[4] = P / (gamma - 1.0) + rho * (u[0] * u[0]) / 2.; in Exact_ShockTube()
114 …void ConvectiveFluxJacobian_Euler(CeedScalar dF[3][5][5], const CeedScalar rho, const CeedScalar u… in ConvectiveFluxJacobian_Euler() argument
123 …dF[i][4][k + 1] = ((i == k) ? (E * gamma / rho - (gamma - 1.) * u_sq / 2.) : 0.) - (gamma - 1.… in ConvectiveFluxJacobian_Euler()
[all …]
|
| H A D | channel.h | 161 const CeedScalar rho = q[0][i]; in Channel_Outflow() local
162 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Channel_Outflow()
174 v[0][i] -= wdetJb * rho * u_normal; in Channel_Outflow()
177 for (CeedInt j = 0; j < 3; j++) v[j + 1][i] -= wdetJb * (rho * u_normal * u[j] + norm[j] * P); in Channel_Outflow()
|
| H A D | stabilization.h | 81 const CeedScalar rho = s.U.density; in Tau_diagPrim() local
98 …tau = Square(rho) * (4. * Square(dts) + velocity_term) + Ctau_v * Square(mu) * DotN((CeedScalar *)… in Tau_diagPrim()
102 Tau_d[0] = Ctau_C * fact / (rho * (gijd_mat[0][0] + gijd_mat[1][1] + gijd_mat[2][2])) * 0.125; in Tau_diagPrim()
|
| H A D | stg_shur14.h | 246 const CeedScalar rho = P0 / (Rd * theta0); in StgShur14Preprocess() local
247 const CeedScalar nu = mu / rho; in StgShur14Preprocess()
293 const CeedScalar rho = P0 / (GasConstant(gas) * theta0); in ICsStg() local
294 const CeedScalar nu = gas->mu / rho; in ICsStg()
352 const CeedScalar rho = prescribe_T ? q[0][i] : P0 / (Rd * theta0); in StgShur14Inflow() local
365 CalcSpectrum(X[1][i], eps, lt, h_node_sep, mu / rho, qn, stg_ctx); in StgShur14Inflow()
371 const CeedScalar E_kinetic = .5 * rho * Dot3(u, u); in StgShur14Inflow()
375 E_internal = rho * cv * theta0; in StgShur14Inflow()
377 P = rho * Rd * theta0; // interior rho with exterior T in StgShur14Inflow()
394 v[0][i] -= wdetJb * rho * u_normal; in StgShur14Inflow()
[all …]
|
| /libCEED/tests/ |
| H A D | t502-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
14 rho[i] = weight[i] * dxdX[i]; in setup()
20 const CeedScalar *rho = in[0], *u = in[1]; in mass() local
23 v[i] = rho[i] * u[i]; in mass()
24 v[Q + i] = rho[i] * u[Q + i]; in mass()
|
| H A D | t507-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
15 rho[i] = weight[i] * dxdX[i]; in setup()
21 const CeedScalar *rho = in[0], (*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in mass() local
25 v[0][i] = rho[i] * u[0][i]; in mass()
26 v[1][i] = rho[i] * u[1][i]; in mass()
|
| H A D | t500-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
14 rho[i] = weight[i] * dxdX[i]; in setup()
20 const CeedScalar *rho = in[0], *u = in[1]; in mass() local
23 v[i] = rho[i] * u[i]; in mass()
|
| H A D | t510-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup()
20 const CeedScalar *rho = in[0], *u = in[1]; in mass() local
23 v[i] = rho[i] * u[i]; in mass()
|
| H A D | t530-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup()
20 const CeedScalar *rho = in[0], *u = in[1]; in mass() local
23 v[i] = rho[i] * u[i]; in mass()
|
| H A D | t537-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup()
20 const CeedScalar *rho = in[0], *u = in[1]; in mass() local
23 for (CeedInt c = 0; c < 2; c++) v[i + Q * c] = rho[i] * u[i + Q * (c ? 0 : 1)]; in mass()
|
| H A D | t596-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
15 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup()
22 const CeedScalar *rho = in[0], *u = in[1]; in mass() local
26 for (CeedInt c = 0; c < num_comp; c++) v[i + c * Q] = rho[i] * c * u[i + c * Q]; in mass()
|
| H A D | t595-operator.h | 11 const CeedScalar *u = in[0], *rho = in[1]; in mass() local
15 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { v[i] = rho[i] * u[i]; } in mass()
|
| H A D | t540-operator.h | 12 CeedScalar *rho = out[0]; in setup_mass() local
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup_mass()
|
| H A D | t591-operator.h | 21 const CeedScalar *u = in[0], *rho = in[1]; in mass() local
25 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { v[i] = rho[i] * u[i]; } in mass()
|
| H A D | t566-operator.h | 12 CeedScalar *rho = out[0]; in setup() local
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup()
|
| H A D | t535-operator.h | 12 CeedScalar *rho = out[0]; in setup_mass() local
14 rho[i] = weight[i] * (J[i + Q * 0] * J[i + Q * 3] - J[i + Q * 1] * J[i + Q * 2]); in setup_mass()
|
| /libCEED/python/tests/output/ |
| H A D | test_523.out | 18 Name: "rho"
38 Name: "rho"
49 Name: "rho"
69 Name: "rho"
|
| H A D | test_504.out | 17 Name: "rho"
28 Name: "rho"
|
| /libCEED/tests/output/ |
| H A D | t523-operator-f.out | 18 Name: "rho"
38 Name: "rho"
48 Name: "rho"
68 Name: "rho"
|
| H A D | t523-operator.out | 21 Name: "rho"
41 Name: "rho"
54 Name: "rho"
74 Name: "rho"
|
| H A D | t504-operator-f.out | 17 Name: "rho"
27 Name: "rho"
|
| H A D | t504-operator.out | 18 Name: "rho"
29 Name: "rho"
|
| /libCEED/examples/fluids/postprocess/ |
| H A D | vortexshedding.py | 8 def coeff(force, rho=1, u=1, D=1, zspan=0.2): argument
10 return 2 * force / (rho * u**2 * S)
|
| /libCEED/examples/petsc/qfunctions/swarm/ |
| H A D | swarmmass.h | 27 const CeedScalar *rho = (const CeedScalar *)in[0]; in Mass() local
34 for (CeedInt j = 0; j < num_comp; j++) v[j][i] = rho[i] * u[j][i]; in Mass()
|