Lines Matching refs:networkdm
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()
125 DM networkdm; in FormFunction() local
131 PetscCall(SNESGetDM(snes, &networkdm)); in FormFunction()
132 PetscCall(DMGetLocalVector(networkdm, &localX)); in FormFunction()
133 PetscCall(DMGetLocalVector(networkdm, &localF)); in FormFunction()
136 PetscCall(DMGlobalToLocalBegin(networkdm, X, INSERT_VALUES, localX)); in FormFunction()
137 PetscCall(DMGlobalToLocalEnd(networkdm, X, INSERT_VALUES, localX)); in FormFunction()
139 PetscCall(DMGlobalToLocalBegin(networkdm, F, INSERT_VALUES, localF)); in FormFunction()
140 PetscCall(DMGlobalToLocalEnd(networkdm, F, INSERT_VALUES, localF)); in FormFunction()
143 PetscCall(DMNetworkGetSubnetwork(networkdm, 0, &nv, &ne, &vtx, &edges)); in FormFunction()
144 PetscCall(FormFunction_Subnet(networkdm, localX, localF, nv, ne, vtx, edges, appctx)); in FormFunction()
147 PetscCall(DMNetworkGetSubnetwork(networkdm, 1, &nv, &ne, &vtx, &edges)); in FormFunction()
148 PetscCall(FormFunction_Subnet(networkdm, localX, localF, nv, ne, vtx, edges, appctx)); in FormFunction()
150 PetscCall(DMRestoreLocalVector(networkdm, &localX)); in FormFunction()
152 PetscCall(DMLocalToGlobalBegin(networkdm, localF, ADD_VALUES, F)); in FormFunction()
153 PetscCall(DMLocalToGlobalEnd(networkdm, localF, ADD_VALUES, F)); in FormFunction()
154 PetscCall(DMRestoreLocalVector(networkdm, &localF)); in FormFunction()
158 PetscErrorCode FormJacobian_Subnet(DM networkdm, Vec localX, Mat J, Mat Jpre, PetscInt nv, PetscInt… in FormJacobian_Subnet() argument
180 PetscCall(DMNetworkIsGhostVertex(networkdm, vtx[v], &ghostvtex)); in FormJacobian_Subnet()
181 PetscCall(DMNetworkGetNumComponents(networkdm, vtx[v], &numComps)); in FormJacobian_Subnet()
183 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vtx[v], ALL_COMPONENTS, &offset)); in FormJacobian_Subnet()
184 PetscCall(DMNetworkGetGlobalVecOffset(networkdm, vtx[v], ALL_COMPONENTS, &goffset)); in FormJacobian_Subnet()
185 PetscCall(DMNetworkGetComponent(networkdm, vtx[v], j, &key, &component, NULL)); in FormJacobian_Subnet()
223 PetscCall(DMNetworkGetSupportingEdges(networkdm, vtx[v], &nconnedges, &connedges)); in FormJacobian_Subnet()
233 PetscCall(DMNetworkGetComponent(networkdm, e, 0, &key, (void **)&branch, NULL)); in FormJacobian_Subnet()
245 PetscCall(DMNetworkGetConnectedVertices(networkdm, e, &cone)); in FormJacobian_Subnet()
249 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vfrom, ALL_COMPONENTS, &offsetfrom)); in FormJacobian_Subnet()
250 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vto, ALL_COMPONENTS, &offsetto)); in FormJacobian_Subnet()
251 PetscCall(DMNetworkGetGlobalVecOffset(networkdm, vfrom, ALL_COMPONENTS, &goffsetfrom)); in FormJacobian_Subnet()
252 PetscCall(DMNetworkGetGlobalVecOffset(networkdm, vto, ALL_COMPONENTS, &goffsetto)); in FormJacobian_Subnet()
263 PetscCall(DMNetworkGetComponent(networkdm, vfrom, 0, &keyf, (void **)&busf, NULL)); in FormJacobian_Subnet()
264 PetscCall(DMNetworkGetComponent(networkdm, vto, 0, &keyt, (void **)&bust, NULL)); in FormJacobian_Subnet()
341 DM networkdm; in FormJacobian() local
349 PetscCall(SNESGetDM(snes, &networkdm)); in FormJacobian()
350 PetscCall(DMGetLocalVector(networkdm, &localX)); in FormJacobian()
352 PetscCall(DMGlobalToLocalBegin(networkdm, X, INSERT_VALUES, localX)); in FormJacobian()
353 PetscCall(DMGlobalToLocalEnd(networkdm, X, INSERT_VALUES, localX)); in FormJacobian()
356 PetscCall(DMNetworkGetSubnetwork(networkdm, 0, &nv, &ne, &vtx, &edges)); in FormJacobian()
357 PetscCall(FormJacobian_Subnet(networkdm, localX, J, Jpre, nv, ne, vtx, edges, appctx)); in FormJacobian()
360 PetscCall(DMNetworkGetSubnetwork(networkdm, 1, &nv, &ne, &vtx, &edges)); in FormJacobian()
361 PetscCall(FormJacobian_Subnet(networkdm, localX, J, Jpre, nv, ne, vtx, edges, appctx)); in FormJacobian()
363 PetscCall(DMRestoreLocalVector(networkdm, &localX)); in FormJacobian()
370 PetscErrorCode SetInitialValues_Subnet(DM networkdm, Vec localX, PetscInt nv, PetscInt ne, const Pe… in SetInitialValues_Subnet() argument
383 PetscCall(DMNetworkIsGhostVertex(networkdm, vtx[i], &ghostvtex)); in SetInitialValues_Subnet()
386 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vtx[i], ALL_COMPONENTS, &offset)); in SetInitialValues_Subnet()
387 PetscCall(DMNetworkGetNumComponents(networkdm, vtx[i], &numComps)); in SetInitialValues_Subnet()
389 PetscCall(DMNetworkGetComponent(networkdm, vtx[i], j, &key, &component, NULL)); in SetInitialValues_Subnet()
406 PetscErrorCode SetInitialValues(DM networkdm, Vec X, void *appctx) in SetInitialValues() argument
413 PetscCall(DMGetLocalVector(networkdm, &localX)); in SetInitialValues()
416 PetscCall(DMGlobalToLocalBegin(networkdm, X, INSERT_VALUES, localX)); in SetInitialValues()
417 PetscCall(DMGlobalToLocalEnd(networkdm, X, INSERT_VALUES, localX)); in SetInitialValues()
420 PetscCall(DMNetworkGetSubnetwork(networkdm, 0, &nv, &ne, &vtx, &edges)); in SetInitialValues()
421 PetscCall(SetInitialValues_Subnet(networkdm, localX, nv, ne, vtx, edges, appctx)); in SetInitialValues()
424 PetscCall(DMNetworkGetSubnetwork(networkdm, 1, &nv, &ne, &vtx, &edges)); in SetInitialValues()
425 PetscCall(SetInitialValues_Subnet(networkdm, localX, nv, ne, vtx, edges, appctx)); in SetInitialValues()
427 PetscCall(DMLocalToGlobalBegin(networkdm, localX, ADD_VALUES, X)); in SetInitialValues()
428 PetscCall(DMLocalToGlobalEnd(networkdm, localX, ADD_VALUES, X)); in SetInitialValues()
429 PetscCall(DMRestoreLocalVector(networkdm, &localX)); in SetInitialValues()
439 DM networkdm; in main() local
458 PetscCall(DMNetworkCreate(PETSC_COMM_WORLD, &networkdm)); in main()
461 …PetscCall(DMNetworkRegisterComponent(networkdm, "branchstruct", sizeof(struct _p_EDGE_Power), &com… in main()
462 …PetscCall(DMNetworkRegisterComponent(networkdm, "busstruct", sizeof(struct _p_VERTEX_Power), &comp… in main()
463 …PetscCall(DMNetworkRegisterComponent(networkdm, "genstruct", sizeof(struct _p_GEN), &componentkey[… in main()
464 …PetscCall(DMNetworkRegisterComponent(networkdm, "loadstruct", sizeof(struct _p_LOAD), &componentke… in main()
499 PetscCall(DMNetworkSetNumSubNetworks(networkdm, PETSC_DECIDE, nsubnet)); in main()
500 PetscCall(DMNetworkAddSubnetwork(networkdm, "", numEdges1, edgelist1, NULL)); in main()
501 PetscCall(DMNetworkAddSubnetwork(networkdm, "", numEdges2, edgelist2, NULL)); in main()
504 PetscCall(DMNetworkLayoutSetUp(networkdm)); in main()
512 PetscCall(DMNetworkGetSubnetwork(networkdm, 0, &nv, &ne, &vtx, &edges)); in main()
514 …for (i = 0; i < ne; i++) PetscCall(DMNetworkAddComponent(networkdm, edges[i], componentkey[0], &pf… in main()
517 PetscCall(DMNetworkAddComponent(networkdm, vtx[i], componentkey[1], &pfdata1->bus[i], 2)); in main()
519 …for (j = 0; j < pfdata1->bus[i].ngen; j++) PetscCall(DMNetworkAddComponent(networkdm, vtx[i], comp… in main()
522 …for (j = 0; j < pfdata1->bus[i].nload; j++) PetscCall(DMNetworkAddComponent(networkdm, vtx[i], com… in main()
530 PetscCall(DMNetworkGetSubnetwork(networkdm, 1, &nv, &ne, &vtx, &edges)); in main()
532 …for (i = 0; i < ne; i++) PetscCall(DMNetworkAddComponent(networkdm, edges[i], componentkey[0], &pf… in main()
535 PetscCall(DMNetworkAddComponent(networkdm, vtx[i], componentkey[1], &pfdata2->bus[i], 2)); in main()
537 …for (j = 0; j < pfdata2->bus[i].ngen; j++) PetscCall(DMNetworkAddComponent(networkdm, vtx[i], comp… in main()
540 …for (j = 0; j < pfdata2->bus[i].nload; j++) PetscCall(DMNetworkAddComponent(networkdm, vtx[i], com… in main()
546 PetscCall(DMSetUp(networkdm)); in main()
568 PetscCall(DMNetworkDistribute(&networkdm, 0)); in main()
575 PetscCall(DMCreateGlobalVector(networkdm, &X)); in main()
578 PetscCall(DMCreateMatrix(networkdm, &J)); in main()
581 PetscCall(SetInitialValues(networkdm, X, &User)); in main()
585 PetscCall(SNESSetDM(snes, networkdm)); in main()
598 PetscCall(DMDestroy(&networkdm)); in main()