Actual source code: pvec2.c
1: /*
2: Code for some of the parallel vector primitives.
3: */
4: #include <../src/vec/vec/impls/mpi/pvecimpl.h>
5: #include <petscblaslapack.h>
7: PetscErrorCode VecDot_MPI(Vec xin, Vec yin, PetscScalar *z)
8: {
9: PetscFunctionBegin;
10: PetscCall(VecXDot_MPI_Default(xin, yin, z, VecDot_Seq));
11: PetscFunctionReturn(PETSC_SUCCESS);
12: }
14: PetscErrorCode VecTDot_MPI(Vec xin, Vec yin, PetscScalar *z)
15: {
16: PetscFunctionBegin;
17: PetscCall(VecXDot_MPI_Default(xin, yin, z, VecTDot_Seq));
18: PetscFunctionReturn(PETSC_SUCCESS);
19: }
21: PetscErrorCode VecMDot_MPI(Vec xin, PetscInt nv, const Vec y[], PetscScalar *z)
22: {
23: PetscFunctionBegin;
24: PetscCall(VecMXDot_MPI_Default(xin, nv, y, z, VecMDot_Seq));
25: PetscFunctionReturn(PETSC_SUCCESS);
26: }
28: PetscErrorCode VecMTDot_MPI(Vec xin, PetscInt nv, const Vec y[], PetscScalar *z)
29: {
30: PetscFunctionBegin;
31: PetscCall(VecMXDot_MPI_Default(xin, nv, y, z, VecMTDot_Seq));
32: PetscFunctionReturn(PETSC_SUCCESS);
33: }
35: PetscErrorCode VecNorm_MPI(Vec xin, NormType type, PetscReal *z)
36: {
37: PetscFunctionBegin;
38: PetscCall(VecNorm_MPI_Default(xin, type, z, VecNorm_Seq));
39: PetscFunctionReturn(PETSC_SUCCESS);
40: }
42: PetscErrorCode VecMax_MPI(Vec xin, PetscInt *idx, PetscReal *z)
43: {
44: const MPI_Op ops[] = {MPIU_MAXLOC, MPIU_MAX};
46: PetscFunctionBegin;
47: PetscCall(VecMinMax_MPI_Default(xin, idx, z, VecMax_Seq, ops));
48: PetscFunctionReturn(PETSC_SUCCESS);
49: }
51: PetscErrorCode VecMin_MPI(Vec xin, PetscInt *idx, PetscReal *z)
52: {
53: const MPI_Op ops[] = {MPIU_MINLOC, MPIU_MIN};
55: PetscFunctionBegin;
56: PetscCall(VecMinMax_MPI_Default(xin, idx, z, VecMin_Seq, ops));
57: PetscFunctionReturn(PETSC_SUCCESS);
58: }