Lines Matching refs:f0
175 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_quadratic_v() argument
179 …f0[0] -= t * (2 * X[0] + 2 * X[1]) + 2 * X[0] * X[0] * X[0] + 4 * X[0] * X[0] * X[1] - 2 * X[0] * … in f0_quadratic_v()
180 …f0[1] -= t * (2 * X[0] - 2 * X[1]) + 4 * X[0] * X[1] * X[1] + 2 * X[0] * X[0] * X[1] - 2 * X[1] * … in f0_quadratic_v()
183 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_quadratic_w() argument
185 f0[0] -= 2 * t + 1 + 3 * X[0] * X[0] - 2 * X[0] * X[1] + X[1] * X[1]; in f0_quadratic_w()
219 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_quadratic_v() argument
229 …f0[0] -= rho * S + rho * (2. * t * (X[0] + X[1]) + 2. * X[0] * X[0] * X[0] + 4. * X[0] * X[0] * X[… in f0_conduct_quadratic_v()
230 …f0[1] -= rho * S + rho * (2. * t * (X[0] - X[1]) + 2. * X[0] * X[0] * X[1] + 4. * X[0] * X[1] * X[… in f0_conduct_quadratic_v()
231 f0[gd] -= rho / PetscSqr(F); in f0_conduct_quadratic_v()
234 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_quadratic_q() argument
239 …f0[0] += p_th * (S + 2. * t + 3. * X[0] * X[0] - 2. * X[0] * X[1] + X[1] * X[1]) / PetscSqr(t + X[… in f0_conduct_quadratic_q()
242 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_quadratic_w() argument
248 …f0[0] -= c_p * p_th * (S + 2. * t + 3. * X[0] * X[0] - 2. * X[0] * X[1] + X[1] * X[1]) / (t + X[0]… in f0_conduct_quadratic_w()
302 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_cubic_v() argument
306 …f0[0] -= (t * (3 * X[0] * X[0] + 3 * X[1] * X[1]) + 3 * X[0] * X[0] * X[0] * X[0] * X[0] + 6 * X[0… in f0_cubic_v()
307 …f0[1] -= (t * (3 * X[0] * X[0] - 6 * X[0] * X[1]) + 3 * X[0] * X[0] * X[0] * X[0] * X[1] + 6 * X[0… in f0_cubic_v()
310 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_cubic_w() argument
314 …f0[0] -= X[0] * X[0] * X[0] * X[0] + 2.0 * X[0] * X[0] * X[0] * X[1] - 3.0 * X[0] * X[0] * X[1] * … in f0_cubic_w()
374 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_cubic_trig_v() argument
378 …f0[0] -= 100. * PetscCosReal(t) * (3 * X[0] * X[0]) + 100. * PetscSinReal(t) * (3 * X[1] * X[1] - … in f0_cubic_trig_v()
379 …f0[1] -= 100. * PetscCosReal(t) * (6 * X[0] * X[0] - 6 * X[0] * X[1]) - 100. * PetscSinReal(t) * (… in f0_cubic_trig_v()
382 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_cubic_trig_w() argument
386 …f0[0] -= 100. * PetscCosReal(t) * X[0] + 100. * PetscSinReal(t) * (X[1] - 1.) + X[0] * X[0] * X[0]… in f0_cubic_trig_w()
478 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_taylor_green_w() argument
483 f0[0] -= 1.0 + vel[0] + vel[1]; in f0_taylor_green_w()
567 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_pipe_v() argument
575 f0[gd] -= rho / PetscSqr(F); in f0_conduct_pipe_v()
578 …l X[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_bd_pipe_v() argument
584 for (e = 0; e < dim; ++e) f0[d] -= sigma[d * dim + e] * n[e]; in f0_conduct_bd_pipe_v()
588 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_pipe_q() argument
590 f0[0] += 0.0; in f0_conduct_pipe_q()
593 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_pipe_w() argument
598 f0[0] -= 2 * k / Pe; in f0_conduct_pipe_w()
683 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_pipe_wiggly_v() argument
692 f0[0] -= eps * 0.5 * PetscSqr(PETSC_PI) * PetscSinReal(PETSC_PI * X[1]); in f0_conduct_pipe_wiggly_v()
693 f0[gd] -= rho / PetscSqr(F); in f0_conduct_pipe_wiggly_v()
696 …l X[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_bd_pipe_wiggly_v() argument
703 for (e = 0; e < dim; ++e) f0[d] -= sigma[d * dim + e] * n[e]; in f0_conduct_bd_pipe_wiggly_v()
707 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_pipe_wiggly_q() argument
709 f0[0] += 0.0; in f0_conduct_pipe_wiggly_q()
712 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_pipe_wiggly_w() argument
718 f0[0] -= 2 * k / Pe * (1.0 + eps * 0.5 * PetscSqr(PETSC_PI) * PetscSinReal(PETSC_PI * X[1])); in f0_conduct_pipe_wiggly_w()
723 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_q() argument
726 for (d = 0, f0[0] = 0.0; d < dim; ++d) f0[0] += u_x[d * dim + d]; in f0_q()
730 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_q() argument
737 f0[0] += -u_t[uOff[TEMP]] * S * p_th / PetscSqr(u[uOff[TEMP]]); in f0_conduct_q()
740 for (d = 0; d < dim; ++d) f0[0] += p_th / u[uOff[TEMP]] * u_x[uOff_x[VEL] + d * dim + d]; in f0_conduct_q()
743 …for (d = 0; d < dim; ++d) f0[0] -= p_th / (u[uOff[TEMP]] * u[uOff[TEMP]]) * u[uOff[VEL] + d] * u_x… in f0_conduct_q()
746 if (NfAux > 0) f0[0] += a[aOff[MASS]]; in f0_conduct_q()
750 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_v() argument
757 f0[c] += u_t[uOff[VEL] + c]; in f0_v()
759 for (d = 0; d < dim; ++d) f0[c] += u[uOff[VEL] + d] * u_x[uOff_x[VEL] + c * dim + d]; in f0_v()
764 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_v() argument
775 for (d = 0; d < dim; ++d) f0[d] = rho * S * u_t[uOff[VEL] + d]; in f0_conduct_v()
779 for (d = 0; d < dim; ++d) f0[c] += rho * u[uOff[VEL] + d] * u_x[uOff_x[VEL] + c * dim + d]; in f0_conduct_v()
783 f0[gdir] += rho / (F * F); in f0_conduct_v()
787 for (d = 0; d < dim; ++d) f0[d] += a[aOff[MOMENTUM] + d]; in f0_conduct_v()
828 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_w() argument
833 f0[0] += u_t[uOff[TEMP]]; in f0_w()
835 for (d = 0; d < dim; ++d) f0[0] += u[uOff[VEL] + d] * u_x[uOff_x[TEMP] + d]; in f0_w()
839 …eal t, const PetscReal X[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[]) in f0_conduct_w() argument
847 f0[0] = c_p * S * p_th / u[uOff[TEMP]] * u_t[uOff[TEMP]]; in f0_conduct_w()
850 …for (d = 0; d < dim; ++d) f0[0] += c_p * p_th / u[uOff[TEMP]] * u[uOff[VEL] + d] * u_x[uOff_x[TEMP… in f0_conduct_w()
853 if (NfAux > 0) f0[0] += a[aOff[ENERGY]]; in f0_conduct_w()