Actual source code: petscvec.h
1: /*
2: Defines the vector component of PETSc. Vectors generally represent
3: degrees of freedom for finite element/finite difference functions
4: on a grid. They have more mathematical structure than simple arrays.
5: */
6: #pragma once
8: #include <petscsys.h>
9: #include <petscsftypes.h>
10: #include <petscis.h>
11: #include <petscdevicetypes.h>
12: #include <petscviewer.h>
14: /* SUBMANSEC = Vec */
16: /*S
17: Vec - Abstract PETSc vector object. Used for holding solutions and right-hand sides for linear systems, nonlinear systems, and time integrators
19: Level: beginner
21: Note:
22: Internally the actual vector representation is generally a simple array but most PETSc code can work on other representations through this abstraction
24: .seealso: [](doc_vector), [](ch_vectors), `VecCreate()`, `VecType`, `VecSetType()`
25: S*/
26: typedef struct _p_Vec *Vec;
28: /*E
29: ScatterMode - Determines the direction of a scatter in `VecScatterBegin()` and `VecScatterEnd()`
31: Values:
32: + `SCATTER_FORWARD` - Scatters the values as dictated by the `VecScatterCreate()` call
33: . `SCATTER_REVERSE` - Moves the values in the opposite direction than the directions indicated in the `VecScatterCreate()` call
34: . `SCATTER_FORWARD_LOCAL` - Scatters the values as dictated by the `VecScatterCreate()` call except NO MPI communication is done
35: - `SCATTER_REVERSE_LOCAL` - Moves the values in the opposite direction than the directions indicated in the `VecScatterCreate()` call
36: except NO MPI communication is done
38: Level: beginner
40: .seealso: [](ch_vectors), `VecScatter`, `VecScatterBegin()`, `VecScatterEnd()`, `SCATTER_FORWARD`, `SCATTER_REVERSE`, `SCATTER_FORWARD_LOCAL`, `SCATTER_REVERSE_LOCAL`
41: E*/
42: typedef enum {
43: SCATTER_FORWARD = 0,
44: SCATTER_REVERSE = 1,
45: SCATTER_FORWARD_LOCAL = 2,
46: SCATTER_REVERSE_LOCAL = 3
47: } ScatterMode;
49: /*MC
50: SCATTER_FORWARD - Scatters the values as dictated by the `VecScatterCreate()` call during `VecScatterBegin()` and `VecScatterEnd()`
52: Level: beginner
54: .seealso: [](ch_vectors), `VecScatter`, `ScatterMode`, `VecScatterCreate()`, `VecScatterBegin()`, `VecScatterEnd()`, `SCATTER_REVERSE`, `SCATTER_FORWARD_LOCAL`,
55: `SCATTER_REVERSE_LOCAL`
56: M*/
58: /*MC
59: SCATTER_REVERSE - Moves the values in the opposite direction than the directions indicated
60: in the `VecScatterCreate()` during `VecScatterBegin()` and `VecScatterEnd()`
62: Level: beginner
64: .seealso: [](ch_vectors), `VecScatter`, `ScatterMode`, `VecScatterCreate()`, `VecScatterBegin()`, `VecScatterEnd()`, `SCATTER_FORWARD`, `SCATTER_FORWARD_LOCAL`,
65: `SCATTER_REVERSE_LOCAL`
66: M*/
68: /*MC
69: SCATTER_FORWARD_LOCAL - Scatters the values as dictated by the `VecScatterCreate()` during `VecScatterBegin()` and `VecScatterEnd()` call except NO parallel communication
70: is done. Any variables that have be moved between processes are ignored
72: Level: developer
74: .seealso: [](ch_vectors), `VecScatter`, `ScatterMode`, `VecScatterCreate()`, `VecScatterBegin()`, `VecScatterEnd()`, `SCATTER_REVERSE`, `SCATTER_FORWARD`,
75: `SCATTER_REVERSE_LOCAL`
76: M*/
78: /*MC
79: SCATTER_REVERSE_LOCAL - Moves the values in the opposite direction than the directions indicated
80: in the `VecScatterCreate()` during `VecScatterBegin()` and `VecScatterEnd()` except NO parallel communication
81: is done. Any variables that have be moved between processes are ignored
83: Level: developer
85: .seealso: [](ch_vectors), `VecScatter`, `ScatterMode`, `VecScatterCreate()`, `VecScatterBegin()`, `VecScatterEnd()`, `SCATTER_FORWARD`, `SCATTER_FORWARD_LOCAL`,
86: `SCATTER_REVERSE`
87: M*/
89: /*J
90: VecType - String with the name of a PETSc vector, `Vec`, type
92: Level: beginner
94: .seealso: [](doc_vector), [](ch_vectors), `VecSetType()`, `Vec`, `VecCreate()`, `VecDestroy()`
95: J*/
96: typedef const char *VecType;
97: #define VECSEQ "seq"
98: #define VECMPI "mpi"
99: #define VECSTANDARD "standard" /* seq on one process and mpi on multiple */
100: #define VECSHARED "shared"
101: #define VECSEQVIENNACL "seqviennacl"
102: #define VECMPIVIENNACL "mpiviennacl"
103: #define VECVIENNACL "viennacl" /* seqviennacl on one process and mpiviennacl on multiple */
104: #define VECSEQCUDA "seqcuda"
105: #define VECMPICUDA "mpicuda"
106: #define VECCUDA "cuda" /* seqcuda on one process and mpicuda on multiple */
107: #define VECSEQHIP "seqhip"
108: #define VECMPIHIP "mpihip"
109: #define VECHIP "hip" /* seqhip on one process and mpihip on multiple */
110: #define VECNEST "nest"
111: #define VECSEQKOKKOS "seqkokkos"
112: #define VECMPIKOKKOS "mpikokkos"
113: #define VECKOKKOS "kokkos" /* seqkokkos on one process and mpikokkos on multiple */
115: /* Dynamic creation and loading functions */
116: PETSC_EXTERN PetscErrorCode VecScatterSetType(VecScatter, VecScatterType);
117: PETSC_EXTERN PetscErrorCode VecScatterGetType(VecScatter, VecScatterType *);
118: PETSC_EXTERN PetscErrorCode VecScatterSetFromOptions(VecScatter);
119: PETSC_EXTERN PetscErrorCode VecScatterRegister(const char[], PetscErrorCode (*)(VecScatter));
120: PETSC_EXTERN PetscErrorCode VecScatterCreate(Vec, IS, Vec, IS, VecScatter *);
122: /* Logging support */
123: #define REAL_FILE_CLASSID 1211213
124: #define VEC_FILE_CLASSID 1211214
125: PETSC_EXTERN PetscClassId VEC_CLASSID;
126: PETSC_EXTERN PetscClassId PETSCSF_CLASSID;
128: PETSC_EXTERN PetscErrorCode VecInitializePackage(void);
129: PETSC_EXTERN PetscErrorCode VecFinalizePackage(void);
131: PETSC_EXTERN PetscErrorCode VecCreate(MPI_Comm, Vec *);
132: PETSC_EXTERN PetscErrorCode VecCreateFromOptions(MPI_Comm, const char *, PetscInt, PetscInt, PetscInt, Vec *);
133: PETSC_EXTERN PetscErrorCode VecCreateSeq(MPI_Comm, PetscInt, Vec *);
134: PETSC_EXTERN PetscErrorCode VecCreateMPI(MPI_Comm, PetscInt, PetscInt, Vec *);
135: PETSC_EXTERN PetscErrorCode VecCreateSeqWithArray(MPI_Comm, PetscInt, PetscInt, const PetscScalar[], Vec *);
136: PETSC_EXTERN PetscErrorCode VecCreateMPIWithArray(MPI_Comm, PetscInt, PetscInt, PetscInt, const PetscScalar[], Vec *);
137: PETSC_EXTERN PetscErrorCode VecCreateShared(MPI_Comm, PetscInt, PetscInt, Vec *);
139: PETSC_EXTERN PetscErrorCode VecSetFromOptions(Vec);
140: PETSC_EXTERN PetscErrorCode VecViewFromOptions(Vec, PetscObject, const char[]);
142: PETSC_EXTERN PetscErrorCode VecSetUp(Vec);
143: PETSC_EXTERN PetscErrorCode VecDestroy(Vec *);
144: PETSC_EXTERN PetscErrorCode VecZeroEntries(Vec);
145: PETSC_EXTERN PetscErrorCode VecSetOptionsPrefix(Vec, const char[]);
146: PETSC_EXTERN PetscErrorCode VecAppendOptionsPrefix(Vec, const char[]);
147: PETSC_EXTERN PetscErrorCode VecGetOptionsPrefix(Vec, const char *[]);
148: PETSC_EXTERN PetscErrorCode VecGetState(Vec, PetscObjectState *);
150: PETSC_EXTERN PetscErrorCode VecSetSizes(Vec, PetscInt, PetscInt);
152: PETSC_EXTERN PetscErrorCode VecDotNorm2(Vec, Vec, PetscScalar *, PetscReal *);
153: PETSC_EXTERN PetscErrorCode VecDot(Vec, Vec, PetscScalar *);
154: PETSC_EXTERN PetscErrorCode VecDotRealPart(Vec, Vec, PetscReal *);
155: PETSC_EXTERN PetscErrorCode VecTDot(Vec, Vec, PetscScalar *);
156: PETSC_EXTERN PetscErrorCode VecMDot(Vec, PetscInt, const Vec[], PetscScalar[]);
157: PETSC_EXTERN PetscErrorCode VecMTDot(Vec, PetscInt, const Vec[], PetscScalar[]);
158: PETSC_EXTERN PetscErrorCode VecGetSubVector(Vec, IS, Vec *);
159: PETSC_EXTERN PetscErrorCode VecRestoreSubVector(Vec, IS, Vec *);
160: PETSC_EXTERN PetscErrorCode VecConcatenate(PetscInt, const Vec[], Vec *, IS *[]);
162: /*E
163: NormType - determines what type of norm to compute with `VecNorm()`, `VecNormBegin()`/`VecNormEnd()` and `MatNorm()`.
165: Values:
166: + `NORM_1` - the one norm, $||v|| = \sum_i | v_i |$. $||A|| = \max_j || A_{*j} ||$, maximum column sum
167: . `NORM_2` - the two norm, $||v|| = sqrt(\sum_i |v_i|^2)$ (vectors only)
168: . `NORM_FROBENIUS` - $||A|| = sqrt(\sum_{ij} |A_{ij}|^2)$, same as `NORM_2` for vectors
169: . `NORM_INFINITY` - $||v|| = \max_i |v_i|$. $||A|| = \max_i || A_{i*} ||_1$, maximum row sum
170: - `NORM_1_AND_2` - computes both the 1 and 2 norm of a vector. The values are stored in two adjacent `PetscReal` memory locations
172: Level: beginner
174: Note:
175: The `v` above represents a `Vec` while the `A` represents a `Mat`
177: .seealso: [](ch_vectors), `Vec`, `Mat`, `VecNorm()`, `VecNormBegin()`, `VecNormEnd()`, `MatNorm()`, `NORM_1`,
178: `NORM_2`, `NORM_FROBENIUS`, `NORM_INFINITY`, `NORM_1_AND_2`, `ReductionType`
179: E*/
180: typedef enum {
181: NORM_1 = 0,
182: NORM_2 = 1,
183: NORM_FROBENIUS = 2,
184: NORM_INFINITY = 3,
185: NORM_1_AND_2 = 4
186: } NormType;
187: PETSC_EXTERN const char *const NormTypes[];
188: #define NORM_MAX NORM_INFINITY
190: /*MC
191: NORM_1 - the one norm, $||v|| = \sum_i | v_i |$. $||A|| = \max_j || A_{*,j} ||$, maximum column sum
193: Level: beginner
195: .seealso: [](ch_vectors), `NormType`, `MatNorm()`, `VecNorm()`, `VecNormBegin()`, `VecNormEnd()`, `NORM_2`, `NORM_FROBENIUS`,
196: `NORM_INFINITY`, `NORM_1_AND_2`
197: M*/
199: /*MC
200: NORM_2 - the two norm, $||v|| = \sqrt{\sum_i |v_i|^2}$ (vectors only)
202: Level: beginner
204: .seealso: [](ch_vectors), `NormType`, `MatNorm()`, `VecNorm()`, `VecNormBegin()`, `VecNormEnd()`, `NORM_1`, `NORM_FROBENIUS`,
205: `NORM_INFINITY`, `NORM_1_AND_2`
206: M*/
208: /*MC
209: NORM_FROBENIUS - $||A|| = \sqrt{\sum_{i,j} |A_{i,j}|^2}$, same as `NORM_2` for vectors
211: Level: beginner
213: .seealso: [](ch_vectors), `NormType`, `MatNorm()`, `VecNorm()`, `VecNormBegin()`, `VecNormEnd()`, `NORM_1`, `NORM_2`,
214: `NORM_INFINITY`, `NORM_1_AND_2`
215: M*/
217: /*MC
218: NORM_INFINITY - $||v|| = \max_i |v_i|$. $||A|| = \max_i || A_{i,*} ||_1$, maximum row sum
220: Level: beginner
222: .seealso: [](ch_vectors), `NormType`, `MatNorm()`, `VecNorm()`, `VecNormBegin()`, `VecNormEnd()`, `NORM_1`, `NORM_2`,
223: `NORM_FROBENIUS`, `NORM_1_AND_2`
224: M*/
226: /*MC
227: NORM_1_AND_2 - computes both the 1 and 2 norm of a vector. The values are stored in two adjacent `PetscReal` memory locations
229: Level: beginner
231: .seealso: [](ch_vectors), `NormType`, `MatNorm()`, `VecNorm()`, `VecNormBegin()`, `VecNormEnd()`, `NORM_1`, `NORM_2`,
232: `NORM_FROBENIUS`, `NORM_INFINITY`
233: M*/
235: /*MC
236: NORM_MAX - see `NORM_INFINITY`
238: Level: beginner
239: M*/
241: /*E
242: ReductionType - determines what type of column reduction (one that is not a type of norm defined in `NormType`) to obtain with `MatGetColumnReductions()`
244: Values:
245: + `REDUCTION_SUM_REALPART` - sum of real part of each matrix column
246: . `REDUCTION_SUM_IMAGINARYPART` - sum of imaginary part of each matrix column
247: . `REDUCTION_MEAN_REALPART` - arithmetic mean of real part of each matrix column
248: - `REDUCTION_MEAN_IMAGINARYPART` - arithmetic mean of imaginary part of each matrix column
250: Level: beginner
252: Developer Note:
253: The constants defined in `ReductionType` MUST BE DISTINCT from those defined in `NormType`.
254: This is because `MatGetColumnReductions()` is used to compute both norms and other types of reductions,
255: and the constants defined in both `NormType` and `ReductionType` are used to designate the desired operation.
257: .seealso: [](ch_vectors), `MatGetColumnReductions()`, `MatGetColumnNorms()`, `NormType`, `REDUCTION_SUM_REALPART`,
258: `REDUCTION_SUM_IMAGINARYPART`, `REDUCTION_MEAN_REALPART`, `REDUCTION_MEAN_IMAGINARYPART`
259: E*/
260: typedef enum {
261: REDUCTION_SUM_REALPART = 10,
262: REDUCTION_MEAN_REALPART = 11,
263: REDUCTION_SUM_IMAGINARYPART = 12,
264: REDUCTION_MEAN_IMAGINARYPART = 13
265: } ReductionType;
267: /*E
268: VecSignMode - How `VecPointwiseSign()` should handle zero value
270: Values:
271: + `VEC_SIGN_ZERO_TO_ZERO` - `-0.0` and `0.0` map to `0.0`
272: . `VEC_SIGN_ZERO_TO_SIGNED_ZERO` - `-0.0` maps to `-0.0` and `0.0` maps to `0.0`
273: - `VEC_SIGN_ZERO_TO_SIGNED_UNIT` - `-0.0` maps to `-1.0` and `0.0` maps to `1.0`
275: Level: advanced
277: .seealso: [](ch_vectors), `VecPointwiseSign()`
278: E*/
279: typedef enum {
280: VEC_SIGN_ZERO_TO_ZERO,
281: VEC_SIGN_ZERO_TO_SIGNED_ZERO,
282: VEC_SIGN_ZERO_TO_SIGNED_UNIT,
283: } VecSignMode;
285: /*MC
286: REDUCTION_SUM_REALPART - sum of real part of a matrix column to obtain with `MatGetColumnReductions()`
288: Level: beginner
290: .seealso: [](ch_vectors), `ReductionType`, `MatGetColumnReductions()`, `REDUCTION_SUM_IMAGINARYPART`, `REDUCTION_MEAN_REALPART`
291: M*/
293: /*MC
294: REDUCTION_SUM_IMAGINARYPART - sum of imaginary part of matrix column to obtain with `MatGetColumnReductions()`
296: Level: beginner
298: .seealso: [](ch_vectors), `ReductionType`, `MatGetColumnReductions()`, `REDUCTION_SUM_REALPART`, `REDUCTION_MEAN_IMAGINARYPART`
299: M*/
301: /*MC
302: REDUCTION_MEAN_REALPART - arithmetic mean of real part of matrix column to obtain with `MatGetColumnReductions()`
304: Level: beginner
306: .seealso: [](ch_vectors), `ReductionType`, `MatGetColumnReductions()`, `REDUCTION_MEAN_IMAGINARYPART`, `REDUCTION_SUM_REALPART`
307: M*/
309: /*MC
310: REDUCTION_MEAN_IMAGINARYPART - arithmetic mean of imaginary part of matrix column to obtain with `MatGetColumnReductions()`
312: Level: beginner
314: .seealso: [](ch_vectors), `ReductionType`, `MatGetColumnReductions()`, `REDUCTION_MEAN_REALPART`, `REDUCTION_SUM_IMAGINARYPART`
315: M*/
317: PETSC_EXTERN PetscErrorCode VecNorm(Vec, NormType, PetscReal *);
318: PETSC_EXTERN PetscErrorCode VecNormAvailable(Vec, NormType, PetscBool *, PetscReal *);
319: PETSC_EXTERN PetscErrorCode VecFlag(Vec, PetscInt);
320: PETSC_EXTERN PetscErrorCode VecNormalize(Vec, PetscReal *);
321: PETSC_EXTERN PetscErrorCode VecSum(Vec, PetscScalar *);
322: PETSC_EXTERN PetscErrorCode VecMean(Vec, PetscScalar *);
323: PETSC_EXTERN PetscErrorCode VecMax(Vec, PetscInt *, PetscReal *);
324: PETSC_EXTERN PetscErrorCode VecMin(Vec, PetscInt *, PetscReal *);
325: PETSC_EXTERN PetscErrorCode VecScale(Vec, PetscScalar);
326: PETSC_EXTERN PetscErrorCode VecCopy(Vec, Vec);
327: PETSC_EXTERN PetscErrorCode VecSetRandom(Vec, PetscRandom);
328: PETSC_EXTERN PetscErrorCode VecSetRandomGaussian(Vec, PetscRandom, PetscReal, PetscReal);
329: PETSC_EXTERN PetscErrorCode VecSet(Vec, PetscScalar);
330: PETSC_EXTERN PetscErrorCode VecSetStdBasis(Vec, PetscInt);
331: PETSC_DEPRECATED_FUNCTION(3, 22, 0, "VecFlag()", ) PetscErrorCode VecSetInf(Vec);
332: PETSC_EXTERN PetscErrorCode VecSwap(Vec, Vec);
333: PETSC_EXTERN PetscErrorCode VecAXPY(Vec, PetscScalar, Vec);
334: PETSC_EXTERN PetscErrorCode VecAXPBY(Vec, PetscScalar, PetscScalar, Vec);
335: PETSC_EXTERN PetscErrorCode VecMAXPY(Vec, PetscInt, const PetscScalar[], Vec[]);
336: PETSC_EXTERN PetscErrorCode VecMAXPBY(Vec, PetscInt, const PetscScalar[], PetscScalar, Vec[]);
337: PETSC_EXTERN PetscErrorCode VecAYPX(Vec, PetscScalar, Vec);
338: PETSC_EXTERN PetscErrorCode VecWAXPY(Vec, PetscScalar, Vec, Vec);
339: PETSC_EXTERN PetscErrorCode VecAXPBYPCZ(Vec, PetscScalar, PetscScalar, PetscScalar, Vec, Vec);
340: PETSC_EXTERN PetscErrorCode VecPointwiseMax(Vec, Vec, Vec);
341: PETSC_EXTERN PetscErrorCode VecPointwiseMaxAbs(Vec, Vec, Vec);
342: PETSC_EXTERN PetscErrorCode VecPointwiseMin(Vec, Vec, Vec);
343: PETSC_EXTERN PetscErrorCode VecPointwiseMult(Vec, Vec, Vec);
344: PETSC_EXTERN PetscErrorCode VecPointwiseDivide(Vec, Vec, Vec);
345: PETSC_EXTERN PetscErrorCode VecPointwiseSign(Vec, Vec, VecSignMode);
346: PETSC_EXTERN PetscErrorCode VecMaxPointwiseDivide(Vec, Vec, PetscReal *);
347: PETSC_EXTERN PetscErrorCode VecShift(Vec, PetscScalar);
348: PETSC_EXTERN PetscErrorCode VecReciprocal(Vec);
349: PETSC_EXTERN PetscErrorCode VecPermute(Vec, IS, PetscBool);
350: PETSC_EXTERN PetscErrorCode VecSqrtAbs(Vec);
351: PETSC_EXTERN PetscErrorCode VecLog(Vec);
352: PETSC_EXTERN PetscErrorCode VecExp(Vec);
353: PETSC_EXTERN PetscErrorCode VecAbs(Vec);
354: PETSC_EXTERN PetscErrorCode VecDuplicate(Vec, Vec *);
355: PETSC_EXTERN PetscErrorCode VecDuplicateVecs(Vec, PetscInt, Vec *[]);
356: PETSC_EXTERN PetscErrorCode VecDestroyVecs(PetscInt, Vec *[]);
357: PETSC_EXTERN PetscErrorCode VecStrideNormAll(Vec, NormType, PetscReal[]);
358: PETSC_EXTERN PetscErrorCode VecStrideMaxAll(Vec, PetscInt[], PetscReal[]);
359: PETSC_EXTERN PetscErrorCode VecStrideMinAll(Vec, PetscInt[], PetscReal[]);
360: PETSC_EXTERN PetscErrorCode VecStrideScaleAll(Vec, const PetscScalar[]);
361: PETSC_EXTERN PetscErrorCode VecStrideSumAll(Vec, PetscScalar *);
362: PETSC_EXTERN PetscErrorCode VecUniqueEntries(Vec, PetscInt *, PetscScalar *[]);
364: PETSC_EXTERN PetscErrorCode VecStrideNorm(Vec, PetscInt, NormType, PetscReal *);
365: PETSC_EXTERN PetscErrorCode VecStrideMax(Vec, PetscInt, PetscInt *, PetscReal *);
366: PETSC_EXTERN PetscErrorCode VecStrideMin(Vec, PetscInt, PetscInt *, PetscReal *);
367: PETSC_EXTERN PetscErrorCode VecStrideScale(Vec, PetscInt, PetscScalar);
368: PETSC_EXTERN PetscErrorCode VecStrideSum(Vec, PetscInt, PetscScalar *);
369: PETSC_EXTERN PetscErrorCode VecStrideSet(Vec, PetscInt, PetscScalar);
371: PETSC_EXTERN PetscErrorCode VecStrideGather(Vec, PetscInt, Vec, InsertMode);
372: PETSC_EXTERN PetscErrorCode VecStrideScatter(Vec, PetscInt, Vec, InsertMode);
373: PETSC_EXTERN PetscErrorCode VecStrideGatherAll(Vec, Vec[], InsertMode);
374: PETSC_EXTERN PetscErrorCode VecStrideScatterAll(Vec[], Vec, InsertMode);
376: PETSC_EXTERN PetscErrorCode VecStrideSubSetScatter(Vec, PetscInt, const PetscInt[], const PetscInt[], Vec, InsertMode);
377: PETSC_EXTERN PetscErrorCode VecStrideSubSetGather(Vec, PetscInt, const PetscInt[], const PetscInt[], Vec, InsertMode);
379: PETSC_EXTERN PetscErrorCode VecSetValues(Vec, PetscInt, const PetscInt[], const PetscScalar[], InsertMode);
380: PETSC_EXTERN PetscErrorCode VecGetValues(Vec, PetscInt, const PetscInt[], PetscScalar[]);
381: PETSC_EXTERN PetscErrorCode VecAssemblyBegin(Vec);
382: PETSC_EXTERN PetscErrorCode VecAssemblyEnd(Vec);
383: PETSC_EXTERN PetscErrorCode VecStashSetInitialSize(Vec, PetscInt, PetscInt);
384: PETSC_EXTERN PetscErrorCode VecStashView(Vec, PetscViewer);
385: PETSC_EXTERN PetscErrorCode VecStashViewFromOptions(Vec, PetscObject, const char[]);
386: PETSC_EXTERN PetscErrorCode VecStashGetInfo(Vec, PetscInt *, PetscInt *, PetscInt *, PetscInt *);
388: PETSC_EXTERN PetscErrorCode VecSetPreallocationCOO(Vec, PetscCount, const PetscInt[]);
389: PETSC_EXTERN PetscErrorCode VecSetPreallocationCOOLocal(Vec, PetscCount, PetscInt[]);
390: PETSC_EXTERN PetscErrorCode VecSetValuesCOO(Vec, const PetscScalar[], InsertMode);
392: /*@C
393: VecSetValue - Set a single entry into a PETSc vector, `Vec`.
395: Not Collective
397: Input Parameters:
398: + v - the vector
399: . row - the row location of the entry
400: . value - the value to insert
401: - mode - either `INSERT_VALUES` or `ADD_VALUES`
403: Level: beginner
405: Notes:
406: For efficiency one should use `VecSetValues()` and set several or
407: many values simultaneously if possible.
409: These values may be cached, so `VecAssemblyBegin()` and `VecAssemblyEnd()`
410: MUST be called after all calls to `VecSetValue()` have been completed.
412: `VecSetValue()` uses 0-based indices in Python, C, and Fortran
414: .seealso: [](ch_vectors), `VecSetValues()`, `VecAssemblyBegin()`, `VecAssemblyEnd()`, `VecSetValuesBlockedLocal()`, `VecSetValueLocal()`
415: @*/
416: static inline PetscErrorCode VecSetValue(Vec v, PetscInt i, PetscScalar va, InsertMode mode)
417: {
418: return VecSetValues(v, 1, &i, &va, mode);
419: }
421: PETSC_EXTERN PetscErrorCode VecSetBlockSize(Vec, PetscInt);
422: PETSC_EXTERN PetscErrorCode VecGetBlockSize(Vec, PetscInt *);
423: PETSC_EXTERN PetscErrorCode VecSetValuesBlocked(Vec, PetscInt, const PetscInt[], const PetscScalar[], InsertMode);
425: /* Dynamic creation and loading functions */
426: PETSC_EXTERN PetscFunctionList VecList;
427: PETSC_EXTERN PetscErrorCode VecSetType(Vec, VecType);
428: PETSC_EXTERN PetscErrorCode VecGetType(Vec, VecType *);
429: PETSC_EXTERN PetscErrorCode VecRegister(const char[], PetscErrorCode (*)(Vec));
430: PETSC_EXTERN PetscErrorCode VecRegisterAll(void);
432: PETSC_EXTERN PetscErrorCode VecScatterBegin(VecScatter, Vec, Vec, InsertMode, ScatterMode);
433: PETSC_EXTERN PetscErrorCode VecScatterEnd(VecScatter, Vec, Vec, InsertMode, ScatterMode);
434: PETSC_EXTERN PetscErrorCode VecScatterDestroy(VecScatter *);
435: PETSC_EXTERN PetscErrorCode VecScatterSetUp(VecScatter);
436: PETSC_EXTERN PetscErrorCode VecScatterCopy(VecScatter, VecScatter *);
437: PETSC_EXTERN PetscErrorCode VecScatterView(VecScatter, PetscViewer);
438: PETSC_EXTERN PetscErrorCode VecScatterViewFromOptions(VecScatter, PetscObject, const char[]);
439: PETSC_EXTERN PetscErrorCode VecScatterRemap(VecScatter, PetscInt[], PetscInt[]);
440: PETSC_EXTERN PetscErrorCode VecScatterGetMerged(VecScatter, PetscBool *);
442: PETSC_EXTERN PetscErrorCode VecGetArray4d(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ****[]);
443: PETSC_EXTERN PetscErrorCode VecRestoreArray4d(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ****[]);
444: PETSC_EXTERN PetscErrorCode VecGetArray3d(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ***[]);
445: PETSC_EXTERN PetscErrorCode VecRestoreArray3d(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ***[]);
446: PETSC_EXTERN PetscErrorCode VecGetArray2d(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar **[]);
447: PETSC_EXTERN PetscErrorCode VecRestoreArray2d(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar **[]);
448: PETSC_EXTERN PetscErrorCode VecGetArray1d(Vec, PetscInt, PetscInt, PetscScalar *[]);
449: PETSC_EXTERN PetscErrorCode VecRestoreArray1d(Vec, PetscInt, PetscInt, PetscScalar *[]);
451: PETSC_EXTERN PetscErrorCode VecGetArray4dWrite(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ****[]);
452: PETSC_EXTERN PetscErrorCode VecGetArray4dWrite(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ****[]);
453: PETSC_EXTERN PetscErrorCode VecRestoreArray4dWrite(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ****[]);
454: PETSC_EXTERN PetscErrorCode VecGetArray3dWrite(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ***[]);
455: PETSC_EXTERN PetscErrorCode VecRestoreArray3dWrite(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ***[]);
456: PETSC_EXTERN PetscErrorCode VecGetArray2dWrite(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar **[]);
457: PETSC_EXTERN PetscErrorCode VecRestoreArray2dWrite(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar **[]);
458: PETSC_EXTERN PetscErrorCode VecGetArray1dWrite(Vec, PetscInt, PetscInt, PetscScalar *[]);
459: PETSC_EXTERN PetscErrorCode VecRestoreArray1dWrite(Vec, PetscInt, PetscInt, PetscScalar *[]);
461: PETSC_EXTERN PetscErrorCode VecGetArray4dRead(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ****[]);
462: PETSC_EXTERN PetscErrorCode VecRestoreArray4dRead(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ****[]);
463: PETSC_EXTERN PetscErrorCode VecGetArray3dRead(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ***[]);
464: PETSC_EXTERN PetscErrorCode VecRestoreArray3dRead(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar ***[]);
465: PETSC_EXTERN PetscErrorCode VecGetArray2dRead(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar **[]);
466: PETSC_EXTERN PetscErrorCode VecRestoreArray2dRead(Vec, PetscInt, PetscInt, PetscInt, PetscInt, PetscScalar **[]);
467: PETSC_EXTERN PetscErrorCode VecGetArray1dRead(Vec, PetscInt, PetscInt, PetscScalar *[]);
468: PETSC_EXTERN PetscErrorCode VecRestoreArray1dRead(Vec, PetscInt, PetscInt, PetscScalar *[]);
470: PETSC_EXTERN PetscErrorCode VecPlaceArray(Vec, const PetscScalar[]);
471: PETSC_EXTERN PetscErrorCode VecResetArray(Vec);
472: PETSC_EXTERN PetscErrorCode VecReplaceArray(Vec, const PetscScalar[]);
474: PETSC_EXTERN PetscErrorCode VecGetArrays(const Vec[], PetscInt, PetscScalar **[]);
475: PETSC_EXTERN PetscErrorCode VecRestoreArrays(const Vec[], PetscInt, PetscScalar **[]);
477: PETSC_EXTERN PetscErrorCode VecView(Vec, PetscViewer);
478: PETSC_EXTERN PetscErrorCode VecViewNative(Vec, PetscViewer);
479: PETSC_EXTERN PetscErrorCode VecEqual(Vec, Vec, PetscBool *);
480: PETSC_EXTERN PetscErrorCode VecLoad(Vec, PetscViewer);
482: PETSC_EXTERN PetscErrorCode VecGetSize(Vec, PetscInt *);
483: PETSC_EXTERN PetscErrorCode VecGetLocalSize(Vec, PetscInt *);
484: PETSC_EXTERN PetscErrorCode VecGetOwnershipRange(Vec, PetscInt *, PetscInt *);
485: PETSC_EXTERN PetscErrorCode VecGetOwnershipRanges(Vec, const PetscInt *[]);
487: PETSC_EXTERN PetscErrorCode VecSetLocalToGlobalMapping(Vec, ISLocalToGlobalMapping);
488: PETSC_EXTERN PetscErrorCode VecSetValuesLocal(Vec, PetscInt, const PetscInt[], const PetscScalar[], InsertMode);
490: PETSC_EXTERN PetscErrorCode VecViennaCLGetCLContext(Vec, PETSC_UINTPTR_T *);
491: PETSC_EXTERN PetscErrorCode VecViennaCLGetCLQueue(Vec, PETSC_UINTPTR_T *);
492: PETSC_EXTERN PetscErrorCode VecViennaCLGetCLMemRead(Vec, PETSC_UINTPTR_T *);
493: PETSC_EXTERN PetscErrorCode VecViennaCLGetCLMemWrite(Vec, PETSC_UINTPTR_T *);
494: PETSC_EXTERN PetscErrorCode VecViennaCLRestoreCLMemWrite(Vec);
495: PETSC_EXTERN PetscErrorCode VecViennaCLGetCLMem(Vec, PETSC_UINTPTR_T *);
496: PETSC_EXTERN PetscErrorCode VecViennaCLRestoreCLMem(Vec);
498: /*@C
499: VecSetValueLocal - Set a single entry into a vector using the local numbering of the vector, see `VecSetValuesLocal()`
501: Not Collective
503: Input Parameters:
504: + v - the vector
505: . row - the local row location of the entry
506: . value - the value to insert
507: - mode - either `INSERT_VALUES` or `ADD_VALUES`
509: Level: beginner
511: Notes:
512: For efficiency one should use `VecSetValuesLocal()` and set several or
513: many values simultaneously if possible.
515: These values may be cached, so `VecAssemblyBegin()` and `VecAssemblyEnd()`
516: MUST be called after all calls to `VecSetValueLocal()` have been completed.
518: See `VecSetLocalToGlobalMapping()` for how the local numbering is defined
520: `VecSetValueLocal()` uses 0-based indices in Fortran as well as in C.
522: .seealso: [](ch_vectors), `VecSetValuesLocal()`, `VecSetValues()`, `VecAssemblyBegin()`, `VecAssemblyEnd()`, `VecSetValuesBlockedLocal()`, `VecSetValue()`,
523: `VecSetLocalToGlobalMapping()`
524: @*/
525: static inline PetscErrorCode VecSetValueLocal(Vec v, PetscInt i, PetscScalar va, InsertMode mode)
526: {
527: return VecSetValuesLocal(v, 1, &i, &va, mode);
528: }
530: /*MC
531: VecCheckAssembled - checks if values have been changed in the vector, by `VecSetValues()` or related routines, but it has not been assembled
533: Synopsis:
534: #include <petscvec.h>
535: VecCheckAssembled(Vec v);
537: Not Collective
539: Input Parameter:
540: . v - the vector to check
542: Level: developer
544: Note:
545: After calls to `VecSetValues()` and related routines one must call `VecAssemblyBegin()` and `VecAssemblyEnd()` before using the vector
547: .seealso: [](ch_vectors), `Vec`, `VecSetValues()`, `VecAssemblyBegin()`, `VecAssemblyEnd()`, `MatAssemblyBegin()`, `MatAssemblyEnd()`
548: M*/
549: #define VecCheckAssembled(v) PetscCheck(v->stash.insertmode == NOT_SET_VALUES, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONGSTATE, "Not for unassembled vector, did you call VecAssemblyBegin()/VecAssemblyEnd()?")
551: PETSC_EXTERN PetscErrorCode VecSetValuesBlockedLocal(Vec, PetscInt, const PetscInt[], const PetscScalar[], InsertMode);
552: PETSC_EXTERN PetscErrorCode VecGetLocalToGlobalMapping(Vec, ISLocalToGlobalMapping *);
554: PETSC_EXTERN PetscErrorCode VecDotBegin(Vec, Vec, PetscScalar *);
555: PETSC_EXTERN PetscErrorCode VecDotEnd(Vec, Vec, PetscScalar *);
556: PETSC_EXTERN PetscErrorCode VecTDotBegin(Vec, Vec, PetscScalar *);
557: PETSC_EXTERN PetscErrorCode VecTDotEnd(Vec, Vec, PetscScalar *);
558: PETSC_EXTERN PetscErrorCode VecNormBegin(Vec, NormType, PetscReal *);
559: PETSC_EXTERN PetscErrorCode VecNormEnd(Vec, NormType, PetscReal *);
560: PETSC_EXTERN PetscErrorCode VecErrorWeightedNorms(Vec, Vec, Vec, NormType, PetscReal, Vec, PetscReal, Vec, PetscReal, PetscReal *, PetscInt *, PetscReal *, PetscInt *, PetscReal *, PetscInt *);
562: PETSC_EXTERN PetscErrorCode VecMDotBegin(Vec, PetscInt, const Vec[], PetscScalar[]);
563: PETSC_EXTERN PetscErrorCode VecMDotEnd(Vec, PetscInt, const Vec[], PetscScalar[]);
564: PETSC_EXTERN PetscErrorCode VecMTDotBegin(Vec, PetscInt, const Vec[], PetscScalar[]);
565: PETSC_EXTERN PetscErrorCode VecMTDotEnd(Vec, PetscInt, const Vec[], PetscScalar[]);
566: PETSC_EXTERN PetscErrorCode PetscCommSplitReductionBegin(MPI_Comm);
568: PETSC_EXTERN PetscErrorCode VecBindToCPU(Vec, PetscBool);
569: PETSC_DEPRECATED_FUNCTION(3, 13, 0, "VecBindToCPU()", ) static inline PetscErrorCode VecPinToCPU(Vec v, PetscBool flg)
570: {
571: return VecBindToCPU(v, flg);
572: }
573: PETSC_EXTERN PetscErrorCode VecBoundToCPU(Vec, PetscBool *);
574: PETSC_EXTERN PetscErrorCode VecSetBindingPropagates(Vec, PetscBool);
575: PETSC_EXTERN PetscErrorCode VecGetBindingPropagates(Vec, PetscBool *);
576: PETSC_EXTERN PetscErrorCode VecSetPinnedMemoryMin(Vec, size_t);
577: PETSC_EXTERN PetscErrorCode VecGetPinnedMemoryMin(Vec, size_t *);
579: PETSC_EXTERN PetscErrorCode VecGetOffloadMask(Vec, PetscOffloadMask *);
581: /*E
582: VecOption - Options that may be set for a vector regarding entries passed to `VecSetValues()` and related routines
584: Values:
585: + `VEC_IGNORE_OFF_PROC_ENTRIES` - causes `VecSetValues()` to ignore entries destined to be stored on a separate processor.
586: This can be used to eliminate the global reduction in `VecAssemblyBegin()` if you know that
587: you have only used `VecSetValues()` to set local elements
588: . `VEC_IGNORE_NEGATIVE_INDICES` - means you can pass negative indices in `ix` in calls to `VecSetValues()` or `VecGetValues()`.
589: These rows are simply ignored.
590: - `VEC_SUBSET_OFF_PROC_ENTRIES` - causes `VecAssemblyBegin()` to assume that the off-process entries will always be a subset
591: (possibly equal) of the off-process entries set on the first assembly which had a true
592: `VEC_SUBSET_OFF_PROC_ENTRIES` and the vector has not changed this flag afterwards. If this
593: assembly is not such first assembly, then this assembly can reuse the communication pattern
594: setup in that first assembly, thus avoiding a global reduction. Subsequent assemblies setting
595: off-process values should use the same `InsertMode` as the first assembly.
597: Level: beginner
599: .seealso: [](ch_matrices), `Vec`, `MatSetOption()`, `VecSetOption()`, `VecSetValues()`, `VecAssemblyBegin()`
600: E*/
601: typedef enum {
602: VEC_IGNORE_OFF_PROC_ENTRIES,
603: VEC_IGNORE_NEGATIVE_INDICES,
604: VEC_SUBSET_OFF_PROC_ENTRIES
605: } VecOption;
607: PETSC_EXTERN PetscErrorCode VecSetOption(Vec, VecOption, PetscBool);
609: PETSC_EXTERN PetscErrorCode VecGetArray(Vec, PetscScalar *[]);
610: PETSC_EXTERN PetscErrorCode VecGetArrayWrite(Vec, PetscScalar *[]);
611: PETSC_EXTERN PetscErrorCode VecGetArrayRead(Vec, const PetscScalar *[]);
612: PETSC_EXTERN PetscErrorCode VecRestoreArray(Vec, PetscScalar *[]);
613: PETSC_EXTERN PetscErrorCode VecRestoreArrayWrite(Vec, PetscScalar *[]);
614: PETSC_EXTERN PetscErrorCode VecRestoreArrayRead(Vec, const PetscScalar *[]);
615: PETSC_EXTERN PetscErrorCode VecCreateLocalVector(Vec, Vec *);
616: PETSC_EXTERN PetscErrorCode VecGetLocalVector(Vec, Vec);
617: PETSC_EXTERN PetscErrorCode VecRestoreLocalVector(Vec, Vec);
618: PETSC_EXTERN PetscErrorCode VecGetLocalVectorRead(Vec, Vec);
619: PETSC_EXTERN PetscErrorCode VecRestoreLocalVectorRead(Vec, Vec);
620: PETSC_EXTERN PetscErrorCode VecGetArrayAndMemType(Vec, PetscScalar *[], PetscMemType *);
621: PETSC_EXTERN PetscErrorCode VecRestoreArrayAndMemType(Vec, PetscScalar *[]);
622: PETSC_EXTERN PetscErrorCode VecGetArrayReadAndMemType(Vec, const PetscScalar *[], PetscMemType *);
623: PETSC_EXTERN PetscErrorCode VecRestoreArrayReadAndMemType(Vec, const PetscScalar *[]);
624: PETSC_EXTERN PetscErrorCode VecGetArrayWriteAndMemType(Vec, PetscScalar *[], PetscMemType *);
625: PETSC_EXTERN PetscErrorCode VecRestoreArrayWriteAndMemType(Vec, PetscScalar *[]);
627: /*@C
628: VecGetArrayPair - Accesses a pair of pointers for two vectors that may be common. When the vectors are not the same the first pointer is read only
630: Logically Collective; No Fortran Support
632: Input Parameters:
633: + x - the vector
634: - y - the second vector
636: Output Parameters:
637: + xv - location to put pointer to the first array
638: - yv - location to put pointer to the second array
640: Level: developer
642: .seealso: [](ch_vectors), `VecGetArray()`, `VecGetArrayRead()`, `VecRestoreArrayPair()`
643: @*/
644: static inline PetscErrorCode VecGetArrayPair(Vec x, Vec y, PetscScalar *xv[], PetscScalar *yv[])
645: {
646: PetscFunctionBegin;
647: PetscCall(VecGetArray(y, yv));
648: if (x == y) *xv = *yv;
649: else PetscCall(VecGetArrayRead(x, (const PetscScalar **)xv));
650: PetscFunctionReturn(PETSC_SUCCESS);
651: }
653: /*@C
654: VecRestoreArrayPair - Returns a pair of pointers for two vectors that may be common obtained with `VecGetArrayPair()`
656: Logically Collective; No Fortran Support
658: Input Parameters:
659: + x - the vector
660: . y - the second vector
661: . xv - location to put pointer to the first array
662: - yv - location to put pointer to the second array
664: Level: developer
666: .seealso: [](ch_vectors), `VecGetArray()`, `VecGetArrayRead()`, `VecGetArrayPair()`
667: @*/
668: static inline PetscErrorCode VecRestoreArrayPair(Vec x, Vec y, PetscScalar *xv[], PetscScalar *yv[])
669: {
670: PetscFunctionBegin;
671: PetscCall(VecRestoreArray(y, yv));
672: if (x != y) PetscCall(VecRestoreArrayRead(x, (const PetscScalar **)xv));
673: PetscFunctionReturn(PETSC_SUCCESS);
674: }
676: PETSC_EXTERN PetscErrorCode VecLockReadPush(Vec);
677: PETSC_EXTERN PetscErrorCode VecLockReadPop(Vec);
678: PETSC_EXTERN PetscErrorCode VecLockWriteSet(Vec, PetscBool);
679: PETSC_EXTERN PetscErrorCode VecLockGet(Vec, PetscInt *);
680: PETSC_EXTERN PetscErrorCode VecLockGetLocation(Vec, const char *[], const char *[], int *);
682: static inline PetscErrorCode VecSetErrorIfLocked(Vec x, PetscInt arg)
683: {
684: PetscInt state;
686: PetscFunctionBegin;
687: PetscCall(VecLockGet(x, &state));
688: if (PetscUnlikely(state != 0)) {
689: const char *file, *func, *name;
690: int line;
692: PetscCall(VecLockGetLocation(x, &file, &func, &line));
693: PetscCall(PetscObjectGetName((PetscObject)x, &name));
694: SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONGSTATE, "Vector '%s' (argument #%" PetscInt_FMT ") was locked for %s access in %s() at %s:%d (line numbers only accurate to function begin)", name, arg, state > 0 ? "read-only" : "write-only", func ? func : "unknown_function", file ? file : "unknown file", line);
695: }
696: PetscFunctionReturn(PETSC_SUCCESS);
697: }
699: /* The three are deprecated */
700: PETSC_DEPRECATED_FUNCTION(3, 11, 0, "VecLockReadPush()", ) static inline PetscErrorCode VecLockPush(Vec v)
701: {
702: return VecLockReadPush(v);
703: }
705: PETSC_DEPRECATED_FUNCTION(3, 11, 0, "VecLockReadPop()", ) static inline PetscErrorCode VecLockPop(Vec v)
706: {
707: return VecLockReadPop(v);
708: }
710: /*MC
711: VecLocked - Deprecated alias for `VecSetErrorIfLocked()`; raises an error if the `Vec` is currently locked for read
713: Synopsis:
714: #include <petscvec.h>
715: PetscErrorCode VecLocked(Vec x, int arg)
717: Not Collective; No Fortran Support
719: Input Parameters:
720: + x - the `Vec` to test
721: - arg - the argument position of `x` in the caller (used in the error message)
723: Level: deprecated
725: Note:
726: Use `VecSetErrorIfLocked()` in new code; this macro is retained only for backwards compatibility.
728: .seealso: `Vec`, `VecSetErrorIfLocked()`, `VecLockReadPush()`, `VecLockReadPop()`
729: M*/
730: #define VecLocked(x, arg) VecSetErrorIfLocked(x, arg) PETSC_DEPRECATED_MACRO(3, 11, 0, "VecSetErrorIfLocked()", )
732: /*E
733: VecOperation - Enumeration of overide-able methods in the `Vec` implementation function-table.
735: Values:
736: + `VECOP_DUPLICATE` - `VecDuplicate()`
737: . `VECOP_SET` - `VecSet()`
738: . `VECOP_VIEW` - `VecView()`
739: . `VECOP_LOAD` - `VecLoad()`
740: . `VECOP_VIEWNATIVE` - `VecViewNative()`
741: - `VECOP_LOADNATIVE` - `VecLoadNative()`
743: Level: advanced
745: Notes:
746: Some operations may serve as the implementation for other routines not listed above. For
747: example `VECOP_SET` can be used to simultaneously overriding the implementation used in
748: `VecSet()`, `VecSetInf()`, and `VecZeroEntries()`.
750: Entries to `VecOperation` are added as needed so if you do not see the operation listed which
751: you'd like to replace, please send mail to `petsc-maint@mcs.anl.gov`!
753: .seealso: [](ch_vectors), `Vec`, `VecSetOperation()`
754: E*/
755: typedef enum {
756: VECOP_DUPLICATE = 0,
757: VECOP_SET = 10,
758: VECOP_VIEW = 33,
759: VECOP_LOAD = 41,
760: VECOP_VIEWNATIVE = 68,
761: VECOP_LOADNATIVE = 69
762: } VecOperation;
763: PETSC_EXTERN PetscErrorCode VecSetOperation(Vec, VecOperation, PetscErrorCodeFn *);
765: /*
766: Routines for dealing with ghosted vectors:
767: vectors with ghost elements at the end of the array.
768: */
769: PETSC_EXTERN PetscErrorCode VecMPISetGhost(Vec, PetscInt, const PetscInt[]);
770: PETSC_EXTERN PetscErrorCode VecCreateGhost(MPI_Comm, PetscInt, PetscInt, PetscInt, const PetscInt[], Vec *);
771: PETSC_EXTERN PetscErrorCode VecCreateGhostWithArray(MPI_Comm, PetscInt, PetscInt, PetscInt, const PetscInt[], const PetscScalar[], Vec *);
772: PETSC_EXTERN PetscErrorCode VecCreateGhostBlock(MPI_Comm, PetscInt, PetscInt, PetscInt, PetscInt, const PetscInt[], Vec *);
773: PETSC_EXTERN PetscErrorCode VecCreateGhostBlockWithArray(MPI_Comm, PetscInt, PetscInt, PetscInt, PetscInt, const PetscInt[], const PetscScalar[], Vec *);
774: PETSC_EXTERN PetscErrorCode VecGhostGetGhostIS(Vec, IS *);
775: PETSC_EXTERN PetscErrorCode VecGhostGetLocalForm(Vec, Vec *);
776: PETSC_EXTERN PetscErrorCode VecGhostRestoreLocalForm(Vec, Vec *);
777: PETSC_EXTERN PetscErrorCode VecGhostIsLocalForm(Vec, Vec, PetscBool *);
778: PETSC_EXTERN PetscErrorCode VecGhostUpdateBegin(Vec, InsertMode, ScatterMode);
779: PETSC_EXTERN PetscErrorCode VecGhostUpdateEnd(Vec, InsertMode, ScatterMode);
781: PETSC_EXTERN PetscErrorCode VecConjugate(Vec);
782: PETSC_EXTERN PetscErrorCode VecImaginaryPart(Vec);
783: PETSC_EXTERN PetscErrorCode VecRealPart(Vec);
785: PETSC_EXTERN PetscErrorCode VecScatterCreateToAll(Vec, VecScatter *, Vec *);
786: PETSC_EXTERN PetscErrorCode VecScatterCreateToZero(Vec, VecScatter *, Vec *);
788: /* These vector calls were included from TAO. They miss vtable entries and GPU implementation */
789: PETSC_EXTERN PetscErrorCode ISComplementVec(IS, Vec, IS *);
790: PETSC_EXTERN PetscErrorCode VecPow(Vec, PetscScalar);
791: PETSC_EXTERN PetscErrorCode VecMedian(Vec, Vec, Vec, Vec);
792: PETSC_EXTERN PetscErrorCode VecWhichInactive(Vec, Vec, Vec, Vec, PetscBool, IS *);
793: PETSC_EXTERN PetscErrorCode VecWhichBetween(Vec, Vec, Vec, IS *);
794: PETSC_EXTERN PetscErrorCode VecWhichBetweenOrEqual(Vec, Vec, Vec, IS *);
795: PETSC_EXTERN PetscErrorCode VecWhichGreaterThan(Vec, Vec, IS *);
796: PETSC_EXTERN PetscErrorCode VecWhichLessThan(Vec, Vec, IS *);
797: PETSC_EXTERN PetscErrorCode VecWhichEqual(Vec, Vec, IS *);
798: PETSC_EXTERN PetscErrorCode VecISAXPY(Vec, IS, PetscScalar, Vec);
799: PETSC_EXTERN PetscErrorCode VecISCopy(Vec, IS, ScatterMode, Vec);
800: PETSC_EXTERN PetscErrorCode VecISSet(Vec, IS, PetscScalar);
801: PETSC_EXTERN PetscErrorCode VecISShift(Vec, IS, PetscScalar);
802: PETSC_EXTERN PetscErrorCode VecBoundGradientProjection(Vec, Vec, Vec, Vec, Vec);
803: PETSC_EXTERN PetscErrorCode VecStepBoundInfo(Vec, Vec, Vec, Vec, PetscReal *, PetscReal *, PetscReal *);
804: PETSC_EXTERN PetscErrorCode VecStepMax(Vec, Vec, PetscReal *);
805: PETSC_EXTERN PetscErrorCode VecStepMaxBounded(Vec, Vec, Vec, Vec, PetscReal *);
807: PETSC_EXTERN PetscErrorCode PetscViewerMathematicaGetVector(PetscViewer, Vec);
808: PETSC_EXTERN PetscErrorCode PetscViewerMathematicaPutVector(PetscViewer, Vec);
810: struct _n_Vecs {
811: PetscInt n;
812: Vec v;
813: };
814: /*S
815: Vecs - Collection of `Vec`s where the storage for the vectors is held in a single contiguous block of memory
817: Level: advanced
819: Notes:
820: Temporary construct for handling multiple right-hand side solves.
822: This is faked by storing a single `Vec` whose array is sized to hold `n` vectors back to back.
824: .seealso: `Vec`, `VecsCreateSeq()`, `VecsCreateSeqWithArray()`, `VecsDuplicate()`, `VecsDestroy()`
825: S*/
826: typedef struct _n_Vecs *Vecs;
827: PETSC_EXTERN PetscErrorCode VecsDestroy(Vecs);
828: PETSC_EXTERN PetscErrorCode VecsCreateSeq(MPI_Comm, PetscInt, PetscInt, Vecs *);
829: PETSC_EXTERN PetscErrorCode VecsCreateSeqWithArray(MPI_Comm, PetscInt, PetscInt, PetscScalar *, Vecs *);
830: PETSC_EXTERN PetscErrorCode VecsDuplicate(Vecs, Vecs *);
832: #if PetscDefined(HAVE_VIENNACL)
833: /*S
834: PetscViennaCLIndices - Opaque handle to an index buffer used by PETSc's ViennaCL `VECVIENNACL` vector backend to perform partial scatters between CPU and GPU memory
836: Level: developer
838: .seealso: `Vec`, `VECVIENNACL`, `VecCreateSeqViennaCL()`, `VecCreateMPIViennaCL()`, `VecViennaCLCopyToGPUSome_Public()`, `VecViennaCLCopyFromGPUSome_Public()`
839: S*/
840: typedef struct _p_PetscViennaCLIndices *PetscViennaCLIndices;
841: PETSC_EXTERN PetscErrorCode VecViennaCLCopyToGPUSome_Public(Vec, PetscViennaCLIndices);
842: PETSC_EXTERN PetscErrorCode VecViennaCLCopyFromGPUSome_Public(Vec, PetscViennaCLIndices);
843: PETSC_EXTERN PetscErrorCode VecCreateSeqViennaCL(MPI_Comm, PetscInt, Vec *);
844: PETSC_EXTERN PetscErrorCode VecCreateMPIViennaCL(MPI_Comm, PetscInt, PetscInt, Vec *);
845: #endif
846: #if PetscDefined(HAVE_CUDA) || PetscDefined(HAVE_HIP)
847: PETSC_EXTERN PetscErrorCode VecScatterInitializeForGPU(VecScatter, Vec);
848: PETSC_EXTERN PetscErrorCode VecScatterFinalizeForGPU(VecScatter);
849: #endif
850: #if PetscDefined(HAVE_KOKKOS_KERNELS)
851: PETSC_EXTERN PetscErrorCode VecCreateSeqKokkos(MPI_Comm, PetscInt, Vec *);
852: PETSC_EXTERN PetscErrorCode VecCreateSeqKokkosWithArray(MPI_Comm, PetscInt, PetscInt, const PetscScalar *, Vec *);
853: PETSC_EXTERN PetscErrorCode VecCreateMPIKokkos(MPI_Comm, PetscInt, PetscInt, Vec *);
854: PETSC_EXTERN PetscErrorCode VecCreateMPIKokkosWithArray(MPI_Comm, PetscInt, PetscInt, PetscInt, const PetscScalar *, Vec *);
855: PETSC_EXTERN PetscErrorCode VecKokkosPlaceArray(Vec, PetscScalar *);
856: PETSC_EXTERN PetscErrorCode VecKokkosResetArray(Vec);
857: #endif
859: PETSC_EXTERN PetscErrorCode VecNestGetSubVecs(Vec, PetscInt *, Vec *[]);
860: PETSC_EXTERN PetscErrorCode VecNestGetSubVecsRead(Vec, PetscInt *, Vec *[]);
861: PETSC_EXTERN PetscErrorCode VecNestRestoreSubVecsRead(Vec, PetscInt *, Vec *[]);
862: PETSC_EXTERN PetscErrorCode VecNestGetSubVec(Vec, PetscInt, Vec *);
863: PETSC_EXTERN PetscErrorCode VecNestSetSubVecs(Vec, PetscInt, PetscInt *, Vec *);
864: PETSC_EXTERN PetscErrorCode VecNestSetSubVec(Vec, PetscInt, Vec);
865: PETSC_EXTERN PetscErrorCode VecCreateNest(MPI_Comm, PetscInt, IS *, Vec *, Vec *);
866: PETSC_EXTERN PetscErrorCode VecNestGetSize(Vec, PetscInt *);
868: PETSC_EXTERN PetscErrorCode PetscOptionsGetVec(PetscOptions, const char[], const char[], Vec, PetscBool *);
869: PETSC_EXTERN PetscErrorCode VecFilter(Vec, PetscReal);
870: PETSC_DEPRECATED_FUNCTION(3, 20, 0, "VecFilter()", ) static inline PetscErrorCode VecChop(Vec v, PetscReal tol)
871: {
872: return VecFilter(v, tol);
873: }
875: PETSC_EXTERN PetscErrorCode VecGetLayout(Vec, PetscLayout *);
876: PETSC_EXTERN PetscErrorCode VecSetLayout(Vec, PetscLayout);
878: PETSC_EXTERN PetscErrorCode PetscSectionVecView(PetscSection, Vec, PetscViewer);
879: PETSC_EXTERN PetscErrorCode VecGetValuesSection(Vec, PetscSection, PetscInt, PetscScalar *[]);
880: PETSC_EXTERN PetscErrorCode VecSetValuesSection(Vec, PetscSection, PetscInt, const PetscScalar[], InsertMode);
881: PETSC_EXTERN PetscErrorCode PetscSectionVecNorm(PetscSection, PetscSection, Vec, NormType, PetscReal[]);
883: /*S
884: VecTagger - Object used to manage the tagging of a subset of indices based on the values of a vector. The
885: motivating application is the selection of cells for refinement or coarsening based on vector containing
886: the values in an error indicator metric.
888: Values:
889: + `VECTAGGERABSOLUTE` - "absolute" values are in a interval (box for complex values) of explicitly defined values
890: . `VECTAGGERRELATIVE` - "relative" values are in a interval (box for complex values) of values relative to the set of all values in the vector
891: . `VECTAGGERCDF` - "cdf" values are in a relative range of the *cumulative distribution* of values in the vector
892: . `VECTAGGEROR` - "or" values are in the union of other tags
893: - `VECTAGGERAND` - "and" values are in the intersection of other tags
895: Level: advanced
897: Developer Note:
898: Why not use a `DMLabel` or similar object
900: .seealso: [](ch_vectors), `Vec`, `VecTaggerType`, `VecTaggerCreate()`
901: S*/
902: typedef struct _p_VecTagger *VecTagger;
904: /*J
905: VecTaggerType - String with the name of a `VecTagger` type
907: Level: advanced
909: .seealso: [](ch_vectors), `Vec`, `VecTagger`, `VecTaggerCreate()`
910: J*/
911: typedef const char *VecTaggerType;
912: #define VECTAGGERABSOLUTE "absolute"
913: #define VECTAGGERRELATIVE "relative"
914: #define VECTAGGERCDF "cdf"
915: #define VECTAGGEROR "or"
916: #define VECTAGGERAND "and"
918: PETSC_EXTERN PetscClassId VEC_TAGGER_CLASSID;
919: PETSC_EXTERN PetscFunctionList VecTaggerList;
920: PETSC_EXTERN PetscErrorCode VecTaggerRegister(const char[], PetscErrorCode (*)(VecTagger));
921: PETSC_EXTERN PetscErrorCode VecTaggerRegisterAll(void);
923: PETSC_EXTERN PetscErrorCode VecTaggerCreate(MPI_Comm, VecTagger *);
924: PETSC_EXTERN PetscErrorCode VecTaggerSetBlockSize(VecTagger, PetscInt);
925: PETSC_EXTERN PetscErrorCode VecTaggerGetBlockSize(VecTagger, PetscInt *);
926: PETSC_EXTERN PetscErrorCode VecTaggerSetType(VecTagger, VecTaggerType);
927: PETSC_EXTERN PetscErrorCode VecTaggerGetType(VecTagger, VecTaggerType *);
928: PETSC_EXTERN PetscErrorCode VecTaggerSetInvert(VecTagger, PetscBool);
929: PETSC_EXTERN PetscErrorCode VecTaggerGetInvert(VecTagger, PetscBool *);
930: PETSC_EXTERN PetscErrorCode VecTaggerSetFromOptions(VecTagger);
931: PETSC_EXTERN PetscErrorCode VecTaggerSetUp(VecTagger);
932: PETSC_EXTERN PetscErrorCode VecTaggerView(VecTagger, PetscViewer);
933: PETSC_EXTERN PetscErrorCode VecTaggerComputeIS(VecTagger, Vec, IS *, PetscBool *);
934: PETSC_EXTERN PetscErrorCode VecTaggerDestroy(VecTagger *);
936: /*S
937: VecTaggerBox - A interval (box for complex numbers) range used to tag values. For real scalars, this is just a closed interval; for complex scalars,
938: the box is the closed region in the complex plane such that real(min) <= real(z) <= real(max) and imag(min) <= imag(z) <= imag(max). `INF` is an acceptable endpoint.
940: Level: beginner
942: .seealso: [](ch_vectors), `Vec`, `VecTagger`, `VecTaggerType`, `VecTaggerCreate()`, `VecTaggerComputeIntervals()`
943: S*/
944: typedef struct {
945: PetscScalar min;
946: PetscScalar max;
947: } VecTaggerBox;
948: PETSC_EXTERN PetscErrorCode VecTaggerComputeBoxes(VecTagger, Vec, PetscInt *, VecTaggerBox *[], PetscBool *);
950: PETSC_EXTERN PetscErrorCode VecTaggerAbsoluteSetBox(VecTagger, VecTaggerBox *);
951: PETSC_EXTERN PetscErrorCode VecTaggerAbsoluteGetBox(VecTagger, const VecTaggerBox *[]);
953: PETSC_EXTERN PetscErrorCode VecTaggerRelativeSetBox(VecTagger, VecTaggerBox *);
954: PETSC_EXTERN PetscErrorCode VecTaggerRelativeGetBox(VecTagger, const VecTaggerBox *[]);
956: PETSC_EXTERN PetscErrorCode VecTaggerCDFSetBox(VecTagger, VecTaggerBox *);
957: PETSC_EXTERN PetscErrorCode VecTaggerCDFGetBox(VecTagger, const VecTaggerBox *[]);
959: /*E
960: VecTaggerCDFMethod - Determines what method is used to compute absolute values from cumulative distribution values (e.g., what value is the preimage of .95 in the cdf).
962: Values:
963: + `VECTAGGER_CDF_GATHER` - gather the data to MPI rank 0, perform the computation and broadcast the result
964: - `VECTAGGER_CDF_ITERATIVE` - compute the results on all ranks iteratively using `MPI_Allreduce()`
966: Level: advanced
968: Note:
969: Relevant only in parallel: in serial it is directly computed.
971: Developer Note:
972: In PETSc enums of this type are referred to with the term type, not method.
974: .seealso: [](ch_vectors), `Vec`, `VecTagger`, `VecTaggerType`, `VecTaggerCreate()`, `VecTaggerCDFSetMethod()`
975: E*/
976: typedef enum {
977: VECTAGGER_CDF_GATHER,
978: VECTAGGER_CDF_ITERATIVE,
979: VECTAGGER_CDF_NUM_METHODS
980: } VecTaggerCDFMethod;
981: PETSC_EXTERN const char *const VecTaggerCDFMethods[];
983: PETSC_EXTERN PetscErrorCode VecTaggerCDFSetMethod(VecTagger, VecTaggerCDFMethod);
984: PETSC_EXTERN PetscErrorCode VecTaggerCDFGetMethod(VecTagger, VecTaggerCDFMethod *);
985: PETSC_EXTERN PetscErrorCode VecTaggerCDFIterativeSetTolerances(VecTagger, PetscInt, PetscReal, PetscReal);
986: PETSC_EXTERN PetscErrorCode VecTaggerCDFIterativeGetTolerances(VecTagger, PetscInt *, PetscReal *, PetscReal *);
988: PETSC_EXTERN PetscErrorCode VecTaggerOrSetSubs(VecTagger, PetscInt, VecTagger *, PetscCopyMode);
989: PETSC_EXTERN PetscErrorCode VecTaggerOrGetSubs(VecTagger, PetscInt *, VecTagger *[]);
991: PETSC_EXTERN PetscErrorCode VecTaggerAndSetSubs(VecTagger, PetscInt, VecTagger *, PetscCopyMode);
992: PETSC_EXTERN PetscErrorCode VecTaggerAndGetSubs(VecTagger, PetscInt *, VecTagger *[]);
994: PETSC_EXTERN PetscErrorCode VecTaggerInitializePackage(void);
995: PETSC_EXTERN PetscErrorCode VecTaggerFinalizePackage(void);
997: PETSC_EXTERN PetscErrorCode PetscViewerMathematicaSetVector(PetscViewer, Vec);
998: PETSC_EXTERN PetscErrorCode PetscViewerMathematicaGetVector(PetscViewer, Vec);
1000: PETSC_EXTERN PetscErrorCode VecFischer(Vec, Vec, Vec, Vec, Vec);
1001: PETSC_EXTERN PetscErrorCode VecSFischer(Vec, Vec, Vec, Vec, PetscReal, Vec);
1003: #if PetscDefined(USE_DEBUG)
1004: /* This is an internal debug-only routine that should not be used by users */
1005: PETSC_SINGLE_LIBRARY_INTERN PetscErrorCode VecValidValues_Internal(Vec, PetscInt, PetscBool);
1006: #else
1007: #define VecValidValues_Internal(...) PETSC_SUCCESS
1008: #endif /* PETSC_USE_DEBUG */
1010: #define VEC_CUPM_NOT_CONFIGURED(impl) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP_SYS, "Must configure PETSc with " PetscStringize(impl) " support to use %s", PETSC_FUNCTION_NAME)
1012: #if PetscDefined(HAVE_CUDA)
1013: #define VEC_CUDA_DECL_OR_STUB(__decl__, ...) PETSC_EXTERN __decl__;
1014: #else
1015: #define VEC_CUDA_DECL_OR_STUB(__decl__, ...) \
1016: static inline __decl__ \
1017: { \
1018: __VA_ARGS__; \
1019: VEC_CUPM_NOT_CONFIGURED(cuda); \
1020: }
1021: #endif /* PETSC_HAVE_CUDA */
1023: /* extra underscore here to make it line up with the cuda versions */
1024: #if PetscDefined(HAVE_HIP)
1025: #define VEC_HIP__DECL_OR_STUB(__decl__, ...) PETSC_EXTERN __decl__;
1026: #else
1027: #define VEC_HIP__DECL_OR_STUB(__decl__, ...) \
1028: static inline __decl__ \
1029: { \
1030: __VA_ARGS__; \
1031: VEC_CUPM_NOT_CONFIGURED(hip); \
1032: }
1033: #endif /* PETSC_HAVE_HIP */
1035: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCreateSeqCUDA(MPI_Comm a, PetscInt b, Vec *c), (void)a, (void)b, (void)c)
1036: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecCreateSeqHIP(MPI_Comm a, PetscInt b, Vec *c), (void)a, (void)b, (void)c)
1038: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCreateSeqCUDAWithArray(MPI_Comm a, PetscInt b, PetscInt c, const PetscScalar *d, Vec *e), (void)a, (void)b, (void)c, (void)d, (void)e)
1039: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecCreateSeqHIPWithArray(MPI_Comm a, PetscInt b, PetscInt c, const PetscScalar *d, Vec *e), (void)a, (void)b, (void)c, (void)d, (void)e)
1041: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCreateSeqCUDAWithArrays(MPI_Comm a, PetscInt b, PetscInt c, const PetscScalar *d, const PetscScalar *e, Vec *f), (void)a, (void)b, (void)c, (void)d, (void)e, (void)f)
1042: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecCreateSeqHIPWithArrays(MPI_Comm a, PetscInt b, PetscInt c, const PetscScalar *d, const PetscScalar *e, Vec *f), (void)a, (void)b, (void)c, (void)d, (void)e, (void)f)
1044: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCreateMPICUDA(MPI_Comm a, PetscInt b, PetscInt c, Vec *d), (void)a, (void)b, (void)c, (void)d)
1045: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecCreateMPIHIP(MPI_Comm a, PetscInt b, PetscInt c, Vec *d), (void)a, (void)b, (void)c, (void)d)
1047: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCreateMPICUDAWithArray(MPI_Comm a, PetscInt b, PetscInt c, PetscInt d, const PetscScalar *e, Vec *f), (void)a, (void)b, (void)c, (void)d, (void)e, (void)f)
1048: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecCreateMPIHIPWithArray(MPI_Comm a, PetscInt b, PetscInt c, PetscInt d, const PetscScalar *e, Vec *f), (void)a, (void)b, (void)c, (void)d, (void)e, (void)f)
1050: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCreateMPICUDAWithArrays(MPI_Comm a, PetscInt b, PetscInt c, PetscInt d, const PetscScalar *e, const PetscScalar *f, Vec *g), (void)a, (void)b, (void)c, (void)d, (void)e, (void)f, (void)g)
1051: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecCreateMPIHIPWithArrays(MPI_Comm a, PetscInt b, PetscInt c, PetscInt d, const PetscScalar *e, const PetscScalar *f, Vec *g), (void)a, (void)b, (void)c, (void)d, (void)e, (void)f, (void)g)
1053: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDAGetArray(Vec a, PetscScalar **b), (void)a, (void)b)
1054: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPGetArray(Vec a, PetscScalar **b), (void)a, (void)b)
1056: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDARestoreArray(Vec a, PetscScalar **b), (void)a, (void)b)
1057: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPRestoreArray(Vec a, PetscScalar **b), (void)a, (void)b)
1059: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDAGetArrayRead(Vec a, const PetscScalar **b), (void)a, (void)b)
1060: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPGetArrayRead(Vec a, const PetscScalar **b), (void)a, (void)b)
1062: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDARestoreArrayRead(Vec a, const PetscScalar **b), (void)a, (void)b)
1063: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPRestoreArrayRead(Vec a, const PetscScalar **b), (void)a, (void)b)
1065: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDAGetArrayWrite(Vec a, PetscScalar **b), (void)a, (void)b)
1066: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPGetArrayWrite(Vec a, PetscScalar **b), (void)a, (void)b)
1068: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDARestoreArrayWrite(Vec a, PetscScalar **b), (void)a, (void)b)
1069: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPRestoreArrayWrite(Vec a, PetscScalar **b), (void)a, (void)b)
1071: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDAPlaceArray(Vec a, const PetscScalar b[]), (void)a, (void)b)
1072: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPPlaceArray(Vec a, const PetscScalar b[]), (void)a, (void)b)
1074: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDAReplaceArray(Vec a, const PetscScalar b[]), (void)a, (void)b)
1075: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPReplaceArray(Vec a, const PetscScalar b[]), (void)a, (void)b)
1077: VEC_CUDA_DECL_OR_STUB(PetscErrorCode VecCUDAResetArray(Vec a), (void)a)
1078: VEC_HIP__DECL_OR_STUB(PetscErrorCode VecHIPResetArray(Vec a), (void)a)
1080: #undef VEC_CUPM_NOT_CONFIGURED
1081: #undef VEC_CUDA_DECL_OR_STUB
1082: #undef VEC_HIP__DECL_OR_STUB