| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busadj.c | 123 PetscScalar theta, Vd, Vq, SE; in SetInitialGuess() local
158 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess()
161 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-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
423 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess()
429 Eqp = Vq + gen->Rs * Iq + gen->Xdp * Id; 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
160 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess()
163 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-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
197 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess()
200 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-axis transient EMF */ in SetInitialGuess()
250 PetscScalar theta, Vd, Vq, SE; in InitialGuess() local
285 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in InitialGuess()
288 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-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
384 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess()
387 Eqp = Vq + Rs[i] * Iq + Xdp[i] * Id; /* q-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 …]
|