GCC Code Coverage Report

Directory: ./
File: src/sys/classes/bv/impls/svec/svec.c
Date: 2025-11-19 04:19:03
Exec Total Coverage
Lines: 407 418 97.4%
Functions: 24 24 100.0%
Branches: 864 1442 59.9%
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 BV implemented as a single Vec
12 */
13
14 #include <slepc/private/bvimpl.h>
15 #include "svec.h"
16
17 998 static PetscErrorCode BVMult_Svec(BV Y,PetscScalar alpha,PetscScalar beta,BV X,Mat Q)
18 {
19 998 BV_SVEC *y = (BV_SVEC*)Y->data,*x = (BV_SVEC*)X->data;
20 998 const PetscScalar *px,*q;
21 998 PetscScalar *py;
22 998 PetscInt ldq;
23
24
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 998 PetscFunctionBegin;
25
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.
 998 PetscCall(VecGetArrayRead(x->v,&px));
26
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.
 998 PetscCall(VecGetArray(y->v,&py));
27
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 998 if (Q) {
28
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.
 798 PetscCall(MatDenseGetLDA(Q,&ldq));
29
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.
 798 PetscCall(MatDenseGetArrayRead(Q,&q));
30
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.
 798 PetscCall(BVMult_BLAS_Private(Y,Y->n,Y->k-Y->l,X->k-X->l,alpha,px+(X->nc+X->l)*X->ld,X->ld,q+Y->l*ldq+X->l,ldq,beta,py+(Y->nc+Y->l)*Y->ld,Y->ld));
31
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.
 798 PetscCall(MatDenseRestoreArrayRead(Q,&q));
32
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.
 200 } else PetscCall(BVAXPY_BLAS_Private(Y,Y->n,Y->k-Y->l,alpha,px+(X->nc+X->l)*X->ld,X->ld,beta,py+(Y->nc+Y->l)*Y->ld,Y->ld));
33
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.
 998 PetscCall(VecRestoreArrayRead(x->v,&px));
34
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.
 998 PetscCall(VecRestoreArray(y->v,&py));
35
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.
 196 PetscFunctionReturn(PETSC_SUCCESS);
36 }
37
38 3674 static PetscErrorCode BVMultVec_Svec(BV X,PetscScalar alpha,PetscScalar beta,Vec y,PetscScalar *q)
39 {
40 3674 BV_SVEC *x = (BV_SVEC*)X->data;
41 3674 PetscScalar *px,*py,*qq=q;
42
43
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 3674 PetscFunctionBegin;
44
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.
 3674 PetscCall(VecGetArray(x->v,&px));
45
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.
 3674 PetscCall(VecGetArray(y,&py));
46
6/8
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 3674 if (!q) PetscCall(VecGetArray(X->buffer,&qq));
47
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.
 3674 PetscCall(BVMultVec_BLAS_Private(X,X->n,X->k-X->l,alpha,px+(X->nc+X->l)*X->ld,X->ld,qq,beta,py));
48
6/8
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 3674 if (!q) PetscCall(VecRestoreArray(X->buffer,&qq));
49
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.
 3674 PetscCall(VecRestoreArray(x->v,&px));
50
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.
 3674 PetscCall(VecRestoreArray(y,&py));
51
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.
 713 PetscFunctionReturn(PETSC_SUCCESS);
52 }
53
54 478 static PetscErrorCode BVMultInPlace_Svec(BV V,Mat Q,PetscInt s,PetscInt e)
55 {
56 478 BV_SVEC *ctx = (BV_SVEC*)V->data;
57 478 PetscScalar *pv;
58 478 const PetscScalar *q;
59 478 PetscInt ldq;
60
61
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 478 PetscFunctionBegin;
62
3/16
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 2 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
 478 if (s>=e || !V->n) PetscFunctionReturn(PETSC_SUCCESS);
63
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.
 478 PetscCall(MatDenseGetLDA(Q,&ldq));
64
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.
 478 PetscCall(VecGetArray(ctx->v,&pv));
65
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.
 478 PetscCall(MatDenseGetArrayRead(Q,&q));
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.
 478 PetscCall(BVMultInPlace_BLAS_Private(V,V->n,V->k-V->l,s-V->l,e-V->l,pv+(V->nc+V->l)*V->ld,V->ld,q+V->l*ldq+V->l,ldq,PETSC_FALSE));
67
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.
 478 PetscCall(MatDenseRestoreArrayRead(Q,&q));
68
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.
 478 PetscCall(VecRestoreArray(ctx->v,&pv));
69
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.
 92 PetscFunctionReturn(PETSC_SUCCESS);
70 }
71
72 10 static PetscErrorCode BVMultInPlaceHermitianTranspose_Svec(BV V,Mat Q,PetscInt s,PetscInt e)
73 {
74 10 BV_SVEC *ctx = (BV_SVEC*)V->data;
75 10 PetscScalar *pv;
76 10 const PetscScalar *q;
77 10 PetscInt ldq;
78
79
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 10 PetscFunctionBegin;
80
3/16
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 2 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
 10 if (s>=e || !V->n) PetscFunctionReturn(PETSC_SUCCESS);
81
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.
 10 PetscCall(MatDenseGetLDA(Q,&ldq));
82
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.
 10 PetscCall(VecGetArray(ctx->v,&pv));
83
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.
 10 PetscCall(MatDenseGetArrayRead(Q,&q));
84
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.
 10 PetscCall(BVMultInPlace_BLAS_Private(V,V->n,V->k-V->l,s-V->l,e-V->l,pv+(V->nc+V->l)*V->ld,V->ld,q+V->l*ldq+V->l,ldq,PETSC_TRUE));
85
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.
 10 PetscCall(MatDenseRestoreArrayRead(Q,&q));
86
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.
 10 PetscCall(VecRestoreArray(ctx->v,&pv));
87
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.
 2 PetscFunctionReturn(PETSC_SUCCESS);
88 }
89
90 2436 static PetscErrorCode BVDot_Svec(BV X,BV Y,Mat M)
91 {
92 2436 BV_SVEC *x = (BV_SVEC*)X->data,*y = (BV_SVEC*)Y->data;
93 2436 const PetscScalar *px,*py;
94 2436 PetscScalar *m;
95 2436 PetscInt ldm;
96
97
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2436 PetscFunctionBegin;
98
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.
 2436 PetscCall(MatDenseGetLDA(M,&ldm));
99
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.
 2436 PetscCall(VecGetArrayRead(x->v,&px));
100
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.
 2436 PetscCall(VecGetArrayRead(y->v,&py));
101
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.
 2436 PetscCall(MatDenseGetArray(M,&m));
102
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.
 2436 PetscCall(BVDot_BLAS_Private(X,Y->k-Y->l,X->k-X->l,X->n,py+(Y->nc+Y->l)*Y->ld,Y->ld,px+(X->nc+X->l)*X->ld,X->ld,m+X->l*ldm+Y->l,ldm,x->mpi));
103
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.
 2436 PetscCall(MatDenseRestoreArray(M,&m));
104
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.
 2436 PetscCall(VecRestoreArrayRead(x->v,&px));
105
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.
 2436 PetscCall(VecRestoreArrayRead(y->v,&py));
106
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.
 482 PetscFunctionReturn(PETSC_SUCCESS);
107 }
108
109 4244 static PetscErrorCode BVDotVec_Svec(BV X,Vec y,PetscScalar *q)
110 {
111 4244 BV_SVEC *x = (BV_SVEC*)X->data;
112 4244 const PetscScalar *px,*py;
113 4244 PetscScalar *qq=q;
114 4244 Vec z = y;
115
116
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 4244 PetscFunctionBegin;
117
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 4244 if (PetscUnlikely(X->matrix)) {
118
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.
 620 PetscCall(BV_IPMatMult(X,y));
119 620 z = X->Bx;
120 }
121
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.
 4244 PetscCall(VecGetArrayRead(x->v,&px));
122
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.
 4244 PetscCall(VecGetArrayRead(z,&py));
123
6/8
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 4244 if (!q) PetscCall(VecGetArray(X->buffer,&qq));
124
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.
 4244 PetscCall(BVDotVec_BLAS_Private(X,X->n,X->k-X->l,px+(X->nc+X->l)*X->ld,X->ld,py,qq,x->mpi));
125
6/8
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 4244 if (!q) PetscCall(VecRestoreArray(X->buffer,&qq));
126
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.
 4244 PetscCall(VecRestoreArrayRead(z,&py));
127
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.
 4244 PetscCall(VecRestoreArrayRead(x->v,&px));
128
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.
 827 PetscFunctionReturn(PETSC_SUCCESS);
129 }
130
131 40 static PetscErrorCode BVDotVec_Local_Svec(BV X,Vec y,PetscScalar *m)
132 {
133 40 BV_SVEC *x = (BV_SVEC*)X->data;
134 40 PetscScalar *px,*py;
135 40 Vec z = y;
136
137
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 40 PetscFunctionBegin;
138
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 40 if (PetscUnlikely(X->matrix)) {
139 ✗ PetscCall(BV_IPMatMult(X,y));
140 ✗ z = X->Bx;
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.
 40 PetscCall(VecGetArray(x->v,&px));
143
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.
 40 PetscCall(VecGetArray(z,&py));
144
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.
 40 PetscCall(BVDotVec_BLAS_Private(X,X->n,X->k-X->l,px+(X->nc+X->l)*X->ld,X->ld,py,m,PETSC_FALSE));
145
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.
 40 PetscCall(VecRestoreArray(z,&py));
146
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.
 40 PetscCall(VecRestoreArray(x->v,&px));
147
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.
 8 PetscFunctionReturn(PETSC_SUCCESS);
148 }
149
150 3338 static PetscErrorCode BVScale_Svec(BV bv,PetscInt j,PetscScalar alpha)
151 {
152 3338 BV_SVEC *ctx = (BV_SVEC*)bv->data;
153 3338 PetscScalar *array;
154
155
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 3338 PetscFunctionBegin;
156
2/14
✓ Branch 0 taken 8 times.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
 3338 if (!bv->n) PetscFunctionReturn(PETSC_SUCCESS);
157
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.
 3338 PetscCall(VecGetArray(ctx->v,&array));
158
6/8
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 3338 if (PetscUnlikely(j<0)) PetscCall(BVScale_BLAS_Private(bv,(bv->k-bv->l)*bv->ld,array+(bv->nc+bv->l)*bv->ld,alpha));
159
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.
 3280 else PetscCall(BVScale_BLAS_Private(bv,bv->n,array+(bv->nc+j)*bv->ld,alpha));
160
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.
 3338 PetscCall(VecRestoreArray(ctx->v,&array));
161
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.
 648 PetscFunctionReturn(PETSC_SUCCESS);
162 }
163
164 1123 static PetscErrorCode BVNorm_Svec(BV bv,PetscInt j,NormType type,PetscReal *val)
165 {
166 1123 BV_SVEC *ctx = (BV_SVEC*)bv->data;
167 1123 const PetscScalar *array;
168
169
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 1123 PetscFunctionBegin;
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.
 1123 PetscCall(VecGetArrayRead(ctx->v,&array));
171
6/8
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 1123 if (PetscUnlikely(j<0)) PetscCall(BVNorm_LAPACK_Private(bv,bv->n,bv->k-bv->l,array+(bv->nc+bv->l)*bv->ld,bv->ld,type,val,ctx->mpi));
172
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.
 425 else PetscCall(BVNorm_LAPACK_Private(bv,bv->n,1,array+(bv->nc+j)*bv->ld,bv->ld,type,val,ctx->mpi));
173
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.
 1123 PetscCall(VecRestoreArrayRead(ctx->v,&array));
174
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.
 217 PetscFunctionReturn(PETSC_SUCCESS);
175 }
176
177 20 static PetscErrorCode BVNorm_Local_Svec(BV bv,PetscInt j,NormType type,PetscReal *val)
178 {
179 20 BV_SVEC *ctx = (BV_SVEC*)bv->data;
180 20 const PetscScalar *array;
181
182
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 20 PetscFunctionBegin;
183
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.
 20 PetscCall(VecGetArrayRead(ctx->v,&array));
184
1/8
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
 20 if (PetscUnlikely(j<0)) PetscCall(BVNorm_LAPACK_Private(bv,bv->n,bv->k-bv->l,array+(bv->nc+bv->l)*bv->ld,bv->ld,type,val,PETSC_FALSE));
185
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.
 20 else PetscCall(BVNorm_LAPACK_Private(bv,bv->n,1,array+(bv->nc+j)*bv->ld,bv->ld,type,val,PETSC_FALSE));
186
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.
 20 PetscCall(VecRestoreArrayRead(ctx->v,&array));
187
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.
 4 PetscFunctionReturn(PETSC_SUCCESS);
188 }
189
190 45 static PetscErrorCode BVNormalize_Svec(BV bv,PetscScalar *eigi)
191 {
192 45 BV_SVEC *ctx = (BV_SVEC*)bv->data;
193 45 PetscScalar *array,*wi=NULL;
194
195
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 45 PetscFunctionBegin;
196
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.
 45 PetscCall(VecGetArray(ctx->v,&array));
197
2/2
✓ Branch 0 taken 5 times.
✓ Branch 1 taken 10 times.
 45 if (eigi) wi = eigi+bv->l;
198
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.
 45 PetscCall(BVNormalize_LAPACK_Private(bv,bv->n,bv->k-bv->l,array+(bv->nc+bv->l)*bv->ld,bv->ld,wi,ctx->mpi));
199
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.
 45 PetscCall(VecRestoreArray(ctx->v,&array));
200
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.
 9 PetscFunctionReturn(PETSC_SUCCESS);
201 }
202
203 1070 static PetscErrorCode BVMatMult_Svec(BV V,Mat A,BV W)
204 {
205 1070 PetscInt j;
206 1070 Mat Vmat,Wmat;
207 1070 Vec vv,ww;
208
209
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 1070 PetscFunctionBegin;
210
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 1070 if (V->vmm) {
211
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.
 980 PetscCall(BVGetMat(V,&Vmat));
212
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.
 980 PetscCall(BVGetMat(W,&Wmat));
213
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.
 980 PetscCall(MatProductCreateWithMat(A,Vmat,NULL,Wmat));
214
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.
 980 PetscCall(MatProductSetType(Wmat,MATPRODUCT_AB));
215
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.
 980 PetscCall(MatProductSetFromOptions(Wmat));
216
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.
 980 PetscCall(MatProductSymbolic(Wmat));
217
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.
 980 PetscCall(MatProductNumeric(Wmat));
218
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.
 980 PetscCall(MatProductClear(Wmat));
219
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.
 980 PetscCall(BVRestoreMat(V,&Vmat));
220
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.
 980 PetscCall(BVRestoreMat(W,&Wmat));
221 } else {
222
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 650 for (j=0;j<V->k-V->l;j++) {
223
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.
 560 PetscCall(BVGetColumn(V,V->l+j,&vv));
224
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.
 560 PetscCall(BVGetColumn(W,W->l+j,&ww));
225
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.
 560 PetscCall(MatMult(A,vv,ww));
226
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.
 560 PetscCall(BVRestoreColumn(V,V->l+j,&vv));
227
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.
 560 PetscCall(BVRestoreColumn(W,W->l+j,&ww));
228 }
229 }
230
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.
 214 PetscFunctionReturn(PETSC_SUCCESS);
231 }
232
233 970 static PetscErrorCode BVCopy_Svec(BV V,BV W)
234 {
235 970 BV_SVEC *v = (BV_SVEC*)V->data,*w = (BV_SVEC*)W->data;
236 970 const PetscScalar *pv;
237 970 PetscScalar *pw;
238 970 PetscInt j;
239
240
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 970 PetscFunctionBegin;
241
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.
 970 PetscCall(VecGetArrayRead(v->v,&pv));
242
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.
 970 PetscCall(VecGetArray(w->v,&pw));
243
7/8
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
 7260 for (j=0;j<V->k-V->l;j++) PetscCall(PetscArraycpy(pw+(W->nc+W->l+j)*W->ld,pv+(V->nc+V->l+j)*V->ld,V->n));
244
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.
 970 PetscCall(VecRestoreArrayRead(v->v,&pv));
245
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.
 970 PetscCall(VecRestoreArray(w->v,&pw));
246
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.
 194 PetscFunctionReturn(PETSC_SUCCESS);
247 }
248
249 20 static PetscErrorCode BVCopyColumn_Svec(BV V,PetscInt j,PetscInt i)
250 {
251 20 BV_SVEC *v = (BV_SVEC*)V->data;
252 20 PetscScalar *pv;
253
254
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 20 PetscFunctionBegin;
255
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.
 20 PetscCall(VecGetArray(v->v,&pv));
256
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.
 20 PetscCall(PetscArraycpy(pv+(V->nc+i)*V->ld,pv+(V->nc+j)*V->ld,V->n));
257
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.
 20 PetscCall(VecRestoreArray(v->v,&pv));
258
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.
 4 PetscFunctionReturn(PETSC_SUCCESS);
259 }
260
261 130 static PetscErrorCode BVResize_Svec(BV bv,PetscInt m,PetscBool copy)
262 {
263 130 BV_SVEC *ctx = (BV_SVEC*)bv->data;
264 130 PetscScalar *pnew;
265 130 const PetscScalar *pv;
266 130 PetscInt bs;
267 130 Vec vnew;
268 130 char str[50];
269
270
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 130 PetscFunctionBegin;
271
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.
 130 PetscCall(PetscLayoutGetBlockSize(bv->map,&bs));
272
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.
 130 PetscCall(VecCreate(PetscObjectComm((PetscObject)bv),&vnew));
273
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.
 130 PetscCall(VecSetType(vnew,bv->vtype));
274
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.
 130 PetscCall(VecSetSizes(vnew,m*bv->ld,PETSC_DECIDE));
275
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.
 130 PetscCall(VecSetBlockSize(vnew,bs));
276
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 130 if (((PetscObject)bv)->name) {
277
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.
 130 PetscCall(PetscSNPrintf(str,sizeof(str),"%s_0",((PetscObject)bv)->name));
278
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.
 130 PetscCall(PetscObjectSetName((PetscObject)vnew,str));
279 }
280
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 130 if (copy) {
281
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.
 20 PetscCall(VecGetArrayRead(ctx->v,&pv));
282
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.
 20 PetscCall(VecGetArray(vnew,&pnew));
283
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.
 20 PetscCall(PetscArraycpy(pnew,pv,PetscMin(m,bv->m)*bv->ld));
284
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.
 20 PetscCall(VecRestoreArrayRead(ctx->v,&pv));
285
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.
 20 PetscCall(VecRestoreArray(vnew,&pnew));
286 }
287
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.
 130 PetscCall(VecDestroy(&ctx->v));
288 130 ctx->v = vnew;
289
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.
 130 PetscFunctionReturn(PETSC_SUCCESS);
290 }
291
292 27979 static PetscErrorCode BVGetColumn_Svec(BV bv,PetscInt j,Vec *v)
293 {
294 27979 BV_SVEC *ctx = (BV_SVEC*)bv->data;
295 27979 PetscScalar *pv;
296 27979 PetscInt l;
297
298
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 27979 PetscFunctionBegin;
299
3/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
 27979 l = BVAvailableVec;
300
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.
 27979 PetscCall(VecGetArray(ctx->v,&pv));
301
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.
 27979 PetscCall(VecPlaceArray(bv->cv[l],pv+(bv->nc+j)*bv->ld));
302
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.
 5381 PetscFunctionReturn(PETSC_SUCCESS);
303 }
304
305 27979 static PetscErrorCode BVRestoreColumn_Svec(BV bv,PetscInt j,Vec *v)
306 {
307 27979 BV_SVEC *ctx = (BV_SVEC*)bv->data;
308 27979 PetscInt l;
309
310
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 27979 PetscFunctionBegin;
311 27979 l = (j==bv->ci[0])? 0: 1;
312
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.
 27979 PetscCall(VecResetArray(bv->cv[l]));
313
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.
 27979 PetscCall(VecRestoreArray(ctx->v,NULL));
314
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.
 5381 PetscFunctionReturn(PETSC_SUCCESS);
315 }
316
317 2518 static PetscErrorCode BVGetArray_Svec(BV bv,PetscScalar **a)
318 {
319 2518 BV_SVEC *ctx = (BV_SVEC*)bv->data;
320
321
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2518 PetscFunctionBegin;
322
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.
 2518 PetscCall(VecGetArray(ctx->v,a));
323
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.
 490 PetscFunctionReturn(PETSC_SUCCESS);
324 }
325
326 2518 static PetscErrorCode BVRestoreArray_Svec(BV bv,PetscScalar **a)
327 {
328 2518 BV_SVEC *ctx = (BV_SVEC*)bv->data;
329
330
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 2518 PetscFunctionBegin;
331
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.
 2518 PetscCall(VecRestoreArray(ctx->v,a));
332
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.
 490 PetscFunctionReturn(PETSC_SUCCESS);
333 }
334
335 120 static PetscErrorCode BVGetArrayRead_Svec(BV bv,const PetscScalar **a)
336 {
337 120 BV_SVEC *ctx = (BV_SVEC*)bv->data;
338
339
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 120 PetscFunctionBegin;
340
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.
 120 PetscCall(VecGetArrayRead(ctx->v,a));
341
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.
 16 PetscFunctionReturn(PETSC_SUCCESS);
342 }
343
344 120 static PetscErrorCode BVRestoreArrayRead_Svec(BV bv,const PetscScalar **a)
345 {
346 120 BV_SVEC *ctx = (BV_SVEC*)bv->data;
347
348
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 120 PetscFunctionBegin;
349
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.
 120 PetscCall(VecRestoreArrayRead(ctx->v,a));
350
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.
 16 PetscFunctionReturn(PETSC_SUCCESS);
351 }
352
353 216 static PetscErrorCode BVView_Svec(BV bv,PetscViewer viewer)
354 {
355 216 PetscInt j;
356 216 Vec v;
357 216 PetscViewerFormat format;
358 216 PetscBool isascii,ismatlab=PETSC_FALSE;
359 216 const char *bvname,*name;
360
361
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 216 PetscFunctionBegin;
362
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.
 216 PetscCall(PetscObjectTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii));
363
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 216 if (isascii) {
364
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.
 216 PetscCall(PetscViewerGetFormat(viewer,&format));
365
8/14
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ 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.
 216 if (format == PETSC_VIEWER_ASCII_INFO || format == PETSC_VIEWER_ASCII_INFO_DETAIL) PetscFunctionReturn(PETSC_SUCCESS);
366
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 180 if (format == PETSC_VIEWER_ASCII_MATLAB) ismatlab = PETSC_TRUE;
367 }
368 ✗ if (ismatlab) {
369 ✗ PetscCall(PetscObjectGetName((PetscObject)bv,&bvname));
370 ✗ PetscCall(PetscViewerASCIIPrintf(viewer,"%s=[];\n",bvname));
371 }
372
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 864 for (j=0;j<bv->m;j++) {
373
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.
 684 PetscCall(BVGetColumn(bv,j,&v));
374
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.
 684 PetscCall(VecView(v,viewer));
375
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
 684 if (ismatlab) {
376 ✗ PetscCall(PetscObjectGetName((PetscObject)v,&name));
377 ✗ PetscCall(PetscViewerASCIIPrintf(viewer,"%s=[%s,%s];clear %s\n",bvname,bvname,name,name));
378 }
379
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.
 684 PetscCall(BVRestoreColumn(bv,j,&v));
380 }
381
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.
 20 PetscFunctionReturn(PETSC_SUCCESS);
382 }
383
384 4606 static PetscErrorCode BVDestroy_Svec(BV bv)
385 {
386 4606 BV_SVEC *ctx = (BV_SVEC*)bv->data;
387
388
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 4606 PetscFunctionBegin;
389
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.
 4606 PetscCall(VecDestroy(&ctx->v));
390
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.
 4606 PetscCall(VecDestroy(&bv->cv[0]));
391
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.
 4606 PetscCall(VecDestroy(&bv->cv[1]));
392
5/8
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 4606 PetscCall(PetscFree(bv->data));
393 4606 bv->cuda = PETSC_FALSE;
394
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.
 4606 PetscFunctionReturn(PETSC_SUCCESS);
395 }
396
397 /*MC
398 BVSVEC - BVSVEC = "svec" - A basis vectors object represented by a single
399 `Vec` object. The local array of the `Vec` stores the local part of the
400 basis vectors column-wise.
401
402 Level: beginner
403
404 .seealso: [](sec:bv), `BV`, `BVType`, `BVSetType()`
405 M*/
406
407 4606 SLEPC_EXTERN PetscErrorCode BVCreate_Svec(BV bv)
408 {
409 4606 BV_SVEC *ctx;
410 4606 PetscInt nloc,N,bs,tglobal=0,tlocal,lsplit,j,lda;
411 4606 PetscBool seq,isdense;
412 4606 PetscScalar *vv;
413 4606 const PetscScalar *aa,*array,*ptr;
414 4606 char str[50];
415 4606 BV parent;
416 4606 Vec vpar;
417 4606 MatType mtype;
418
419
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 4606 PetscFunctionBegin;
420
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.
 4606 PetscCall(PetscNew(&ctx));
421 4606 bv->data = (void*)ctx;
422
423
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.
 4606 PetscCall(PetscStrcmpAny(bv->vtype,&bv->cuda,VECSEQCUDA,VECMPICUDA,""));
424
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.
 4606 PetscCall(PetscStrcmpAny(bv->vtype,&bv->hip,VECSEQHIP,VECMPIHIP,""));
425
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.
 4606 PetscCall(PetscStrcmpAny(bv->vtype,&ctx->mpi,VECMPI,VECMPICUDA,VECMPIHIP,""));
426
427
3/6
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
 4606 PetscCall(PetscStrcmp(bv->vtype,VECSEQ,&seq));
428
7/10
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 4 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 4606 PetscCheck(seq || ctx->mpi || bv->cuda || bv->hip,PetscObjectComm((PetscObject)bv),PETSC_ERR_SUP,"BVSVEC does not support the requested vector type: %s",bv->vtype);
429
430
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.
 4606 PetscCall(PetscLayoutGetLocalSize(bv->map,&nloc));
431
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.
 4606 PetscCall(PetscLayoutGetSize(bv->map,&N));
432
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.
 4606 PetscCall(PetscLayoutGetBlockSize(bv->map,&bs));
433
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.
 4606 PetscCall(BV_SetDefaultLD(bv,nloc));
434 4606 tlocal = bv->m*bv->ld;
435
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.
 4606 PetscCall(PetscIntMultError(bv->m,N,&tglobal)); /* just to check integer overflow */
436
437
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 4606 if (PetscUnlikely(bv->issplit)) {
438 /* split BV: create Vec sharing the memory of the parent BV */
439 1232 parent = bv->splitparent;
440 1232 lsplit = parent->lsplit;
441 1232 vpar = ((BV_SVEC*)parent->data)->v;
442
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 2 times.
 1232 if (bv->cuda) {
443 #if defined(PETSC_HAVE_CUDA)
444
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 132 PetscCall(VecCUDAGetArrayRead(vpar,&array));
445
3/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 2 times.
 132 if (bv->issplit>0) ptr = (bv->issplit==1)? array: array+lsplit*bv->ld;
446 ✗ else ptr = (bv->issplit==1)? array: array-lsplit;
447
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 132 PetscCall(VecCUDARestoreArrayRead(vpar,&array));
448
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
 132 if (ctx->mpi) PetscCall(VecCreateMPICUDAWithArray(PetscObjectComm((PetscObject)bv),bs,tlocal,PETSC_DECIDE,NULL,&ctx->v));
449 ✗ else PetscCall(VecCreateSeqCUDAWithArray(PetscObjectComm((PetscObject)bv),bs,tlocal,NULL,&ctx->v));
450
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 132 PetscCall(VecCUDAPlaceArray(ctx->v,ptr));
451 #endif
452
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 2 times.
 1100 } else if (bv->hip) {
453 #if defined(PETSC_HAVE_HIP)
454
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 132 PetscCall(VecHIPGetArrayRead(vpar,&array));
455
3/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 2 times.
 132 if (bv->issplit>0) ptr = (bv->issplit==1)? array: array+lsplit*bv->ld;
456 ✗ else ptr = (bv->issplit==1)? array: array-lsplit;
457
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 132 PetscCall(VecHIPRestoreArrayRead(vpar,&array));
458
2/4
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
 132 if (ctx->mpi) PetscCall(VecCreateMPIHIPWithArray(PetscObjectComm((PetscObject)bv),bs,tlocal,PETSC_DECIDE,NULL,&ctx->v));
459 ✗ else PetscCall(VecCreateSeqHIPWithArray(PetscObjectComm((PetscObject)bv),bs,tlocal,NULL,&ctx->v));
460
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
 132 PetscCall(VecHIPPlaceArray(ctx->v,ptr));
461 #endif
462 } else {
463
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.
 968 PetscCall(VecGetArrayRead(vpar,&array));
464
4/4
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
 968 if (bv->issplit>0) ptr = (bv->issplit==1)? array: array+lsplit*bv->ld;
465
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 168 else ptr = (bv->issplit==1)? array: array-lsplit;
466
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.
 968 PetscCall(VecRestoreArrayRead(vpar,&array));
467
6/8
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
 968 if (ctx->mpi) PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)bv),bs,tlocal,PETSC_DECIDE,NULL,&ctx->v));
468
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.
 56 else PetscCall(VecCreateSeqWithArray(PetscObjectComm((PetscObject)bv),bs,tlocal,NULL,&ctx->v));
469
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.
 968 PetscCall(VecPlaceArray(ctx->v,ptr));
470 }
471 } else {
472 /* regular BV: create Vec to store the BV entries */
473
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.
 3374 PetscCall(VecCreate(PetscObjectComm((PetscObject)bv),&ctx->v));
474
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.
 3374 PetscCall(VecSetType(ctx->v,bv->vtype));
475
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.
 3374 PetscCall(VecSetSizes(ctx->v,tlocal,PETSC_DECIDE));
476
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.
 3374 PetscCall(VecSetBlockSize(ctx->v,bs));
477 }
478
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 4606 if (((PetscObject)bv)->name) {
479
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.
 1522 PetscCall(PetscSNPrintf(str,sizeof(str),"%s_0",((PetscObject)bv)->name));
480
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.
 1522 PetscCall(PetscObjectSetName((PetscObject)ctx->v,str));
481 }
482
483
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 4606 if (PetscUnlikely(bv->Acreate)) {
484
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.
 20 PetscCall(MatGetType(bv->Acreate,&mtype));
485
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.
 20 PetscCall(PetscStrcmpAny(mtype,&isdense,MATSEQDENSE,MATMPIDENSE,""));
486
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 20 PetscCheck(isdense,PetscObjectComm((PetscObject)bv->Acreate),PETSC_ERR_SUP,"BVSVEC requires a dense matrix in BVCreateFromMat()");
487
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.
 20 PetscCall(MatDenseGetArrayRead(bv->Acreate,&aa));
488
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.
 20 PetscCall(MatDenseGetLDA(bv->Acreate,&lda));
489
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.
 20 PetscCall(VecGetArray(ctx->v,&vv));
490
7/8
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
 180 for (j=0;j<bv->m;j++) PetscCall(PetscArraycpy(vv+j*bv->ld,aa+j*lda,bv->n));
491
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.
 20 PetscCall(VecRestoreArray(ctx->v,&vv));
492
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.
 20 PetscCall(MatDenseRestoreArrayRead(bv->Acreate,&aa));
493
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.
 20 PetscCall(MatDestroy(&bv->Acreate));
494 }
495
496
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.
 4606 PetscCall(BVCreateVecEmpty(bv,&bv->cv[0]));
497
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.
 4606 PetscCall(BVCreateVecEmpty(bv,&bv->cv[1]));
498
499
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 2 times.
 4606 if (bv->cuda) {
500 #if defined(PETSC_HAVE_CUDA)
501 470 bv->ops->mult = BVMult_Svec_CUDA;
502 470 bv->ops->multvec = BVMultVec_Svec_CUDA;
503 470 bv->ops->multinplace = BVMultInPlace_Svec_CUDA;
504 470 bv->ops->multinplacetrans = BVMultInPlaceHermitianTranspose_Svec_CUDA;
505 470 bv->ops->dot = BVDot_Svec_CUDA;
506 470 bv->ops->dotvec = BVDotVec_Svec_CUDA;
507 470 bv->ops->dotvec_local = BVDotVec_Local_Svec_CUDA;
508 470 bv->ops->scale = BVScale_Svec_CUDA;
509 470 bv->ops->norm = BVNorm_Svec_CUDA;
510 470 bv->ops->norm_local = BVNorm_Local_Svec_CUDA;
511 470 bv->ops->normalize = BVNormalize_Svec_CUDA;
512 470 bv->ops->matmult = BVMatMult_Svec_CUDA;
513 470 bv->ops->copy = BVCopy_Svec_CUDA;
514 470 bv->ops->copycolumn = BVCopyColumn_Svec_CUDA;
515 470 bv->ops->resize = BVResize_Svec_CUDA;
516 470 bv->ops->getcolumn = BVGetColumn_Svec_CUDA;
517 470 bv->ops->restorecolumn = BVRestoreColumn_Svec_CUDA;
518 470 bv->ops->restoresplit = BVRestoreSplit_Svec_CUDA;
519 470 bv->ops->restoresplitrows = BVRestoreSplitRows_Svec_CUDA;
520 470 bv->ops->getmat = BVGetMat_Svec_CUDA;
521 470 bv->ops->restoremat = BVRestoreMat_Svec_CUDA;
522 #endif
523
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 2 times.
 4136 } else if (bv->hip) {
524 #if defined(PETSC_HAVE_HIP)
525 470 bv->ops->mult = BVMult_Svec_HIP;
526 470 bv->ops->multvec = BVMultVec_Svec_HIP;
527 470 bv->ops->multinplace = BVMultInPlace_Svec_HIP;
528 470 bv->ops->multinplacetrans = BVMultInPlaceHermitianTranspose_Svec_HIP;
529 470 bv->ops->dot = BVDot_Svec_HIP;
530 470 bv->ops->dotvec = BVDotVec_Svec_HIP;
531 470 bv->ops->dotvec_local = BVDotVec_Local_Svec_HIP;
532 470 bv->ops->scale = BVScale_Svec_HIP;
533 470 bv->ops->norm = BVNorm_Svec_HIP;
534 470 bv->ops->norm_local = BVNorm_Local_Svec_HIP;
535 470 bv->ops->normalize = BVNormalize_Svec_HIP;
536 470 bv->ops->matmult = BVMatMult_Svec_HIP;
537 470 bv->ops->copy = BVCopy_Svec_HIP;
538 470 bv->ops->copycolumn = BVCopyColumn_Svec_HIP;
539 470 bv->ops->resize = BVResize_Svec_HIP;
540 470 bv->ops->getcolumn = BVGetColumn_Svec_HIP;
541 470 bv->ops->restorecolumn = BVRestoreColumn_Svec_HIP;
542 470 bv->ops->restoresplit = BVRestoreSplit_Svec_HIP;
543 470 bv->ops->restoresplitrows = BVRestoreSplitRows_Svec_HIP;
544 470 bv->ops->getmat = BVGetMat_Svec_HIP;
545 470 bv->ops->restoremat = BVRestoreMat_Svec_HIP;
546 #endif
547 } else {
548 3666 bv->ops->mult = BVMult_Svec;
549 3666 bv->ops->multvec = BVMultVec_Svec;
550 3666 bv->ops->multinplace = BVMultInPlace_Svec;
551 3666 bv->ops->multinplacetrans = BVMultInPlaceHermitianTranspose_Svec;
552 3666 bv->ops->dot = BVDot_Svec;
553 3666 bv->ops->dotvec = BVDotVec_Svec;
554 3666 bv->ops->dotvec_local = BVDotVec_Local_Svec;
555 3666 bv->ops->scale = BVScale_Svec;
556 3666 bv->ops->norm = BVNorm_Svec;
557 3666 bv->ops->norm_local = BVNorm_Local_Svec;
558 3666 bv->ops->normalize = BVNormalize_Svec;
559 3666 bv->ops->matmult = BVMatMult_Svec;
560 3666 bv->ops->copy = BVCopy_Svec;
561 3666 bv->ops->copycolumn = BVCopyColumn_Svec;
562 3666 bv->ops->resize = BVResize_Svec;
563 3666 bv->ops->getcolumn = BVGetColumn_Svec;
564 3666 bv->ops->restorecolumn = BVRestoreColumn_Svec;
565 3666 bv->ops->getmat = BVGetMat_Default;
566 3666 bv->ops->restoremat = BVRestoreMat_Default;
567 }
568 4606 bv->ops->getarray = BVGetArray_Svec;
569 4606 bv->ops->restorearray = BVRestoreArray_Svec;
570 4606 bv->ops->getarrayread = BVGetArrayRead_Svec;
571 4606 bv->ops->restorearrayread = BVRestoreArrayRead_Svec;
572 4606 bv->ops->destroy = BVDestroy_Svec;
573
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 4606 if (!ctx->mpi) bv->ops->view = BVView_Svec;
574
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.
 710 PetscFunctionReturn(PETSC_SUCCESS);
575 }
576