11   Mat             A, B, C, D;  in main()  local
103     PetscCall(MatDuplicate(C, MAT_COPY_VALUES, &D));  in main()
105     PetscCall(MatDestroy(&D));  in main()
114     PetscCall(MatTransposeMatMult(A, A, MAT_INITIAL_MATRIX, fill, &D)); /* D = A^T*A */  in main()
115     PetscCall(MatTransposeMatMult(A, A, MAT_REUSE_MATRIX, fill, &D));  in main()
116     PetscCall(MatTransposeMatMultEqual(A, A, D, 10, &equal));  in main()
120     PetscCall(MatDuplicate(D, MAT_COPY_VALUES, &C));  in main()
124     PetscCall(MatSetRandom(D, NULL));  in main()
125     PetscCall(MatMatMult(A, D, MAT_INITIAL_MATRIX, fill, &E));  in main()
126     PetscCall(MatMatMult(A, D, MAT_REUSE_MATRIX, fill, &E));  in main()
127     PetscCall(MatMatMultEqual(A, D, E, 10, &equal));  in main()
131     PetscCall(MatDestroy(&D));  in main()
151     PetscCall(MatTransposeMatMult(A, B, MAT_INITIAL_MATRIX, fill, &D));  in main()
152     PetscCall(MatTransposeMatMultEqual(A, B, D, 10, &equal));  in main()
155     PetscCall(MatDestroy(&D));  in main()
169     PetscCall(MatMatTransposeMult(A, B, MAT_INITIAL_MATRIX, fill, &D)); /* D = A*A^T */  in main()
172     PetscCall(MatDuplicate(D, MAT_COPY_VALUES, &C));  in main()
174     PetscCall(MatMatTransposeMult(A, B, MAT_REUSE_MATRIX, fill, &D));  in main()
175     PetscCall(MatAXPY(C, -1., D, SAME_NONZERO_PATTERN));  in main()
178     PetscCall(MatNorm(D, NORM_FROBENIUS, &scale));  in main()
179 …PetscCheck(diff <= PETSC_SMALL * scale, PetscObjectComm((PetscObject)D), PETSC_ERR_PLIB, "MatMatTr…  in main()
182     PetscCall(MatMatTransposeMultEqual(A, B, D, 10, &equal));  in main()
184     PetscCall(MatDestroy(&D));  in main()