| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busadj.c | 119 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local 143 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in SetInitialGuess() 145 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess() 147 IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2; in SetInitialGuess() 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() 157 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 158 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 208 PetscScalar Vr, Vi, Vm, Vm2; in ResidualFunction() local 258 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in ResidualFunction() [all …]
|
| H A D | ex9bus.c | 142 PetscScalar Efd, RF, VR, Vr, Vi, Vm; in EventFunction() local 161 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in EventFunction() 163 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in EventFunction() 194 PetscScalar Vr, Vi, Vm; in PostEventFunction() local 245 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in PostEventFunction() 247 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in PostEventFunction() 345 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local 369 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in SetInitialGuess() 371 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess() 373 IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2; in SetInitialGuess() [all …]
|
| H A D | ex9busoptfd.c | 121 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local 145 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in SetInitialGuess() 147 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess() 149 IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2; in SetInitialGuess() 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() 159 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 160 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 210 PetscScalar Vr, Vi, Vm, Vm2; in ResidualFunction() local 260 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in ResidualFunction() [all …]
|
| H A D | ex9busopt.c | 158 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in SetInitialGuess() local 182 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in SetInitialGuess() 184 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess() 186 IGr = (Vr * PG[i] + Vi * QG[i]) / Vm2; in SetInitialGuess() 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() 196 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 197 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 246 PetscScalar Vr, Vi, IGr, IGi, Vm, Vm2; in InitialGuess() local 270 Vr = xnet[2 * gbus[i]]; /* Real part of generator terminal voltage */ in InitialGuess() [all …]
|
| H A D | ex9busdmnetwork.c | 366 PetscScalar Vr = 0, Vi = 0, Vm = 0, Vm2; /* Terminal voltage variables */ in SetInitialGuess() local 400 Vr = bus->vr; in SetInitialGuess() 405 Vm = PetscSqrtScalar(Vr * Vr + Vi * Vi); in SetInitialGuess() 408 IGr = (Vr * gen->PG + Vi * gen->QG) / Vm2; in SetInitialGuess() 410 IGi = (Vi * gen->PG - Vr * gen->QG) / Vm2; in SetInitialGuess() 413 delta = atan2(Vi + gen->Xq * IGr, Vr - gen->Xq * IGi); in SetInitialGuess() 422 Vd = Vr * PetscCosScalar(theta) - Vi * PetscSinScalar(theta); in SetInitialGuess() 423 Vq = Vr * PetscSinScalar(theta) + Vi * PetscCosScalar(theta); in SetInitialGuess() 494 PetscScalar Vr = 0, Vi = 0; in FormIFunction() local 543 Vr = xarr[offsetbus]; in FormIFunction() [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex20opt_p.c | 86 static PetscErrorCode RHSHessianProductUU(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductUU() argument 97 PetscCall(VecGetArrayRead(Vr, &vr)); in RHSHessianProductUU() 111 PetscCall(VecRestoreArrayRead(Vr, &vr)); in RHSHessianProductUU() 116 static PetscErrorCode RHSHessianProductUP(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductUP() argument 126 PetscCall(VecGetArrayRead(Vr, &vr)); in RHSHessianProductUP() 139 PetscCall(VecRestoreArrayRead(Vr, &vr)); in RHSHessianProductUP() 144 static PetscErrorCode RHSHessianProductPU(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductPU() argument 154 PetscCall(VecGetArrayRead(Vr, &vr)); in RHSHessianProductPU() 167 PetscCall(VecRestoreArrayRead(Vr, &vr)); in RHSHessianProductPU() 172 static PetscErrorCode RHSHessianProductPP(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductPP() argument [all …]
|
| H A D | ex20opt_ic.c | 78 static PetscErrorCode RHSHessianProductUU(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductUU() argument 89 PetscCall(VecGetArrayRead(Vr, &vr)); in RHSHessianProductUU() 102 PetscCall(VecRestoreArrayRead(Vr, &vr)); in RHSHessianProductUU() 175 static PetscErrorCode IHessianProductUU(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Petsc… in IHessianProductUU() argument 186 PetscCall(VecGetArrayRead(Vr, &vr)); in IHessianProductUU() 198 PetscCall(VecRestoreArrayRead(Vr, &vr)); in IHessianProductUU()
|
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 469 …cErrorCode (*ihessianproductfunc1)(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, PetscCtx … in TSSetIHessianProduct() argument 510 PetscErrorCode TSComputeIHessianProductFunctionUU(TS ts, PetscReal t, Vec U, Vec Vl[], Vec Vr, Vec … in TSComputeIHessianProductFunctionUU() argument 517 …duct 1 for sensitivity analysis", (*ts->ihessianproduct_fuu)(ts, t, U, Vl, Vr, VHV, ts->ihessianpr… in TSComputeIHessianProductFunctionUU() 522 PetscCall(TSComputeRHSHessianProductFunctionUU(ts, t, U, Vl, Vr, VHV)); in TSComputeIHessianProductFunctionUU() 551 PetscErrorCode TSComputeIHessianProductFunctionUP(TS ts, PetscReal t, Vec U, Vec Vl[], Vec Vr, Vec … in TSComputeIHessianProductFunctionUP() argument 558 …duct 2 for sensitivity analysis", (*ts->ihessianproduct_fup)(ts, t, U, Vl, Vr, VHV, ts->ihessianpr… in TSComputeIHessianProductFunctionUP() 563 PetscCall(TSComputeRHSHessianProductFunctionUP(ts, t, U, Vl, Vr, VHV)); in TSComputeIHessianProductFunctionUP() 592 PetscErrorCode TSComputeIHessianProductFunctionPU(TS ts, PetscReal t, Vec U, Vec Vl[], Vec Vr, Vec … in TSComputeIHessianProductFunctionPU() argument 599 …duct 3 for sensitivity analysis", (*ts->ihessianproduct_fpu)(ts, t, U, Vl, Vr, VHV, ts->ihessianpr… in TSComputeIHessianProductFunctionPU() 604 PetscCall(TSComputeRHSHessianProductFunctionPU(ts, t, U, Vl, Vr, VHV)); in TSComputeIHessianProductFunctionPU() [all …]
|
| /petsc/src/ts/tutorials/optimal_control/ |
| H A D | ex1.c | 78 static PetscErrorCode RHSHessianProductUU(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductUU() argument 84 static PetscErrorCode RHSHessianProductUP(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductUP() argument 90 static PetscErrorCode RHSHessianProductPU(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductPU() argument 96 static PetscErrorCode RHSHessianProductPP(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VHV, Pet… in RHSHessianProductPP() argument 110 PetscCall(VecGetArrayRead(Vr, &vr)); in RHSHessianProductPP() 128 PetscCall(VecRestoreArrayRead(Vr, &vr)); in RHSHessianProductPP() 134 static PetscErrorCode IntegrandHessianProductUU(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VH… in IntegrandHessianProductUU() argument 147 PetscCall(VecGetArrayRead(Vr, &vr)); in IntegrandHessianProductUU() 158 PetscCall(VecRestoreArrayRead(Vr, &vr)); in IntegrandHessianProductUU() 163 static PetscErrorCode IntegrandHessianProductUP(TS ts, PetscReal t, Vec U, Vec *Vl, Vec Vr, Vec *VH… in IntegrandHessianProductUP() argument [all …]
|