190   Mat            *Jv;  in JunctionCreateJacobian()  local
202   PetscCall(PetscCalloc1(2 * nedges + 1, &Jv));  in JunctionCreateJacobian()
218         Jv[2 * e + 1] = Jin[0];  in JunctionCreateJacobian()
220         Jv[2 * e + 1] = Jin[1];  in JunctionCreateJacobian()
222       Jv[2 * e + 2] = Jin[2];  in JunctionCreateJacobian()
223       PetscCall(PetscObjectReference((PetscObject)Jv[2 * e + 1]));  in JunctionCreateJacobian()
224       PetscCall(PetscObjectReference((PetscObject)Jv[2 * e + 2]));  in JunctionCreateJacobian()
227       PetscCall(MatCreate(PETSC_COMM_SELF, &Jv[2 * e + 1]));  in JunctionCreateJacobian()
229       PetscCall(MatSetSizes(Jv[2 * e + 1], PETSC_DECIDE, PETSC_DECIDE, M, N));  in JunctionCreateJacobian()
230       PetscCall(MatSetFromOptions(Jv[2 * e + 1]));  in JunctionCreateJacobian()
231       PetscCall(MatSetOption(Jv[2 * e + 1], MAT_STRUCTURE_ONLY, PETSC_TRUE));  in JunctionCreateJacobian()
232       PetscCall(MatSeqAIJSetPreallocation(Jv[2 * e + 1], 2, NULL));  in JunctionCreateJacobian()
240         PetscCall(MatSetValues(Jv[2 * e + 1], 2, rows, 2, cols, zeros, INSERT_VALUES));  in JunctionCreateJacobian()
242       PetscCall(MatAssemblyBegin(Jv[2 * e + 1], MAT_FINAL_ASSEMBLY));  in JunctionCreateJacobian()
243       PetscCall(MatAssemblyEnd(Jv[2 * e + 1], MAT_FINAL_ASSEMBLY));  in JunctionCreateJacobian()
247       PetscCall(MatCreate(PETSC_COMM_SELF, &Jv[2 * e + 2]));  in JunctionCreateJacobian()
248 …PetscCall(MatSetSizes(Jv[2 * e + 2], PETSC_DECIDE, PETSC_DECIDE, M, M)); /* empty matrix, sizes ca…  in JunctionCreateJacobian()
249       PetscCall(MatSetFromOptions(Jv[2 * e + 2]));  in JunctionCreateJacobian()
250       PetscCall(MatSetOption(Jv[2 * e + 2], MAT_STRUCTURE_ONLY, PETSC_TRUE));  in JunctionCreateJacobian()
251       PetscCall(MatSeqAIJSetPreallocation(Jv[2 * e + 2], 1, NULL));  in JunctionCreateJacobian()
252       PetscCall(MatAssemblyBegin(Jv[2 * e + 2], MAT_FINAL_ASSEMBLY));  in JunctionCreateJacobian()
253       PetscCall(MatAssemblyEnd(Jv[2 * e + 2], MAT_FINAL_ASSEMBLY));  in JunctionCreateJacobian()
258   *J = Jv;  in JunctionCreateJacobian()
264   Mat            *Jv = junc->jacobian;  in JunctionDestroyJacobian()  local
269   if (!Jv) PetscFunctionReturn(PETSC_SUCCESS);  in JunctionDestroyJacobian()
273     PetscCall(MatDestroy(&Jv[2 * e + 1]));  in JunctionDestroyJacobian()
274     PetscCall(MatDestroy(&Jv[2 * e + 2]));  in JunctionDestroyJacobian()
276   PetscCall(PetscFree(Jv));  in JunctionDestroyJacobian()