101   PetscInt     M, rows[2], cols[2], *nz;  in PipeCreateJacobian()  local
119   PetscCall(MatGetSize(Jpipe[0], &M, NULL));  in PipeCreateJacobian()
120   PetscCall(PetscCalloc2(M, &nz, 4, &aa));  in PipeCreateJacobian()
123   PetscCall(MatSetSizes(Jpipe[1], PETSC_DECIDE, PETSC_DECIDE, M, 2));  in PipeCreateJacobian()
139   PetscCall(MatSetSizes(Jpipe[2], PETSC_DECIDE, PETSC_DECIDE, M, 2));  in PipeCreateJacobian()
144   nz[M - 2] = 2;  in PipeCreateJacobian()
145   nz[M - 1] = 2;  in PipeCreateJacobian()
146   rows[0]   = M - 2;  in PipeCreateJacobian()
147   rows[1]   = M - 1;  in PipeCreateJacobian()
191   PetscInt        nedges, e, i, M, N, *rows, *cols;  in JunctionCreateJacobian()  local
205   PetscCall(DMNetworkGetComponent(dm, v, -1, NULL, NULL, &M));  in JunctionCreateJacobian()
206   PetscCheck(M == 2, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "%" PetscInt_FMT " != 2", M);  in JunctionCreateJacobian()
207   PetscCall(PetscMalloc3(M, &rows, M, &cols, M * M, &zeros));  in JunctionCreateJacobian()
208   PetscCall(PetscArrayzero(zeros, M * M));  in JunctionCreateJacobian()
209   for (i = 0; i < M; i++) rows[i] = i;  in JunctionCreateJacobian()
229       PetscCall(MatSetSizes(Jv[2 * e + 1], PETSC_DECIDE, PETSC_DECIDE, M, N));  in JunctionCreateJacobian()
248 …PetscCall(MatSetSizes(Jv[2 * e + 2], PETSC_DECIDE, PETSC_DECIDE, M, M)); /* empty matrix, sizes ca…  in JunctionCreateJacobian()