19 PetscErrorCode WashNetworkDistribute(MPI_Comm comm, Wash wash)  in WashNetworkDistribute()  argument
23   PetscInt   *edgelist = wash->edgelist, *nvtx = NULL, *vtxDone = NULL;  in WashNetworkDistribute()
30   numEdges    = wash->nedge;  in WashNetworkDistribute()
31   numVertices = wash->nvertex;  in WashNetworkDistribute()
37   wash->Nedge = numEdges;  in WashNetworkDistribute()
38   wash->nedge = nedges;  in WashNetworkDistribute()
52     vtype = wash->vtype;  in WashNetworkDistribute()
69   wash->edgelist = edgelist;  in WashNetworkDistribute()
92   wash->Nvertex = numVertices;  in WashNetworkDistribute()
93   wash->nvertex = nvertices;  in WashNetworkDistribute()
94   wash->vtype   = vtype;  in WashNetworkDistribute()
100   Wash               wash = (Wash)ctx;  in WASHIFunction()  local
117   localX    = wash->localX;  in WASHIFunction()
118   localXdot = wash->localXdot;  in WASHIFunction()
163       juncf[0] = juncx[0].h - wash->H0;  in WASHIFunction()
226 PetscErrorCode WASHSetInitialSolution(DM networkdm, Vec X, Wash wash)  in WASHSetInitialSolution()  argument
250     PetscCall(PipeComputeSteadyState(pipe, wash->Q0, wash->H0));  in WASHSetInitialSolution()
266     if (junction->type == RESERVOIR) (xarr + offsetfrom)[1] = wash->H0; /* 1st H */  in WASHSetInitialSolution()
270     if (junction->type == VALVE) (xarr + offsetto)[0] = wash->QL; /* last Q */  in WASHSetInitialSolution()
407 PetscErrorCode WashNetworkCleanUp(Wash wash)  in WashNetworkCleanUp()  argument
412   PetscCallMPI(MPI_Comm_rank(wash->comm, &rank));  in WashNetworkCleanUp()
413   PetscCall(PetscFree(wash->edgelist));  in WashNetworkCleanUp()
414   PetscCall(PetscFree(wash->vtype));  in WashNetworkCleanUp()
415   if (rank == 0) PetscCall(PetscFree2(wash->junction, wash->pipe));  in WashNetworkCleanUp()
423   Wash        wash = NULL;  in WashNetworkCreate()  local
433   PetscCall(PetscCalloc1(1, &wash));  in WashNetworkCreate()
434   wash->comm       = comm;  in WashNetworkCreate()
435   *wash_ptr        = wash;  in WashNetworkCreate()
436   wash->Q0         = 0.477432; /* RESERVOIR */  in WashNetworkCreate()
437   wash->H0         = 150.0;  in WashNetworkCreate()
438   wash->HL         = 143.488; /* VALVE */  in WashNetworkCreate()
439   wash->QL         = 0.0;  in WashNetworkCreate()
440   wash->nnodes_loc = 0;  in WashNetworkCreate()
459     wash->nedge   = npipes;  in WashNetworkCreate()
460     wash->nvertex = npipes + 1;  in WashNetworkCreate()
467       numVertices = wash->nvertex;  in WashNetworkCreate()
468       numEdges    = wash->nedge;  in WashNetworkCreate()
502     wash->nedge   = npipes;  in WashNetworkCreate()
503     wash->nvertex = npipes + 1;  in WashNetworkCreate()
507       numVertices = wash->nvertex;  in WashNetworkCreate()
508       numEdges    = wash->nedge;  in WashNetworkCreate()
543     wash->nedge   = npipes;  in WashNetworkCreate()
544     wash->nvertex = npipes + 1;  in WashNetworkCreate()
548       numVertices = wash->nvertex;  in WashNetworkCreate()
549       numEdges    = wash->nedge;  in WashNetworkCreate()
587     wash->vtype = vtype;  in WashNetworkCreate()
590   *wash_ptr      = wash;  in WashNetworkCreate()
591   wash->nedge    = numEdges;  in WashNetworkCreate()
592   wash->nvertex  = numVertices;  in WashNetworkCreate()
593   wash->edgelist = edgelist;  in WashNetworkCreate()
594   wash->junction = junctions;  in WashNetworkCreate()
595   wash->pipe     = pipes;  in WashNetworkCreate()
603     PetscCall(WashNetworkDistribute(comm, wash));  in WashNetworkCreate()
611   Wash              wash;  in main()  local
654   PetscCall(WashNetworkCreate(comm, pipesCase, &wash));  in main()
655   nedges    = wash->nedge;  in main()
656   edgelist  = wash->edgelist;  in main()
657   vtype     = wash->vtype;  in main()
658   junctions = wash->junction;  in main()
659   pipes     = wash->pipe;  in main()
719   PetscCall(WashNetworkCleanUp(wash));  in main()
730   PetscCall(DMCreateLocalVector(networkdm, &wash->localX));  in main()
731   PetscCall(DMCreateLocalVector(networkdm, &wash->localXdot));  in main()
743 …wash->nnodes_loc += pipe->nnodes;        /* local total number of nodes, will be used by PipesView…  in main()
775   PetscCall(TSSetIFunction(ts, NULL, WASHIFunction, wash));  in main()
784   PetscCall(WASHSetInitialSolution(networkdm, X, wash));  in main()
819   PetscCall(VecDestroy(&wash->localX));  in main()
820   PetscCall(VecDestroy(&wash->localXdot));  in main()
837   PetscCall(PetscFree(wash));  in main()