9   Mat         A, B, C, D;  in main()  local
27   PetscCall(MatCreate(PETSC_COMM_WORLD, &A));  in main()
28   PetscCall(MatSetSizes(A, PETSC_DECIDE, PETSC_DECIDE, 3, 3));  in main()
29   PetscCall(MatSetType(A, MATAIJ));  in main()
30   PetscCall(MatSetFromOptions(A));  in main()
31   PetscCall(MatSetUp(A));  in main()
32   PetscCall(MatSetOption(A, MAT_IGNORE_ZERO_ENTRIES, PETSC_TRUE));  in main()
34   if (rank == 0) PetscCall(MatSetValues(A, 3, ij, 3, ij, a, ADD_VALUES));  in main()
35   PetscCall(MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY));  in main()
36   PetscCall(MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY));  in main()
39   PetscCall(MatTranspose(A, MAT_INITIAL_MATRIX, &B));        /* B = A^T */  in main()
40   PetscCall(MatMatMult(B, A, MAT_INITIAL_MATRIX, fill, &C)); /* C = B*A */  in main()
41   PetscCall(MatMatMult(B, A, MAT_REUSE_MATRIX, fill, &C));   /* recompute C=B*A */  in main()
43   PetscCall(MatMatMultEqual(B, A, C, 10, &isequal));  in main()
46   PetscCall(MatMatMult(C, A, MAT_INITIAL_MATRIX, fill, &D)); /* D = C*A = (A^T*A)*A */  in main()
47   PetscCall(MatMatMult(C, A, MAT_REUSE_MATRIX, fill, &D));  in main()
48   PetscCall(MatMatMultEqual(C, A, D, 10, &isequal));  in main()
56   PetscCall(MatDuplicate(A, MAT_COPY_VALUES, &B));        /* B = A */  in main()
57   PetscCall(MatPtAP(A, B, MAT_INITIAL_MATRIX, fill, &C)); /* C = B^T*A*B */  in main()
58   PetscCall(MatPtAPMultEqual(A, B, C, 10, &isequal));  in main()
62   PetscCall(MatPtAP(A, B, MAT_REUSE_MATRIX, fill, &C));  in main()
63   PetscCall(MatPtAPMultEqual(A, B, C, 10, &isequal));  in main()
71     PetscCall(MatConvert(A, MATDENSE, MAT_INITIAL_MATRIX, &Adense));  in main()
84     PetscCall(MatMatTransposeMult(A, A, MAT_INITIAL_MATRIX, fill, &D)); /* D = A*A^T */  in main()
85     PetscCall(MatScale(A, 2.0));  in main()
86     PetscCall(MatMatTransposeMult(A, A, MAT_REUSE_MATRIX, fill, &D));  in main()
87     PetscCall(MatMatTransposeMultEqual(A, A, D, 10, &isequal));  in main()
91   PetscCall(MatDestroy(&A));  in main()
104   Mat       A, P, C;  in testPTAPRectangular()  local
108   PetscCall(MatCreateSeqAIJ(PETSC_COMM_WORLD, rows, rows, 1, NULL, &A));  in testPTAPRectangular()
109   for (i = 0; i < rows; i++) PetscCall(MatSetValue(A, i, i, 1.0, INSERT_VALUES));  in testPTAPRectangular()
110   PetscCall(MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY));  in testPTAPRectangular()
111   PetscCall(MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY));  in testPTAPRectangular()
133   PetscCall(MatPtAP(A, P, MAT_INITIAL_MATRIX, 1.0, &C));  in testPTAPRectangular()
169   PetscCall(MatDestroy(&A));  in testPTAPRectangular()