Lines Matching refs:networkdm
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()
155 PetscCall(DMNetworkGetLocalVecOffset(networkdm, v, ALL_COMPONENTS, &varoffset)); in WASHIFunction()
168 PetscCall(DMNetworkGetEdgeRange(networkdm, &eStart, &eEnd)); in WASHIFunction()
170 PetscCall(DMNetworkGetComponent(networkdm, e, 0, &type, (void **)&pipe, NULL)); in WASHIFunction()
171 PetscCall(DMNetworkGetLocalVecOffset(networkdm, e, ALL_COMPONENTS, &varoffset)); in WASHIFunction()
186 PetscCall(DMNetworkGetConnectedVertices(networkdm, e, &cone)); in WASHIFunction()
189 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vfrom, ALL_COMPONENTS, &offsetfrom)); in WASHIFunction()
190 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vto, ALL_COMPONENTS, &offsetto)); in WASHIFunction()
193 PetscCall(DMNetworkGetComponent(networkdm, vfrom, 0, &type, (void **)&junction, NULL)); in WASHIFunction()
202 PetscCall(DMNetworkGetComponent(networkdm, vto, 0, &type, (void **)&junction, NULL)); in WASHIFunction()
216 PetscCall(DMLocalToGlobalBegin(networkdm, localF, ADD_VALUES, F)); in WASHIFunction()
217 PetscCall(DMLocalToGlobalEnd(networkdm, localF, ADD_VALUES, F)); in WASHIFunction()
218 PetscCall(DMRestoreLocalVector(networkdm, &localF)); in WASHIFunction()
226 PetscErrorCode WASHSetInitialSolution(DM networkdm, Vec X, Wash wash) in WASHSetInitialSolution() argument
240 PetscCall(DMGetLocalVector(networkdm, &localX)); in WASHSetInitialSolution()
244 PetscCall(DMNetworkGetEdgeRange(networkdm, &eStart, &eEnd)); in WASHSetInitialSolution()
246 PetscCall(DMNetworkGetLocalVecOffset(networkdm, e, ALL_COMPONENTS, &varoffset)); in WASHSetInitialSolution()
247 PetscCall(DMNetworkGetComponent(networkdm, e, 0, &type, (void **)&pipe, NULL)); in WASHSetInitialSolution()
258 PetscCall(DMNetworkGetConnectedVertices(networkdm, e, &cone)); in WASHSetInitialSolution()
261 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vfrom, ALL_COMPONENTS, &offsetfrom)); in WASHSetInitialSolution()
262 PetscCall(DMNetworkGetLocalVecOffset(networkdm, vto, ALL_COMPONENTS, &offsetto)); in WASHSetInitialSolution()
265 PetscCall(DMNetworkGetComponent(networkdm, vfrom, 0, &vkey, (void **)&junction, NULL)); in WASHSetInitialSolution()
269 PetscCall(DMNetworkGetComponent(networkdm, vto, 0, &vkey, (void **)&junction, NULL)); in WASHSetInitialSolution()
275 PetscCall(DMLocalToGlobalBegin(networkdm, localX, ADD_VALUES, X)); in WASHSetInitialSolution()
276 PetscCall(DMLocalToGlobalEnd(networkdm, localX, ADD_VALUES, X)); in WASHSetInitialSolution()
277 PetscCall(DMRestoreLocalVector(networkdm, &localX)); in WASHSetInitialSolution()
298 PetscErrorCode PipesView(DM networkdm, PetscInt KeyPipe, Vec X) in PipesView() argument
307 PetscCall(PetscObjectGetComm((PetscObject)networkdm, &comm)); in PipesView()
317 …PetscCall(DMNetworkCreateIS(networkdm, numkeys, &KeyPipe, blocksize, numselectedvariable, selected… in PipesView()
340 …PetscCall(DMNetworkCreateIS(networkdm, numkeys, &KeyPipe, blocksize, numselectedvariable, selected… in PipesView()
354 PetscCall(DMNetworkCreateIS(networkdm, numkeys, &KeyPipe, blocksize, NULL, NULL, &isfrom)); in PipesView()
356 PetscCall(DMNetworkCreateIS(networkdm, numkeys, &KeyPipe, NULL, NULL, NULL, &isfrom)); in PipesView()
380 PetscErrorCode ISJunctionsView(DM networkdm, PetscInt KeyJunc) in ISJunctionsView() argument
388 PetscCall(PetscObjectGetComm((PetscObject)networkdm, &comm)); in ISJunctionsView()
392 PetscCall(DMNetworkCreateIS(networkdm, numkeys, &KeyJunc, NULL, NULL, NULL, &isfrom)); in ISJunctionsView()
398 PetscCall(DMNetworkCreateLocalIS(networkdm, numkeys, &KeyJunc, NULL, NULL, NULL, &isfrom)); in ISJunctionsView()
618 DM networkdm; in main() local
642 PetscCall(DMNetworkCreate(PETSC_COMM_WORLD, &networkdm)); in main()
643 PetscCall(PetscObjectGetComm((PetscObject)networkdm, &comm)); in main()
647 if (size == 1 && monipipes) PetscCall(DMNetworkMonitorCreate(networkdm, &monitor)); in main()
650 …PetscCall(DMNetworkRegisterComponent(networkdm, "junctionstruct", sizeof(struct _p_Junction), &Key… in main()
651 PetscCall(DMNetworkRegisterComponent(networkdm, "pipestruct", sizeof(struct _p_Pipe), &KeyPipe)); in main()
662 PetscCall(DMNetworkSetNumSubNetworks(networkdm, PETSC_DECIDE, 1)); in main()
663 PetscCall(DMNetworkAddSubnetwork(networkdm, NULL, nedges, edgelist, NULL)); in main()
665 PetscCall(DMNetworkLayoutSetUp(networkdm)); in main()
667 PetscCall(DMNetworkGetEdgeRange(networkdm, &eStart, &eEnd)); in main()
668 PetscCall(DMNetworkGetVertexRange(networkdm, &vStart, &vEnd)); in main()
679 …PetscCall(DMNetworkAddComponent(networkdm, e, KeyPipe, &pipes[e - eStart], 2 * pipes[e - eStart].n… in main()
689 …for (v = vStart; v < vEnd; v++) PetscCall(DMNetworkAddComponent(networkdm, v, KeyJunction, &juncti… in main()
694 PetscCall(DMNetworkGetPlex(networkdm, &plexdm)); in main()
701 PetscCall(DMSetUp(networkdm)); in main()
704 PetscCall(DMView(networkdm, PETSC_VIEWER_STDOUT_WORLD)); in main()
709 PetscCall(DMNetworkGetConnectedVertices(networkdm, e, &cone)); in main()
711 PetscCall(DMNetworkGetComponent(networkdm, cone[0], 0, &vkey, (void **)&junction, NULL)); in main()
715 PetscCall(DMNetworkGetComponent(networkdm, cone[1], 0, &vkey, (void **)&junction, NULL)); in main()
722 PetscCall(DMNetworkDistribute(&networkdm, 0)); in main()
725 PetscCall(DMView(networkdm, PETSC_VIEWER_STDOUT_WORLD)); in main()
729 PetscCall(DMCreateGlobalVector(networkdm, &X)); in main()
730 PetscCall(DMCreateLocalVector(networkdm, &wash->localX)); in main()
731 PetscCall(DMCreateLocalVector(networkdm, &wash->localXdot)); in main()
735 PetscCall(DMNetworkGetVertexRange(networkdm, &vStart, &vEnd)); in main()
738 PetscCall(DMNetworkHasJacobian(networkdm, userJac, userJac)); in main()
739 PetscCall(DMNetworkGetEdgeRange(networkdm, &eStart, &eEnd)); in main()
741 PetscCall(DMNetworkGetComponent(networkdm, e, 0, &type, (void **)&pipe, NULL)); in main()
754 PetscCall(DMNetworkEdgeSetMatrix(networkdm, e, J)); in main()
759 PetscCall(DMNetworkGetVertexRange(networkdm, &vStart, &vEnd)); in main()
762 PetscCall(JunctionCreateJacobian(networkdm, v, NULL, &J)); in main()
763 PetscCall(DMNetworkVertexSetMatrix(networkdm, v, J)); in main()
765 PetscCall(DMNetworkGetComponent(networkdm, v, 0, &vkey, (void **)&junction, NULL)); in main()
774 PetscCall(TSSetDM(ts, (DM)networkdm)); in main()
784 PetscCall(WASHSetInitialSolution(networkdm, X, wash)); in main()
808 if (viewpipes) PetscCall(PipesView(networkdm, KeyPipe, X)); in main()
813 if (viewjuncs) PetscCall(ISJunctionsView(networkdm, KeyJunction)); in main()
823 PetscCall(DMNetworkGetEdgeRange(networkdm, &eStart, &eEnd)); in main()
825 PetscCall(DMNetworkGetComponent(networkdm, i, 0, &key, (void **)&pipe, NULL)); in main()
830 PetscCall(DMNetworkGetComponent(networkdm, v, 0, &vkey, (void **)&junction, NULL)); in main()
831 PetscCall(JunctionDestroyJacobian(networkdm, v, junction)); in main()
836 PetscCall(DMDestroy(&networkdm)); in main()