| /petsc/src/mat/tests/ |
| H A D | ex70.c | 576 Mat TE; in main() local
579 PetscCall(MatComputeOperator(T, MATDENSE, &TE)); in main()
580 PetscCall(MatView(TE, NULL)); in main()
582 PetscCall(MatAXPY(TE, -1.0, PtAP, SAME_NONZERO_PATTERN)); in main()
583 PetscCall(MatView(TE, NULL)); in main()
584 PetscCall(MatDestroy(&TE)); in main()
597 Mat TE; in main() local
600 PetscCall(MatComputeOperator(T, MATDENSE, &TE)); in main()
601 PetscCall(MatView(TE, NULL)); in main()
604 PetscCall(MatAXPY(TE, -1.0, RARt, SAME_NONZERO_PATTERN)); in main()
[all …]
|
| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busdmnetwork.c | 48 PetscScalar TE; /* Exciter time constant */ member
139 static const PetscScalar TE[3] = {0.314, 0.314, 0.314}; /* Exciter time constant */ in read_data() local
229 exc[i * 3 + j].TE = TE[j]; in read_data()
660 PetscScalar k1, k2, KE, TE, TF, KA, KF, Vref, TA; /* Generator parameters */ in FormIFunction() local
673 TE = exc->TE; in FormIFunction()
684 farr[offsetexc] = (KE * Efd + SE - VR) / TE + xdotarr[offsetexc]; in FormIFunction()
|
| H A D | ex9busadj.c | 49 const PetscScalar TE[3] = {0.314, 0.314, 0.314}; /* Exciter time constant */ variable
285 fgen[idx + 6] = (KE[i] * Efd + SE - VR) / TE[i]; in ResidualFunction()
599 val[0] = (KE[i] + dSE_dEfd) / TE[i]; in ResidualJacobian()
600 val[1] = -1 / TE[i]; in ResidualJacobian()
|
| H A D | ex9busoptfd.c | 46 const PetscScalar TE[3] = {0.314, 0.314, 0.314}; /* Exciter time constant */ variable
286 fgen[idx + 6] = (KE[i] * Efd + SE - VR) / TE[i]; in ResidualFunction()
601 val[0] = (KE[i] + dSE_dEfd) / TE[i]; in ResidualJacobian()
602 val[1] = -1 / TE[i]; in ResidualJacobian()
|
| H A D | ex9bus.c | 76 const PetscScalar TE[3] = {0.314, 0.314, 0.314}; /* Exciter time constant */ variable
512 fgen[idx + 6] = (-KE[i] * Efd - SE + VR) / TE[i]; in ResidualFunction()
859 val[0] = (-KE[i] - dSE_dEfd) / TE[i]; in ResidualJacobian()
860 val[1] = 1 / TE[i]; in ResidualJacobian()
|
| H A D | ex9busopt.c | 56 const PetscScalar TE[3] = {0.314, 0.314, 0.314}; /* Exciter time constant */ variable
431 fgen[idx + 6] = (KE[i] * Efd + SE - VR) / TE[i]; in ResidualFunction()
743 val[0] = (KE[i] + dSE_dEfd) / TE[i]; in ResidualJacobian()
744 val[1] = -1 / TE[i]; in ResidualJacobian()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex30.c | 468 PetscScalar TE = 0.0, TN = 0.0, TS = 0.0, TW = 0.0, dPdx, residual, z_scale; in XMomentumResidual() local
479 TE = param->potentialT * x[j][i + 1].T * PetscExpScalar((j - 0.5) * dz * z_scale); in XMomentumResidual()
490 etaE = Viscosity(TE, epsE, dz * j, param); in XMomentumResidual()
519 PetscScalar TE = 0.0, TN = 0.0, TS = 0.0, TW = 0.0, dPdz, residual, z_scale; in ZMomentumResidual() local
531 if (i == ilim) TE = TW; in ZMomentumResidual()
532 else TE = param->potentialT * TInterp(x, i, j) * PetscExpScalar(j * dz * z_scale); in ZMomentumResidual()
543 etaE = Viscosity(TE, epsE, dz * (j + 0.5), param); in ZMomentumResidual()
589 PetscScalar TE, TN, TS, TW, residual; in EnergyResidual() local
621 TE = (x[j][i].T + x[j][i + 1].T) / 2.0; in EnergyResidual()
625 fE = uE * TE * dz; in EnergyResidual()
|