2186   if (vflg && !network->Jv && nVertices) {  in DMNetworkHasJacobian()
2202     PetscCall(PetscCalloc1(2 * nedges_total + nVertices, &network->Jv));  in DMNetworkHasJacobian()
2263 …PetscCheck(network->Jv, PetscObjectComm((PetscObject)dm), PETSC_ERR_ORDER, "Must call DMNetworkHas…  in DMNetworkVertexSetMatrix()
2267     network->Jv[vptr[p - vStart]] = J[0]; /* Set Jacobian for this vertex */  in DMNetworkVertexSetMatrix()
2271     for (i = 1; i <= 2 * nedges; i++) network->Jv[vptr[p - vStart] + i] = J[i];  in DMNetworkVertexSetMatrix()
2569       if (network->Jv) {  in DMCreateMatrix_Network()
2570         Juser = network->Jv[vptr[v - vStart] + 2 * e + 1]; /* Jacobian(v,e) */  in DMCreateMatrix_Network()
2581       if (network->Jv) {  in DMCreateMatrix_Network()
2582         Juser = network->Jv[vptr[v - vStart] + 2 * e + 2]; /* Jacobian(v,vc) */  in DMCreateMatrix_Network()
2596       if (network->Jv) {  in DMCreateMatrix_Network()
2597         Juser = network->Jv[vptr[v - vStart]]; /* Jacobian(v,v) */  in DMCreateMatrix_Network()
2681       if (network->Jv) {  in DMCreateMatrix_Network()
2682         Juser = network->Jv[vptr[v - vStart] + 2 * e + 1]; /* Jacobian(v,e) */  in DMCreateMatrix_Network()
2693       if (network->Jv) {  in DMCreateMatrix_Network()
2694         Juser = network->Jv[vptr[v - vStart] + 2 * e + 2]; /* Jacobian(v,vc) */  in DMCreateMatrix_Network()
2702       if (network->Jv) {  in DMCreateMatrix_Network()
2703         Juser = network->Jv[vptr[v - vStart]]; /* Jacobian(v,v) */  in DMCreateMatrix_Network()
2751   if (network->Jv) {  in DMDestroy_Network()
2753     PetscCall(PetscFree(network->Jv));  in DMDestroy_Network()