| 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 an array of Vecs sharing a contiguous array for elements | ||
| 12 | */ | ||
| 13 | |||
| 14 | #include <slepc/private/bvimpl.h> | ||
| 15 | |||
| 16 | typedef struct { | ||
| 17 | Vec *V; | ||
| 18 | PetscScalar *array; | ||
| 19 | PetscBool mpi; | ||
| 20 | } BV_CONTIGUOUS; | ||
| 21 | |||
| 22 | 799 | static PetscErrorCode BVMult_Contiguous(BV Y,PetscScalar alpha,PetscScalar beta,BV X,Mat Q) | |
| 23 | { | ||
| 24 | 799 | BV_CONTIGUOUS *y = (BV_CONTIGUOUS*)Y->data,*x = (BV_CONTIGUOUS*)X->data; | |
| 25 | 799 | const PetscScalar *q; | |
| 26 | 799 | PetscInt ldq; | |
| 27 | |||
| 28 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
799 | PetscFunctionBegin; |
| 29 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
799 | if (Q) { |
| 30 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
639 | PetscCall(MatDenseGetLDA(Q,&ldq)); |
| 31 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
639 | PetscCall(MatDenseGetArrayRead(Q,&q)); |
| 32 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
639 | PetscCall(BVMult_BLAS_Private(Y,Y->n,Y->k-Y->l,X->k-X->l,alpha,x->array+(X->nc+X->l)*X->ld,X->ld,q+Y->l*ldq+X->l,ldq,beta,y->array+(Y->nc+Y->l)*Y->ld,Y->ld)); |
| 33 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
639 | PetscCall(MatDenseRestoreArrayRead(Q,&q)); |
| 34 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
160 | } else PetscCall(BVAXPY_BLAS_Private(Y,Y->n,Y->k-Y->l,alpha,x->array+(X->nc+X->l)*X->ld,X->ld,beta,y->array+(Y->nc+Y->l)*Y->ld,Y->ld)); |
| 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 | 2784 | static PetscErrorCode BVMultVec_Contiguous(BV X,PetscScalar alpha,PetscScalar beta,Vec y,PetscScalar *q) | |
| 39 | { | ||
| 40 | 2784 | BV_CONTIGUOUS *x = (BV_CONTIGUOUS*)X->data; | |
| 41 | 2784 | PetscScalar *py,*qq=q; | |
| 42 | |||
| 43 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
2784 | PetscFunctionBegin; |
| 44 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2784 | PetscCall(VecGetArray(y,&py)); |
| 45 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
2784 | if (!q) PetscCall(VecGetArray(X->buffer,&qq)); |
| 46 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2784 | PetscCall(BVMultVec_BLAS_Private(X,X->n,X->k-X->l,alpha,x->array+(X->nc+X->l)*X->ld,X->ld,qq,beta,py)); |
| 47 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
2784 | if (!q) PetscCall(VecRestoreArray(X->buffer,&qq)); |
| 48 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2784 | PetscCall(VecRestoreArray(y,&py)); |
| 49 |
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.
|
693 | PetscFunctionReturn(PETSC_SUCCESS); |
| 50 | } | ||
| 51 | |||
| 52 | 383 | static PetscErrorCode BVMultInPlace_Contiguous(BV V,Mat Q,PetscInt s,PetscInt e) | |
| 53 | { | ||
| 54 | 383 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)V->data; | |
| 55 | 383 | const PetscScalar *q; | |
| 56 | 383 | PetscInt ldq; | |
| 57 | |||
| 58 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
383 | PetscFunctionBegin; |
| 59 |
3/16✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 6 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.
|
383 | if (s>=e || !V->n) PetscFunctionReturn(PETSC_SUCCESS); |
| 60 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
383 | PetscCall(MatDenseGetLDA(Q,&ldq)); |
| 61 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
383 | PetscCall(MatDenseGetArrayRead(Q,&q)); |
| 62 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
383 | PetscCall(BVMultInPlace_BLAS_Private(V,V->n,V->k-V->l,s-V->l,e-V->l,ctx->array+(V->nc+V->l)*V->ld,V->ld,q+V->l*ldq+V->l,ldq,PETSC_FALSE)); |
| 63 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
383 | PetscCall(MatDenseRestoreArrayRead(Q,&q)); |
| 64 |
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); |
| 65 | } | ||
| 66 | |||
| 67 | 8 | static PetscErrorCode BVMultInPlaceHermitianTranspose_Contiguous(BV V,Mat Q,PetscInt s,PetscInt e) | |
| 68 | { | ||
| 69 | 8 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)V->data; | |
| 70 | 8 | const PetscScalar *q; | |
| 71 | 8 | PetscInt ldq; | |
| 72 | |||
| 73 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
8 | PetscFunctionBegin; |
| 74 |
3/16✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 6 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.
|
8 | if (s>=e || !V->n) PetscFunctionReturn(PETSC_SUCCESS); |
| 75 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
8 | PetscCall(MatDenseGetLDA(Q,&ldq)); |
| 76 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
8 | PetscCall(MatDenseGetArrayRead(Q,&q)); |
| 77 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
8 | PetscCall(BVMultInPlace_BLAS_Private(V,V->n,V->k-V->l,s-V->l,e-V->l,ctx->array+(V->nc+V->l)*V->ld,V->ld,q+V->l*ldq+V->l,ldq,PETSC_TRUE)); |
| 78 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
8 | PetscCall(MatDenseRestoreArrayRead(Q,&q)); |
| 79 |
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); |
| 80 | } | ||
| 81 | |||
| 82 | 1895 | static PetscErrorCode BVDot_Contiguous(BV X,BV Y,Mat M) | |
| 83 | { | ||
| 84 | 1895 | BV_CONTIGUOUS *x = (BV_CONTIGUOUS*)X->data,*y = (BV_CONTIGUOUS*)Y->data; | |
| 85 | 1895 | PetscScalar *m; | |
| 86 | 1895 | PetscInt ldm; | |
| 87 | |||
| 88 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1895 | PetscFunctionBegin; |
| 89 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1895 | PetscCall(MatDenseGetLDA(M,&ldm)); |
| 90 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1895 | PetscCall(MatDenseGetArray(M,&m)); |
| 91 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1895 | PetscCall(BVDot_BLAS_Private(X,Y->k-Y->l,X->k-X->l,X->n,y->array+(Y->nc+Y->l)*Y->ld,Y->ld,x->array+(X->nc+X->l)*X->ld,X->ld,m+X->l*ldm+Y->l,ldm,x->mpi)); |
| 92 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1895 | PetscCall(MatDenseRestoreArray(M,&m)); |
| 93 |
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.
|
470 | PetscFunctionReturn(PETSC_SUCCESS); |
| 94 | } | ||
| 95 | |||
| 96 | 2840 | static PetscErrorCode BVDotVec_Contiguous(BV X,Vec y,PetscScalar *q) | |
| 97 | { | ||
| 98 | 2840 | BV_CONTIGUOUS *x = (BV_CONTIGUOUS*)X->data; | |
| 99 | 2840 | const PetscScalar *py; | |
| 100 | 2840 | PetscScalar *qq=q; | |
| 101 | 2840 | Vec z = y; | |
| 102 | |||
| 103 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
2840 | PetscFunctionBegin; |
| 104 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
2840 | if (PetscUnlikely(X->matrix)) { |
| 105 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
456 | PetscCall(BV_IPMatMult(X,y)); |
| 106 | 456 | z = X->Bx; | |
| 107 | } | ||
| 108 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2840 | PetscCall(VecGetArrayRead(z,&py)); |
| 109 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
2840 | if (!q) PetscCall(VecGetArray(X->buffer,&qq)); |
| 110 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2840 | PetscCall(BVDotVec_BLAS_Private(X,X->n,X->k-X->l,x->array+(X->nc+X->l)*X->ld,X->ld,py,qq,x->mpi)); |
| 111 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
2840 | if (!q) PetscCall(VecRestoreArray(X->buffer,&qq)); |
| 112 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2840 | PetscCall(VecRestoreArrayRead(z,&py)); |
| 113 |
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.
|
707 | PetscFunctionReturn(PETSC_SUCCESS); |
| 114 | } | ||
| 115 | |||
| 116 | 32 | static PetscErrorCode BVDotVec_Local_Contiguous(BV X,Vec y,PetscScalar *m) | |
| 117 | { | ||
| 118 | 32 | BV_CONTIGUOUS *x = (BV_CONTIGUOUS*)X->data; | |
| 119 | 32 | PetscScalar *py; | |
| 120 | 32 | Vec z = y; | |
| 121 | |||
| 122 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
32 | PetscFunctionBegin; |
| 123 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
|
32 | if (PetscUnlikely(X->matrix)) { |
| 124 | ✗ | PetscCall(BV_IPMatMult(X,y)); | |
| 125 | ✗ | z = X->Bx; | |
| 126 | } | ||
| 127 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
32 | PetscCall(VecGetArray(z,&py)); |
| 128 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
32 | PetscCall(BVDotVec_BLAS_Private(X,X->n,X->k-X->l,x->array+(X->nc+X->l)*X->ld,X->ld,py,m,PETSC_FALSE)); |
| 129 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
32 | PetscCall(VecRestoreArray(z,&py)); |
| 130 |
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); |
| 131 | } | ||
| 132 | |||
| 133 | 2549 | static PetscErrorCode BVScale_Contiguous(BV bv,PetscInt j,PetscScalar alpha) | |
| 134 | { | ||
| 135 | 2549 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 136 | |||
| 137 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
2549 | PetscFunctionBegin; |
| 138 |
2/14✓ Branch 0 taken 6 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.
|
2549 | if (!bv->n) PetscFunctionReturn(PETSC_SUCCESS); |
| 139 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
2549 | if (PetscUnlikely(j<0)) PetscCall(BVScale_BLAS_Private(bv,(bv->k-bv->l)*bv->ld,ctx->array+(bv->nc+bv->l)*bv->ld,alpha)); |
| 140 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2502 | else PetscCall(BVScale_BLAS_Private(bv,bv->n,ctx->array+(bv->nc+j)*bv->ld,alpha)); |
| 141 |
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.
|
634 | PetscFunctionReturn(PETSC_SUCCESS); |
| 142 | } | ||
| 143 | |||
| 144 | 900 | static PetscErrorCode BVNorm_Contiguous(BV bv,PetscInt j,NormType type,PetscReal *val) | |
| 145 | { | ||
| 146 | 900 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 147 | |||
| 148 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
900 | PetscFunctionBegin; |
| 149 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
900 | if (PetscUnlikely(j<0)) PetscCall(BVNorm_LAPACK_Private(bv,bv->n,bv->k-bv->l,ctx->array+(bv->nc+bv->l)*bv->ld,bv->ld,type,val,ctx->mpi)); |
| 150 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
341 | else PetscCall(BVNorm_LAPACK_Private(bv,bv->n,1,ctx->array+(bv->nc+j)*bv->ld,bv->ld,type,val,ctx->mpi)); |
| 151 |
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.
|
218 | PetscFunctionReturn(PETSC_SUCCESS); |
| 152 | } | ||
| 153 | |||
| 154 | 16 | static PetscErrorCode BVNorm_Local_Contiguous(BV bv,PetscInt j,NormType type,PetscReal *val) | |
| 155 | { | ||
| 156 | 16 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 157 | |||
| 158 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
16 | PetscFunctionBegin; |
| 159 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 8 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.
|
16 | if (PetscUnlikely(j<0)) PetscCall(BVNorm_LAPACK_Private(bv,bv->n,bv->k-bv->l,ctx->array+(bv->nc+bv->l)*bv->ld,bv->ld,type,val,PETSC_FALSE)); |
| 160 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | else PetscCall(BVNorm_LAPACK_Private(bv,bv->n,1,ctx->array+(bv->nc+j)*bv->ld,bv->ld,type,val,PETSC_FALSE)); |
| 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.
|
4 | PetscFunctionReturn(PETSC_SUCCESS); |
| 162 | } | ||
| 163 | |||
| 164 | 36 | static PetscErrorCode BVNormalize_Contiguous(BV bv,PetscScalar *eigi) | |
| 165 | { | ||
| 166 | 36 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 167 | 36 | PetscScalar *wi=NULL; | |
| 168 | |||
| 169 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
36 | PetscFunctionBegin; |
| 170 |
2/2✓ Branch 0 taken 4 times.
✓ Branch 1 taken 8 times.
|
36 | if (eigi) wi = eigi+bv->l; |
| 171 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(BVNormalize_LAPACK_Private(bv,bv->n,bv->k-bv->l,ctx->array+(bv->nc+bv->l)*bv->ld,bv->ld,wi,ctx->mpi)); |
| 172 |
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); |
| 173 | } | ||
| 174 | |||
| 175 | 824 | static PetscErrorCode BVMatMult_Contiguous(BV V,Mat A,BV W) | |
| 176 | { | ||
| 177 | 824 | BV_CONTIGUOUS *v = (BV_CONTIGUOUS*)V->data,*w = (BV_CONTIGUOUS*)W->data; | |
| 178 | 824 | PetscInt j; | |
| 179 | 824 | Mat Vmat,Wmat; | |
| 180 | |||
| 181 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
824 | PetscFunctionBegin; |
| 182 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
824 | if (V->vmm) { |
| 183 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(BVGetMat(V,&Vmat)); |
| 184 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(BVGetMat(W,&Wmat)); |
| 185 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(MatProductCreateWithMat(A,Vmat,NULL,Wmat)); |
| 186 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(MatProductSetType(Wmat,MATPRODUCT_AB)); |
| 187 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(MatProductSetFromOptions(Wmat)); |
| 188 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(MatProductSymbolic(Wmat)); |
| 189 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(MatProductNumeric(Wmat)); |
| 190 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(MatProductClear(Wmat)); |
| 191 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(BVRestoreMat(V,&Vmat)); |
| 192 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
784 | PetscCall(BVRestoreMat(W,&Wmat)); |
| 193 | } else { | ||
| 194 |
7/8✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 6 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
|
304 | for (j=0;j<V->k-V->l;j++) PetscCall(MatMult(A,v->V[V->nc+V->l+j],w->V[W->nc+W->l+j])); |
| 195 | } | ||
| 196 |
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.
|
206 | PetscFunctionReturn(PETSC_SUCCESS); |
| 197 | } | ||
| 198 | |||
| 199 | 776 | static PetscErrorCode BVCopy_Contiguous(BV V,BV W) | |
| 200 | { | ||
| 201 | 776 | BV_CONTIGUOUS *v = (BV_CONTIGUOUS*)V->data,*w = (BV_CONTIGUOUS*)W->data; | |
| 202 | 776 | PetscInt j; | |
| 203 | |||
| 204 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
776 | PetscFunctionBegin; |
| 205 |
7/8✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 6 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
|
5808 | for (j=0;j<V->k-V->l;j++) PetscCall(PetscArraycpy(w->array+(W->nc+W->l+j)*W->ld,v->array+(V->nc+V->l+j)*V->ld,V->n)); |
| 206 |
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); |
| 207 | } | ||
| 208 | |||
| 209 | 16 | static PetscErrorCode BVCopyColumn_Contiguous(BV V,PetscInt j,PetscInt i) | |
| 210 | { | ||
| 211 | 16 | BV_CONTIGUOUS *v = (BV_CONTIGUOUS*)V->data; | |
| 212 | |||
| 213 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
16 | PetscFunctionBegin; |
| 214 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(PetscArraycpy(v->array+(V->nc+i)*V->ld,v->array+(V->nc+j)*V->ld,V->n)); |
| 215 |
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); |
| 216 | } | ||
| 217 | |||
| 218 | 64 | static PetscErrorCode BVResize_Contiguous(BV bv,PetscInt m,PetscBool copy) | |
| 219 | { | ||
| 220 | 64 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 221 | 64 | PetscInt j,bs; | |
| 222 | 64 | PetscScalar *newarray; | |
| 223 | 64 | Vec *newV; | |
| 224 | 64 | char str[50]; | |
| 225 | |||
| 226 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
64 | PetscFunctionBegin; |
| 227 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
64 | PetscCall(PetscLayoutGetBlockSize(bv->map,&bs)); |
| 228 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
64 | PetscCall(PetscCalloc1(m*bv->ld,&newarray)); |
| 229 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
64 | PetscCall(PetscMalloc1(m,&newV)); |
| 230 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
752 | for (j=0;j<m;j++) { |
| 231 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
688 | if (ctx->mpi) PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)bv),bs,bv->n,PETSC_DECIDE,newarray+j*bv->ld,newV+j)); |
| 232 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
688 | else PetscCall(VecCreateSeqWithArray(PetscObjectComm((PetscObject)bv),bs,bv->n,newarray+j*bv->ld,newV+j)); |
| 233 | } | ||
| 234 |
1/2✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
|
64 | if (((PetscObject)bv)->name) { |
| 235 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
752 | for (j=0;j<m;j++) { |
| 236 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
688 | PetscCall(PetscSNPrintf(str,sizeof(str),"%s_%" PetscInt_FMT,((PetscObject)bv)->name,j)); |
| 237 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
688 | PetscCall(PetscObjectSetName((PetscObject)newV[j],str)); |
| 238 | } | ||
| 239 | } | ||
| 240 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
64 | if (copy) PetscCall(PetscArraycpy(newarray,ctx->array,PetscMin(m,bv->m)*bv->ld)); |
| 241 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
64 | PetscCall(VecDestroyVecs(bv->m,&ctx->V)); |
| 242 | 64 | ctx->V = newV; | |
| 243 |
5/8✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
64 | PetscCall(PetscFree(ctx->array)); |
| 244 | 64 | ctx->array = newarray; | |
| 245 |
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.
|
64 | PetscFunctionReturn(PETSC_SUCCESS); |
| 246 | } | ||
| 247 | |||
| 248 | 20126 | static PetscErrorCode BVGetColumn_Contiguous(BV bv,PetscInt j,Vec *v) | |
| 249 | { | ||
| 250 | 20126 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 251 | 20126 | PetscInt l; | |
| 252 | |||
| 253 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
20126 | PetscFunctionBegin; |
| 254 |
3/4✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 8 times.
|
20126 | l = BVAvailableVec; |
| 255 | 20126 | bv->cv[l] = ctx->V[bv->nc+j]; | |
| 256 |
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.
|
20126 | PetscFunctionReturn(PETSC_SUCCESS); |
| 257 | } | ||
| 258 | |||
| 259 | 20126 | static PetscErrorCode BVRestoreColumn_Contiguous(BV bv,PetscInt j,Vec *v) | |
| 260 | { | ||
| 261 | 20126 | PetscInt l; | |
| 262 | |||
| 263 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
20126 | PetscFunctionBegin; |
| 264 | 20126 | l = (j==bv->ci[0])? 0: 1; | |
| 265 | 20126 | bv->cv[l] = NULL; | |
| 266 |
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.
|
20126 | PetscFunctionReturn(PETSC_SUCCESS); |
| 267 | } | ||
| 268 | |||
| 269 | 1960 | static PetscErrorCode BVGetArray_Contiguous(BV bv,PetscScalar **a) | |
| 270 | { | ||
| 271 | 1960 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 272 | |||
| 273 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1960 | PetscFunctionBegin; |
| 274 | 1960 | *a = ctx->array; | |
| 275 |
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.
|
1960 | PetscFunctionReturn(PETSC_SUCCESS); |
| 276 | } | ||
| 277 | |||
| 278 | 64 | static PetscErrorCode BVGetArrayRead_Contiguous(BV bv,const PetscScalar **a) | |
| 279 | { | ||
| 280 | 64 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 281 | |||
| 282 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
64 | PetscFunctionBegin; |
| 283 | 64 | *a = ctx->array; | |
| 284 |
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.
|
64 | PetscFunctionReturn(PETSC_SUCCESS); |
| 285 | } | ||
| 286 | |||
| 287 | 2674 | static PetscErrorCode BVDestroy_Contiguous(BV bv) | |
| 288 | { | ||
| 289 | 2674 | BV_CONTIGUOUS *ctx = (BV_CONTIGUOUS*)bv->data; | |
| 290 | |||
| 291 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
2674 | PetscFunctionBegin; |
| 292 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
2674 | if (!bv->issplit) { |
| 293 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2034 | PetscCall(VecDestroyVecs(bv->nc+bv->m,&ctx->V)); |
| 294 |
5/8✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
2034 | PetscCall(PetscFree(ctx->array)); |
| 295 | } | ||
| 296 |
5/8✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
2674 | PetscCall(PetscFree(bv->data)); |
| 297 |
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.
|
664 | PetscFunctionReturn(PETSC_SUCCESS); |
| 298 | } | ||
| 299 | |||
| 300 | 96 | static PetscErrorCode BVView_Contiguous(BV bv,PetscViewer viewer) | |
| 301 | { | ||
| 302 | 96 | PetscInt j; | |
| 303 | 96 | Vec v; | |
| 304 | 96 | PetscViewerFormat format; | |
| 305 | 96 | PetscBool isascii,ismatlab=PETSC_FALSE; | |
| 306 | 96 | const char *bvname,*name; | |
| 307 | |||
| 308 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
96 | PetscFunctionBegin; |
| 309 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
96 | PetscCall(PetscObjectTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii)); |
| 310 |
1/2✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
|
96 | if (isascii) { |
| 311 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
96 | PetscCall(PetscViewerGetFormat(viewer,&format)); |
| 312 |
8/14✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 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.
|
96 | if (format == PETSC_VIEWER_ASCII_INFO || format == PETSC_VIEWER_ASCII_INFO_DETAIL) PetscFunctionReturn(PETSC_SUCCESS); |
| 313 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
|
80 | if (format == PETSC_VIEWER_ASCII_MATLAB) ismatlab = PETSC_TRUE; |
| 314 | } | ||
| 315 | ✗ | if (ismatlab) { | |
| 316 | ✗ | PetscCall(PetscObjectGetName((PetscObject)bv,&bvname)); | |
| 317 | ✗ | PetscCall(PetscViewerASCIIPrintf(viewer,"%s=[];\n",bvname)); | |
| 318 | } | ||
| 319 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
384 | for (j=0;j<bv->m;j++) { |
| 320 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
304 | PetscCall(BVGetColumn(bv,j,&v)); |
| 321 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
304 | PetscCall(VecView(v,viewer)); |
| 322 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
|
304 | if (ismatlab) { |
| 323 | ✗ | PetscCall(PetscObjectGetName((PetscObject)v,&name)); | |
| 324 | ✗ | PetscCall(PetscViewerASCIIPrintf(viewer,"%s=[%s,%s];clear %s\n",bvname,bvname,name,name)); | |
| 325 | } | ||
| 326 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
304 | PetscCall(BVRestoreColumn(bv,j,&v)); |
| 327 | } | ||
| 328 |
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); |
| 329 | } | ||
| 330 | |||
| 331 | /*MC | ||
| 332 | BVCONTIGUOUS - BVCONTIGUOUS = "contiguous" - A basis vectors object represented | ||
| 333 | by an array of `Vec`s sharing a contiguous array of `PetscScalar` storing the | ||
| 334 | elements column-wise. | ||
| 335 | |||
| 336 | Level: beginner | ||
| 337 | |||
| 338 | .seealso: [](sec:bv), `BV`, `BVType`, `BVSetType()` | ||
| 339 | M*/ | ||
| 340 | |||
| 341 | 2674 | SLEPC_EXTERN PetscErrorCode BVCreate_Contiguous(BV bv) | |
| 342 | { | ||
| 343 | 2674 | BV_CONTIGUOUS *ctx; | |
| 344 | 2674 | PetscInt j,nloc,bs,lsplit,lda; | |
| 345 | 2674 | PetscBool seq,isdense; | |
| 346 | 2674 | PetscScalar *aa; | |
| 347 | 2674 | char str[50]; | |
| 348 | 2674 | PetscScalar *array; | |
| 349 | 2674 | BV parent; | |
| 350 | 2674 | Vec *Vpar; | |
| 351 | 2674 | MatType mtype; | |
| 352 | |||
| 353 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
2674 | PetscFunctionBegin; |
| 354 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2674 | PetscCall(PetscNew(&ctx)); |
| 355 | 2674 | bv->data = (void*)ctx; | |
| 356 | |||
| 357 |
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.
|
2674 | PetscCall(PetscStrcmp(bv->vtype,VECMPI,&ctx->mpi)); |
| 358 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
2674 | if (!ctx->mpi) { |
| 359 |
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.
|
478 | PetscCall(PetscStrcmp(bv->vtype,VECSEQ,&seq)); |
| 360 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
478 | PetscCheck(seq,PetscObjectComm((PetscObject)bv),PETSC_ERR_SUP,"Cannot create a contiguous BV from a non-standard vector type: %s",bv->vtype); |
| 361 | } | ||
| 362 | |||
| 363 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2674 | PetscCall(PetscLayoutGetLocalSize(bv->map,&nloc)); |
| 364 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2674 | PetscCall(PetscLayoutGetBlockSize(bv->map,&bs)); |
| 365 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2674 | PetscCall(BV_SetDefaultLD(bv,nloc)); |
| 366 | |||
| 367 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
2674 | if (PetscUnlikely(bv->issplit)) { |
| 368 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
640 | PetscCheck(bv->issplit>0,PetscObjectComm((PetscObject)bv),PETSC_ERR_SUP,"BVCONTIGUOUS does not support BVGetSplitRows()"); |
| 369 | /* split BV: share memory and Vecs of the parent BV */ | ||
| 370 | 640 | parent = bv->splitparent; | |
| 371 | 640 | lsplit = parent->lsplit; | |
| 372 | 640 | Vpar = ((BV_CONTIGUOUS*)parent->data)->V; | |
| 373 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
640 | ctx->V = (bv->issplit==1)? Vpar: Vpar+lsplit; |
| 374 | 640 | array = ((BV_CONTIGUOUS*)parent->data)->array; | |
| 375 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
640 | ctx->array = (bv->issplit==1)? array: array+lsplit*bv->ld; |
| 376 | } else { | ||
| 377 | /* regular BV: allocate memory and Vecs for the BV entries */ | ||
| 378 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2034 | PetscCall(PetscCalloc1(bv->m*bv->ld,&ctx->array)); |
| 379 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2034 | PetscCall(PetscMalloc1(bv->m,&ctx->V)); |
| 380 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
17516 | for (j=0;j<bv->m;j++) { |
| 381 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
15482 | if (ctx->mpi) PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)bv),bs,nloc,PETSC_DECIDE,ctx->array+j*bv->ld,ctx->V+j)); |
| 382 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
15482 | else PetscCall(VecCreateSeqWithArray(PetscObjectComm((PetscObject)bv),bs,nloc,ctx->array+j*bv->ld,ctx->V+j)); |
| 383 | } | ||
| 384 | } | ||
| 385 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
2674 | if (((PetscObject)bv)->name) { |
| 386 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
7600 | for (j=0;j<bv->m;j++) { |
| 387 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6736 | PetscCall(PetscSNPrintf(str,sizeof(str),"%s_%" PetscInt_FMT,((PetscObject)bv)->name,j)); |
| 388 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6736 | PetscCall(PetscObjectSetName((PetscObject)ctx->V[j],str)); |
| 389 | } | ||
| 390 | } | ||
| 391 | |||
| 392 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
2674 | if (PetscUnlikely(bv->Acreate)) { |
| 393 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(MatGetType(bv->Acreate,&mtype)); |
| 394 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(PetscStrcmpAny(mtype,&isdense,MATSEQDENSE,MATMPIDENSE,"")); |
| 395 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
16 | PetscCheck(isdense,PetscObjectComm((PetscObject)bv->Acreate),PETSC_ERR_SUP,"BVCONTIGUOUS requires a dense matrix in BVCreateFromMat()"); |
| 396 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(MatDenseGetArray(bv->Acreate,&aa)); |
| 397 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(MatDenseGetLDA(bv->Acreate,&lda)); |
| 398 |
7/8✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 6 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
|
144 | for (j=0;j<bv->m;j++) PetscCall(PetscArraycpy(ctx->array+j*bv->ld,aa+j*lda,bv->n)); |
| 399 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(MatDenseRestoreArray(bv->Acreate,&aa)); |
| 400 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(MatDestroy(&bv->Acreate)); |
| 401 | } | ||
| 402 | |||
| 403 | 2674 | bv->ops->mult = BVMult_Contiguous; | |
| 404 | 2674 | bv->ops->multvec = BVMultVec_Contiguous; | |
| 405 | 2674 | bv->ops->multinplace = BVMultInPlace_Contiguous; | |
| 406 | 2674 | bv->ops->multinplacetrans = BVMultInPlaceHermitianTranspose_Contiguous; | |
| 407 | 2674 | bv->ops->dot = BVDot_Contiguous; | |
| 408 | 2674 | bv->ops->dotvec = BVDotVec_Contiguous; | |
| 409 | 2674 | bv->ops->dotvec_local = BVDotVec_Local_Contiguous; | |
| 410 | 2674 | bv->ops->scale = BVScale_Contiguous; | |
| 411 | 2674 | bv->ops->norm = BVNorm_Contiguous; | |
| 412 | 2674 | bv->ops->norm_local = BVNorm_Local_Contiguous; | |
| 413 | 2674 | bv->ops->normalize = BVNormalize_Contiguous; | |
| 414 | 2674 | bv->ops->matmult = BVMatMult_Contiguous; | |
| 415 | 2674 | bv->ops->copy = BVCopy_Contiguous; | |
| 416 | 2674 | bv->ops->copycolumn = BVCopyColumn_Contiguous; | |
| 417 | 2674 | bv->ops->resize = BVResize_Contiguous; | |
| 418 | 2674 | bv->ops->getcolumn = BVGetColumn_Contiguous; | |
| 419 | 2674 | bv->ops->restorecolumn = BVRestoreColumn_Contiguous; | |
| 420 | 2674 | bv->ops->getarray = BVGetArray_Contiguous; | |
| 421 | 2674 | bv->ops->getarrayread = BVGetArrayRead_Contiguous; | |
| 422 | 2674 | bv->ops->getmat = BVGetMat_Default; | |
| 423 | 2674 | bv->ops->restoremat = BVRestoreMat_Default; | |
| 424 | 2674 | bv->ops->destroy = BVDestroy_Contiguous; | |
| 425 | 2674 | bv->ops->view = BVView_Contiguous; | |
| 426 |
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.
|
2674 | PetscFunctionReturn(PETSC_SUCCESS); |
| 427 | } | ||
| 428 |