| /petsc/src/snes/tutorials/network/power/ |
| H A D | pffunctions.c | 34 const PetscScalar *xarr; in FormJacobian_Power_private() local
47 PetscCall(VecGetArrayRead(localX, &xarr)); in FormJacobian_Power_private()
80 Vm = xarr[offset + 1]; in FormJacobian_Power_private()
129 thetaf = xarr[offsetfrom]; in FormJacobian_Power_private()
130 Vmf = xarr[offsetfrom + 1]; in FormJacobian_Power_private()
131 thetat = xarr[offsetto]; in FormJacobian_Power_private()
132 Vmt = xarr[offsetto + 1]; in FormJacobian_Power_private()
210 PetscCall(VecRestoreArrayRead(localX, &xarr)); in FormJacobian_Power_private()
244 const PetscScalar *xarr; in FormFunction_Power() local
256 PetscCall(VecGetArrayRead(localX, &xarr)); in FormFunction_Power()
[all …]
|
| H A D | power2.c | 15 const PetscScalar *xarr; in FormFunction_Subnet() local
20 PetscCall(VecGetArrayRead(localX, &xarr)); in FormFunction_Subnet()
46 farr[offset] = xarr[offset] - bus->va * PETSC_PI / 180.0; in FormFunction_Subnet()
47 farr[offset + 1] = xarr[offset + 1] - bus->vm; in FormFunction_Subnet()
52 Vm = xarr[offset + 1]; in FormFunction_Subnet()
86 thetaf = xarr[offsetfrom]; in FormFunction_Subnet()
87 Vmf = xarr[offsetfrom + 1]; in FormFunction_Subnet()
88 thetat = xarr[offsetto]; in FormFunction_Subnet()
89 Vmt = xarr[offsetto + 1]; in FormFunction_Subnet()
116 if (bus && bus->ide == PV_BUS) farr[offset + 1] = xarr[offset + 1] - bus->vm; in FormFunction_Subnet()
[all …]
|
| /petsc/src/snes/tutorials/network/water/ |
| H A D | waterfunctions.c | 23 const PetscScalar *xarr; in FormFunction_Water() local
37 PetscCall(VecGetArrayRead(localX, &xarr)); in FormFunction_Water()
59 hf = xarr[offsetnode1]; in FormFunction_Water()
60 ht = xarr[offsetnode2]; in FormFunction_Water()
88 farr[offset] = xarr[offset] - reservoir->head; in FormFunction_Water()
91 farr[offset] = xarr[offset] - (tank->elev + tank->initlvl); in FormFunction_Water()
95 PetscCall(VecRestoreArrayRead(localX, &xarr)); in FormFunction_Water()
203 PetscScalar *xarr; in SetInitialGuess_Water() local
206 PetscCall(VecGetArray(localX, &xarr)); in SetInitialGuess_Water()
215 xarr[offset] = 100; in SetInitialGuess_Water()
[all …]
|
| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busdmnetwork.c | 365 PetscScalar *xarr; in SetInitialGuess() local
384 PetscCall(VecGetArray(localX, &xarr)); in SetInitialGuess()
397 xarr[offset] = bus->vr; in SetInitialGuess()
398 xarr[offset + 1] = bus->vi; in SetInitialGuess()
433 xarr[offset] = Eqp; in SetInitialGuess()
434 xarr[offset + 1] = Edp; in SetInitialGuess()
435 xarr[offset + 2] = delta; in SetInitialGuess()
436 xarr[offset + 3] = W_S; in SetInitialGuess()
437 xarr[offset + 4] = Id; in SetInitialGuess()
438 xarr[offset + 5] = Iq; in SetInitialGuess()
[all …]
|
| H A D | ex9busopt.c | 970 PetscScalar *xarr; in DRDUJacobianTranspose() local
974 PetscCall(MatDenseGetColumn(DRDU, 0, &xarr)); in DRDUJacobianTranspose()
975 PetscCall(VecPlaceArray(drdu_col, xarr)); in DRDUJacobianTranspose()
998 PetscCall(MatDenseRestoreColumn(DRDU, &xarr)); in DRDUJacobianTranspose()
|
| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 280 PetscScalar *xarr; in TSAdjointStepBEuler_Private() local
310 PetscCall(MatDenseGetColumn(quadJ, nadj, &xarr)); in TSAdjointStepBEuler_Private()
311 PetscCall(VecPlaceArray(ts->vec_drdu_col, xarr)); in TSAdjointStepBEuler_Private()
314 PetscCall(MatDenseRestoreColumn(quadJ, &xarr)); in TSAdjointStepBEuler_Private()
336 PetscCall(MatDenseGetColumn(ts->mat_sensip, 0, &xarr)); in TSAdjointStepBEuler_Private()
337 PetscCall(VecPlaceArray(ts->vec_sensip_col, xarr)); in TSAdjointStepBEuler_Private()
393 PetscCall(MatDenseGetColumn(quadJp, nadj, &xarr)); in TSAdjointStepBEuler_Private()
394 PetscCall(VecPlaceArray(ts->vec_drdp_col, xarr)); in TSAdjointStepBEuler_Private()
397 PetscCall(MatDenseRestoreColumn(quadJp, &xarr)); in TSAdjointStepBEuler_Private()
410 PetscCall(MatDenseRestoreColumn(ts->mat_sensip, &xarr)); in TSAdjointStepBEuler_Private()
[all …]
|
| /petsc/src/snes/tutorials/network/ |
| H A D | ex1.c | 137 const PetscScalar *xarr, *xoldarr; in FormFunction_Dummy() local
145 PetscCall(VecGetArrayRead(localX, &xarr)); in FormFunction_Dummy()
155 for (j = 0; j < nvar; j++) farr[offset + j] = xarr[offset + j] - xoldarr[offset + j]; in FormFunction_Dummy()
158 PetscCall(VecRestoreArrayRead(localX, &xarr)); in FormFunction_Dummy()
270 PetscScalar *xarr; in SetInitialGuess() local
289 PetscCall(VecGetArray(localX, &xarr)); in SetInitialGuess()
301 xarr[offset] = bus->va * PETSC_PI / 180.0; in SetInitialGuess()
302 xarr[offset + 1] = bus->vm; in SetInitialGuess()
306 xarr[offset + 1] = gen->vs; in SetInitialGuess()
310 xarr[offset] = 100; in SetInitialGuess()
[all …]
|
| /petsc/src/ts/tutorials/network/ |
| H A D | pipes.c | 114 const PetscScalar *xarr, *xdotarr, *xoldarr; in WASHIFunction() local
136 PetscCall(VecGetArrayRead(localX, &xarr)); in WASHIFunction()
156 juncx = (PipeField *)(xarr + varoffset); in WASHIFunction()
172 pipex = (PipeField *)(xarr + varoffset); in WASHIFunction()
195 juncx = (PipeField *)(xarr + offsetfrom); in WASHIFunction()
204 juncx = (PipeField *)(xarr + offsetto); in WASHIFunction()
212 PetscCall(VecRestoreArrayRead(localX, &xarr)); in WASHIFunction()
232 PetscScalar *xarr; in WASHSetInitialSolution() local
241 PetscCall(VecGetArray(localX, &xarr)); in WASHSetInitialSolution()
255 for (k = 0; k < nx; k++) (xarr + varoffset)[k] = xarray[k]; in WASHSetInitialSolution()
[all …]
|
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | rk.c | 707 PetscScalar *xarr; in TSForwardStep_RK() local
708 PetscCall(MatDenseGetColumn(MatsFwdSensipTemp[i], 0, &xarr)); in TSForwardStep_RK()
709 PetscCall(VecPlaceArray(rk->VecDeltaFwdSensipCol, xarr)); in TSForwardStep_RK()
712 PetscCall(MatDenseRestoreColumn(MatsFwdSensipTemp[i], &xarr)); in TSForwardStep_RK()
879 PetscScalar *w = rk->work, *xarr; in TSAdjointStep_RK() local
920 PetscCall(MatDenseGetColumn(Jquad, nadj, &xarr)); in TSAdjointStep_RK()
921 PetscCall(VecPlaceArray(VecDRDUTransCol, xarr)); in TSAdjointStep_RK()
924 PetscCall(MatDenseRestoreColumn(Jquad, &xarr)); in TSAdjointStep_RK()
940 PetscCall(MatDenseGetColumn(quadts->Jacprhs, nadj, &xarr)); in TSAdjointStep_RK()
941 PetscCall(VecPlaceArray(VecDRDPTransCol, xarr)); in TSAdjointStep_RK()
[all …]
|
| /petsc/src/tao/constrained/impls/ipm/ |
| H A D | pdipm.c | 57 const PetscScalar *xarr, *xuarr, *xlarr, *garr, *harr; in TaoPDIPMUpdateConstraints() local
61 PetscCall(VecGetArrayRead(x, &xarr)); in TaoPDIPMUpdateConstraints()
81 carr[offset + k] = xarr[i] - xuarr[i]; in TaoPDIPMUpdateConstraints()
102 carr[offset + k] = xuarr[i] - xarr[i]; in TaoPDIPMUpdateConstraints()
112 carr[offset + k] = xarr[i] - xlarr[i]; in TaoPDIPMUpdateConstraints()
123 carr[offset + k] = xuarr[i] - xarr[i]; in TaoPDIPMUpdateConstraints()
124 carr[offset1 + k] = xarr[i] - xlarr[i]; in TaoPDIPMUpdateConstraints()
130 PetscCall(VecRestoreArrayRead(x, &xarr)); in TaoPDIPMUpdateConstraints()
225 const PetscScalar *xarr, *h; in TaoPDIPMInitializeSolution() local
231 PetscCall(VecGetArrayRead(tao->solution, &xarr)); in TaoPDIPMInitializeSolution()
[all …]
|
| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 1004 PetscScalar *xarr; in TSAdjointSetForward() local
1016 PetscCall(MatDenseGetColumn(A, 0, &xarr)); in TSAdjointSetForward()
1017 PetscCall(VecPlaceArray(sp, xarr)); in TSAdjointSetForward()
1029 PetscCall(MatDenseRestoreColumn(A, &xarr)); in TSAdjointSetForward()
|
| /petsc/src/mat/impls/hypre/ |
| H A D | mhypre.c | 2247 HYPRE_Complex *xarr, *yarr; local
2258 xarr = NULL;
2261 xarr = hypre_CSRMatrixData(xloc);
2270 …PetscCallBLAS("BLASaxpy", BLASaxpy_(&bnz, &a, (PetscScalar *)xarr, &one, (PetscScalar *)yarr, &one…
2277 xarr = NULL;
2280 xarr = hypre_CSRMatrixData(xloc);
2289 …PetscCallBLAS("BLASaxpy", BLASaxpy_(&bnz, &a, (PetscScalar *)xarr, &one, (PetscScalar *)yarr, &one…
|
| /petsc/src/vec/vec/impls/seq/cupm/ |
| H A D | vecseqcupm_impl.hpp | 1152 …ispatch_(PetscDeviceContext dctx, cupmStream_t stream, const PetscScalar *xarr, const Vec yin[], P… in MDot_kernel_dispatch_() argument
1162 xarr, size, results, DeviceArrayRead(dctx, yin[Idx]).data()... in MDot_kernel_dispatch_()
1171 …ispatch_(PetscDeviceContext dctx, cupmStream_t stream, const PetscScalar *xarr, const Vec yin[], P… in MDot_kernel_dispatch_() argument
1174 …PetscCall(MDot_kernel_dispatch_(dctx, stream, xarr, yin + yidx, size, results + yidx * MDOT_WORKGR… in MDot_kernel_dispatch_()
|
| /petsc/src/mat/impls/aij/seq/seqhipsparse/ |
| H A D | aijhipsparse.hip.cxx | 3085 thrust::device_ptr<PetscScalar> xarr = thrust::device_pointer_cast(xarray); in MatMultAddKernel_SeqAIJHIPSPARSE() local
3090 …e(hipsparsestruct->workVector->begin(), thrust::make_permutation_iterator(xarr, matstruct->cprowIn… in MatMultAddKernel_SeqAIJHIPSPARSE()
3091 …e(hipsparsestruct->workVector->begin(), thrust::make_permutation_iterator(xarr, matstruct->cprowIn… in MatMultAddKernel_SeqAIJHIPSPARSE()
|
| /petsc/src/mat/impls/aij/seq/seqcusparse/ |
| H A D | aijcusparse.cu | 3624 thrust::device_ptr<PetscScalar> xarr = thrust::device_pointer_cast(xarray); in MatMultAddKernel_SeqAIJCUSPARSE() local
3630 …le(cusparsestruct->workVector->begin(), thrust::make_permutation_iterator(xarr, matstruct->cprowIn… in MatMultAddKernel_SeqAIJCUSPARSE()
3631 …le(cusparsestruct->workVector->begin(), thrust::make_permutation_iterator(xarr, matstruct->cprowIn… in MatMultAddKernel_SeqAIJCUSPARSE()
|