| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9bus.c | 142 PetscScalar Efd, RF, VR, Vr, Vi, Vm; in EventFunction() local
163 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in EventFunction()
168 … = PetscRealPart((VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]); in EventFunction()
173 … = PetscRealPart((VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]); in EventFunction()
194 PetscScalar Vr, Vi, Vm; in PostEventFunction() local
247 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in PostEventFunction()
250 fvalue = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]; in PostEventFunction()
257 fvalue = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]; in PostEventFunction()
345 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
371 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess()
[all …]
|
| H A D | ex9busadj.c | 119 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
145 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess()
146 Vm2 = Vm * Vm; in SetInitialGuess()
188 Vref[i] = Vm + (VR / KA[i]); in SetInitialGuess()
208 PetscScalar Vr, Vi, Vm, Vm2; in ResidualFunction() local
280 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in ResidualFunction()
287 … fgen[idx + 8] = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]; in ResidualFunction()
296 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in ResidualFunction()
297 Vm2 = Vm * Vm; in ResidualFunction()
300 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap… in ResidualFunction()
[all …]
|
| H A D | ex9busoptfd.c | 121 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
147 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess()
148 Vm2 = Vm * Vm; in SetInitialGuess()
190 Vref[i] = Vm + (VR / KA[i]); in SetInitialGuess()
210 PetscScalar Vr, Vi, Vm, Vm2; in ResidualFunction() local
281 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in ResidualFunction()
288 … fgen[idx + 8] = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]; in ResidualFunction()
297 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in ResidualFunction()
298 Vm2 = Vm * Vm; in ResidualFunction()
301 …for (k = 0; k < ld_nsegsp[i]; k++) PD += ld_alphap[k] * PD0[i] * PetscPowScalar(Vm / Vm0, ld_betap… in ResidualFunction()
[all …]
|
| H A D | ex9busopt.c | 158 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
184 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess()
185 Vm2 = Vm * Vm; in SetInitialGuess()
227 Vref[i] = Vm + (VR / KA[i]); in SetInitialGuess()
246 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in InitialGuess() local
272 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in InitialGuess()
273 Vm2 = Vm * Vm; in InitialGuess()
355 PetscScalar Vr, Vi, Vm, Vm2; in ResidualFunction() local
426 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in ResidualFunction()
433 … fgen[idx + 8] = (VR - KA[i] * RF + KA[i] * KF[i] * Efd / TF[i] - KA[i] * (Vref[i] - Vm)) / TA[i]; in ResidualFunction()
[all …]
|
| H A D | ex9busdmnetwork.c | 366 PetscScalar Vr = 0, Vi = 0, Vm = 0, Vm2; /* Terminal voltage variables */ in SetInitialGuess() local
405 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess()
406 Vm2 = Vm * Vm; in SetInitialGuess()
454 exc->Vref = Vm + (VR / exc->KA); in SetInitialGuess()
527 PetscScalar Vm, Vm2, Vm0; in FormIFunction() local
680 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in FormIFunction()
686 …farr[offsetexc + 2] = (VR - KA * RF + KA * KF * Efd / TF - KA * (Vref - Vm)) / TA + xdotarr[offset… in FormIFunction()
710 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in FormIFunction()
711 Vm2 = Vm * Vm; in FormIFunction()
714 … for (k = 0; k < ld_nsegsp; k++) PD += ld_alphap[k] * PD0 * PetscPowScalar(Vm / Vm0, ld_betap[k]); in FormIFunction()
[all …]
|
| /petsc/src/snes/tutorials/network/power/ |
| H A D | pffunctions.c | 39 PetscScalar Vm; in FormJacobian_Power_private() local
80 Vm = xarr[offset + 1]; in FormJacobian_Power_private()
89 values[1] = 2.0 * Vm * bus->gl / Sbase; in FormJacobian_Power_private()
90 values[3] = -2.0 * Vm * bus->bl / Sbase; in FormJacobian_Power_private()
247 PetscScalar Vm; in FormFunction_Power() local
279 Vm = xarr[offset + 1]; in FormFunction_Power()
282 farr[offset] += Vm * Vm * bus->gl / Sbase; in FormFunction_Power()
283 if (bus->ide != PV_BUS) farr[offset + 1] += -Vm * Vm * bus->bl / Sbase; in FormFunction_Power()
|
| H A D | power2.c | 25 PetscScalar Vm; in FormFunction_Subnet() local
52 Vm = xarr[offset + 1]; in FormFunction_Subnet()
55 farr[offset] += Vm * Vm * bus->gl / Sbase; in FormFunction_Subnet()
56 if (bus->ide != PV_BUS) farr[offset + 1] += -Vm * Vm * bus->bl / Sbase; in FormFunction_Subnet()
173 PetscScalar Vm; in FormJacobian_Subnet() local
208 Vm = xarr[offset + 1]; in FormJacobian_Subnet()
217 values[1] = 2.0 * Vm * bus->gl / Sbase; in FormJacobian_Subnet()
218 values[3] = -2.0 * Vm * bus->bl / Sbase; in FormJacobian_Subnet()
|
| H A D | case9.m | 16 % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
|