341       Mat      sA, sS, dA, dS, At, St;  in main()  local
349       PetscCall(MatCreateSubMatrix(A, r, c, MAT_INITIAL_MATRIX, &sA));  in main()
351       PetscCall(MatMultAddEqual(sA, sS, 10, &flg));  in main()
353       PetscCall(MatMultTransposeAddEqual(sA, sS, 10, &flg));  in main()
355       PetscCall(MatConvert(sA, MATDENSE, MAT_INITIAL_MATRIX, &dA));  in main()
360       PetscCall(MatDestroy(&sA));  in main()
366       PetscCall(MatCreateSubMatrix(At, c, r, MAT_INITIAL_MATRIX, &sA));  in main()
368       PetscCall(MatMultAddEqual(sA, sS, 10, &flg));  in main()
370       PetscCall(MatMultTransposeAddEqual(sA, sS, 10, &flg));  in main()
372       PetscCall(MatConvert(sA, MATDENSE, MAT_INITIAL_MATRIX, &dA));  in main()
377       PetscCall(MatDestroy(&sA));  in main()