Home
last modified time | relevance | path

Searched refs:omega (Results 1 – 25 of 97) sorted by relevance

1234

/petsc/src/ksp/pc/impls/sor/
H A Dsor.tex18 …1/2} = x_i^n + \omega A_{ii}^{-1}( b_i - \sum_{j < i} A_{ij} x_j^{n+1/2} - \sum_{j \ge i} A_{ij} …
22 …^{n+1/2} + \omega A_{ii}^{-1}( b_i - \sum_{j \le i} A_{ij} x_j^{n+1/2} - \sum_{j > i} A_{ij} x_j^…
24 … over-relaxation because generally $ \omega $ is greater than one, though on occasion underrelaxat…
28 x_i^{1/2} = \omega A_{ii}^{-1}( b_i - \sum_{j < i} A_{ij} x_j^{1/2})
32 x_i = (1 - \omega) x_i^{1/2} + \omega A_{ii}^{-1}( b_i - \sum_{j < i} A_{ij} x_j^{1/2} - \sum_{j > …
38 x^{1/2} = \omega (L + D)^{-1} b
41 x = (1 - \omega) x^{1/2} + \omega (U + D)^{-1}(b - L x^{1/2}) = x^{1/2} + \omega (U+D)^{-1}(b - A x…
50 tmp = - (x[i] = omega*t[i]*aidiag[i]);
71 x[i] = (1-omega)*x[i] + omega*sum*aidiag[i];
86 x = \omega (L + D)^{-1}b
[all …]
H A Dsor.c10 PetscReal omega; member
33 PetscCall(MatSOR(pc->pmat, x, jac->omega, (MatSORType)flag, jac->fshift, jac->its, jac->lits, y)); in PCApply_SOR()
47 PetscCall(MatSOR(pc->pmat, x, jac->omega, (MatSORType)flag, jac->fshift, jac->its, jac->lits, y)); in PCApplyTranspose_SOR()
59 PetscCall(MatSOR(pc->pmat, b, jac->omega, stype, jac->fshift, its * jac->its, jac->lits, y)); in PCApplyRichardson_SOR()
70 PetscReal omega; in PCSetFromOptions_SOR() local
74 …("-pc_sor_omega", "relaxation factor (0 < omega < 2)", "PCSORSetOmega", jac->omega, &omega, &flg)); in PCSetFromOptions_SOR()
75 if (flg) PetscCall(PCSORSetOmega(pc, omega)); in PCSetFromOptions_SOR()
116 … iterations = %" PetscInt_FMT ", omega = %g\n", sortype, jac->its, jac->lits, (double)jac->omega)); in PCView_SOR()
130 static PetscErrorCode PCSORSetOmega_SOR(PC pc, PetscReal omega) in PCSORSetOmega_SOR() argument
135 …PetscCheck(omega > 0.0 && omega < 2.0, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_OUTOFRANGE,… in PCSORSetOmega_SOR()
[all …]
/petsc/src/ts/tutorials/hamiltonian/
H A Dex1.c31 PetscReal omega; member
49 vres[0] = -user->omega * user->omega * x[0]; in RHSFunction2()
79 r[1] = -user->omega * user->omega * u[0]; in RHSFunction()
97 …menergy = (u[1] * u[1] + user->omega * user->omega * u[0] * u[0] - user->omega * user->omega * dt … in Monitor()
98 energy = (u[1] * u[1] + user->omega * user->omega * u[0] * u[0]) / 2.; in Monitor()
128 user.omega = 64.; in main()
132 PetscCall(PetscOptionsReal("-omega", "parameter", "<64>", user.omega, &user.omega, NULL)); in main()
184 …e)ftime, (double)(0.2 * PetscCosReal(user.omega * ftime)), (double)(-0.2 * user.omega * PetscSinRe… in main()
/petsc/src/snes/tutorials/
H A Dex19.c77 PetscScalar u, v, omega, temp; member
221 x[j][i].omega = 0.0; in FormInitialGuess()
273 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormFunctionLocal()
286 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormFunctionLocal()
299 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormFunctionLocal()
312 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormFunctionLocal()
336 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal()
342 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal()
345 u = x[j][i].omega; in FormFunctionLocal()
346 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormFunctionLocal()
[all …]
/petsc/src/ksp/pc/impls/eisens/
H A Deisen.c12 PetscReal omega; member
24 PetscCall(MatSOR(eis->A, b, eis->omega, SOR_EISENSTAT, 0.0, 1, 1, x)); in PCMult_Eisenstat()
96 …PetscCall(MatSOR(eis->A, eis->b[pc->presolvedone - 1], eis->omega, SOR_APPLY_UPPER, 0.0, 1, 1, x)); in PCPreSolve_Eisenstat()
104 …PetscCall(MatSOR(eis->A, eis->b[pc->presolvedone - 1], eis->omega, (MatSORType)(SOR_ZERO_INITIAL_G… in PCPreSolve_Eisenstat()
119 …PetscCall(MatSOR(eis->A, eis->b[pc->presolvedone], eis->omega, (MatSORType)(SOR_ZERO_INITIAL_GUESS… in PCPostSolve_Eisenstat()
154 PetscReal omega; in PCSetFromOptions_Eisenstat() local
158 …nstat_omega", "Relaxation factor 0 < omega < 2", "PCEisenstatSetOmega", eis->omega, &omega, &flg)); in PCSetFromOptions_Eisenstat()
159 if (flg) PetscCall(PCEisenstatSetOmega(pc, omega)); in PCSetFromOptions_Eisenstat()
174 PetscCall(PetscViewerASCIIPrintf(viewer, " omega = %g\n", (double)eis->omega)); in PCView_Eisenstat()
210 static PetscErrorCode PCEisenstatSetOmega_Eisenstat(PC pc, PetscReal omega) in PCEisenstatSetOmega_Eisenstat() argument
[all …]
/petsc/src/ts/tutorials/
H A Dex26.c98 PetscScalar u, v, omega, temp; member
245 x[j][i].omega = 0.0; in FormInitialSolution()
297 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormIFunctionLocal()
310 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormIFunctionLocal()
323 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormIFunctionLocal()
336 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormIFunctionLocal()
361 f[j][i].u = udot + uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormIFunctionLocal()
368 f[j][i].v = udot + uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormIFunctionLocal()
371 u = x[j][i].omega; in FormIFunctionLocal()
372 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormIFunctionLocal()
[all …]
H A Dex43.c12 static void Exact(PetscReal t, PetscReal omega, PetscReal xi, PetscReal u0, PetscReal v0, PetscReal… in Exact() argument
16 PetscReal a = xi * omega; in Exact()
17 PetscReal w = PetscSqrtReal(1 - xi * xi) * omega; in Exact()
23 PetscReal w = PetscSqrtReal(xi * xi - 1) * omega; in Exact()
29 PetscReal a = xi * omega; in Exact()
H A Dex77.c51 PetscReal omega; /* Rotation speed in MMS benchmark */ member
124 x[0] = R0 * PetscCosReal(p->omega * time + theta0); in trig_trig_x()
125 x[1] = R0 * PetscSinReal(p->omega * time + theta0); in trig_trig_x()
132 u[0] = -p->omega * X[1]; in trig_trig_u()
133 u[1] = p->omega * X[0]; in trig_trig_u()
162 const PetscReal omega = PetscRealPart(constants[3]); in f0_trig_trig_v() local
171 f0[0] -= 1.0 - omega * omega * X[0]; in f0_trig_trig_v()
172 f0[1] -= 1.0 - omega * omega * X[1]; in f0_trig_trig_v()
177 const PetscReal omega = PetscRealPart(constants[3]); in f0_trig_trig_w() local
181 f0[0] += u_t[uOff[2]] - (1.0 + omega * (X[0] - X[1])); in f0_trig_trig_w()
[all …]
/petsc/src/dm/impls/swarm/tests/
H A Dex4.c59 PetscReal omega; /* Oscillation frequency omega */ member
67 options->omega = 64.0; in ProcessOptions()
72 …etscOptionsReal("-omega", "Oscillator frequency", "ex4.c", options->omega, &options->omega, NULL)); in ProcessOptions()
131 const PetscReal omega = ((AppCtx *)ctx)->omega; in RHSFunction() local
147 g[(p * 2 + 1) * dim + d] = -PetscSqr(omega) * u[(p * 2 + 0) * dim + d]; in RHSFunction()
161 PetscScalar vals[4] = {0., 1., -PetscSqr(((AppCtx *)ctx)->omega), 0.}; in RHSJacobian()
200 const PetscReal omega = ((AppCtx *)ctx)->omega; in RHSFunctionV() local
209 for (p = 0; p < Np; ++p) vres[p] = -PetscSqr(omega) * x[p]; in RHSFunctionV()
240 const PetscReal omega = ((AppCtx *)ctx)->omega; in RHSObjectiveF() local
254 const PetscReal E = 0.5 * (v2 + PetscSqr(omega) * x2); in RHSObjectiveF()
[all …]
H A Dex5.c113 const PetscReal omega = vy0 / x0; in RHSFunction() local
117 g[(p * 2 + 1) * dim + d] = -PetscSqr(omega) * u[(p * 2 + 0) * dim + d]; in RHSFunction()
148 const PetscReal omega = vy0 / x0; in RHSJacobian() local
149 PetscScalar vals[4] = {0., 1., -PetscSqr(omega), 0.}; in RHSJacobian()
205 const PetscReal omega = vy0 / x0; in RHSFunctionV() local
207 for (d = 0; d < dim; ++d) vres[p * dim + d] = -PetscSqr(omega) * x[p * dim + d]; in RHSFunctionV()
357 const PetscReal omega = PetscSqrtReal(1000. / r0) / r0; in circleMultipleV() local
359 v[0] = -r0 * omega * PetscSinReal(th0); in circleMultipleV()
360 v[1] = r0 * omega * PetscCosReal(th0); in circleMultipleV()
450 const PetscReal omega = v0 / r0; in ComputeError() local
[all …]
H A Dex6.c546 const PetscReal omega = vy0 / x0; in RHSFunction() local
550 g[(p * 2 + 1) * dim + d] = E[p * dim + d] - PetscSqr(omega) * u[(p * 2 + 0) * dim + d]; in RHSFunction()
584 const PetscReal omega = vy0 / x0; in RHSJacobian() local
585 PetscScalar vals[4] = {0., 1., -PetscSqr(omega), 0.}; in RHSJacobian()
648 const PetscReal omega = vy0 / x0; in RHSFunctionV() local
650 for (d = 0; d < dim; ++d) vres[p * dim + d] = E[p * dim + d] - PetscSqr(omega) * x[p * dim + d]; in RHSFunctionV()
721 const PetscReal omega = vy0 / x0; in RHSObjectiveF() local
724 E += 0.5 * q_p * m_p * (v2) + 0.5 * PetscSqr(omega) * (x2); in RHSObjectiveF()
769 const PetscReal omega = vy0 / x0; in RHSFunctionG() local
771 …g[(p * 2 + 0) * dim + d] = -(q_p / m_p) * E[p * dim + d] + PetscSqr(omega) * u[(p * 2 + 0) * dim +… in RHSFunctionG()
[all …]
/petsc/src/snes/tests/
H A Dex69.c15 PetscScalar u, v, omega, temp; member
196 x[j][i].omega = 0.0; in FormInitialGuess()
254 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormFunctionLocal()
267 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormFunctionLocal()
280 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormFunctionLocal()
293 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormFunctionLocal()
317 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal()
323 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal()
326 u = x[j][i].omega; in FormFunctionLocal()
327 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormFunctionLocal()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/pipebcgs/
H A Dpipebcgs.c18 PetscScalar rho, rhoold, alpha, beta, omega = 0.0, d1, d2, d3; in KSPSolve_PIPEBCGS() local
99 …PetscCall(VecAXPBYPCZ(P2, 1.0, -beta * omega, beta, R2, S2)); /* p2 <- beta * p2 + r2 - beta * ome… in KSPSolve_PIPEBCGS()
100 …PetscCall(VecAXPBYPCZ(S, 1.0, -beta * omega, beta, W, Z)); /* s <- beta * s + w - beta * ome… in KSPSolve_PIPEBCGS()
101 …PetscCall(VecAXPBYPCZ(S2, 1.0, -beta * omega, beta, W2, Z2)); /* s2 <- beta * s2 + w2 - beta * ome… in KSPSolve_PIPEBCGS()
102 …PetscCall(VecAXPBYPCZ(Z, 1.0, -beta * omega, beta, T, V)); /* z <- beta * z + t - beta * ome… in KSPSolve_PIPEBCGS()
135 omega = d1 / d2; /* omega <- (y'q) / (y'y) */ in KSPSolve_PIPEBCGS()
136 PetscCall(VecAXPBYPCZ(X, alpha, omega, 1.0, P2, Q2)); /* x <- alpha * p2 + omega * q2 + x */ in KSPSolve_PIPEBCGS()
137 PetscCall(VecWAXPY(R, -omega, Y, Q)); /* r <- q - omega y */ in KSPSolve_PIPEBCGS()
139 PetscCall(VecAYPX(R2, -omega, Q2)); /* r2 <- q2 - omega r2 */ in KSPSolve_PIPEBCGS()
141 PetscCall(VecAYPX(W, -omega, Y)); /* w <- y - omega w */ in KSPSolve_PIPEBCGS()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Drelax.h122 PetscErrorCode MatSOR_SeqSBAIJ_ushort(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fs… in MatSOR_SeqSBAIJ_ushort() argument
124 PetscErrorCode MatSOR_SeqSBAIJ(Mat A, Vec bb, PetscReal omega, MatSORType flag, PetscReal fshift, P… in MatSOR_SeqSBAIJ_ushort()
172 sum = b[i] * d / omega; in MatSOR_SeqSBAIJ_ushort()
191 tmp = -(x[i] = omega * t[i] * aidiag[i]); in MatSOR_SeqSBAIJ_ushort()
211 x[i] = omega * sum * aidiag[i]; in MatSOR_SeqSBAIJ_ushort()
226 x[i] = (1 - omega) * x[i] + omega * sum * aidiag[i]; in MatSOR_SeqSBAIJ_ushort()
258 x[i] = (1 - omega) * x[i] + omega * sum * aidiag[i]; in MatSOR_SeqSBAIJ_ushort()
288 x[i] = (1 - omega) * x[i] + omega * sum * aidiag[i]; in MatSOR_SeqSBAIJ_ushort()
/petsc/src/ksp/ksp/impls/bcgs/qmrcgs/
H A Dqmrcgs.c18 PetscScalar eta, rho1, rho2, alpha, eta2, omega, beta, cf, cf1, uu; in KSPSolve_QMRCGS() local
156 omega = uu / vv; /* omega <- uu/vv; */ in KSPSolve_QMRCGS()
159 PetscCall(VecWAXPY(R, -omega, T, S)); /* r <- s - omega t */ in KSPSolve_QMRCGS()
177 eta = c * c * omega; in KSPSolve_QMRCGS()
179 cf1 = eta2 * theta2 * theta2 / omega; in KSPSolve_QMRCGS()
191 beta = (alpha * rho2) / (omega * rho1); in KSPSolve_QMRCGS()
192 …PetscCall(VecAXPBYPCZ(P, 1.0, -omega * beta, beta, R, V)); /* p <- r - omega * beta* v + beta * p … in KSPSolve_QMRCGS()
/petsc/src/mat/tests/
H A Dex3.c11 PetscReal omega = 1.0, norm; in main() local
15 PetscCall(PetscOptionsGetReal(NULL, NULL, "-omega", &omega, NULL)); in main()
57 PetscCall(MatSOR(C, b, omega, SOR_FORWARD_SWEEP, 0.0, 1, 1, x)); in main()
/petsc/src/ksp/ksp/impls/bcgs/fbcgs/
H A Dfbcgs.c18 PetscScalar rho, rhoold, alpha, beta, omega, omegaold, d1; in KSPSolve_FBCGS() local
117 omega = d1 / d2; /* omega <- (t's) / (t't) */ in KSPSolve_FBCGS()
118 PetscCall(VecAXPBYPCZ(X, alpha, omega, 1.0, P2, S2)); /* x <- alpha * p2 + omega * s2 + x */ in KSPSolve_FBCGS()
120 PetscCall(VecWAXPY(R, -omega, T, S)); /* r <- s - omega t */ in KSPSolve_FBCGS()
124 omegaold = omega; in KSPSolve_FBCGS()
/petsc/src/ksp/ksp/impls/bcgs/fbcgsr/
H A Dfbcgsr.c20 PetscScalar tau, sigma, alpha, omega, beta; in KSPSolve_FBCGSR() local
170 omega = xi2 / xi3; in KSPSolve_FBCGSR()
172 rho = PetscSqrtReal(PetscAbsScalar(xi1 - omega * xi2)); /* residual norm */ in KSPSolve_FBCGSR()
175 PetscCall(VecAXPBYPCZ(X, alpha, omega, 1.0, P2, S2)); /* x <- alpha * p2 + omega * s2 + x */ in KSPSolve_FBCGSR()
191 r[j] = s[j] - omega * t[j]; /* r <- s - omega t */ in KSPSolve_FBCGSR()
192 p[j] = r[j] + beta * (p[j] - omega * v[j]); /* p <- r + beta * (p - omega v) */ in KSPSolve_FBCGSR()
/petsc/src/ksp/ksp/impls/bcgs/
H A Dbcgs.c21 PetscScalar rho, rhoold, alpha, beta, omega, omegaold, d1; in KSPSolve_BCGS() local
110 omega = d1 / d2; /* w <- (t's) / (t't) */ in KSPSolve_BCGS()
111 PetscCall(VecAXPBYPCZ(X, alpha, omega, 1.0, P, S)); /* x <- alpha * p + omega * s + x */ in KSPSolve_BCGS()
112 PetscCall(VecWAXPY(R, -omega, T, S)); /* r <- s - w t */ in KSPSolve_BCGS()
119 omegaold = omega; in KSPSolve_BCGS()
/petsc/src/ts/tutorials/power_grid/output/
H A Dex3sa_2.out2 trajectory sensitivity: d[phi(tf)]/d[pm] d[omega(tf)]/d[pm]
/petsc/src/ksp/ksp/tutorials/output/
H A Dex1_2_aijcusparse.out20 type = symmetric, iterations = 1, local iterations = 1, omega = 1.
44 type = symmetric, iterations = 1, local iterations = 1, omega = 1.
H A Dex1_3_aijcusparse.out20 omega = 1.
45 omega = 1.
H A Dex29_4.out64 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
94 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
124 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
154 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.
/petsc/src/ksp/ksp/impls/bcgsl/
H A Dbcgsl.c14 PetscScalar alpha, beta, omega, sigma; in KSPSolve_BCGSL() local
63 rho0 = omega = 1; in KSPSolve_BCGSL()
97 rho0 = -omega * rho0; in KSPSolve_BCGSL()
252 omega = AY0c[bcgsl->ell]; in KSPSolve_BCGSL()
253 for (h = bcgsl->ell; h > 0 && omega == 0.0; h--) omega = AY0c[h]; in KSPSolve_BCGSL()
254 if (omega == 0.0) { in KSPSolve_BCGSL()
/petsc/src/ksp/ksp/tests/output/
H A Dex34_1.out11 type = local_symmetric, iterations = 1, local iterations = 1, omega = 1.1

1234