3 static PetscErrorCode MatTransposeAXPY_Private(Mat Y, PetscScalar a, Mat X, MatStructure str, Mat T)  in MatTransposeAXPY_Private()  argument
10 …if (T == X) PetscCall(MatShellGetScalingShifts(T, &vshift, &vscale, (Vec *)MAT_SHELL_NOT_ALLOWED, …  in MatTransposeAXPY_Private()
18     if (T == X) {  in MatTransposeAXPY_Private()
24       PetscCall(MatTranspose(X, MAT_INITIAL_MATRIX, &F));  in MatTransposeAXPY_Private()
28     if (T == X) {  in MatTransposeAXPY_Private()
34       PetscCall(MatHermitianTranspose(X, MAT_INITIAL_MATRIX, &F));  in MatTransposeAXPY_Private()
43 static PetscErrorCode MatAXPY_BasicWithTypeCompare(Mat Y, PetscScalar a, Mat X, MatStructure str)  in MatAXPY_BasicWithTypeCompare()  argument
54       PetscCall((*f)(Y, a, X, str));  in MatAXPY_BasicWithTypeCompare()
61       PetscCall(PetscObjectTypeCompare((PetscObject)X, MATNEST, &flg));  in MatAXPY_BasicWithTypeCompare()
63         PetscCall(MatAXPY_Dense_Nest(Y, a, X));  in MatAXPY_BasicWithTypeCompare()
67     PetscCall(PetscObjectTypeCompareAny((PetscObject)X, &flg, MATSCALAPACK, MATELEMENTAL, ""));  in MatAXPY_BasicWithTypeCompare()
71       PetscCall(MatConvert(X, ((PetscObject)Y)->type_name, MAT_INITIAL_MATRIX, &C));  in MatAXPY_BasicWithTypeCompare()
77   PetscCall(MatAXPY_Basic(Y, a, X, str));  in MatAXPY_BasicWithTypeCompare()
96 PetscErrorCode MatAXPY(Mat Y, PetscScalar a, Mat X, MatStructure str)  in MatAXPY()  argument
105   PetscValidHeaderSpecific(X, MAT_CLASSID, 3);  in MatAXPY()
106   PetscCheckSameComm(Y, 1, X, 3);  in MatAXPY()
107   PetscCall(MatGetSize(X, &M1, &N1));  in MatAXPY()
109   PetscCall(MatGetLocalSize(X, &m1, &n1));  in MatAXPY()
114 …PetscCheck(X->assembled, PetscObjectComm((PetscObject)X), PETSC_ERR_ARG_WRONGSTATE, "Not for unass…  in MatAXPY()
116   if (Y == X) {  in MatAXPY()
120   PetscCall(PetscObjectObjectTypeCompare((PetscObject)X, (PetscObject)Y, &sametype));  in MatAXPY()
122   if (Y->ops->axpy && (sametype || X->ops->axpy == Y->ops->axpy)) {  in MatAXPY()
123     PetscUseTypeMethod(Y, axpy, a, X, str);  in MatAXPY()
125 …PetscCall(PetscObjectTypeCompareAny((PetscObject)X, &transpose, MATTRANSPOSEVIRTUAL, MATHERMITIANT…  in MatAXPY()
127       PetscCall(MatTransposeAXPY_Private(Y, a, X, str, X));  in MatAXPY()
131         PetscCall(MatTransposeAXPY_Private(Y, a, X, str, Y));  in MatAXPY()
133         PetscCall(MatAXPY_BasicWithTypeCompare(Y, a, X, str));  in MatAXPY()
141 PetscErrorCode MatAXPY_Basic_Preallocate(Mat Y, Mat X, Mat *B)  in MatAXPY_Basic_Preallocate()  argument
149     PetscCall((*preall)(Y, X, B));  in MatAXPY_Basic_Preallocate()
156     PetscCall(MatGetRowUpperTriangular(X));  in MatAXPY_Basic_Preallocate()
172       PetscCall(MatGetRow(X, r, &ncols, &row, &vals));  in MatAXPY_Basic_Preallocate()
174       PetscCall(MatRestoreRow(X, r, &ncols, &row, &vals));  in MatAXPY_Basic_Preallocate()
180     PetscCall(MatRestoreRowUpperTriangular(X));  in MatAXPY_Basic_Preallocate()
191 PetscErrorCode MatAXPY_Basic(Mat Y, PetscScalar a, Mat X, MatStructure str)  in MatAXPY_Basic()  argument
208     PetscCall(MatGetSize(X, &m, &n));  in MatAXPY_Basic()
209     PetscCall(MatGetOwnershipRange(X, &start, &end));  in MatAXPY_Basic()
210     PetscCall(MatGetRowUpperTriangular(X));  in MatAXPY_Basic()
213         PetscCall(MatGetRow(X, i, &ncols, &row, &vals));  in MatAXPY_Basic()
215         PetscCall(MatRestoreRow(X, i, &ncols, &row, &vals));  in MatAXPY_Basic()
222         PetscCall(MatGetRow(X, i, &ncols, &row, &vals));  in MatAXPY_Basic()
229         PetscCall(MatRestoreRow(X, i, &ncols, &row, &vals));  in MatAXPY_Basic()
233     PetscCall(MatRestoreRowUpperTriangular(X));  in MatAXPY_Basic()
242     PetscCall(MatAXPY_Basic_Preallocate(Y, X, &B));  in MatAXPY_Basic()
243     PetscCall(MatAXPY_BasicWithPreallocation(B, Y, a, X, str));  in MatAXPY_Basic()
249 PetscErrorCode MatAXPY_BasicWithPreallocation(Mat B, Mat Y, PetscScalar a, Mat X, MatStructure str)  in MatAXPY_BasicWithPreallocation()  argument
258   PetscCall(MatGetSize(X, &m, &n));  in MatAXPY_BasicWithPreallocation()
259   PetscCall(MatGetOwnershipRange(X, &start, &end));  in MatAXPY_BasicWithPreallocation()
261   PetscCall(MatGetRowUpperTriangular(X));  in MatAXPY_BasicWithPreallocation()
268       PetscCall(MatGetRow(X, i, &ncols, &row, &vals));  in MatAXPY_BasicWithPreallocation()
270       PetscCall(MatRestoreRow(X, i, &ncols, &row, &vals));  in MatAXPY_BasicWithPreallocation()
281       PetscCall(MatGetRow(X, i, &ncols, &row, &vals));  in MatAXPY_BasicWithPreallocation()
288       PetscCall(MatRestoreRow(X, i, &ncols, &row, &vals));  in MatAXPY_BasicWithPreallocation()
293   PetscCall(MatRestoreRowUpperTriangular(X));  in MatAXPY_BasicWithPreallocation()
408 PetscErrorCode MatAYPX(Mat Y, PetscScalar a, Mat X, MatStructure str)  in MatAYPX()  argument
412   PetscCall(MatAXPY(Y, 1.0, X, str));  in MatAYPX()