| /petsc/src/ts/tutorials/ |
| H A D | ex53.c | 225 u[0] = -sum * PetscSinReal(time); in linear_trig_p_t() 254 f0[0] += PetscSinReal(t) * sum / M; in f0_quadratic_trig_p() 308 …for (d = 0; d < dim; ++d) u[d] = PetscSinReal(2. * PETSC_PI * x[d]) - (d > 0 ? 2.0 * x[d - 1] * x[… in trig_u() 352 …r(2. * PETSC_PI) * PetscSinReal(2. * PETSC_PI * x[d]) * (2. * G + lambda) + 2.0 * (G + lambda) - 2… in f0_trig_linear_u() 353 … PetscSqr(2. * PETSC_PI) * PetscSinReal(2. * PETSC_PI * x[dim - 1]) * (2. * G + lambda) - 2. * PET… in f0_trig_linear_u() 499 …if (m % 2 == 1) F2_z += (-4.0 / (m * PETSC_PI * L)) * PetscSinReal(0.5 * m * PETSC_PI * zstar) * (… in terzaghi_2d_eps() 537 …if (m % 2 == 1) F1 += (4.0 / (m * PETSC_PI)) * PetscSinReal(0.5 * m * PETSC_PI * zstar) * PetscExp… in terzaghi_2d_p() 611 …if (m % 2 == 1) F2_zt += ((-m * PETSC_PI * c) / (L * L * L)) * PetscSinReal(0.5 * m * PETSC_PI * z… in terzaghi_2d_eps_t() 661 …if (m % 2 == 1) F1_t += ((-m * PETSC_PI * c) / PetscSqr(L)) * PetscSinReal(0.5 * m * PETSC_PI * zs… in terzaghi_2d_p_t() 698 …p += (PetscSinReal(aa) / (aa - PetscSinReal(aa) * PetscCosReal(aa))) * (PetscCosReal((aa * x[0]) /… in mandel_drainage_pressure() [all …]
|
| H A D | ex46.c | 62 u[0] = PetscSinReal(time + x[0]) * PetscSinReal(time + x[1]); in mms2_u_2d() 69 *p = PetscSinReal(time + x[0] - x[1]); in mms2_p_2d() 101 …) * PetscSinReal(2 * t + 2 * x[0]) + PetscSinReal(2 * t + x[0] + x[1]) + PetscCosReal(t + x[0] - x… in f0_mms2_u() 102 …f0[1] -= (-Re * ((1.0L / 2.0L) * PetscSinReal(2 * t + 2 * x[1]) + PetscSinReal(2 * t + x[0] + x[1]… in f0_mms2_u()
|
| H A D | ex6.c | 274 … (i = 0; i < appctx->m; i++) u_localptr[i] = PetscSinReal(PETSC_PI * i * 6. * appctx->h) + 3. * Pe… in InitialConditions() 318 … 0; i < appctx->m; i++) s_localptr[i] = PetscSinReal(PetscRealPart(sc1) * (PetscReal)i) * ex1 + 3.… in ExactSolution()
|
| H A D | ex36.c | 40 *U = 0.4 * PetscSinReal(200 * PETSC_PI * t); in Ue() 58 f[0] = ydot[0] / 1.e6 - ydot[1] / 1.e6 - PetscSinReal(200 * PETSC_PI * t) / 2500. + y[0] / 1000.; in IFunctionImplicit()
|
| H A D | ex43.c | 20 u = PetscExpReal(-a * t) * (C1 * PetscSinReal(w * t) + C2 * PetscCosReal(w * t)); in Exact() 21 …a * t) * (C1 * PetscSinReal(w * t) + C2 * PetscCosReal(w * t)) + w * PetscExpReal(-a * t) * (C1 * … in Exact()
|
| H A D | ex76.c | 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() 352 u[0] = -100. * PetscSinReal(time); in cubic_trig_u_t() 370 T[0] = -100. * PetscSinReal(time); in cubic_trig_T_t() 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() 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() 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() 450 …PetscSinReal(PETSC_PI * (X[0] - time)) * PetscSinReal(PETSC_PI * (X[1] - time)) - PetscCosReal(PET… in taylor_green_u_t() 451 …PetscSinReal(PETSC_PI * (X[0] - time)) * PetscSinReal(PETSC_PI * (X[1] - time)) - PetscCosReal(PET… in taylor_green_u_t() [all …]
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | ex5opt_ic.c | 101 …(y <= 1.5)) u[j][i].v = .25 * PetscPowReal(PetscSinReal(4.0 * PETSC_PI * x), 2.0) * PetscPowReal(P… in InitialConditions() 133 …<= 1.5)) u[j][i].v = .25 * PetscPowReal(PetscSinReal(4.0 * PETSC_PI * 0.5 * x), 2.0) * PetscPowRea… in PerturbedInitialConditions() 166 …u[j][i].v = PetscPowReal(PetscSinReal(4.0 * PETSC_PI * (x - 0.5)), 2.0) * PetscPowReal(PetscSinRea… in PerturbedInitialConditions2() 198 …(y <= 2.0)) u[j][i].v = .25 * PetscPowReal(PetscSinReal(4.0 * PETSC_PI * x), 2.0) * PetscPowReal(P… in PerturbedInitialConditions3()
|
| H A D | ex6.c | 141 …for (i = mstart; i < mend; i++) u[i] = PetscSinReal(PETSC_PI * i * 6. * h) + 3. * PetscSinReal(PET… in InitialConditions() 180 …for (i = mstart; i < mend; i++) u[i] = PetscSinReal(PI6 * (i * h - a * t)) + 3. * PetscSinReal(PI2… in Solution()
|
| H A D | ex5.c | 72 …u[j][i].v = PetscPowReal(PetscSinReal(4.0 * PETSC_PI * x), 2.0) * PetscPowReal(PetscSinReal(4.0 * … in InitialConditions()
|
| H A D | ex5adj.c | 183 …u[j][i].v = PetscPowReal(PetscSinReal(4.0 * PETSC_PI * x), 2.0) * PetscPowReal(PetscSinReal(4.0 * … in InitialConditions()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex5.c | 193 u[0] = PetscSinReal(PETSC_PI * x) * PetscSinReal(PETSC_PI * y); in MMSSolution2() 202 …(PETSC_PI) * PetscSinReal(PETSC_PI * x) * PetscSinReal(PETSC_PI * y) - user->param * PetscExpReal(… in MMSForcing2() 211 …u[0] = PetscSinReal(user->m * PETSC_PI * x * (1 - y)) * PetscSinReal(user->n * PETSC_PI * y * (1 -… in MMSSolution3() 221 …PetscSinReal(m * PETSC_PI * x * (1 - y)) * PetscSinReal(n * PETSC_PI * (1 - x) * y)) * lambda) + P… in MMSForcing3()
|
| H A D | ex55.c | 146 u[0] = PetscSinReal(PETSC_PI * x) * PetscSinReal(PETSC_PI * y); in MMSSolution2() 153 …(PETSC_PI) * PetscSinReal(PETSC_PI * x) * PetscSinReal(PETSC_PI * y) - user->param * PetscExpReal(… in MMSForcing2() 160 …u[0] = PetscSinReal(user->m * PETSC_PI * x * (1 - y)) * PetscSinReal(user->n * PETSC_PI * y * (1 -… in MMSSolution3() 168 …PetscSinReal(m * PETSC_PI * x * (1 - y)) * PetscSinReal(n * PETSC_PI * (1 - x) * y)) * lambda) + P… in MMSForcing3()
|
| H A D | ex26.c | 42 for (d = 0; d < dim; ++d) *u += PetscSinReal(2.0 * PETSC_PI * x[d]); in trig_u() 58 … d < dim; ++d) f0[0] -= 4.0 * PetscSqr(PETSC_PI) * PetscSinReal(2.0 * PETSC_PI * x[d]) + PetscSinR… in f0_trig_u()
|
| H A D | ex62.c | 200 u[0] = (dim - 1) * PetscSinReal(PETSC_PI * x[0]); in trig_u() 202 u[0] += PetscSinReal(PETSC_PI * x[c]); in trig_u() 212 for (d = 0, u[0] = 0.0; d < dim; ++d) u[0] += PetscSinReal(2.0 * PETSC_PI * x[d]); in trig_p() 221 …osReal(2.0 * PETSC_PI * x[0]) - (dim - 1) * mu * PetscSqr(PETSC_PI) * PetscSinReal(PETSC_PI * x[0]… in f0_trig_u() 223 f0[0] -= mu * PetscSqr(PETSC_PI) * PetscSinReal(PETSC_PI * x[d]); in f0_trig_u()
|
| H A D | ex17.c | 149 u[0] = PetscSinReal(2.0 * PETSC_PI * x[0]); in trig_2d_u() 150 u[1] = PetscSinReal(2.0 * PETSC_PI * x[1]) - 2.0 * x[0] * x[1]; in trig_2d_u() 156 u[0] = PetscSinReal(2.0 * PETSC_PI * x[0]); in trig_3d_u() 157 u[1] = PetscSinReal(2.0 * PETSC_PI * x[1]) - 2.0 * x[0] * x[1]; in trig_3d_u() 158 u[2] = PetscSinReal(2.0 * PETSC_PI * x[2]) - 2.0 * x[1] * x[2]; in trig_3d_u() 175 …for (d = 0; d < dim; ++d) f0[d] += -4.0 * PetscSqr(PETSC_PI) * PetscSinReal(2.0 * PETSC_PI * x[d]); in f0_vlap_trig_u() 205 …for (PetscInt d = 0; d < dim; ++d) f0[d] += -(2.0 * mu + lambda) * fact * PetscSinReal(2.0 * PETSC… in f0_elas_trig_u()
|
| H A D | ex20.c | 38 for (d = 0; d < dim; ++d) u[0] += PetscSinReal(2.0 * PETSC_PI * x[d]); in trig_u() 45 …for (d = 0; d < dim; ++d) f0[0] += -4.0 * PetscSqr(PETSC_PI) * PetscSinReal(2.0 * PETSC_PI * x[d]); in f0_trig_u()
|
| H A D | ex27.c | 143 for (PetscInt d = 0; d < dim; ++d) u[0] *= PetscSinReal(n * PETSC_PI * x[d]); in trig_u() 150 … < Nc; c++) u[c] = n * PETSC_PI * PetscCosReal(n * PETSC_PI * x[c]) * PetscSinReal(n * PETSC_PI * … in trig_q() 168 u[0] = PetscSinReal(50. * xy); in sensor_u() 200 …f0[0] = -2.0 * PetscSqr(n * PETSC_PI) * PetscSinReal(n * PETSC_PI * x[0]) * PetscSinReal(n * PETSC… in f0_trig_primal() 209 f0[0] = -2500.0 * PetscSinReal(50. * xy) * (xref * xref + yref * yref) * 4.0; in f0_sensor_primal() 264 …f0[0] = -2.0 * PetscSqr(n * PETSC_PI) * PetscSinReal(n * PETSC_PI * x[0]) * PetscSinReal(n * PETSC… in f0_mixed_trig_u() 273 f0[0] = -2500.0 * PetscSinReal(50. * xy) * (xref * xref + yref * yref) * 4.0; in f0_mixed_sensor_u()
|
| H A D | ex8.c | 58 *u = (y - 0.5) * PetscSinReal(pi2 * x) + (x - 0.5) * PetscSinReal(pi2 * y); in sins_2d()
|
| /petsc/src/dm/tests/ |
| H A D | ex12.c | 46 …globalptr[i] = PetscSinReal((PETSC_PI * j * 6) / ((PetscReal)M) + 1.2 * PetscSinReal((PETSC_PI * j… in main()
|
| H A D | ex3.c | 66 …globalptr[i - 1] = PetscSinReal((PETSC_PI * j * 6) / ((PetscReal)M) + 1.2 * PetscSinReal((PETSC_PI… in main()
|
| /petsc/src/ts/event/tests/ |
| H A D | ex4.c | 198 …if (Ctx->F == 0 || Ctx->F == -1) gval[n++] = PetscSinReal(Ctx->pi * t) / Ctx->pi; // FUNC-0, roots… in EventFunction() 204 …if (Ctx->F == 3 || Ctx->F == -1) gval[n++] = t - PetscExpReal(PetscSinReal(t)) + 1; … in EventFunction() 210 …tx->F == -1) gval[n++] = (t * t * (t * t / 3.0 + PetscSqrtReal(2.0) * PetscSinReal(t)) - PetscSqrt… in EventFunction() 211 …if (Ctx->F == 10 || Ctx->F == -1) gval[n++] = ((t * t + 1) * PetscSinReal(t) - PetscExpReal(PetscS… in EventFunction()
|
| /petsc/src/ts/tutorials/optimal_control/ |
| H A D | ex1.c | 38 f[1] = v[step] * PetscSinReal(w[step]); in RHSFunction() 61 Jp[0][1] = -v[step] * PetscSinReal(w[step]); in RHSJacobianP() 62 Jp[1][0] = PetscSinReal(w[step]); in RHSJacobianP() 114 dJpdP[0][0][1] = -PetscSinReal(w[step]); in RHSHessianProductPP() 115 dJpdP[0][1][0] = -PetscSinReal(w[step]); in RHSHessianProductPP() 119 dJpdP[1][1][1] = -v[step] * PetscSinReal(w[step]); in RHSHessianProductPP() 192 dy = u[1] - actx->lv * t * PetscSinReal(actx->lw); in CostIntegrand() 210 dy = u[1] - actx->lv * t * PetscSinReal(actx->lw); in DRDUJacobianTranspose()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex42.c | 1170 v[0] = PetscPowRealInt(z, 3) * PetscExpReal(y) * PetscSinReal(PETSC_PI * x); in evaluate_MS_FrankKamentski() 1171 v[1] = z * PetscCosReal(2.0 * PETSC_PI * x) * PetscExpReal(-y) * PetscSinReal(PETSC_PI * y); in evaluate_MS_FrankKamentski() 1172 …ealInt(z, 2) * (PetscCosReal(2.0 * PETSC_PI * x) * PetscExpReal(-y) * PetscSinReal(PETSC_PI * y) /… in evaluate_MS_FrankKamentski() 1195 …PetscSinReal(2.0 * PETSC_PI * x) - 1.0 * PETSC_PI * PetscExpReal(-y) * PetscSinReal(PETSC_PI * y) … in evaluate_MS_FrankKamentski() 1196 …PetscSinReal(PETSC_PI * y) / 2 + PetscPowRealInt(PETSC_PI, 2) * PetscCosReal(2.0 * PETSC_PI * x) *… in evaluate_MS_FrankKamentski() 1197 …PetscSinReal(PETSC_PI * y) + 2.0 * PetscPowRealInt(PETSC_PI, 3) * PetscCosReal(PETSC_PI * y) * Pet… in evaluate_MS_FrankKamentski() 1876 … cell->fy[p] = -PetscSinReal(2.2 * PETSC_PI * coord_y) * PetscCosReal(1.0 * PETSC_PI * coord_x); in solve_stokes_3d_coupled() 1930 … cell->fy[p] = -PetscSinReal(3.0 * PETSC_PI * coord_y) * PetscCosReal(1.0 * PETSC_PI * coord_x); in solve_stokes_3d_coupled() 1995 cell->fy[p] = -PetscSinReal(2.2 * PETSC_PI * yp) * PetscCosReal(1.0 * PETSC_PI * xp); in solve_stokes_3d_coupled()
|
| /petsc/src/ksp/pc/impls/mg/ |
| H A D | mgadapt.c | 29 for (c = 0; c < Nc; ++c) u[c] = PetscSinReal(PETSC_PI * (k + 1) * coords[0]); in xsin() 36 for (c = 0; c < Nc; ++c) u[c] = PetscSinReal(PETSC_PI * (k + 1) * coords[1]); in ysin() 43 for (c = 0; c < Nc; ++c) u[c] = PetscSinReal(PETSC_PI * (k + 1) * coords[2]); in zsin()
|
| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | ex19.c | 280 …ordData[faceCPStartRow + (3 * ii) + 0] += -0.05 * deltaR * PetscCosReal(phi) * PetscSinReal(theta); in main() 281 …oordData[faceCPStartRow + (3 * ii) + 1] += -0.05 * deltaR * PetscSinReal(phi) * PetscSinReal(theta… in main()
|