| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busadj.c | 123 PetscScalar theta, Vd, Vq, SE; in SetInitialGuess() local
157 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess()
160 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
212 PetscScalar Vd, Vq, SE; in ResidualFunction() local
261 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in ResidualFunction()
271 fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id; in ResidualFunction()
272 fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq; in ResidualFunction()
280 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in ResidualFunction()
438 PetscScalar Vd, Vq; in ResidualJacobian() local
515 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in ResidualJacobian()
[all …]
|
| H A D | ex9busdmnetwork.c | 370 PetscScalar Vd, Vq; /* Generator dq axis voltages */ in SetInitialGuess() local
422 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess()
426 Edp = Vd + gen->Rs * Id - gen->Xqp * Iq; in SetInitialGuess()
491 PetscScalar Vd = 0, Vq = 0, SE; in FormIFunction() local
639 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in FormIFunction()
649 farr[offsetgen + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id; in FormIFunction()
650 farr[offsetgen + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq; in FormIFunction()
680 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in FormIFunction()
848 PetscScalar Vd, Vq, IGr, IGi, Zdq_inv[4], det; in AlgFunction() local
873 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in AlgFunction()
[all …]
|
| H A D | ex9busoptfd.c | 125 PetscScalar theta, Vd, Vq, SE; in SetInitialGuess() local
159 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess()
162 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
214 PetscScalar Vd, Vq, SE; in ResidualFunction() local
263 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in ResidualFunction()
272 fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id; in ResidualFunction()
273 fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq; in ResidualFunction()
281 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in ResidualFunction()
439 PetscScalar Vd, Vq; in ResidualJacobian() local
516 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in ResidualJacobian()
[all …]
|
| H A D | ex9busopt.c | 162 PetscScalar theta, Vd, Vq, SE; in SetInitialGuess() local
196 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess()
199 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
250 PetscScalar theta, Vd, Vq, SE; in InitialGuess() local
284 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in InitialGuess()
287 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in InitialGuess()
359 PetscScalar Vd, Vq, SE; in ResidualFunction() local
408 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in ResidualFunction()
417 fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id; in ResidualFunction()
418 fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq; in ResidualFunction()
[all …]
|
| H A D | ex9bus.c | 349 PetscScalar theta, Vd, Vq, SE; in SetInitialGuess() local
383 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess()
386 Edp = Vd + Rs[i] * Id - Xqp[i] * Iq; /* d-axis transient EMF */ in SetInitialGuess()
440 PetscScalar Vd, Vq, SE; in ResidualFunction() local
489 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in ResidualFunction()
498 fgen[idx + 4] = Zdq_inv[0] * (-Edp + Vd) + Zdq_inv[1] * (-Eqp + Vq) + Id; in ResidualFunction()
499 fgen[idx + 5] = Zdq_inv[2] * (-Edp + Vd) + Zdq_inv[3] * (-Eqp + Vq) + Iq; in ResidualFunction()
507 Vm = PetscSqrtScalar(Vd * Vd + Vq * Vq); in ResidualFunction()
696 PetscScalar Vd, Vq; in ResidualJacobian() local
774 PetscCall(ri2dq(Vr, Vi, delta, &Vd, &Vq)); in ResidualJacobian()
[all …]
|