impls/network/network.c

143 - edgelist - list of edges for this subnetwork, this is a one dimensional array with pairs of entries being the two vertices (in global numbering
152   There is no copy involved in this operation, only the pointer is referenced. The `edgelist` should
174   edgelist = [0 2 | 1 2]
178   edgelist = [3 5 | 4 5]
194   edgelist = [0 2 | 2 1 | 1 3]
198   edgelist = NULL
204   edgelist = NULL
207   edgelist = [0 3 | 3 2 | 2 1]
212 PetscErrorCode DMNetworkAddSubnetwork(DM dm, const char *name, PetscInt ne, PetscInt edgelist[], PetscInt *netnum)
219   for (i = 0; i < ne; i++) PetscCheck(edgelist[2 * i] != edgelist[2 * i + 1], PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Edge %" PetscInt_FMT " has the same vertex %" PetscInt_FMT " at each endpoint", i, edgelist[2 * i]);
222   if (ne) nvtx_min = nvtx_max = edgelist[0];
224     nvtx_min = PetscMin(nvtx_min, edgelist[i]);
225     nvtx_max = PetscMax(nvtx_max, edgelist[i]);
227     nvtx_min = PetscMin(nvtx_min, edgelist[i]);
228     nvtx_max = PetscMax(nvtx_max, edgelist[i]);
238     PetscCall(PetscBTSet(table, edgelist[i++] - nvtx_min));
239     PetscCall(PetscBTSet(table, edgelist[i++] - nvtx_min));
246   /* Get global total Nvtx = max(edgelist[])+1 for this subnet */
258   network->cloneshared->subnet[i].edgelist = edgelist;
518   GetEdgelist_Coupling - Get an integrated edgelist for dmplex from user-provided subnet[].edgelist when subnets are coupled by shared vertices
524 +  edges - the integrated edgelist for dmplex
594       i              = network->cloneshared->subnet[net].edgelist[2 * j] + (network->cloneshared->subnet[net].vStart - vrange);
598       i                  = network->cloneshared->subnet[net].edgelist[2 * j + 1] + (network->cloneshared->subnet[net].vStart - vrange);
677   /* Create LOCAL edgelist in global vertex ordering for the network by concatenating local input edgelists of the subnetworks */
688         edges[2 * ctr]     = network->cloneshared->subnet[i].vStart + network->cloneshared->subnet[i].edgelist[2 * j];
689         edges[2 * ctr + 1] = network->cloneshared->subnet[i].vStart + network->cloneshared->subnet[i].edgelist[2 * j + 1];
751         vfrom                                           = network->cloneshared->subnet[i].edgelist[2 * ctr]; /* =subnet[i].idx, Global index! */
762         vto                                           = network->cloneshared->subnet[i].edgelist[2 * ctr + 1]; /* =subnet[i].idx, Global index! */
866   PetscCheck(Nsvtx <= 2 * nsubnet, PETSC_COMM_SELF, PETSC_ERR_SUP, "allocate more space for coupling edgelist");