Lines Matching refs:P
104 Mat A, P, C; in testPTAPRectangular() local
114 PetscCall(MatCreateSeqAIJ(PETSC_COMM_WORLD, rows, cols, 5, NULL, &P)); in testPTAPRectangular()
115 PetscCall(MatSetValue(P, 0, 0, 1.0, INSERT_VALUES)); in testPTAPRectangular()
116 PetscCall(MatSetValue(P, 0, 1, 2.0, INSERT_VALUES)); in testPTAPRectangular()
117 PetscCall(MatSetValue(P, 0, 2, 0.0, INSERT_VALUES)); in testPTAPRectangular()
119 PetscCall(MatSetValue(P, 0, 3, -1.0, INSERT_VALUES)); in testPTAPRectangular()
121 PetscCall(MatSetValue(P, 1, 0, 0.0, INSERT_VALUES)); in testPTAPRectangular()
122 PetscCall(MatSetValue(P, 1, 1, -1.0, INSERT_VALUES)); in testPTAPRectangular()
123 PetscCall(MatSetValue(P, 1, 2, 1.0, INSERT_VALUES)); in testPTAPRectangular()
125 PetscCall(MatSetValue(P, 2, 0, 3.0, INSERT_VALUES)); in testPTAPRectangular()
126 PetscCall(MatSetValue(P, 2, 1, 0.0, INSERT_VALUES)); in testPTAPRectangular()
127 PetscCall(MatSetValue(P, 2, 2, -3.0, INSERT_VALUES)); in testPTAPRectangular()
129 PetscCall(MatAssemblyBegin(P, MAT_FINAL_ASSEMBLY)); in testPTAPRectangular()
130 PetscCall(MatAssemblyEnd(P, MAT_FINAL_ASSEMBLY)); in testPTAPRectangular()
133 PetscCall(MatPtAP(A, P, MAT_INITIAL_MATRIX, 1.0, &C)); in testPTAPRectangular()
170 PetscCall(MatDestroy(&P)); in testPTAPRectangular()