| /petsc/src/snes/tutorials/network/power/ |
| H A D | pffunctions.c | 147 …values[0] = Vmf * Vmt * (Gft * -PetscSinScalar(thetaft) + Bft * PetscCosScalar(thetaft)); … in FormJacobian_Power_private() 148 …[1] = 2.0 * Gff * Vmf + Vmt * (Gft * PetscCosScalar(thetaft) + Bft * PetscSinScalar(thetaft)); /* … in FormJacobian_Power_private() 149 …values[2] = Vmf * Vmt * (Gft * PetscSinScalar(thetaft) + Bft * -PetscCosScalar(thetaft)); … in FormJacobian_Power_private() 150 …values[3] = Vmf * (Gft * PetscCosScalar(thetaft) + Bft * PetscSinScalar(thetaft)); … in FormJacobian_Power_private() 161 … values[0] = Vmf * Vmt * (Bft * PetscSinScalar(thetaft) + Gft * PetscCosScalar(thetaft)); in FormJacobian_Power_private() 162 …ues[1] = -2.0 * Bff * Vmf + Vmt * (-Bft * PetscCosScalar(thetaft) + Gft * PetscSinScalar(thetaft)); in FormJacobian_Power_private() 163 … values[2] = Vmf * Vmt * (-Bft * PetscSinScalar(thetaft) + Gft * -PetscCosScalar(thetaft)); in FormJacobian_Power_private() 164 values[3] = Vmf * (-Bft * PetscCosScalar(thetaft) + Gft * PetscSinScalar(thetaft)); in FormJacobian_Power_private() 176 …values[0] = Vmt * Vmf * (Gtf * -PetscSinScalar(thetatf) + Btf * PetscCosScalar(thetaft)); … in FormJacobian_Power_private() 177 …[1] = 2.0 * Gtt * Vmt + Vmf * (Gtf * PetscCosScalar(thetatf) + Btf * PetscSinScalar(thetatf)); /* … in FormJacobian_Power_private() [all …]
|
| H A D | power2.c | 94 …] += Gff * Vmf * Vmf + Vmf * Vmt * (Gft * PetscCosScalar(thetaft) + Bft * PetscSinScalar(thetaft)); in FormFunction_Subnet() 95 …+= -Bff * Vmf * Vmf + Vmf * Vmt * (-Bft * PetscCosScalar(thetaft) + Gft * PetscSinScalar(thetaft)); in FormFunction_Subnet() 97 …] += Gtt * Vmt * Vmt + Vmt * Vmf * (Gtf * PetscCosScalar(thetatf) + Btf * PetscSinScalar(thetatf)); in FormFunction_Subnet() 98 …+= -Btt * Vmt * Vmt + Vmt * Vmf * (-Btf * PetscCosScalar(thetatf) + Gtf * PetscSinScalar(thetatf)); in FormFunction_Subnet() 274 …values[0] = Vmf * Vmt * (Gft * -PetscSinScalar(thetaft) + Bft * PetscCosScalar(thetaft)); … in FormJacobian_Subnet() 275 …[1] = 2.0 * Gff * Vmf + Vmt * (Gft * PetscCosScalar(thetaft) + Bft * PetscSinScalar(thetaft)); /* … in FormJacobian_Subnet() 276 …values[2] = Vmf * Vmt * (Gft * PetscSinScalar(thetaft) + Bft * -PetscCosScalar(thetaft)); … in FormJacobian_Subnet() 277 …values[3] = Vmf * (Gft * PetscCosScalar(thetaft) + Bft * PetscSinScalar(thetaft)); … in FormJacobian_Subnet() 288 … values[0] = Vmf * Vmt * (Bft * PetscSinScalar(thetaft) + Gft * PetscCosScalar(thetaft)); in FormJacobian_Subnet() 289 …ues[1] = -2.0 * Bff * Vmf + Vmt * (-Bft * PetscCosScalar(thetaft) + Gft * PetscSinScalar(thetaft)); in FormJacobian_Subnet() [all …]
|
| /petsc/src/ts/tutorials/multirate/ |
| H A D | ex2.c | 26 f[0] = -PetscSinScalar(ctx->a * t); in RHSFunction() 27 …f[1] = ctx->b * PetscCosScalar(ctx->b * t) * u[0] - PetscSinScalar(ctx->b * t) * PetscSinScalar(ct… in RHSFunction() 41 f[0] = -PetscSinScalar(ctx->a * t); in RHSFunctionslow() 55 …f[0] = ctx->b * PetscCosScalar(ctx->b * t) * u[0] - PetscSinScalar(ctx->b * t) * PetscSinScalar(ct… in RHSFunctionfast() 71 u[1] = PetscSinScalar(ctx->b * t) * PetscCosScalar(ctx->a * t) / ctx->a; in sol_true()
|
| H A D | ex3.c | 26 …05 * (-2.0 + u[1] * u[1] - PetscCosScalar(5.0 * t)) / (2.0 * u[1]) - PetscSinScalar(t) / (2.0 * u[… in RHSFunction() 27 …-2.0 + u[1] * u[1] - PetscCosScalar(5.0 * t)) / (2.0 * u[1]) - 5.0 * PetscSinScalar(5.0 * t) / (2.… in RHSFunction() 41 …05 * (-2.0 + u[1] * u[1] - PetscCosScalar(5.0 * t)) / (2.0 * u[1]) - PetscSinScalar(t) / (2.0 * u[… in RHSFunctionslow() 55 …-2.0 + u[1] * u[1] - PetscCosScalar(5.0 * t)) / (2.0 * u[1]) - 5.0 * PetscSinScalar(5.0 * t) / (2.… in RHSFunctionfast()
|
| H A D | ex3fastslowsplit.c | 31 …05 * (-2.0 + u[1] * u[1] - PetscCosScalar(5.0 * t)) / (2.0 * u[1]) - PetscSinScalar(t) / (2.0 * u[… in RHSFunctionslow() 45 …f[0] = 0.05 * (-1.0 + u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) - 5.0 * PetscSinScalar(5.0 *… in RHSFunctionfast()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex25.c | 114 v[0] = (-1.0 - user->e * PetscSinScalar(2.0 * PETSC_PI * user->k * xlow)) / h; in ComputeMatrix() 116 …v[1] = (2.0 + user->e * PetscSinScalar(2.0 * PETSC_PI * user->k * xlow) + user->e * PetscSinSc… in ComputeMatrix() 118 v[2] = (-1.0 - user->e * PetscSinScalar(2.0 * PETSC_PI * user->k * xhigh)) / h; in ComputeMatrix()
|
| H A D | ex13.c | 83 solution[Ii] = PetscSinScalar(2. * PETSC_PI * x) * PetscSinScalar(2. * PETSC_PI * y); in main() 84 … * PETSC_PI * x) * PetscSinScalar(2 * PETSC_PI * y) + 8 * PETSC_PI * PETSC_PI * x * PetscSinScalar… in main()
|
| H A D | ex74.c | 266 uarr[i - is] = PetscExpScalar(-4.0 * pi * pi * a * t) * PetscSinScalar(2 * pi * x); in ExactSolution() 269 uarr[i - is] = PetscSinScalar(2 * pi * (x - a * t)); in ExactSolution()
|
| H A D | ex28.c | 58 val[0] = val[1] = PetscSinScalar((PetscScalar)i * Hx); in ComputeInitialSolution()
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | ex3.c | 183 …for (i = mstart; i < mend; i++) u[i] = PetscSinScalar(PETSC_PI * i * 6. * h) + 3. * PetscSinScalar… in InitialConditions() 228 …start; i < mend; i++) u[i] = PetscSinScalar(sc1 * (PetscReal)i + appctx->a * PETSC_PI * 6. * t) * … in Solution()
|
| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busadj.c | 100 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri() 101 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri() 109 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq() 110 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq() 154 Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */ in SetInitialGuess() 155 Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */ in SetInitialGuess() 157 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 158 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 524 dVd_dVr = PetscSinScalar(delta); in ResidualJacobian() 527 dVq_dVi = PetscSinScalar(delta); in ResidualJacobian() [all …]
|
| H A D | ex9busoptfd.c | 102 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri() 103 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri() 111 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq() 112 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq() 156 Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */ in SetInitialGuess() 157 Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */ in SetInitialGuess() 159 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 160 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 525 dVd_dVr = PetscSinScalar(delta); in ResidualJacobian() 528 dVq_dVi = PetscSinScalar(delta); in ResidualJacobian() [all …]
|
| H A D | ex9busopt.c | 114 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri() 115 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri() 123 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq() 124 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq() 193 Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */ in SetInitialGuess() 194 Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */ in SetInitialGuess() 196 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 197 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 281 Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */ in InitialGuess() 282 Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */ in InitialGuess() [all …]
|
| H A D | ex9bus.c | 300 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri() 301 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri() 309 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq() 310 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq() 380 Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); /* d-axis stator current */ in SetInitialGuess() 381 Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); /* q-axis stator current */ in SetInitialGuess() 383 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 384 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 783 dVd_dVr = PetscSinScalar(delta); in ResidualJacobian() 786 dVq_dVi = PetscSinScalar(delta); in ResidualJacobian() [all …]
|
| H A D | ex9busdmnetwork.c | 417 Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); in SetInitialGuess() 420 Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); in SetInitialGuess() 422 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 423 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 469 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri() 470 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri() 478 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq() 479 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex4.c | 317 …base; i < myend; i++) u_localptr[i - mybase] = PetscSinScalar(PETSC_PI * i * 6. * h) + 3. * PetscS… in InitialConditions() 370 …se; i < myend; i++) s_localptr[i - mybase] = PetscSinScalar(sc1 * (PetscReal)i) * ex1 + 3. * Petsc… in ExactSolution()
|
| H A D | ex3.c | 308 …or (i = 0; i < appctx->m; i++) u_localptr[i] = PetscSinScalar(PETSC_PI * i * 6. * h) + 3. * PetscS… in InitialConditions() 355 … (i = 0; i < appctx->m; i++) s_localptr[i] = PetscSinScalar(sc1 * (PetscReal)i) * ex1 + 3. * Petsc… in ExactSolution()
|
| H A D | ex74.c | 131 uarr[i - is] = PetscExpScalar(-4.0 * pi * pi * a * t) * PetscSinScalar(2 * pi * x); in ExactSolution() 134 uarr[i - is] = PetscSinScalar(2 * pi * (x - a * t)); in ExactSolution()
|
| H A D | ex51.c | 28 f[1] = PetscSinScalar(u[1]); in RHSFunction()
|
| /petsc/src/ts/tutorials/power_grid/ |
| H A D | ex3.h | 65 …f[1] = ctx->omega_s / (2.0 * ctx->H) * (ctx->Pm - Pmax * PetscSinScalar(u[0]) - ctx->D * (u[1] - c… in RHSFunction() 115 …f[1] = 2.0 * ctx->H / ctx->omega_s * udot[1] + Pmax * PetscSinScalar(u[0]) + ctx->D * (u[1] - ctx-… in IFunction()
|
| H A D | ex1.c | 47 …f[0] = 2.0 * ctx->H * udot[0] / ctx->omega_s + ctx->E * ctx->V * PetscSinScalar(u[1]) / ctx->X - R; in IFunction()
|
| H A D | ex6.c | 206 f = (user->ws / (2 * user->H)) * (user->PM_min - user->Pmax * PetscSinScalar(x)); in adv2() 228 fthetac = user->ws / (2 * user->H) * (user->PM_min - user->Pmax * PetscSinScalar(theta)); in BoundaryConditions() 337 fthetac = user->ws / (2 * user->H) * (user->PM_min - user->Pmax * PetscSinScalar(xi)); in IJacobian() 435 …c3 = (user->ws / (2.0 * user->H)) * (user->PM_min - user->Pmax * PetscSinScalar(xi)) / (2 *… in IJacobian()
|
| H A D | ex8.c | 160 user->PM_min = user->Pmax * PetscSinScalar(user->mux); in ini_bou() 193 …f = (user->ws / (2 * user->H)) * (user->PM_min - user->Pmax * PetscSinScalar(x) - user->D * (y … in adv2() 304 …c3 = (user->ws / (2.0 * user->H)) * (user->PM_min - user->Pmax * PetscSinScalar(xi) - use… in IJacobian()
|
| /petsc/src/dm/tests/ |
| H A D | ex11.c | 53 …for (i = xs; i < xs + m; i++) aglobal[j][i][k] = PetscSinScalar(2.0 * PETSC_PI * (k + 1) * xy[j][i… in main()
|
| /petsc/src/snes/tests/ |
| H A D | ex4.c | 182 ff[0] = PetscSinScalar(3.0 * xx[0]) + xx[0]; in FormFunction2()
|