| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | /* | ||
| 2 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 3 | SLEPc - Scalable Library for Eigenvalue Problem Computations | ||
| 4 | Copyright (c) 2002-, Universitat Politecnica de Valencia, Spain | ||
| 5 | |||
| 6 | This file is part of SLEPc. | ||
| 7 | SLEPc is distributed under a 2-clause BSD license (see LICENSE). | ||
| 8 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 9 | */ | ||
| 10 | |||
| 11 | #include <petsc/private/vecimpl.h> | ||
| 12 | |||
| 13 | #if defined(__WITH_MPI__) | ||
| 14 | #define __SUF__(A) A##_MPI | ||
| 15 | #else | ||
| 16 | #define __SUF__(A) A##_Seq | ||
| 17 | #endif | ||
| 18 | #define __QUOTEME_(x) #x | ||
| 19 | #define __QUOTEME(x) __QUOTEME_(x) | ||
| 20 | #define __SUF_C__(A) __QUOTEME(__SUF__(A)) | ||
| 21 | |||
| 22 | 1078350 | static PetscErrorCode __SUF__(VecDot_Comp)(Vec a,Vec b,PetscScalar *z) | |
| 23 | { | ||
| 24 | 1078350 | PetscScalar sum = 0.0,work; | |
| 25 | 1078350 | PetscInt i; | |
| 26 | 1078350 | Vec_Comp *as = (Vec_Comp*)a->data,*bs = (Vec_Comp*)b->data; | |
| 27 | |||
| 28 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1078350 | PetscFunctionBegin; |
| 29 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1078350 | SlepcValidVecComp(a,1); |
| 30 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1078350 | SlepcValidVecComp(b,2); |
| 31 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
1078350 | if (as->x[0]->ops->dot_local) { |
| 32 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
2159456 | for (i=0,sum=0.0;i<as->n->n;i++) { |
| 33 |
5/10✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
1081106 | PetscUseTypeMethod(as->x[i],dot_local,bs->x[i],&work); |
| 34 | 1081106 | sum += work; | |
| 35 | } | ||
| 36 | #if defined(__WITH_MPI__) | ||
| 37 | 1078350 | work = sum; | |
| 38 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
1078350 | PetscCallMPI(MPIU_Allreduce(&work,&sum,1,MPIU_SCALAR,MPIU_SUM,PetscObjectComm((PetscObject)a))); |
| 39 | #endif | ||
| 40 | } else { | ||
| 41 | ✗ | for (i=0,sum=0.0;i<as->n->n;i++) { | |
| 42 | ✗ | PetscCall(VecDot(as->x[i],bs->x[i],&work)); | |
| 43 | ✗ | sum += work; | |
| 44 | } | ||
| 45 | } | ||
| 46 | 1078350 | *z = sum; | |
| 47 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
1078350 | PetscFunctionReturn(PETSC_SUCCESS); |
| 48 | } | ||
| 49 | |||
| 50 | 798034 | static PetscErrorCode __SUF__(VecMDot_Comp)(Vec a,PetscInt n,const Vec b[],PetscScalar *z) | |
| 51 | { | ||
| 52 | 798034 | PetscScalar *work,*work0,*r; | |
| 53 | 798034 | Vec_Comp *as = (Vec_Comp*)a->data; | |
| 54 | 798034 | Vec *bx; | |
| 55 | 798034 | PetscInt i,j; | |
| 56 | |||
| 57 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
798034 | PetscFunctionBegin; |
| 58 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
798034 | SlepcValidVecComp(a,1); |
| 59 |
4/6✓ Branch 0 taken 8 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 2 times.
|
2761400 | for (i=0;i<n;i++) SlepcValidVecComp(b[i],3); |
| 60 | |||
| 61 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
798034 | if (as->n->n == 0) { |
| 62 | ✗ | *z = 0; | |
| 63 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 64 | } | ||
| 65 | |||
| 66 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
798034 | PetscCall(PetscMalloc2(n,&work0,n,&bx)); |
| 67 | |||
| 68 | #if defined(__WITH_MPI__) | ||
| 69 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
798034 | if (as->x[0]->ops->mdot_local) { |
| 70 | 798034 | r = work0; work = z; | |
| 71 | } else | ||
| 72 | #endif | ||
| 73 | { | ||
| 74 | ✗ | r = z; work = work0; | |
| 75 | } | ||
| 76 | |||
| 77 | /* z[i] <- a.x' * b[i].x */ | ||
| 78 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
2761400 | for (i=0;i<n;i++) r[i] = 0.0; |
| 79 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1597820 | for (j=0;j<as->n->n;j++) { |
| 80 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
2766656 | for (i=0;i<n;i++) bx[i] = ((Vec_Comp*)b[i]->data)->x[j]; |
| 81 |
6/12✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 8 times.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
799786 | if (as->x[0]->ops->mdot_local) PetscUseTypeMethod(as->x[j],mdot_local,n,bx,work); |
| 82 | ✗ | else PetscCall(VecMDot(as->x[j],n,bx,work)); | |
| 83 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
2766656 | for (i=0;i<n;i++) r[i] += work[i]; |
| 84 | } | ||
| 85 | |||
| 86 | /* If def(__WITH_MPI__) and exists mdot_local */ | ||
| 87 | #if defined(__WITH_MPI__) | ||
| 88 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
798034 | if (as->x[0]->ops->mdot_local) { |
| 89 | /* z[i] <- Allreduce(work[i]) */ | ||
| 90 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
798034 | PetscCallMPI(MPIU_Allreduce(r,z,n,MPIU_SCALAR,MPIU_SUM,PetscObjectComm((PetscObject)a))); |
| 91 | } | ||
| 92 | #endif | ||
| 93 | |||
| 94 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
798034 | PetscCall(PetscFree2(work0,bx)); |
| 95 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
142278 | PetscFunctionReturn(PETSC_SUCCESS); |
| 96 | } | ||
| 97 | |||
| 98 | 60 | static PetscErrorCode __SUF__(VecTDot_Comp)(Vec a,Vec b,PetscScalar *z) | |
| 99 | { | ||
| 100 | 60 | PetscScalar sum = 0.0,work; | |
| 101 | 60 | PetscInt i; | |
| 102 | 60 | Vec_Comp *as = (Vec_Comp*)a->data,*bs = (Vec_Comp*)b->data; | |
| 103 | |||
| 104 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
60 | PetscFunctionBegin; |
| 105 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
60 | SlepcValidVecComp(a,1); |
| 106 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
60 | SlepcValidVecComp(b,2); |
| 107 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
60 | if (as->x[0]->ops->tdot_local) { |
| 108 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
180 | for (i=0,sum=0.0;i<as->n->n;i++) { |
| 109 |
5/10✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
120 | PetscUseTypeMethod(as->x[i],tdot_local,bs->x[i],&work); |
| 110 | 120 | sum += work; | |
| 111 | } | ||
| 112 | #if defined(__WITH_MPI__) | ||
| 113 | 60 | work = sum; | |
| 114 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
60 | PetscCallMPI(MPIU_Allreduce(&work,&sum,1,MPIU_SCALAR,MPIU_SUM,PetscObjectComm((PetscObject)a))); |
| 115 | #endif | ||
| 116 | } else { | ||
| 117 | ✗ | for (i=0,sum=0.0;i<as->n->n;i++) { | |
| 118 | ✗ | PetscCall(VecTDot(as->x[i],bs->x[i],&work)); | |
| 119 | ✗ | sum += work; | |
| 120 | } | ||
| 121 | } | ||
| 122 | 60 | *z = sum; | |
| 123 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
60 | PetscFunctionReturn(PETSC_SUCCESS); |
| 124 | } | ||
| 125 | |||
| 126 | 60 | static PetscErrorCode __SUF__(VecMTDot_Comp)(Vec a,PetscInt n,const Vec b[],PetscScalar *z) | |
| 127 | { | ||
| 128 | 60 | PetscScalar *work,*work0,*r; | |
| 129 | 60 | Vec_Comp *as = (Vec_Comp*)a->data; | |
| 130 | 60 | Vec *bx; | |
| 131 | 60 | PetscInt i,j; | |
| 132 | |||
| 133 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
60 | PetscFunctionBegin; |
| 134 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
60 | SlepcValidVecComp(a,1); |
| 135 |
4/6✓ Branch 0 taken 8 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 2 times.
|
180 | for (i=0;i<n;i++) SlepcValidVecComp(b[i],3); |
| 136 | |||
| 137 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
60 | if (as->n->n == 0) { |
| 138 | ✗ | *z = 0; | |
| 139 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 140 | } | ||
| 141 | |||
| 142 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
60 | PetscCall(PetscMalloc2(n,&work0,n,&bx)); |
| 143 | |||
| 144 | #if defined(__WITH_MPI__) | ||
| 145 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
60 | if (as->x[0]->ops->mtdot_local) { |
| 146 | 60 | r = work0; work = z; | |
| 147 | } else | ||
| 148 | #endif | ||
| 149 | { | ||
| 150 | ✗ | r = z; work = work0; | |
| 151 | } | ||
| 152 | |||
| 153 | /* z[i] <- a.x' * b[i].x */ | ||
| 154 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
180 | for (i=0;i<n;i++) r[i] = 0.0; |
| 155 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
180 | for (j=0;j<as->n->n;j++) { |
| 156 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
360 | for (i=0;i<n;i++) bx[i] = ((Vec_Comp*)b[i]->data)->x[j]; |
| 157 |
6/12✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 8 times.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
120 | if (as->x[0]->ops->mtdot_local) PetscUseTypeMethod(as->x[j],mtdot_local,n,bx,work); |
| 158 | ✗ | else PetscCall(VecMTDot(as->x[j],n,bx,work)); | |
| 159 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
360 | for (i=0;i<n;i++) r[i] += work[i]; |
| 160 | } | ||
| 161 | |||
| 162 | /* If def(__WITH_MPI__) and exists mtdot_local */ | ||
| 163 | #if defined(__WITH_MPI__) | ||
| 164 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
60 | if (as->x[0]->ops->mtdot_local) { |
| 165 | /* z[i] <- Allreduce(work[i]) */ | ||
| 166 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
60 | PetscCallMPI(MPIU_Allreduce(r,z,n,MPIU_SCALAR,MPIU_SUM,PetscObjectComm((PetscObject)a))); |
| 167 | } | ||
| 168 | #endif | ||
| 169 | |||
| 170 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
60 | PetscCall(PetscFree2(work0,bx)); |
| 171 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
12 | PetscFunctionReturn(PETSC_SUCCESS); |
| 172 | } | ||
| 173 | |||
| 174 | 842988 | static PetscErrorCode __SUF__(VecNorm_Comp)(Vec a,NormType t,PetscReal *norm) | |
| 175 | { | ||
| 176 | 842988 | PetscReal work[3],s=0.0; | |
| 177 | 842988 | Vec_Comp *as = (Vec_Comp*)a->data; | |
| 178 | 842988 | PetscInt i; | |
| 179 | |||
| 180 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
842988 | PetscFunctionBegin; |
| 181 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
842988 | SlepcValidVecComp(a,1); |
| 182 | /* Initialize norm */ | ||
| 183 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 9 times.
|
842988 | switch (t) { |
| 184 | 240 | case NORM_1: case NORM_INFINITY: *norm = 0.0; break; | |
| 185 | 842688 | case NORM_2: case NORM_FROBENIUS: *norm = 1.0; s = 0.0; break; | |
| 186 | 60 | case NORM_1_AND_2: norm[0] = 0.0; norm[1] = 1.0; s = 0.0; break; | |
| 187 | } | ||
| 188 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1688640 | for (i=0;i<as->n->n;i++) { |
| 189 |
6/12✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 8 times.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
845652 | if (as->x[0]->ops->norm_local) PetscUseTypeMethod(as->x[i],norm_local,t,work); |
| 190 | ✗ | else PetscCall(VecNorm(as->x[i],t,work)); | |
| 191 | /* norm+= work */ | ||
| 192 |
4/5✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
|
845652 | switch (t) { |
| 193 | 120 | case NORM_1: *norm+= *work; break; | |
| 194 | 845052 | case NORM_2: case NORM_FROBENIUS: AddNorm2(norm,&s,*work); break; | |
| 195 | 120 | case NORM_1_AND_2: norm[0]+= work[0]; AddNorm2(&norm[1],&s,work[1]); break; | |
| 196 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
480 | case NORM_INFINITY: *norm = PetscMax(*norm,*work); break; |
| 197 | } | ||
| 198 | 303400 | } | |
| 199 | |||
| 200 | /* If def(__WITH_MPI__) and exists norm_local */ | ||
| 201 | #if defined(__WITH_MPI__) | ||
| 202 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
842988 | if (as->x[0]->ops->norm_local) { |
| 203 | 842988 | PetscReal work0[3]; | |
| 204 | /* norm <- Allreduce(work) */ | ||
| 205 |
4/5✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
|
842988 | switch (t) { |
| 206 | 60 | case NORM_1: | |
| 207 | 60 | work[0] = *norm; | |
| 208 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
60 | PetscCallMPI(MPIU_Allreduce(work,norm,1,MPIU_REAL,MPIU_SUM,PetscObjectComm((PetscObject)a))); |
| 209 | break; | ||
| 210 | 842688 | case NORM_2: case NORM_FROBENIUS: | |
| 211 | 842688 | work[0] = *norm; work[1] = s; | |
| 212 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
842688 | PetscCallMPI(MPIU_Allreduce(work,work0,1,MPIU_NORM2,MPIU_NORM2_SUM,PetscObjectComm((PetscObject)a))); |
| 213 | 842688 | *norm = GetNorm2(work0[0],work0[1]); | |
| 214 | 842688 | break; | |
| 215 | 60 | case NORM_1_AND_2: | |
| 216 | 60 | work[0] = norm[0]; work[1] = norm[1]; work[2] = s; | |
| 217 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
60 | PetscCallMPI(MPIU_Allreduce(work,work0,1,MPIU_NORM1_AND_2,MPIU_NORM2_SUM,PetscObjectComm((PetscObject)a))); |
| 218 | 60 | norm[0] = work0[0]; | |
| 219 | 60 | norm[1] = GetNorm2(work0[1],work0[2]); | |
| 220 | 60 | break; | |
| 221 | 180 | case NORM_INFINITY: | |
| 222 | 180 | work[0] = *norm; | |
| 223 |
28/58✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
|
180 | PetscCallMPI(MPIU_Allreduce(work,norm,1,MPIU_REAL,MPIU_MAX,PetscObjectComm((PetscObject)a))); |
| 224 | break; | ||
| 225 | } | ||
| 226 | 303304 | } | |
| 227 | #else | ||
| 228 | /* Norm correction */ | ||
| 229 | ✗ | switch (t) { | |
| 230 | ✗ | case NORM_2: case NORM_FROBENIUS: *norm = GetNorm2(*norm,s); break; | |
| 231 | ✗ | case NORM_1_AND_2: norm[1] = GetNorm2(norm[1],s); break; | |
| 232 | ✗ | default: ; | |
| 233 | } | ||
| 234 | #endif | ||
| 235 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
151408 | PetscFunctionReturn(PETSC_SUCCESS); |
| 236 | } | ||
| 237 | |||
| 238 | 8210 | PetscErrorCode __SUF__(VecDotNorm2_Comp)(Vec v,Vec w,PetscScalar *dp,PetscScalar *nm) | |
| 239 | { | ||
| 240 | 8210 | PetscScalar dp0=0.0,nm0=0.0,dp1=0.0,nm1=0.0; | |
| 241 | 8210 | const PetscScalar *vx,*wx; | |
| 242 | 8210 | Vec_Comp *vs = (Vec_Comp*)v->data,*ws = (Vec_Comp*)w->data; | |
| 243 | 8210 | PetscInt i,n; | |
| 244 | 8210 | PetscBool t0,t1; | |
| 245 | #if defined(__WITH_MPI__) | ||
| 246 | 4090 | PetscScalar work[4]; | |
| 247 | #endif | ||
| 248 | |||
| 249 | 8210 | PetscFunctionBegin; | |
| 250 | /* Compute recursively the local part */ | ||
| 251 | 8210 | PetscCall(PetscObjectTypeCompare((PetscObject)v,VECCOMP,&t0)); | |
| 252 | 8210 | PetscCall(PetscObjectTypeCompare((PetscObject)w,VECCOMP,&t1)); | |
| 253 | 8210 | if (t0 && t1) { | |
| 254 | 823 | SlepcValidVecComp(v,1); | |
| 255 | 823 | SlepcValidVecComp(w,2); | |
| 256 | 8210 | for (i=0;i<vs->n->n;i++) { | |
| 257 | 4120 | PetscCall(VecDotNorm2_Comp_Seq(vs->x[i],ws->x[i],&dp1,&nm1)); | |
| 258 | 4120 | dp0 += dp1; | |
| 259 | 4120 | nm0 += nm1; | |
| 260 | } | ||
| 261 | 4120 | } else if (!t0 && !t1) { | |
| 262 | 4120 | PetscCall(VecGetLocalSize(v,&n)); | |
| 263 | 4120 | PetscCall(VecGetArrayRead(v,&vx)); | |
| 264 | 4120 | PetscCall(VecGetArrayRead(w,&wx)); | |
| 265 | 450880 | for (i=0;i<n;i++) { | |
| 266 | 446760 | dp0 += vx[i]*PetscConj(wx[i]); | |
| 267 | 446760 | nm0 += wx[i]*PetscConj(wx[i]); | |
| 268 | } | ||
| 269 | 4120 | PetscCall(VecRestoreArrayRead(v,&vx)); | |
| 270 | 4120 | PetscCall(VecRestoreArrayRead(w,&wx)); | |
| 271 | ✗ | } else SETERRQ(PetscObjectComm((PetscObject)v),PETSC_ERR_ARG_INCOMP,"Incompatible vector types"); | |
| 272 | |||
| 273 | #if defined(__WITH_MPI__) | ||
| 274 | /* [dp, nm] <- Allreduce([dp0, nm0]) */ | ||
| 275 | 4090 | work[0] = dp0; work[1] = nm0; | |
| 276 | 4090 | PetscCallMPI(MPIU_Allreduce(work,&work[2],2,MPIU_SCALAR,MPIU_SUM,PetscObjectComm((PetscObject)v))); | |
| 277 | 4090 | *dp = work[2]; *nm = work[3]; | |
| 278 | #else | ||
| 279 | 4120 | *dp = dp0; *nm = nm0; | |
| 280 | #endif | ||
| 281 | 8210 | PetscFunctionReturn(PETSC_SUCCESS); | |
| 282 | } | ||
| 283 | |||
| 284 | #undef __SUF__ | ||
| 285 | #undef __QUOTEME | ||
| 286 | #undef __SUF_C__ | ||
| 287 |