| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busdmnetwork.c | 367 PetscScalar IGr, IGi; /* Generator real and imaginary current */ in SetInitialGuess() local
410 IGi = (Vi * gen->PG - Vr * gen->QG) / Vm2; in SetInitialGuess()
413 delta = atan2(Vi + gen->Xq * IGr, Vr - gen->Xq * IGi); in SetInitialGuess()
417 Id = IGr * PetscCosScalar(theta) - IGi * PetscSinScalar(theta); in SetInitialGuess()
420 Iq = IGr * PetscSinScalar(theta) + IGi * PetscCosScalar(theta); in SetInitialGuess()
604 PetscScalar IGr, IGi, Zdq_inv[4], det; in FormIFunction() local
652 PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi)); in FormIFunction()
655 farr[offsetbus] -= IGi; in FormIFunction()
848 PetscScalar Vd, Vq, IGr, IGi, Zdq_inv[4], det; in AlgFunction() local
887 PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi)); in AlgFunction()
[all …]
|
| H A D | ex9busopt.c | 158 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
187 IGi = (Vi * PG[i] - Vr * QG[i]) / Vm2; in SetInitialGuess()
189 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
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()
246 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in InitialGuess() local
275 IGi = (Vi * PGv[i] - Vr * QG[i]) / Vm2; in InitialGuess()
277 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in InitialGuess()
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 | ex9busadj.c | 119 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
148 IGi = (Vi * PG[i] - Vr * QG[i]) / Vm2; in SetInitialGuess()
150 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
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()
213 PetscScalar IGr, IGi, IDr, IDi; in ResidualFunction() local
275 PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi)); in ResidualFunction()
277 fnet[2 * gbus[i]] -= IGi; in ResidualFunction()
|
| H A D | ex9busoptfd.c | 121 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
150 IGi = (Vi * PG[i] - Vr * QG[i]) / Vm2; in SetInitialGuess()
152 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
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()
215 PetscScalar IGr, IGi, IDr, IDi; in ResidualFunction() local
276 PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi)); in ResidualFunction()
278 fnet[2 * gbus[i]] -= IGi; in ResidualFunction()
|
| H A D | ex9bus.c | 345 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local
374 IGi = (Vi * PG[i] - Vr * QG[i]) / Vm2; in SetInitialGuess()
376 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
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()
441 PetscScalar IGr, IGi, IDr, IDi; in ResidualFunction() local
502 PetscCall(dq2ri(Id, Iq, delta, &IGr, &IGi)); in ResidualFunction()
504 fnet[2 * gbus[i]] -= IGi; in ResidualFunction()
|