341       Mat      sA, sS, dA, dS, At, St;  in main()  local
350       PetscCall(MatCreateSubMatrix(S, r, c, MAT_INITIAL_MATRIX, &sS));  in main()
351       PetscCall(MatMultAddEqual(sA, sS, 10, &flg));  in main()
353       PetscCall(MatMultTransposeAddEqual(sA, sS, 10, &flg));  in main()
356       PetscCall(MatConvert(sS, MATDENSE, MAT_INITIAL_MATRIX, &dS));  in main()
361       PetscCall(MatDestroy(&sS));  in main()
367       PetscCall(MatCreateSubMatrix(St, c, r, MAT_INITIAL_MATRIX, &sS));  in main()
368       PetscCall(MatMultAddEqual(sA, sS, 10, &flg));  in main()
370       PetscCall(MatMultTransposeAddEqual(sA, sS, 10, &flg));  in main()
373       PetscCall(MatConvert(sS, MATDENSE, MAT_INITIAL_MATRIX, &dS));  in main()
378       PetscCall(MatDestroy(&sS));  in main()