Lines Matching refs:CeedScalar
72 …R CeedInt Exact_AdvectionGeneric(CeedInt dim, CeedScalar time, const CeedScalar X[], CeedInt Nf, C… in Exact_AdvectionGeneric()
74 const CeedScalar rc = context->rc; in Exact_AdvectionGeneric()
75 const CeedScalar lx = context->lx; in Exact_AdvectionGeneric()
76 const CeedScalar ly = context->ly; in Exact_AdvectionGeneric()
77 const CeedScalar lz = dim == 2 ? 0. : context->lz; in Exact_AdvectionGeneric()
78 const CeedScalar *wind = context->wind; in Exact_AdvectionGeneric()
80 const CeedScalar center[3] = {0.5 * lx, 0.5 * ly, 0.5 * lz}; in Exact_AdvectionGeneric()
81 const CeedScalar theta = dim == 2 ? M_PI / 3 : M_PI; in Exact_AdvectionGeneric()
82 …const CeedScalar x0[3] = {center[0] + .25 * lx * cos(theta + time), center[1] + .25 * ly * sin… in Exact_AdvectionGeneric()
84 const CeedScalar x = X[0], y = X[1], z = dim == 2 ? 0. : X[2]; in Exact_AdvectionGeneric()
86 CeedScalar r = 0.; in Exact_AdvectionGeneric()
138 CeedScalar half_width = context->lx / 2; in Exact_AdvectionGeneric()
142 CeedScalar skewed_barrier[3] = {wind[0], wind[1], 0}; in Exact_AdvectionGeneric()
143 CeedScalar inflow_to_point[3] = {x - context->lx / 2, y, 0}; in Exact_AdvectionGeneric()
144 CeedScalar cross_product[3] = {0}; in Exact_AdvectionGeneric()
145 const CeedScalar boundary_threshold = 20 * CEED_EPSILON; in Exact_AdvectionGeneric()
164 CEED_QFUNCTION(ICsAdvection)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *… in ICsAdvection()
165 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsAdvection()
166 CeedScalar(*q0)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ICsAdvection()
169 const CeedScalar x[] = {X[0][i], X[1][i], X[2][i]}; in ICsAdvection()
170 CeedScalar q[5] = {0.}; in ICsAdvection()
181 CEED_QFUNCTION(ICsAdvection2d)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const… in ICsAdvection2d()
182 const CeedScalar(*X)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in ICsAdvection2d()
183 CeedScalar(*q0)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ICsAdvection2d()
187 const CeedScalar x[] = {X[0][i], X[1][i]}; in ICsAdvection2d()
188 CeedScalar q[5] = {0.}; in ICsAdvection2d()
196 …ataUnpack_ND(CeedInt N, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedSca… in QdataUnpack_ND()
210 …aryUnpack_ND(CeedInt N, CeedInt Q, CeedInt i, const CeedScalar *q_data, CeedScalar *wdetJ, CeedSca… in QdataBoundaryUnpack_ND()
211 CeedScalar *normal) { in QdataBoundaryUnpack_ND()
221 if (dXdx) StoredValuesUnpack(Q, i, 4, 6, q_data, (CeedScalar *)dXdx); in QdataBoundaryUnpack_ND()
228 … StateVariable state_var, const CeedScalar *grad_q, const CeedScalar *dXdx, in StatePhysicalGradientFromReference_ND()
233 CeedScalar dqi[5]; in StatePhysicalGradientFromReference_ND()
239 CeedScalar U[5] = {0.}; in StatePhysicalGradientFromReference_ND()
245 CeedScalar dqi[5]; in StatePhysicalGradientFromReference_ND()
257 CEED_QFUNCTION_HELPER CeedScalar Tau(AdvectionContext context, const State s, const CeedScalar *dXd… in Tau()
260 CeedScalar uX[3] = {0.}; in Tau()
266 CeedScalar gijd_mat[9] = {0.}, gij_uj[3] = {0.}; in Tau()
280 …void IFunction_AdvectionGeneric(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in IFunction_AdvectionGeneric()
281 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in IFunction_AdvectionGeneric()
282 const CeedScalar(*grad_q) = in[1]; in IFunction_AdvectionGeneric()
283 const CeedScalar(*q_dot)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[2]; in IFunction_AdvectionGeneric()
284 const CeedScalar(*q_data) = in[3]; in IFunction_AdvectionGeneric()
286 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in IFunction_AdvectionGeneric()
287 CeedScalar(*grad_v)[5][CEED_Q_VLA] = (CeedScalar(*)[5][CEED_Q_VLA])out[1]; in IFunction_AdvectionGeneric()
288 CeedScalar *jac_data = out[2]; in IFunction_AdvectionGeneric()
291 const CeedScalar zeros[14] = {0.}; in IFunction_AdvectionGeneric()
297 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in IFunction_AdvectionGeneric()
300 CeedScalar wdetJ, dXdx[9]; in IFunction_AdvectionGeneric()
305 const CeedScalar Grad_E[3] = {grad_s[0].U.E_total, grad_s[1].U.E_total, grad_s[2].U.E_total}; in IFunction_AdvectionGeneric()
312 CeedScalar div_u = 0; in IFunction_AdvectionGeneric()
318 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in IFunction_AdvectionGeneric()
319 CeedScalar strong_res = q_dot[4][i] + strong_conv; in IFunction_AdvectionGeneric()
323 CeedScalar uX[3] = {0.}; in IFunction_AdvectionGeneric()
333 CeedScalar Fe[3], Fe_dXdx[3] = {0.}; in IFunction_AdvectionGeneric()
340 const CeedScalar TauS = Tau(context, s, dXdx, dim); in IFunction_AdvectionGeneric()
355 CEED_QFUNCTION(IFunction_Advection)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *… in IFunction_Advection()
360 CEED_QFUNCTION(IFunction_Advection2d)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar… in IFunction_Advection2d()
365 …d MassFunction_AdvectionGeneric(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in MassFunction_AdvectionGeneric()
366 const CeedScalar(*q_dot)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in MassFunction_AdvectionGeneric()
367 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in MassFunction_AdvectionGeneric()
368 const CeedScalar(*q_data) = in[2]; in MassFunction_AdvectionGeneric()
370 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in MassFunction_AdvectionGeneric()
371 CeedScalar(*grad_v)[5][CEED_Q_VLA] = (CeedScalar(*)[5][CEED_Q_VLA])out[1]; in MassFunction_AdvectionGeneric()
378 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in MassFunction_AdvectionGeneric()
380 CeedScalar wdetJ, dXdx[9]; in MassFunction_AdvectionGeneric()
392 CeedScalar uX[3] = {0.}; in MassFunction_AdvectionGeneric()
394 const CeedScalar TauS = Tau(context, s, dXdx, dim); in MassFunction_AdvectionGeneric()
407 …UNCTION(MassFunction_Advection)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in MassFunction_Advection()
412 …CTION(MassFunction_Advection2D)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in MassFunction_Advection2D()
420 …id RHSFunction_AdvectionGeneric(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in RHSFunction_AdvectionGeneric()
421 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in RHSFunction_AdvectionGeneric()
422 const CeedScalar(*grad_q) = in[1]; in RHSFunction_AdvectionGeneric()
423 const CeedScalar(*q_data) = in[2]; in RHSFunction_AdvectionGeneric()
425 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in RHSFunction_AdvectionGeneric()
426 CeedScalar(*grad_v)[5][CEED_Q_VLA] = (CeedScalar(*)[5][CEED_Q_VLA])out[1]; in RHSFunction_AdvectionGeneric()
433 const CeedScalar qi[5] = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]}; in RHSFunction_AdvectionGeneric()
436 CeedScalar wdetJ, dXdx[9]; in RHSFunction_AdvectionGeneric()
441 const CeedScalar Grad_E[3] = {grad_s[0].U.E_total, grad_s[1].U.E_total, grad_s[2].U.E_total}; in RHSFunction_AdvectionGeneric()
448 CeedScalar div_u = 0; in RHSFunction_AdvectionGeneric()
454 CeedScalar strong_conv = s.U.E_total * div_u + DotN(s.Y.velocity, Grad_E, dim); in RHSFunction_AdvectionGeneric()
456 CeedScalar uX[3] = {0.}; in RHSFunction_AdvectionGeneric()
468 CeedScalar Fe[3], Fe_dXdx[3] = {0.}; in RHSFunction_AdvectionGeneric()
475 const CeedScalar TauS = Tau(context, s, dXdx, dim); in RHSFunction_AdvectionGeneric()
487 CEED_QFUNCTION(RHS_Advection)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const … in RHS_Advection()
492 CEED_QFUNCTION(RHS_Advection2d)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *cons… in RHS_Advection2d()
511 …ION(Advection_InOutFlowGeneric)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *con… in Advection_InOutFlowGeneric()
512 const CeedScalar(*q)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0]; in Advection_InOutFlowGeneric()
513 const CeedScalar(*q_data_sur) = in[2]; in Advection_InOutFlowGeneric()
515 CeedScalar(*v)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in Advection_InOutFlowGeneric()
517 const CeedScalar E_wind = context->E_wind; in Advection_InOutFlowGeneric()
518 const CeedScalar strong_form = context->strong_form; in Advection_InOutFlowGeneric()
522 const CeedScalar rho = q[0][i]; in Advection_InOutFlowGeneric()
523 const CeedScalar u[3] = {q[1][i] / rho, q[2][i] / rho, q[3][i] / rho}; in Advection_InOutFlowGeneric()
524 const CeedScalar E = q[4][i]; in Advection_InOutFlowGeneric()
526 CeedScalar wdetJb, norm[3]; in Advection_InOutFlowGeneric()
530 const CeedScalar u_normal = DotN(norm, u, dim); in Advection_InOutFlowGeneric()
546 CEED_QFUNCTION(Advection_InOutFlow)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *… in Advection_InOutFlow()
551 CEED_QFUNCTION(Advection2d_InOutFlow)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar… in Advection2d_InOutFlow()