159   Mat               A, B, Q, Acondensed, Bcondensed, BT, J, AplusJ, QInv, D, AplusD, JplusD, U;  in main()  local
243   PetscCall(MatDuplicate(Q, MAT_DO_NOT_COPY_VALUES, &QInv));  in main()
244   PetscCall(MatDiagonalSet(QInv, Qdiag, INSERT_VALUES));  in main()
245   PetscCall(MatAssemblyBegin(QInv, MAT_FINAL_ASSEMBLY));  in main()
246   PetscCall(MatAssemblyEnd(QInv, MAT_FINAL_ASSEMBLY));  in main()
250   PetscCall(MatMatMatMult(Bcondensed, QInv, BT, MAT_INITIAL_MATRIX, PETSC_CURRENT, &J));  in main()
263   PetscCall(MatDiagonalSet(QInv, Qdiag, INSERT_VALUES));  in main()
264   PetscCall(MatAssemblyBegin(QInv, MAT_FINAL_ASSEMBLY));  in main()
265   PetscCall(MatAssemblyEnd(QInv, MAT_FINAL_ASSEMBLY));  in main()
267   PetscCall(MatMatMult(Bcondensed, QInv, MAT_INITIAL_MATRIX, PETSC_CURRENT, &U));  in main()
332   PetscCall(MatDestroy(&QInv));  in main()