100   Mat         *Jpipe;  in PipeCreateJacobian()  local
111   PetscCall(PetscMalloc1(3, &Jpipe));  in PipeCreateJacobian()
115   PetscCall(DMCreateMatrix(pipe->da, &Jpipe[0]));  in PipeCreateJacobian()
119   PetscCall(MatGetSize(Jpipe[0], &M, NULL));  in PipeCreateJacobian()
122   PetscCall(MatCreate(PETSC_COMM_SELF, &Jpipe[1]));  in PipeCreateJacobian()
123   PetscCall(MatSetSizes(Jpipe[1], PETSC_DECIDE, PETSC_DECIDE, M, 2));  in PipeCreateJacobian()
124   PetscCall(MatSetFromOptions(Jpipe[1]));  in PipeCreateJacobian()
125   PetscCall(MatSetOption(Jpipe[1], MAT_STRUCTURE_ONLY, PETSC_TRUE));  in PipeCreateJacobian()
132   PetscCall(MatSeqAIJSetPreallocation(Jpipe[1], 0, nz));  in PipeCreateJacobian()
133   PetscCall(MatSetValues(Jpipe[1], 2, rows, 2, cols, aa, INSERT_VALUES));  in PipeCreateJacobian()
134   PetscCall(MatAssemblyBegin(Jpipe[1], MAT_FINAL_ASSEMBLY));  in PipeCreateJacobian()
135   PetscCall(MatAssemblyEnd(Jpipe[1], MAT_FINAL_ASSEMBLY));  in PipeCreateJacobian()
138   PetscCall(MatCreate(PETSC_COMM_SELF, &Jpipe[2]));  in PipeCreateJacobian()
139   PetscCall(MatSetSizes(Jpipe[2], PETSC_DECIDE, PETSC_DECIDE, M, 2));  in PipeCreateJacobian()
140   PetscCall(MatSetFromOptions(Jpipe[2]));  in PipeCreateJacobian()
141   PetscCall(MatSetOption(Jpipe[2], MAT_STRUCTURE_ONLY, PETSC_TRUE));  in PipeCreateJacobian()
148   PetscCall(MatSeqAIJSetPreallocation(Jpipe[2], 0, nz));  in PipeCreateJacobian()
149   PetscCall(MatSetValues(Jpipe[2], 2, rows, 2, cols, aa, INSERT_VALUES));  in PipeCreateJacobian()
150   PetscCall(MatAssemblyBegin(Jpipe[2], MAT_FINAL_ASSEMBLY));  in PipeCreateJacobian()
151   PetscCall(MatAssemblyEnd(Jpipe[2], MAT_FINAL_ASSEMBLY));  in PipeCreateJacobian()
155   *J             = Jpipe;  in PipeCreateJacobian()
156   pipe->jacobian = Jpipe;  in PipeCreateJacobian()
162   Mat     *Jpipe = pipe->jacobian;  in PipeDestroyJacobian()  local
166   if (Jpipe) {  in PipeDestroyJacobian()
167     for (i = 0; i < 3; i++) PetscCall(MatDestroy(&Jpipe[i]));  in PipeDestroyJacobian()
169   PetscCall(PetscFree(Jpipe));  in PipeDestroyJacobian()