| /petsc/src/snes/tutorials/network/power/ |
| H A D | power.c | 14 DM networkdm; in FormFunction() local
21 PetscCall(SNESGetDM(snes, &networkdm)); in FormFunction()
22 PetscCall(DMGetLocalVector(networkdm, &localX)); in FormFunction()
23 PetscCall(DMGetLocalVector(networkdm, &localF)); in FormFunction()
27 PetscCall(DMGlobalToLocalBegin(networkdm, X, INSERT_VALUES, localX)); in FormFunction()
28 PetscCall(DMGlobalToLocalEnd(networkdm, X, INSERT_VALUES, localX)); in FormFunction()
30 PetscCall(DMNetworkGetSubnetwork(networkdm, 0, &nv, &ne, &vtx, &edges)); in FormFunction()
31 PetscCall(FormFunction_Power(networkdm, localX, localF, nv, ne, vtx, edges, User)); in FormFunction()
33 PetscCall(DMRestoreLocalVector(networkdm, &localX)); in FormFunction()
35 PetscCall(DMLocalToGlobalBegin(networkdm, localF, ADD_VALUES, F)); in FormFunction()
[all …]
|
| H A D | power2.c | 11 PetscErrorCode FormFunction_Subnet(DM networkdm, Vec localX, Vec localF, PetscInt nv, PetscInt ne, … in FormFunction_Subnet() argument
34 PetscCall(DMNetworkIsGhostVertex(networkdm, vtx[v], &ghostvtex)); in FormFunction_Subnet()
35 PetscCall(DMNetworkGetNumComponents(networkdm, vtx[v], &numComps)); in FormFunction_Subnet()
36 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vtx[v], ALL_COMPONENTS, &offset)); in FormFunction_Subnet()
38 PetscCall(DMNetworkGetComponent(networkdm, vtx[v], j, &key, &component, NULL)); in FormFunction_Subnet()
59 PetscCall(DMNetworkGetSupportingEdges(networkdm, vtx[v], &nconnedges, &connedges)); in FormFunction_Subnet()
68 PetscCall(DMNetworkGetComponent(networkdm, e, 0, &keye, (void **)&branch, NULL)); in FormFunction_Subnet()
79 PetscCall(DMNetworkGetConnectedVertices(networkdm, e, &cone)); in FormFunction_Subnet()
83 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vfrom, ALL_COMPONENTS, &offsetfrom)); in FormFunction_Subnet()
84 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vto, ALL_COMPONENTS, &offsetto)); in FormFunction_Subnet()
[all …]
|
| H A D | pffunctions.c | 32 PetscErrorCode FormJacobian_Power_private(DM networkdm, Vec localX, Mat J, PetscInt nv, PetscInt ne… in FormJacobian_Power_private() argument
50 PetscCall(DMNetworkIsGhostVertex(networkdm, vtx[v], &ghostvtex)); in FormJacobian_Power_private()
52 PetscCall(DMNetworkGetNumComponents(networkdm, vtx[v], &numComps)); in FormJacobian_Power_private()
54 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vtx[v], ALL_COMPONENTS, &offset)); in FormJacobian_Power_private()
55 PetscCall(DMNetworkGetGlobalVecOffset(networkdm, vtx[v], ALL_COMPONENTS, &goffset)); in FormJacobian_Power_private()
56 PetscCall(DMNetworkGetComponent(networkdm, vtx[v], j, &key, &component, NULL)); in FormJacobian_Power_private()
95 PetscCall(DMNetworkGetSupportingEdges(networkdm, vtx[v], &nconnedges, &connedges)); in FormJacobian_Power_private()
106 PetscCall(DMNetworkGetComponent(networkdm, e, 0, &key, (void **)&branch, NULL)); in FormJacobian_Power_private()
118 PetscCall(DMNetworkGetConnectedVertices(networkdm, edges[e], &cone)); in FormJacobian_Power_private()
122 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vfrom, ALL_COMPONENTS, &offsetfrom)); in FormJacobian_Power_private()
[all …]
|
| /petsc/src/ksp/ksp/tutorials/network/ |
| H A D | ex1_nest.c | 138 PetscErrorCode FormOperator(DM networkdm, Mat A, Vec b) in FormOperator() argument
157 PetscCall(DMGetLocalVector(networkdm, &localb)); in FormOperator()
171 PetscCall(DMNetworkGetEdgeRange(networkdm, &eStart, &eEnd)); in FormOperator()
172 PetscCall(DMNetworkGetVertexRange(networkdm, &vStart, &vEnd)); in FormOperator()
175 PetscCall(DMNetworkGetComponent(networkdm, e, 0, &key, (void **)&branch, NULL)); in FormOperator()
176 PetscCall(DMNetworkGetEdgeOffset(networkdm, e, &lofst)); in FormOperator()
178 PetscCall(DMNetworkGetConnectedVertices(networkdm, e, &cone)); in FormOperator()
179 PetscCall(DMNetworkGetVertexOffset(networkdm, cone[0], &lofst_fr)); in FormOperator()
180 PetscCall(DMNetworkGetVertexOffset(networkdm, cone[1], &lofst_to)); in FormOperator()
197 PetscCall(DMNetworkGetComponent(networkdm, cone[0], 0, &key, (void **)&node, NULL)); in FormOperator()
[all …]
|
| H A D | ex2.c | 167 PetscErrorCode FormOperator(DM networkdm, Mat A, Vec b) in FormOperator() argument
179 PetscCall(DMGetLocalVector(networkdm, &localb)); in FormOperator()
194 PetscCall(DMNetworkGetEdgeRange(networkdm, &eStart, &eEnd)); in FormOperator()
196 PetscCall(DMNetworkGetComponent(networkdm, e, 0, NULL, (void **)&branch, NULL)); in FormOperator()
197 PetscCall(DMNetworkGetLocalVecOffset(networkdm, e, ALL_COMPONENTS, &lofst)); in FormOperator()
199 PetscCall(DMNetworkGetConnectedVertices(networkdm, e, &cone)); in FormOperator()
200 PetscCall(DMNetworkGetLocalVecOffset(networkdm, cone[0], ALL_COMPONENTS, &lofst_fr)); in FormOperator()
201 PetscCall(DMNetworkGetLocalVecOffset(networkdm, cone[1], ALL_COMPONENTS, &lofst_to)); in FormOperator()
215 PetscCall(DMNetworkGetComponent(networkdm, cone[0], 0, NULL, (void **)&node, NULL)); in FormOperator()
225 PetscCall(DMNetworkGetComponent(networkdm, cone[1], 0, NULL, (void **)&node, NULL)); in FormOperator()
[all …]
|
| /petsc/src/snes/tutorials/network/ |
| H A D | ex1.c | 50 DM networkdm, dmcoords; in UserMonitor() local
57 PetscCall(SNESGetDM(snes, &networkdm)); in UserMonitor()
58 PetscCall(DMGetCoordinateDM(networkdm, &dmcoords)); in UserMonitor()
68 PetscCall(DMNetworkIsGhostVertex(networkdm, v, &ghost)); in UserMonitor()
70 PetscCall(DMNetworkGetGlobalVertexIndex(networkdm, v, &gidx)); in UserMonitor()
72 PetscCall(DMNetworkGetGlobalVecOffset(networkdm, v, 0, &offset)); in UserMonitor()
81 PetscCall(DMView(networkdm, PETSC_VIEWER_DRAW_WORLD)); in UserMonitor()
89 DM networkdm; in FormJacobian_subPower() local
99 PetscCall(SNESGetDM(snes, &networkdm)); in FormJacobian_subPower()
100 PetscCall(DMGetLocalVector(networkdm, &localX)); in FormJacobian_subPower()
[all …]
|
| /petsc/src/snes/tutorials/network/water/ |
| H A D | waterfunctions.c | 21 PetscErrorCode FormFunction_Water(DM networkdm, Vec localX, Vec localF, PetscInt nv, PetscInt ne, c… in FormFunction_Water() argument
42 PetscCall(DMNetworkGetComponent(networkdm, edges[i], 0, &key, (void **)&edge, NULL)); in FormFunction_Water()
45 PetscCall(DMNetworkGetConnectedVertices(networkdm, edges[i], &cone)); in FormFunction_Water()
50 PetscCall(DMNetworkGetNumComponents(networkdm, vnode1, &ncomp)); in FormFunction_Water()
51 … PetscCall(DMNetworkGetComponent(networkdm, vnode1, ncomp - 1, &key, (void **)&vertexnode1, NULL)); in FormFunction_Water()
52 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vnode1, ncomp - 1, &offsetnode1)); in FormFunction_Water()
54 PetscCall(DMNetworkGetNumComponents(networkdm, vnode2, &ncomp)); in FormFunction_Water()
55 … PetscCall(DMNetworkGetComponent(networkdm, vnode2, ncomp - 1, &key, (void **)&vertexnode2, NULL)); in FormFunction_Water()
56 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vnode2, ncomp - 1, &offsetnode2)); in FormFunction_Water()
77 PetscCall(DMNetworkIsGhostVertex(networkdm, vtx[i], &ghostvtex)); in FormFunction_Water()
[all …]
|
| H A D | water.c | 16 DM networkdm; in main() local
29 PetscCall(DMNetworkCreate(PETSC_COMM_WORLD, &networkdm)); in main()
32 …PetscCall(DMNetworkRegisterComponent(networkdm, "edgestruct", sizeof(struct _p_EDGE_Water), &appct… in main()
33 …PetscCall(DMNetworkRegisterComponent(networkdm, "busstruct", sizeof(struct _p_VERTEX_Water), &appc… in main()
54 PetscCall(DMNetworkSetNumSubNetworks(networkdm, PETSC_DECIDE, 1)); in main()
55 PetscCall(DMNetworkAddSubnetwork(networkdm, "", waterdata->nedge, edgelist, NULL)); in main()
59 PetscCall(DMNetworkLayoutSetUp(networkdm)); in main()
64 PetscCall(DMNetworkGetSubnetwork(networkdm, 0, &nv, &ne, &vtx, &edges)); in main()
66 …for (i = 0; i < ne; i++) PetscCall(DMNetworkAddComponent(networkdm, edges[i], appctx.compkey_edge,… in main()
68 …for (i = 0; i < nv; i++) PetscCall(DMNetworkAddComponent(networkdm, vtx[i], appctx.compkey_vtx, &w… in main()
[all …]
|
| /petsc/src/ts/tutorials/network/ |
| H A D | pipes.c | 101 DM networkdm; in WASHIFunction() local
121 PetscCall(TSGetDM(ts, &networkdm)); in WASHIFunction()
123 PetscCall(DMGetLocalVector(networkdm, &localF)); in WASHIFunction()
130 PetscCall(DMGlobalToLocalBegin(networkdm, X, INSERT_VALUES, localX)); in WASHIFunction()
131 PetscCall(DMGlobalToLocalEnd(networkdm, X, INSERT_VALUES, localX)); in WASHIFunction()
133 PetscCall(DMGlobalToLocalBegin(networkdm, Xdot, INSERT_VALUES, localXdot)); in WASHIFunction()
134 PetscCall(DMGlobalToLocalEnd(networkdm, Xdot, INSERT_VALUES, localXdot)); in WASHIFunction()
149 PetscCall(DMNetworkGetVertexRange(networkdm, &vStart, &vEnd)); in WASHIFunction()
151 PetscCall(DMNetworkIsGhostVertex(networkdm, v, &ghost)); in WASHIFunction()
154 PetscCall(DMNetworkGetComponent(networkdm, v, 0, &type, (void **)&junction, NULL)); in WASHIFunction()
[all …]
|
| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busdmnetwork.c | 356 PetscErrorCode SetInitialGuess(DM networkdm, Vec X) in SetInitialGuess() argument
377 PetscCall(DMNetworkGetVertexRange(networkdm, &vStart, &vEnd)); in SetInitialGuess()
378 PetscCall(DMGetLocalVector(networkdm, &localX)); in SetInitialGuess()
381 PetscCall(DMGlobalToLocalBegin(networkdm, X, INSERT_VALUES, localX)); in SetInitialGuess()
382 PetscCall(DMGlobalToLocalEnd(networkdm, X, INSERT_VALUES, localX)); in SetInitialGuess()
387 PetscCall(DMNetworkIsGhostVertex(networkdm, v, &ghostvtex)); in SetInitialGuess()
390 PetscCall(DMNetworkGetNumComponents(networkdm, v, &numComps)); in SetInitialGuess()
392 PetscCall(DMNetworkGetComponent(networkdm, v, j, &key, &component, NULL)); in SetInitialGuess()
396 PetscCall(DMNetworkGetLocalVecOffset(networkdm, v, j, &offset)); in SetInitialGuess()
404 PetscCall(DMNetworkGetLocalVecOffset(networkdm, v, j, &offset)); in SetInitialGuess()
[all …]
|
| /petsc/src/dm/impls/network/ |
| H A D | networkcreate.c | 42 static PetscErrorCode VecView_Network_Seq(DM networkdm, Vec X, PetscViewer viewer) in VecView_Network_Seq() argument
51 PetscCall(DMNetworkGetEdgeRange(networkdm, &Start, &End)); in VecView_Network_Seq()
53 PetscCall(DMNetworkGetComponent(networkdm, e, ALL_COMPONENTS, NULL, NULL, &nvar)); in VecView_Network_Seq()
56 PetscCall(DMNetworkGetLocalVecOffset(networkdm, e, ALL_COMPONENTS, &offset)); in VecView_Network_Seq()
57 PetscCall(DMNetworkGetGlobalEdgeIndex(networkdm, e, &id)); in VecView_Network_Seq()
64 PetscCall(DMNetworkGetVertexRange(networkdm, &Start, &End)); in VecView_Network_Seq()
66 PetscCall(DMNetworkGetComponent(networkdm, v, ALL_COMPONENTS, NULL, NULL, &nvar)); in VecView_Network_Seq()
69 PetscCall(DMNetworkGetLocalVecOffset(networkdm, v, ALL_COMPONENTS, &offset)); in VecView_Network_Seq()
70 PetscCall(DMNetworkGetGlobalVertexIndex(networkdm, v, &id)); in VecView_Network_Seq()
80 static PetscErrorCode VecView_Network_MPI(DM networkdm, Vec X, PetscViewer viewer) in VecView_Network_MPI() argument
[all …]
|