Lines Matching refs:X
145 static PetscErrorCode quadratic_u(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, P… in quadratic_u() argument
147 u[0] = time + X[0] * X[0] + X[1] * X[1]; in quadratic_u()
148 u[1] = time + 2.0 * X[0] * X[0] - 2.0 * X[0] * X[1]; in quadratic_u()
151 static PetscErrorCode quadratic_u_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf,… in quadratic_u_t() argument
158 static PetscErrorCode quadratic_p(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, P… in quadratic_p() argument
160 p[0] = X[0] + X[1] - 1.0; in quadratic_p()
164 static PetscErrorCode quadratic_T(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, P… in quadratic_T() argument
166 T[0] = time + X[0] + X[1] + 1.0; in quadratic_T()
169 static PetscErrorCode quadratic_T_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf,… in quadratic_T_t() argument
175 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… 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 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… 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 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_quadratic_v() argument
226 const PetscReal rho = p_th / (t + X[0] + X[1] + 1.); in f0_conduct_quadratic_v()
229 …0] -= rho * S + rho * (2. * t * (X[0] + X[1]) + 2. * X[0] * X[0] * X[0] + 4. * X[0] * X[0] * X[1] … in f0_conduct_quadratic_v()
230 …1] -= rho * S + rho * (2. * t * (X[0] - X[1]) + 2. * X[0] * X[0] * X[1] + 4. * X[0] * X[1] * X[1] … in f0_conduct_quadratic_v()
234 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… 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 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… 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()
272 static PetscErrorCode cubic_u(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, Petsc… in cubic_u() argument
274 u[0] = time + X[0] * X[0] * X[0] + X[1] * X[1] * X[1]; in cubic_u()
275 u[1] = time + 2.0 * X[0] * X[0] * X[0] - 3.0 * X[0] * X[0] * X[1]; in cubic_u()
278 static PetscErrorCode cubic_u_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, Pet… in cubic_u_t() argument
285 static PetscErrorCode cubic_p(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, Petsc… in cubic_p() argument
287 p[0] = 3.0 * X[0] * X[0] / 2.0 + 3.0 * X[1] * X[1] / 2.0 - 1.0; in cubic_p()
291 static PetscErrorCode cubic_T(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, Petsc… in cubic_T() argument
293 T[0] = time + X[0] * X[0] / 2.0 + X[1] * X[1] / 2.0; in cubic_T()
296 static PetscErrorCode cubic_T_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, Pet… in cubic_T_t() argument
302 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_cubic_v() argument
306 …X[0] * X[0] + 3 * X[1] * X[1]) + 3 * X[0] * X[0] * X[0] * X[0] * X[0] + 6 * X[0] * X[0] * X[0] * X… in f0_cubic_v()
307 …X[0] * X[0] - 6 * X[0] * X[1]) + 3 * X[0] * X[0] * X[0] * X[0] * X[1] + 6 * X[0] * X[0] * X[1] * X… in f0_cubic_v()
310 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_cubic_w() argument
314 …0[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] * X… in f0_cubic_w()
344 static PetscErrorCode cubic_trig_u(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, … in cubic_trig_u() argument
346 u[0] = 100. * PetscCosReal(time) + X[0] * X[0] * X[0] + X[1] * X[1] * X[1]; in cubic_trig_u()
347 u[1] = 100. * PetscSinReal(time) + 2.0 * X[0] * X[0] * X[0] - 3.0 * X[0] * X[0] * X[1]; in cubic_trig_u()
350 static PetscErrorCode cubic_trig_u_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf… in cubic_trig_u_t() argument
357 static PetscErrorCode cubic_trig_p(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, … in cubic_trig_p() argument
359 p[0] = 3.0 * X[0] * X[0] / 2.0 + 3.0 * X[1] * X[1] / 2.0 - 1.0; in cubic_trig_p()
363 static PetscErrorCode cubic_trig_T(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf, … in cubic_trig_T() argument
365 T[0] = 100. * PetscCosReal(time) + X[0] * X[0] / 2.0 + X[1] * X[1] / 2.0; in cubic_trig_T()
368 static PetscErrorCode cubic_trig_T_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf… in cubic_trig_T_t() argument
374 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_cubic_trig_v() argument
378 …X[0] * X[0]) + 100. * PetscSinReal(t) * (3 * X[1] * X[1] - 1.) + 3 * X[0] * X[0] * X[0] * X[0] * X… in f0_cubic_trig_v()
379 …X[0] * X[0] - 6 * X[0] * X[1]) - 100. * PetscSinReal(t) * (3 * X[0] * X[0]) + 3 * X[0] * X[0] * X[… in f0_cubic_trig_v()
382 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_cubic_trig_w() argument
386 …X[0] + 100. * PetscSinReal(t) * (X[1] - 1.) + X[0] * X[0] * X[0] * X[0] + 2.0 * X[0] * X[0] * X[0]… in f0_cubic_trig_w()
442 static PetscErrorCode taylor_green_u(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf… in taylor_green_u() argument
444 …u[0] = 1 - PetscCosReal(PETSC_PI * (X[0] - time)) * PetscSinReal(PETSC_PI * (X[1] - time)) * Petsc… in taylor_green_u()
445 …u[1] = 1 + PetscSinReal(PETSC_PI * (X[0] - time)) * PetscCosReal(PETSC_PI * (X[1] - time)) * Petsc… in taylor_green_u()
448 static PetscErrorCode taylor_green_u_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt … in taylor_green_u_t() argument
450 …X[0] - time)) * PetscSinReal(PETSC_PI * (X[1] - time)) - PetscCosReal(PETSC_PI * (X[0] - time)) * … in taylor_green_u_t()
451 …X[0] - time)) * PetscSinReal(PETSC_PI * (X[1] - time)) - PetscCosReal(PETSC_PI * (X[0] - time)) * … in taylor_green_u_t()
455 static PetscErrorCode taylor_green_p(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf… in taylor_green_p() argument
457 …p[0] = -0.25 * (PetscCosReal(2 * PETSC_PI * (X[0] - time)) + PetscCosReal(2 * PETSC_PI * (X[1] - t… in taylor_green_p()
461 static PetscErrorCode taylor_green_p_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt … in taylor_green_p_t() argument
463 … (X[0] - time)) + PetscSinReal(2 * PETSC_PI * (X[1] - time))) + PETSC_PI * (PetscCosReal(2 * PETSC… in taylor_green_p_t()
467 static PetscErrorCode taylor_green_T(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt Nf… in taylor_green_T() argument
469 T[0] = time + X[0] + X[1]; in taylor_green_T()
472 static PetscErrorCode taylor_green_T_t(PetscInt Dim, PetscReal time, const PetscReal X[], PetscInt … in taylor_green_T_t() argument
478 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_taylor_green_w() argument
482 PetscCallAbort(PETSC_COMM_SELF, taylor_green_u(dim, t, X, Nf, vel, NULL)); in f0_taylor_green_w()
528 static PetscErrorCode pipe_u(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf, PetscS… in pipe_u() argument
532 u[0] = (0.5 * param->Reynolds / param->mu) * X[1] * (1.0 - X[1]); in pipe_u()
536 static PetscErrorCode pipe_u_t(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf, Pets… in pipe_u_t() argument
543 static PetscErrorCode pipe_p(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf, PetscS… in pipe_p() argument
545 p[0] = -X[0]; in pipe_p()
548 static PetscErrorCode pipe_p_t(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf, Pets… in pipe_p_t() argument
554 static PetscErrorCode pipe_T(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf, PetscS… in pipe_T() argument
558 T[0] = X[1] * (1.0 - X[1]) + param->T_in; in pipe_T()
561 static PetscErrorCode pipe_T_t(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf, Pets… in pipe_T_t() argument
567 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_pipe_v() argument
572 const PetscReal rho = p_th / (X[1] * (1. - X[1]) + T_in); in f0_conduct_pipe_v()
578 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], const PetscReal n… in f0_conduct_bd_pipe_v() argument
580 …PetscReal sigma[4] = {X[0], (PetscReal)(0.5 * (1. - 2. * X[1])), (PetscReal)(0.5 * (1. - 2. * X[1]… in f0_conduct_bd_pipe_v()
588 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_pipe_q() argument
593 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_pipe_w() argument
644 static PetscErrorCode pipe_wiggly_u(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf,… in pipe_wiggly_u() argument
648 …0] = (0.5 * param->Reynolds / param->mu) * (X[1] * (1.0 - X[1]) + param->epsilon * PetscSinReal(PE… in pipe_wiggly_u()
652 static PetscErrorCode pipe_wiggly_u_t(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt N… in pipe_wiggly_u_t() argument
659 static PetscErrorCode pipe_wiggly_p(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf,… in pipe_wiggly_p() argument
661 p[0] = -X[0]; in pipe_wiggly_p()
664 static PetscErrorCode pipe_wiggly_p_t(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt N… in pipe_wiggly_p_t() argument
670 static PetscErrorCode pipe_wiggly_T(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt Nf,… in pipe_wiggly_T() argument
674 T[0] = X[1] * (1.0 - X[1]) + param->epsilon * PetscSinReal(PETSC_PI * X[1]) + param->T_in; in pipe_wiggly_T()
677 static PetscErrorCode pipe_wiggly_T_t(PetscInt dim, PetscReal time, const PetscReal X[], PetscInt N… in pipe_wiggly_T_t() argument
683 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_pipe_wiggly_v() argument
689 const PetscReal rho = p_th / (X[1] * (1. - X[1]) + T_in); in f0_conduct_pipe_wiggly_v()
692 f0[0] -= eps * 0.5 * PetscSqr(PETSC_PI) * PetscSinReal(PETSC_PI * X[1]); in f0_conduct_pipe_wiggly_v()
696 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], const PetscReal n… in f0_conduct_bd_pipe_wiggly_v() argument
699 …X[0], (PetscReal)(0.5 * (1. - 2. * X[1]) + eps * 0.5 * PETSC_PI * PetscCosReal(PETSC_PI * X[1])), … in f0_conduct_bd_pipe_wiggly_v()
707 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_pipe_wiggly_q() argument
712 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… 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 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_q() argument
730 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_q() argument
750 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_v() argument
764 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_v() argument
792 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f1_v() argument
805 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f1_conduct_v() argument
828 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_w() argument
839 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f0_conduct_w() argument
856 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f1_w() argument
865 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in f1_conduct_w() argument
1556 …st PetscScalar a_t[], const PetscScalar a_x[], PetscReal t, const PetscReal X[], PetscInt numConst… in divergence() argument