52   Vec         W1, W2, W3, diff;  in TestMatrix()  local
71   PetscCall(VecDuplicate(X, &W2));  in TestMatrix()
83   PetscCall(MatMultTranspose(A, W1, W2));  in TestMatrix()
84   PetscCall(VecView(W2, viewer));  in TestMatrix()
87   PetscCall(MatMultHermitianTranspose(A, W1, W2));  in TestMatrix()
88   PetscCall(VecConjugate(W2));  in TestMatrix()
89   PetscCall(VecView(W2, viewer));  in TestMatrix()
94   PetscCall(MatMultAdd(A, W1, W1, W2));  in TestMatrix()
95   PetscCall(VecView(W2, viewer));  in TestMatrix()
96   PetscCall(VecAXPY(W2, -1.0, diff));  in TestMatrix()
97   PetscCall(VecNorm(W2, NORM_2, &nrm));  in TestMatrix()
102   PetscCall(VecSet(W2, -1.0));  in TestMatrix()
103   PetscCall(MatMultAdd(A, W1, W2, W2));  in TestMatrix()
104   PetscCall(VecView(W2, viewer));  in TestMatrix()
105   PetscCall(VecAXPY(W2, -1.0, diff));  in TestMatrix()
106   PetscCall(VecNorm(W2, NORM_2, &nrm));  in TestMatrix()
116   PetscCall(MatMultTransposeAdd(A, W1, W1, W2));  in TestMatrix()
117   PetscCall(VecView(W2, viewer));  in TestMatrix()
118   PetscCall(VecAXPY(W2, -1.0, diff));  in TestMatrix()
119   PetscCall(VecNorm(W2, NORM_2, &nrm));  in TestMatrix()
124   PetscCall(VecSet(W2, -1.0));  in TestMatrix()
125   PetscCall(MatMultTransposeAdd(A, W1, W2, W2));  in TestMatrix()
126   PetscCall(VecView(W2, viewer));  in TestMatrix()
127   PetscCall(VecAXPY(W2, -1.0, diff));  in TestMatrix()
128   PetscCall(VecNorm(W2, NORM_2, &nrm));  in TestMatrix()
159   PetscCall(MatGetDiagonal(A, W2));  in TestMatrix()
160   PetscCall(VecView(W2, viewer));  in TestMatrix()
162   PetscCall(VecAXPY(diff, -1.0, W2));  in TestMatrix()
180   PetscCall(VecDestroy(&W2));  in TestMatrix()