| 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 | Tensor BV that is represented in compact form as V = (I otimes U) S | ||
| 12 | */ | ||
| 13 | |||
| 14 | #include <slepc/private/bvimpl.h> /*I "slepcbv.h" I*/ | ||
| 15 | #include <slepcblaslapack.h> | ||
| 16 | |||
| 17 | typedef struct { | ||
| 18 | BV U; /* first factor */ | ||
| 19 | Mat S; /* second factor */ | ||
| 20 | PetscScalar *qB; /* auxiliary matrix used in non-standard inner products */ | ||
| 21 | PetscScalar *sw; /* work space */ | ||
| 22 | PetscInt d; /* degree of the tensor BV */ | ||
| 23 | PetscInt ld; /* leading dimension of a single block in S */ | ||
| 24 | PetscInt puk; /* copy of the k value */ | ||
| 25 | Vec u; /* auxiliary work vector */ | ||
| 26 | } BV_TENSOR; | ||
| 27 | |||
| 28 | 13162 | static PetscErrorCode BVMultInPlace_Tensor(BV V,Mat Q,PetscInt s,PetscInt e) | |
| 29 | { | ||
| 30 | 13162 | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 31 | 13162 | PetscScalar *pS; | |
| 32 | 13162 | const PetscScalar *q; | |
| 33 | 13162 | PetscInt ldq,lds = ctx->ld*ctx->d; | |
| 34 | |||
| 35 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
13162 | PetscFunctionBegin; |
| 36 |
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.
|
13162 | PetscCall(MatDenseGetLDA(Q,&ldq)); |
| 37 |
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.
|
13162 | PetscCall(MatDenseGetArray(ctx->S,&pS)); |
| 38 |
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.
|
13162 | PetscCall(MatDenseGetArrayRead(Q,&q)); |
| 39 |
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.
|
13162 | PetscCall(BVMultInPlace_BLAS_Private(V,lds,V->k-V->l,s-V->l,e-V->l,pS+(V->nc+V->l)*lds,lds,q+V->l*ldq+V->l,ldq,PETSC_FALSE)); |
| 40 |
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.
|
13162 | PetscCall(MatDenseRestoreArrayRead(Q,&q)); |
| 41 |
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.
|
13162 | PetscCall(MatDenseRestoreArray(ctx->S,&pS)); |
| 42 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
2526 | PetscFunctionReturn(PETSC_SUCCESS); |
| 43 | } | ||
| 44 | |||
| 45 | ✗ | static PetscErrorCode BVMultInPlaceHermitianTranspose_Tensor(BV V,Mat Q,PetscInt s,PetscInt e) | |
| 46 | { | ||
| 47 | ✗ | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 48 | ✗ | PetscScalar *pS; | |
| 49 | ✗ | const PetscScalar *q; | |
| 50 | ✗ | PetscInt ldq,lds = ctx->ld*ctx->d; | |
| 51 | |||
| 52 | ✗ | PetscFunctionBegin; | |
| 53 | ✗ | PetscCall(MatDenseGetLDA(Q,&ldq)); | |
| 54 | ✗ | PetscCall(MatDenseGetArray(ctx->S,&pS)); | |
| 55 | ✗ | PetscCall(MatDenseGetArrayRead(Q,&q)); | |
| 56 | ✗ | PetscCall(BVMultInPlace_BLAS_Private(V,lds,V->k-V->l,s-V->l,e-V->l,pS+(V->nc+V->l)*lds,lds,q+V->l*ldq+V->l,ldq,PETSC_TRUE)); | |
| 57 | ✗ | PetscCall(MatDenseRestoreArrayRead(Q,&q)); | |
| 58 | ✗ | PetscCall(MatDenseRestoreArray(ctx->S,&pS)); | |
| 59 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 60 | } | ||
| 61 | |||
| 62 | 104 | static PetscErrorCode BVDot_Tensor(BV X,BV Y,Mat M) | |
| 63 | { | ||
| 64 | 104 | BV_TENSOR *x = (BV_TENSOR*)X->data,*y = (BV_TENSOR*)Y->data; | |
| 65 | 104 | PetscScalar *m; | |
| 66 | 104 | const PetscScalar *px,*py; | |
| 67 | 104 | PetscInt ldm,lds = x->ld*x->d; | |
| 68 | |||
| 69 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
104 | PetscFunctionBegin; |
| 70 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
104 | PetscCheck(x->U==y->U,PetscObjectComm((PetscObject)X),PETSC_ERR_SUP,"BVDot() in BVTENSOR requires that both operands have the same U factor"); |
| 71 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
104 | PetscCheck(lds==y->ld*y->d,PetscObjectComm((PetscObject)X),PETSC_ERR_ARG_SIZ,"Mismatching dimensions ld*d %" PetscInt_FMT " %" PetscInt_FMT,lds,y->ld*y->d); |
| 72 |
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.
|
104 | PetscCall(MatDenseGetLDA(M,&ldm)); |
| 73 |
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.
|
104 | PetscCall(MatDenseGetArrayRead(x->S,&px)); |
| 74 |
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.
|
104 | PetscCall(MatDenseGetArrayRead(y->S,&py)); |
| 75 |
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.
|
104 | PetscCall(MatDenseGetArray(M,&m)); |
| 76 |
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.
|
104 | PetscCall(BVDot_BLAS_Private(X,Y->k-Y->l,X->k-X->l,lds,py+(Y->nc+Y->l)*lds,lds,px+(X->nc+X->l)*lds,lds,m+X->l*ldm+Y->l,ldm,PETSC_FALSE)); |
| 77 |
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.
|
104 | PetscCall(MatDenseRestoreArray(M,&m)); |
| 78 |
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.
|
104 | PetscCall(MatDenseRestoreArrayRead(x->S,&px)); |
| 79 |
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.
|
104 | PetscCall(MatDenseRestoreArrayRead(y->S,&py)); |
| 80 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
16 | PetscFunctionReturn(PETSC_SUCCESS); |
| 81 | } | ||
| 82 | |||
| 83 | 158806 | static PetscErrorCode BVScale_Tensor(BV bv,PetscInt j,PetscScalar alpha) | |
| 84 | { | ||
| 85 | 158806 | BV_TENSOR *ctx = (BV_TENSOR*)bv->data; | |
| 86 | 158806 | PetscScalar *pS; | |
| 87 | 158806 | PetscInt lds = ctx->ld*ctx->d; | |
| 88 | |||
| 89 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
158806 | PetscFunctionBegin; |
| 90 |
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.
|
158806 | PetscCall(MatDenseGetArray(ctx->S,&pS)); |
| 91 |
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.
|
158806 | if (PetscUnlikely(j<0)) PetscCall(BVScale_BLAS_Private(bv,(bv->k-bv->l)*lds,pS+(bv->nc+bv->l)*lds,alpha)); |
| 92 |
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.
|
158806 | else PetscCall(BVScale_BLAS_Private(bv,lds,pS+(bv->nc+j)*lds,alpha)); |
| 93 |
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.
|
158806 | PetscCall(MatDenseRestoreArray(ctx->S,&pS)); |
| 94 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
30161 | PetscFunctionReturn(PETSC_SUCCESS); |
| 95 | } | ||
| 96 | |||
| 97 | 104 | static PetscErrorCode BVNorm_Tensor(BV bv,PetscInt j,NormType type,PetscReal *val) | |
| 98 | { | ||
| 99 | 104 | BV_TENSOR *ctx = (BV_TENSOR*)bv->data; | |
| 100 | 104 | const PetscScalar *pS; | |
| 101 | 104 | PetscInt lds = ctx->ld*ctx->d; | |
| 102 | |||
| 103 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
104 | PetscFunctionBegin; |
| 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.
|
104 | PetscCall(MatDenseGetArrayRead(ctx->S,&pS)); |
| 105 |
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.
|
104 | if (j<0) PetscCall(BVNorm_LAPACK_Private(bv,lds,bv->k-bv->l,pS+(bv->nc+bv->l)*lds,lds,type,val,PETSC_FALSE)); |
| 106 | ✗ | else PetscCall(BVNorm_LAPACK_Private(bv,lds,1,pS+(bv->nc+j)*lds,lds,type,val,PETSC_FALSE)); | |
| 107 |
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.
|
104 | PetscCall(MatDenseRestoreArrayRead(ctx->S,&pS)); |
| 108 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
16 | PetscFunctionReturn(PETSC_SUCCESS); |
| 109 | } | ||
| 110 | |||
| 111 | 12868 | static PetscErrorCode BVCopyColumn_Tensor(BV V,PetscInt j,PetscInt i) | |
| 112 | { | ||
| 113 | 12868 | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 114 | 12868 | PetscScalar *pS; | |
| 115 | 12868 | PetscInt lds = ctx->ld*ctx->d; | |
| 116 | |||
| 117 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
12868 | PetscFunctionBegin; |
| 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.
|
12868 | PetscCall(MatDenseGetArray(ctx->S,&pS)); |
| 119 |
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.
|
12868 | PetscCall(PetscArraycpy(pS+(V->nc+i)*lds,pS+(V->nc+j)*lds,lds)); |
| 120 |
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.
|
12868 | PetscCall(MatDenseRestoreArray(ctx->S,&pS)); |
| 121 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
2472 | PetscFunctionReturn(PETSC_SUCCESS); |
| 122 | } | ||
| 123 | |||
| 124 | 462914 | static PetscErrorCode BVTensorNormColumn(BV bv,PetscInt j,PetscReal *norm) | |
| 125 | { | ||
| 126 | 462914 | BV_TENSOR *ctx = (BV_TENSOR*)bv->data; | |
| 127 | 462914 | PetscBLASInt one=1,lds_; | |
| 128 | 462914 | PetscScalar sone=1.0,szero=0.0,*x,dot; | |
| 129 | 462914 | const PetscScalar *S; | |
| 130 | 462914 | PetscReal alpha=1.0,scale=0.0,aval; | |
| 131 | 462914 | PetscInt i,lds,ld=ctx->ld; | |
| 132 | |||
| 133 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
462914 | PetscFunctionBegin; |
| 134 | 462914 | lds = ld*ctx->d; | |
| 135 |
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.
|
462914 | PetscCall(MatDenseGetArrayRead(ctx->S,&S)); |
| 136 |
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.
|
462914 | PetscCall(PetscBLASIntCast(lds,&lds_)); |
| 137 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
462914 | if (PetscUnlikely(ctx->qB)) { |
| 138 | 49862 | x = ctx->sw; | |
| 139 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
49862 | PetscCallBLAS("BLASgemv",BLASgemv_("N",&lds_,&lds_,&sone,ctx->qB,&lds_,S+j*lds,&one,&szero,x,&one)); |
| 140 | 49862 | dot = PetscRealPart(BLASdot_(&lds_,S+j*lds,&one,x,&one)); | |
| 141 |
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.
|
49862 | PetscCall(BV_SafeSqrt(bv,dot,norm)); |
| 142 | } else { | ||
| 143 | /* Compute *norm = BLASnrm2_(&lds_,S+j*lds,&one); */ | ||
| 144 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
413052 | if (lds==1) *norm = PetscAbsScalar(S[j*lds]); |
| 145 | else { | ||
| 146 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
76940970 | for (i=0;i<lds;i++) { |
| 147 | 76527918 | aval = PetscAbsScalar(S[i+j*lds]); | |
| 148 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
76527918 | if (aval!=0.0) { |
| 149 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
55383085 | if (PetscUnlikely(scale<aval)) { |
| 150 | 2569987 | alpha = 1.0 + alpha*PetscSqr(scale/aval); | |
| 151 | 2569987 | scale = aval; | |
| 152 | 52813098 | } else alpha += PetscSqr(aval/scale); | |
| 153 | } | ||
| 154 | } | ||
| 155 | 413052 | *norm = scale*PetscSqrtReal(alpha); | |
| 156 | } | ||
| 157 | } | ||
| 158 |
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.
|
462914 | PetscCall(MatDenseRestoreArrayRead(ctx->S,&S)); |
| 159 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
89534 | PetscFunctionReturn(PETSC_SUCCESS); |
| 160 | } | ||
| 161 | |||
| 162 | 305120 | static PetscErrorCode BVOrthogonalizeGS1_Tensor(BV bv,PetscInt k,Vec v,PetscBool *which,PetscScalar *h,PetscScalar *c,PetscReal *onorm,PetscReal *norm) | |
| 163 | { | ||
| 164 | 305120 | BV_TENSOR *ctx = (BV_TENSOR*)bv->data; | |
| 165 | 305120 | PetscScalar *pS,*cc,*x,dot,sonem=-1.0,sone=1.0,szero=0.0; | |
| 166 | 305120 | PetscInt i,lds = ctx->ld*ctx->d; | |
| 167 | 305120 | PetscBLASInt lds_,k_,one=1; | |
| 168 | 305120 | const PetscScalar *omega; | |
| 169 | |||
| 170 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
305120 | PetscFunctionBegin; |
| 171 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
305120 | PetscCheck(!v,PetscObjectComm((PetscObject)bv),PETSC_ERR_SUP,"Orthogonalization against an external vector is not allowed in BVTENSOR"); |
| 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.
|
305120 | PetscCall(MatDenseGetArray(ctx->S,&pS)); |
| 173 |
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.
|
305120 | if (!c) PetscCall(VecGetArray(bv->buffer,&cc)); |
| 174 | ✗ | else cc = c; | |
| 175 |
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.
|
305120 | PetscCall(PetscBLASIntCast(lds,&lds_)); |
| 176 |
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.
|
305120 | PetscCall(PetscBLASIntCast(k,&k_)); |
| 177 | |||
| 178 |
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.
|
305120 | if (onorm) PetscCall(BVTensorNormColumn(bv,k,onorm)); |
| 179 | |||
| 180 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
305120 | if (ctx->qB) x = ctx->sw; |
| 181 | 206130 | else x = pS+k*lds; | |
| 182 | |||
| 183 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
305120 | if (PetscUnlikely(bv->orthog_type==BV_ORTHOG_MGS)) { /* modified Gram-Schmidt */ |
| 184 | |||
| 185 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
1160 | if (PetscUnlikely(bv->indef)) { /* signature */ |
| 186 | ✗ | PetscCall(VecGetArrayRead(bv->omega,&omega)); | |
| 187 | } | ||
| 188 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
16450 | for (i=-bv->nc;i<k;i++) { |
| 189 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
15290 | if (which && i>=0 && !which[i]) continue; |
| 190 |
1/26✗ 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.
✗ 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.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
|
15290 | if (ctx->qB) PetscCallBLAS("BLASgemv",BLASgemv_("N",&lds_,&lds_,&sone,ctx->qB,&lds_,pS+k*lds,&one,&szero,x,&one)); |
| 191 | /* c_i = (s_k, s_i) */ | ||
| 192 | 15290 | dot = PetscRealPart(BLASdot_(&lds_,pS+i*lds,&one,x,&one)); | |
| 193 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
15290 | if (bv->indef) dot /= PetscRealPart(omega[i]); |
| 194 |
5/6✓ Branch 0 taken 6 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
21406 | PetscCall(BV_SetValue(bv,i,0,cc,dot)); |
| 195 | /* s_k = s_k - c_i s_i */ | ||
| 196 | 15290 | dot = -dot; | |
| 197 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
15290 | PetscCallBLAS("BLASaxpy",BLASaxpy_(&lds_,&dot,pS+i*lds,&one,pS+k*lds,&one)); |
| 198 | } | ||
| 199 |
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.
|
1160 | if (PetscUnlikely(bv->indef)) PetscCall(VecRestoreArrayRead(bv->omega,&omega)); |
| 200 | |||
| 201 | } else { /* classical Gram-Schmidt */ | ||
| 202 |
14/26✓ 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.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
|
303960 | if (ctx->qB) PetscCallBLAS("BLASgemv",BLASgemv_("N",&lds_,&lds_,&sone,ctx->qB,&lds_,pS+k*lds,&one,&szero,x,&one)); |
| 203 | |||
| 204 | /* cc = S_{0:k-1}^* s_k */ | ||
| 205 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
303960 | PetscCallBLAS("BLASgemv",BLASgemv_("C",&lds_,&k_,&sone,pS,&lds_,x,&one,&szero,cc,&one)); |
| 206 | |||
| 207 | /* s_k = s_k - S_{0:k-1} cc */ | ||
| 208 |
7/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 6 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 4 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
349916 | if (PetscUnlikely(bv->indef)) PetscCall(BV_ApplySignature(bv,k,cc,PETSC_TRUE)); |
| 209 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
303960 | PetscCallBLAS("BLASgemv",BLASgemv_("N",&lds_,&k_,&sonem,pS,&lds_,cc,&one,&sone,pS+k*lds,&one)); |
| 210 |
7/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 6 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 4 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
349916 | if (PetscUnlikely(bv->indef)) PetscCall(BV_ApplySignature(bv,k,cc,PETSC_FALSE)); |
| 211 | } | ||
| 212 | |||
| 213 |
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.
|
305120 | if (norm) PetscCall(BVTensorNormColumn(bv,k,norm)); |
| 214 |
5/6✓ Branch 0 taken 6 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
433381 | PetscCall(BV_AddCoefficients(bv,k,h,cc)); |
| 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.
|
305120 | PetscCall(MatDenseRestoreArray(ctx->S,&pS)); |
| 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.
|
305120 | PetscCall(VecRestoreArray(bv->buffer,&cc)); |
| 217 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
57946 | PetscFunctionReturn(PETSC_SUCCESS); |
| 218 | } | ||
| 219 | |||
| 220 | 104 | static PetscErrorCode BVView_Tensor(BV bv,PetscViewer viewer) | |
| 221 | { | ||
| 222 | 104 | BV_TENSOR *ctx = (BV_TENSOR*)bv->data; | |
| 223 | 104 | PetscViewerFormat format; | |
| 224 | 104 | PetscBool isascii; | |
| 225 | 104 | const char *bvname,*uname,*sname; | |
| 226 | |||
| 227 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
104 | PetscFunctionBegin; |
| 228 |
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.
|
104 | PetscCall(PetscObjectTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii)); |
| 229 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
104 | if (isascii) { |
| 230 |
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.
|
104 | PetscCall(PetscViewerGetFormat(viewer,&format)); |
| 231 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
104 | if (format == PETSC_VIEWER_ASCII_INFO || format == PETSC_VIEWER_ASCII_INFO_DETAIL) { |
| 232 |
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.
|
104 | PetscCall(PetscViewerASCIIPrintf(viewer,"number of tensor blocks (degree): %" PetscInt_FMT "\n",ctx->d)); |
| 233 |
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.
|
104 | PetscCall(PetscViewerASCIIPrintf(viewer,"number of columns of U factor: %" PetscInt_FMT "\n",ctx->ld)); |
| 234 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
16 | PetscFunctionReturn(PETSC_SUCCESS); |
| 235 | } | ||
| 236 | ✗ | PetscCall(BVView(ctx->U,viewer)); | |
| 237 | ✗ | PetscCall(MatView(ctx->S,viewer)); | |
| 238 | ✗ | if (format == PETSC_VIEWER_ASCII_MATLAB) { | |
| 239 | ✗ | PetscCall(PetscObjectGetName((PetscObject)bv,&bvname)); | |
| 240 | ✗ | PetscCall(PetscObjectGetName((PetscObject)ctx->U,&uname)); | |
| 241 | ✗ | PetscCall(PetscObjectGetName((PetscObject)ctx->S,&sname)); | |
| 242 | ✗ | PetscCall(PetscViewerASCIIPrintf(viewer,"%s=kron(eye(%" PetscInt_FMT "),%s)*%s(:,1:%" PetscInt_FMT ");\n",bvname,ctx->d,uname,sname,bv->k)); | |
| 243 | } | ||
| 244 | } else { | ||
| 245 | ✗ | PetscCall(BVView(ctx->U,viewer)); | |
| 246 | ✗ | PetscCall(MatView(ctx->S,viewer)); | |
| 247 | } | ||
| 248 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 249 | } | ||
| 250 | |||
| 251 | 62080 | static PetscErrorCode BVTensorUpdateMatrix(BV V,PetscInt ini,PetscInt end) | |
| 252 | { | ||
| 253 | 62080 | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 254 | 62080 | PetscInt i,j,r,c,l,k,ld=ctx->ld,lds=ctx->d*ctx->ld; | |
| 255 | 62080 | PetscScalar *qB,*sqB; | |
| 256 | 62080 | Vec u; | |
| 257 | 62080 | Mat A; | |
| 258 | |||
| 259 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
62080 | PetscFunctionBegin; |
| 260 |
9/16✓ 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 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
62080 | if (!V->matrix) PetscFunctionReturn(PETSC_SUCCESS); |
| 261 | 49673 | l = ctx->U->l; k = ctx->U->k; | |
| 262 | /* update inner product matrix */ | ||
| 263 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
49673 | if (!ctx->qB) { |
| 264 |
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.
|
367 | PetscCall(PetscCalloc2(lds*lds,&ctx->qB,lds,&ctx->sw)); |
| 265 |
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.
|
367 | PetscCall(BVCreateVec(ctx->U,&ctx->u)); |
| 266 | } | ||
| 267 | 49673 | ctx->U->l = 0; | |
| 268 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
149019 | for (r=0;r<ctx->d;r++) { |
| 269 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
248365 | for (c=0;c<=r;c++) { |
| 270 |
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.
|
149019 | PetscCall(MatNestGetSubMat(V->matrix,r,c,&A)); |
| 271 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
149019 | if (A) { |
| 272 | 99626 | qB = ctx->qB+c*ld*lds+r*ld; | |
| 273 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
200558 | for (i=ini;i<end;i++) { |
| 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.
|
100932 | PetscCall(BVGetColumn(ctx->U,i,&u)); |
| 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.
|
100932 | PetscCall(MatMult(A,u,ctx->u)); |
| 276 | 100932 | ctx->U->k = i+1; | |
| 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.
|
100932 | PetscCall(BVDotVec(ctx->U,ctx->u,qB+i*lds)); |
| 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.
|
100932 | PetscCall(BVRestoreColumn(ctx->U,i,&u)); |
| 279 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
2435218 | for (j=0;j<i;j++) qB[i+j*lds] = PetscConj(qB[j+i*lds]); |
| 280 | 44980 | qB[i*lds+i] = PetscRealPart(qB[i+i*lds]); | |
| 281 | } | ||
| 282 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
99626 | if (PetscUnlikely(c!=r)) { |
| 283 | 560 | sqB = ctx->qB+r*ld*lds+c*ld; | |
| 284 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
7960 | for (i=ini;i<end;i++) for (j=0;j<=i;j++) { |
| 285 | 6840 | sqB[i+j*lds] = PetscConj(qB[j+i*lds]); | |
| 286 | 6840 | sqB[j+i*lds] = qB[j+i*lds]; | |
| 287 | } | ||
| 288 | } | ||
| 289 | } | ||
| 290 | } | ||
| 291 | } | ||
| 292 | 49673 | ctx->U->l = l; ctx->U->k = k; | |
| 293 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
49673 | PetscFunctionReturn(PETSC_SUCCESS); |
| 294 | } | ||
| 295 | |||
| 296 | 1729 | static PetscErrorCode BVTensorBuildFirstColumn_Tensor(BV V,PetscInt k) | |
| 297 | { | ||
| 298 | 1729 | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 299 | 1729 | PetscInt i,nq=0; | |
| 300 | 1729 | PetscScalar *pS,*omega; | |
| 301 | 1729 | PetscReal norm; | |
| 302 | 1729 | PetscBool breakdown=PETSC_FALSE; | |
| 303 | |||
| 304 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
1729 | PetscFunctionBegin; |
| 305 |
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.
|
1729 | PetscCall(MatDenseGetArray(ctx->S,&pS)); |
| 306 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
9368 | for (i=0;i<ctx->d;i++) { |
| 307 |
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.
|
7639 | if (i>=k) PetscCall(BVSetRandomColumn(ctx->U,nq)); |
| 308 |
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.
|
370 | else PetscCall(BVCopyColumn(ctx->U,i,nq)); |
| 309 |
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.
|
7639 | PetscCall(BVOrthogonalizeColumn(ctx->U,nq,pS+i*ctx->ld,&norm,&breakdown)); |
| 310 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 5 times.
|
7639 | if (!breakdown) { |
| 311 |
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.
|
7634 | PetscCall(BVScaleColumn(ctx->U,nq,1.0/norm)); |
| 312 | 7634 | pS[nq+i*ctx->ld] = norm; | |
| 313 | 7634 | nq++; | |
| 314 | } | ||
| 315 | } | ||
| 316 |
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.
|
1729 | PetscCall(MatDenseRestoreArray(ctx->S,&pS)); |
| 317 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1729 | PetscCheck(nq,PetscObjectComm((PetscObject)V),PETSC_ERR_ARG_SIZ,"Cannot build first column of tensor BV; U should contain k=%" PetscInt_FMT " nonzero columns",k); |
| 318 |
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.
|
1729 | PetscCall(BVTensorUpdateMatrix(V,0,nq)); |
| 319 |
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.
|
1729 | PetscCall(BVTensorNormColumn(V,0,&norm)); |
| 320 |
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.
|
1729 | PetscCall(BVScale_Tensor(V,0,1.0/norm)); |
| 321 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1729 | if (V->indef) { |
| 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.
|
367 | PetscCall(BV_AllocateSignature(V)); |
| 323 |
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.
|
367 | PetscCall(VecGetArray(V->omega,&omega)); |
| 324 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
367 | omega[0] = (norm<0.0)? -1.0: 1.0; |
| 325 |
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.
|
367 | PetscCall(VecRestoreArray(V->omega,&omega)); |
| 326 | } | ||
| 327 | /* set active columns */ | ||
| 328 | 1729 | ctx->U->l = 0; | |
| 329 | 1729 | ctx->U->k = nq; | |
| 330 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
1729 | PetscFunctionReturn(PETSC_SUCCESS); |
| 331 | } | ||
| 332 | |||
| 333 | /*@ | ||
| 334 | BVTensorBuildFirstColumn - Builds the first column of the tensor basis vectors | ||
| 335 | `V` from the data contained in the first `k` columns of `U`. | ||
| 336 | |||
| 337 | Collective | ||
| 338 | |||
| 339 | Input Parameters: | ||
| 340 | + V - the basis vectors context | ||
| 341 | - k - the number of columns of `U` with relevant information | ||
| 342 | |||
| 343 | Notes: | ||
| 344 | At most $d$ columns are considered, where $d$ is the degree of the tensor `BV`. | ||
| 345 | Given $V = (I \otimes U) S$, this function computes the first column of $V$, that | ||
| 346 | is, it computes the coefficients of the first column of $S$ by orthogonalizing | ||
| 347 | the first $d$ columns of $U$. If `k` is less than $d$ (or linearly dependent columns | ||
| 348 | are found) then additional random columns are used. | ||
| 349 | |||
| 350 | The computed column has unit norm. | ||
| 351 | |||
| 352 | Level: advanced | ||
| 353 | |||
| 354 | .seealso: [](sec:bv), `BVTENSOR`, `BVCreateTensor()`, `BVTensorGetDegree()` | ||
| 355 | @*/ | ||
| 356 | 1729 | PetscErrorCode BVTensorBuildFirstColumn(BV V,PetscInt k) | |
| 357 | { | ||
| 358 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
1729 | PetscFunctionBegin; |
| 359 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
|
1729 | PetscValidHeaderSpecific(V,BV_CLASSID,1); |
| 360 |
31/70✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✗ Branch 16 not taken.
✓ Branch 17 taken 2 times.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✓ Branch 35 taken 2 times.
✓ Branch 36 taken 2 times.
✗ Branch 37 not taken.
✓ Branch 38 taken 2 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 2 times.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✓ Branch 47 taken 2 times.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✓ Branch 51 taken 2 times.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✓ Branch 56 taken 2 times.
✗ Branch 57 not taken.
✓ Branch 58 taken 2 times.
✗ Branch 59 not taken.
✓ Branch 60 taken 2 times.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✓ Branch 63 taken 2 times.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✓ Branch 67 taken 2 times.
✗ Branch 68 not taken.
✗ Branch 69 not taken.
|
1729 | PetscValidLogicalCollectiveInt(V,k,2); |
| 361 |
9/16✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 10 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
1729 | PetscUseMethod(V,"BVTensorBuildFirstColumn_C",(BV,PetscInt),(V,k)); |
| 362 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
1729 | PetscFunctionReturn(PETSC_SUCCESS); |
| 363 | } | ||
| 364 | |||
| 365 | 13193 | static PetscErrorCode BVTensorCompress_Tensor(BV V,PetscInt newc) | |
| 366 | { | ||
| 367 | 13193 | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 368 | 13193 | PetscInt nwu=0,nnc,nrow,lwa,r,c; | |
| 369 | 13193 | PetscInt i,j,k,n,lds=ctx->ld*ctx->d,deg=ctx->d,lock,cs1=V->k,rs1=ctx->U->k,rk=0,offu; | |
| 370 | 13193 | PetscScalar *S,*M,*Z,*pQ,*SS,*SS2,t,sone=1.0,zero=0.0,mone=-1.0,*p,*tau,*work,*qB,*sqB; | |
| 371 | 13193 | PetscReal *sg,tol,*rwork; | |
| 372 | 13193 | PetscBLASInt ld_,cs1_,rs1_,cs1tdeg,n_,lw_,newc_,newctdeg,nnc_,nrow_,nnctdeg,lds_,rk_; | |
| 373 | 13193 | Mat Q,A; | |
| 374 | |||
| 375 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
13193 | PetscFunctionBegin; |
| 376 |
9/16✓ 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 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
13193 | if (!cs1) PetscFunctionReturn(PETSC_SUCCESS); |
| 377 | 13141 | lwa = 6*ctx->ld*lds+2*cs1; | |
| 378 | 13141 | n = PetscMin(rs1,deg*cs1); | |
| 379 | 13141 | lock = ctx->U->l; | |
| 380 | 13141 | nnc = cs1-lock-newc; | |
| 381 | 13141 | nrow = rs1-lock; | |
| 382 |
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.
|
13141 | PetscCall(PetscCalloc6(deg*newc*nnc,&SS,newc*nnc,&SS2,(rs1+lock+newc)*n,&pQ,deg*rs1,&tau,lwa,&work,6*n,&rwork)); |
| 383 | 13141 | offu = lock*(rs1+1); | |
| 384 | 13141 | M = work+nwu; | |
| 385 | 13141 | nwu += rs1*cs1*deg; | |
| 386 | 13141 | sg = rwork; | |
| 387 | 13141 | Z = work+nwu; | |
| 388 | 13141 | nwu += deg*cs1*n; | |
| 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.
|
13141 | PetscCall(PetscBLASIntCast(n,&n_)); |
| 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.
|
13141 | PetscCall(PetscBLASIntCast(nnc,&nnc_)); |
| 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.
|
13141 | PetscCall(PetscBLASIntCast(cs1,&cs1_)); |
| 392 |
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.
|
13141 | PetscCall(PetscBLASIntCast(rs1,&rs1_)); |
| 393 |
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.
|
13141 | PetscCall(PetscBLASIntCast(newc,&newc_)); |
| 394 |
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.
|
13141 | PetscCall(PetscBLASIntCast(newc*deg,&newctdeg)); |
| 395 |
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.
|
13141 | PetscCall(PetscBLASIntCast(nnc*deg,&nnctdeg)); |
| 396 |
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.
|
13141 | PetscCall(PetscBLASIntCast(cs1*deg,&cs1tdeg)); |
| 397 |
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.
|
13141 | PetscCall(PetscBLASIntCast(lwa-nwu,&lw_)); |
| 398 |
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.
|
13141 | PetscCall(PetscBLASIntCast(nrow,&nrow_)); |
| 399 |
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.
|
13141 | PetscCall(PetscBLASIntCast(lds,&lds_)); |
| 400 |
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.
|
13141 | PetscCall(MatDenseGetArray(ctx->S,&S)); |
| 401 | |||
| 402 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
13141 | if (newc>0) { |
| 403 | /* truncate columns associated with new converged eigenpairs */ | ||
| 404 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
20648 | for (j=0;j<deg;j++) { |
| 405 |
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.
|
71117 | for (i=lock;i<lock+newc;i++) PetscCall(PetscArraycpy(M+(i-lock+j*newc)*nrow,S+i*lds+j*ctx->ld+lock,nrow)); |
| 406 | } | ||
| 407 |
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.
|
3153 | PetscCall(PetscFPTrapPush(PETSC_FP_TRAP_OFF)); |
| 408 | #if !PetscDefined(USE_COMPLEX) | ||
| 409 |
16/32✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 1 times.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 1 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 1 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
|
1581 | PetscCallLAPACKInfo("LAPACKgesvd",LAPACKgesvd_("S","S",&nrow_,&newctdeg,M,&nrow_,sg,pQ+offu,&rs1_,Z,&n_,work+nwu,&lw_,&info)); |
| 410 | #else | ||
| 411 |
16/32✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 1 times.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 1 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 1 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
|
1572 | PetscCallLAPACKInfo("LAPACKgesvd",LAPACKgesvd_("S","S",&nrow_,&newctdeg,M,&nrow_,sg,pQ+offu,&rs1_,Z,&n_,work+nwu,&lw_,rwork+n,&info)); |
| 412 | #endif | ||
| 413 |
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.
|
3153 | PetscCall(PetscFPTrapPop()); |
| 414 | /* SVD has rank min(newc,nrow) */ | ||
| 415 | 3153 | rk = PetscMin(newc,nrow); | |
| 416 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
19386 | for (i=0;i<rk;i++) { |
| 417 | 16233 | t = sg[i]; | |
| 418 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
16233 | PetscCallBLAS("BLASscal",BLASscal_(&newctdeg,&t,Z+i,&n_)); |
| 419 | } | ||
| 420 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
20648 | for (i=0;i<deg;i++) { |
| 421 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
71117 | for (j=lock;j<lock+newc;j++) { |
| 422 |
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.
|
53622 | PetscCall(PetscArraycpy(S+j*lds+i*ctx->ld+lock,Z+(newc*i+j-lock)*n,rk)); |
| 423 |
4/6✓ Branch 0 taken 3 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
59038 | PetscCall(PetscArrayzero(S+j*lds+i*ctx->ld+lock+rk,(ctx->ld-lock-rk))); |
| 424 | } | ||
| 425 | } | ||
| 426 | /* | ||
| 427 | update columns associated with non-converged vectors, orthogonalize | ||
| 428 | against pQ so that next M has rank nnc+d-1 instead of nrow+d-1 | ||
| 429 | */ | ||
| 430 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
20648 | for (i=0;i<deg;i++) { |
| 431 |
15/26✓ Branch 0 taken 10 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ 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 taken 2 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
|
17495 | PetscCallBLAS("BLASgemm",BLASgemm_("C","N",&newc_,&nnc_,&nrow_,&sone,pQ+offu,&rs1_,S+(lock+newc)*lds+i*ctx->ld+lock,&lds_,&zero,PetscSafePointerPlusOffset(SS,i*nnc*newc),&newc_)); |
| 432 |
15/26✓ Branch 0 taken 10 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ 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 taken 2 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
|
17495 | PetscCallBLAS("BLASgemm",BLASgemm_("N","N",&nrow_,&nnc_,&newc_,&mone,pQ+offu,&rs1_,PetscSafePointerPlusOffset(SS,i*nnc*newc),&newc_,&sone,S+(lock+newc)*lds+i*ctx->ld+lock,&lds_)); |
| 433 | /* repeat orthogonalization step */ | ||
| 434 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
17495 | PetscCallBLAS("BLASgemm",BLASgemm_("C","N",&newc_,&nnc_,&nrow_,&sone,pQ+offu,&rs1_,S+(lock+newc)*lds+i*ctx->ld+lock,&lds_,&zero,SS2,&newc_)); |
| 435 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
17495 | PetscCallBLAS("BLASgemm",BLASgemm_("N","N",&nrow_,&nnc_,&newc_,&mone,pQ+offu,&rs1_,SS2,&newc_,&sone,S+(lock+newc)*lds+i*ctx->ld+lock,&lds_)); |
| 436 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
172560 | for (j=0;j<newc*nnc;j++) *(SS+i*newc*nnc+j) += SS2[j]; |
| 437 | } | ||
| 438 | } | ||
| 439 | |||
| 440 | /* truncate columns associated with non-converged eigenpairs */ | ||
| 441 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
69484 | for (j=0;j<deg;j++) { |
| 442 |
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.
|
572561 | for (i=lock+newc;i<cs1;i++) PetscCall(PetscArraycpy(M+(i-lock-newc+j*nnc)*nrow,S+i*lds+j*ctx->ld+lock,nrow)); |
| 443 | } | ||
| 444 |
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.
|
13141 | PetscCall(PetscFPTrapPush(PETSC_FP_TRAP_OFF)); |
| 445 | #if !PetscDefined(USE_COMPLEX) | ||
| 446 |
16/32✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 1 times.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 1 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 1 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
|
6659 | PetscCallLAPACKInfo("LAPACKgesvd",LAPACKgesvd_("S","S",&nrow_,&nnctdeg,M,&nrow_,sg,pQ+offu+newc*rs1,&rs1_,Z,&n_,work+nwu,&lw_,&info)); |
| 447 | #else | ||
| 448 |
16/32✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 1 times.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 1 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 1 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
|
6482 | PetscCallLAPACKInfo("LAPACKgesvd",LAPACKgesvd_("S","S",&nrow_,&nnctdeg,M,&nrow_,sg,pQ+offu+newc*rs1,&rs1_,Z,&n_,work+nwu,&lw_,rwork+n,&info)); |
| 449 | #endif | ||
| 450 |
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.
|
13141 | PetscCall(PetscFPTrapPop()); |
| 451 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
13141 | tol = PetscMax(rs1,deg*cs1)*PETSC_MACHINE_EPSILON*sg[0]; |
| 452 | 13141 | rk = 0; | |
| 453 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
301923 | for (i=0;i<PetscMin(nrow,nnctdeg);i++) if (sg[i]>tol) rk++; |
| 454 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
13141 | rk = PetscMin(nnc+deg-1,rk); |
| 455 | /* the SVD has rank (at most) nnc+deg-1 */ | ||
| 456 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
186555 | for (i=0;i<rk;i++) { |
| 457 | 173414 | t = sg[i]; | |
| 458 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
173414 | PetscCallBLAS("BLASscal",BLASscal_(&nnctdeg,&t,Z+i,&n_)); |
| 459 | } | ||
| 460 | /* update S */ | ||
| 461 |
4/6✓ Branch 0 taken 3 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13141 | PetscCall(PetscArrayzero(S+cs1*lds,(V->m-cs1)*lds)); |
| 462 | 13141 | k = ctx->ld-lock-newc-rk; | |
| 463 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
69484 | for (i=0;i<deg;i++) { |
| 464 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
572561 | for (j=lock+newc;j<cs1;j++) { |
| 465 |
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.
|
516218 | PetscCall(PetscArraycpy(S+j*lds+i*ctx->ld+lock+newc,Z+(nnc*i+j-lock-newc)*n,rk)); |
| 466 |
4/6✓ Branch 0 taken 3 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
568807 | PetscCall(PetscArrayzero(S+j*lds+i*ctx->ld+lock+newc+rk,k)); |
| 467 | } | ||
| 468 | } | ||
| 469 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
13141 | if (newc>0) { |
| 470 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
20648 | for (i=0;i<deg;i++) { |
| 471 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
17495 | p = PetscSafePointerPlusOffset(SS,i*nnc*newc); |
| 472 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
112772 | for (j=lock+newc;j<cs1;j++) { |
| 473 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
250342 | for (k=0;k<newc;k++) S[j*lds+i*ctx->ld+lock+k] = *(p++); |
| 474 | } | ||
| 475 | } | ||
| 476 | } | ||
| 477 | |||
| 478 | /* orthogonalize pQ */ | ||
| 479 | 13141 | rk = rk+newc; | |
| 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.
|
13141 | PetscCall(PetscBLASIntCast(rk,&rk_)); |
| 481 |
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.
|
13141 | PetscCall(PetscBLASIntCast(cs1-lock,&nnc_)); |
| 482 |
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.
|
13141 | PetscCall(PetscFPTrapPush(PETSC_FP_TRAP_OFF)); |
| 483 |
16/32✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 8 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
|
13141 | PetscCallLAPACKInfo("LAPACKgeqrf",LAPACKgeqrf_(&nrow_,&rk_,pQ+offu,&rs1_,tau,work+nwu,&lw_,&info)); |
| 484 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
69484 | for (i=0;i<deg;i++) { |
| 485 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
56343 | PetscCallBLAS("BLAStrmm",BLAStrmm_("L","U","N","N",&rk_,&nnc_,&sone,pQ+offu,&rs1_,S+lock*lds+lock+i*ctx->ld,&lds_)); |
| 486 | } | ||
| 487 |
16/32✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 8 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
|
13141 | PetscCallLAPACKInfo("LAPACKorgqr",LAPACKorgqr_(&nrow_,&rk_,&rk_,pQ+offu,&rs1_,tau,work+nwu,&lw_,&info)); |
| 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.
|
13141 | PetscCall(PetscFPTrapPop()); |
| 489 | |||
| 490 | /* update vectors U(:,idx) = U*Q(:,idx) */ | ||
| 491 | 13141 | rk = rk+lock; | |
| 492 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
28501 | for (i=0;i<lock;i++) pQ[i*(1+rs1)] = 1.0; |
| 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.
|
13141 | PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,rs1,rk,pQ,&Q)); |
| 494 | 13141 | ctx->U->k = rs1; | |
| 495 |
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.
|
13141 | PetscCall(BVMultInPlace(ctx->U,Q,lock,rk)); |
| 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.
|
13141 | PetscCall(MatDestroy(&Q)); |
| 497 | |||
| 498 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
13141 | if (ctx->qB) { |
| 499 | /* update matrix qB */ | ||
| 500 |
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.
|
3990 | PetscCall(PetscBLASIntCast(ctx->ld,&ld_)); |
| 501 |
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.
|
3990 | PetscCall(PetscBLASIntCast(rk,&rk_)); |
| 502 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11970 | for (r=0;r<ctx->d;r++) { |
| 503 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
19950 | for (c=0;c<=r;c++) { |
| 504 |
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.
|
11970 | PetscCall(MatNestGetSubMat(V->matrix,r,c,&A)); |
| 505 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11970 | if (A) { |
| 506 | 8010 | qB = ctx->qB+r*ctx->ld+c*ctx->ld*lds; | |
| 507 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
8010 | PetscCallBLAS("BLASgemm",BLASgemm_("N","N",&rs1_,&rk_,&rs1_,&sone,qB,&lds_,pQ,&rs1_,&zero,work+nwu,&rs1_)); |
| 508 |
12/24✓ 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.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
|
8010 | PetscCallBLAS("BLASgemm",BLASgemm_("C","N",&rk_,&rk_,&rs1_,&sone,pQ,&rs1_,work+nwu,&rs1_,&zero,qB,&lds_)); |
| 509 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
141734 | for (i=0;i<rk;i++) { |
| 510 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1417694 | for (j=0;j<i;j++) qB[i+j*lds] = PetscConj(qB[j+i*lds]); |
| 511 | 133724 | qB[i+i*lds] = PetscRealPart(qB[i+i*lds]); | |
| 512 | } | ||
| 513 |
7/8✓ Branch 0 taken 3 times.
✓ Branch 1 taken 7 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.
|
133060 | for (i=rk;i<ctx->ld;i++) PetscCall(PetscArrayzero(qB+i*lds,ctx->ld)); |
| 514 |
7/8✓ Branch 0 taken 3 times.
✓ Branch 1 taken 7 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.
|
141734 | for (i=0;i<rk;i++) PetscCall(PetscArrayzero(qB+i*lds+rk,(ctx->ld-rk))); |
| 515 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
8010 | if (c!=r) { |
| 516 | 60 | sqB = ctx->qB+r*ctx->ld*lds+c*ctx->ld; | |
| 517 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
25260 | for (i=0;i<ctx->ld;i++) for (j=0;j<ctx->ld;j++) sqB[i+j*lds] = PetscConj(qB[j+i*lds]); |
| 518 | } | ||
| 519 | } | ||
| 520 | } | ||
| 521 | } | ||
| 522 | } | ||
| 523 | |||
| 524 | /* free work space */ | ||
| 525 |
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.
|
13141 | PetscCall(PetscFree6(SS,SS2,pQ,tau,work,rwork)); |
| 526 |
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.
|
13141 | PetscCall(MatDenseRestoreArray(ctx->S,&S)); |
| 527 | |||
| 528 | /* set active columns */ | ||
| 529 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
13141 | if (newc) ctx->U->l += newc; |
| 530 | 13141 | ctx->U->k = rk; | |
| 531 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
13141 | PetscFunctionReturn(PETSC_SUCCESS); |
| 532 | } | ||
| 533 | |||
| 534 | /*@ | ||
| 535 | BVTensorCompress - Updates the $U$ and $S$ factors of the tensor basis vectors | ||
| 536 | object `V` by means of an SVD, removing redundant information. | ||
| 537 | |||
| 538 | Collective | ||
| 539 | |||
| 540 | Input Parameters: | ||
| 541 | + V - the tensor basis vectors context | ||
| 542 | - newc - additional columns to be locked | ||
| 543 | |||
| 544 | Notes: | ||
| 545 | This function is typically used when restarting Krylov solvers. Truncating a | ||
| 546 | tensor `BV` $V = (I \otimes U) S$ to its leading columns amounts to keeping the | ||
| 547 | leading columns of $S$. However, to effectively reduce the size of the | ||
| 548 | decomposition, it is necessary to compress it in a way that fewer columns of | ||
| 549 | $U$ are employed. This can be achieved by means of an update that involves the | ||
| 550 | SVD of the low-rank matrix $[S_0, S_1, \dots, S_{d-1}]$, where $S_i$ are the pieces of $S$. | ||
| 551 | |||
| 552 | If `newc` is nonzero, then `newc` columns are added to the leading columns of `V`. | ||
| 553 | This means that the corresponding columns of the $U$ and $S$ factors will remain | ||
| 554 | invariant in subsequent operations. | ||
| 555 | |||
| 556 | Level: advanced | ||
| 557 | |||
| 558 | .seealso: [](sec:bv), `BVTENSOR`, `BVCreateTensor()`, `BVSetActiveColumns()` | ||
| 559 | @*/ | ||
| 560 | 13193 | PetscErrorCode BVTensorCompress(BV V,PetscInt newc) | |
| 561 | { | ||
| 562 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
13193 | PetscFunctionBegin; |
| 563 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
|
13193 | PetscValidHeaderSpecific(V,BV_CLASSID,1); |
| 564 |
31/70✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✗ Branch 16 not taken.
✓ Branch 17 taken 2 times.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✓ Branch 35 taken 2 times.
✓ Branch 36 taken 2 times.
✗ Branch 37 not taken.
✓ Branch 38 taken 2 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 2 times.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✓ Branch 47 taken 2 times.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✓ Branch 51 taken 2 times.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✓ Branch 56 taken 2 times.
✗ Branch 57 not taken.
✓ Branch 58 taken 2 times.
✗ Branch 59 not taken.
✓ Branch 60 taken 2 times.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✓ Branch 63 taken 2 times.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✓ Branch 67 taken 2 times.
✗ Branch 68 not taken.
✗ Branch 69 not taken.
|
13193 | PetscValidLogicalCollectiveInt(V,newc,2); |
| 565 |
9/16✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 10 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
13193 | PetscUseMethod(V,"BVTensorCompress_C",(BV,PetscInt),(V,newc)); |
| 566 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
13193 | PetscFunctionReturn(PETSC_SUCCESS); |
| 567 | } | ||
| 568 | |||
| 569 | 3442 | static PetscErrorCode BVTensorGetDegree_Tensor(BV bv,PetscInt *d) | |
| 570 | { | ||
| 571 | 3442 | BV_TENSOR *ctx = (BV_TENSOR*)bv->data; | |
| 572 | |||
| 573 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
3442 | PetscFunctionBegin; |
| 574 | 3442 | *d = ctx->d; | |
| 575 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
3442 | PetscFunctionReturn(PETSC_SUCCESS); |
| 576 | } | ||
| 577 | |||
| 578 | /*@ | ||
| 579 | BVTensorGetDegree - Returns the number of blocks (degree) of the tensor `BV`. | ||
| 580 | |||
| 581 | Not Collective | ||
| 582 | |||
| 583 | Input Parameter: | ||
| 584 | . bv - the basis vectors context | ||
| 585 | |||
| 586 | Output Parameter: | ||
| 587 | . d - the degree | ||
| 588 | |||
| 589 | Level: advanced | ||
| 590 | |||
| 591 | .seealso: [](sec:bv), `BVTENSOR`, `BVCreateTensor()` | ||
| 592 | @*/ | ||
| 593 | 3442 | PetscErrorCode BVTensorGetDegree(BV bv,PetscInt *d) | |
| 594 | { | ||
| 595 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
3442 | PetscFunctionBegin; |
| 596 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
|
3442 | PetscValidHeaderSpecific(bv,BV_CLASSID,1); |
| 597 |
2/8✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
3442 | PetscAssertPointer(d,2); |
| 598 |
9/16✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 10 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
3442 | PetscUseMethod(bv,"BVTensorGetDegree_C",(BV,PetscInt*),(bv,d)); |
| 599 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
3442 | PetscFunctionReturn(PETSC_SUCCESS); |
| 600 | } | ||
| 601 | |||
| 602 | 60351 | static PetscErrorCode BVTensorGetFactors_Tensor(BV V,BV *U,Mat *S) | |
| 603 | { | ||
| 604 | 60351 | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 605 | |||
| 606 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
60351 | PetscFunctionBegin; |
| 607 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
60351 | PetscCheck(ctx->puk==-1,PetscObjectComm((PetscObject)V),PETSC_ERR_ORDER,"Previous call to BVTensonGetFactors without a BVTensorRestoreFactors call"); |
| 608 | 60351 | ctx->puk = ctx->U->k; | |
| 609 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
60351 | if (U) *U = ctx->U; |
| 610 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
60351 | if (S) *S = ctx->S; |
| 611 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
10973 | PetscFunctionReturn(PETSC_SUCCESS); |
| 612 | } | ||
| 613 | |||
| 614 | /*@ | ||
| 615 | BVTensorGetFactors - Returns the two factors involved in the definition of the | ||
| 616 | tensor basis vectors object, $V = (I \otimes U) S$. | ||
| 617 | |||
| 618 | Logically Collective | ||
| 619 | |||
| 620 | Input Parameter: | ||
| 621 | . V - the basis vectors context | ||
| 622 | |||
| 623 | Output Parameters: | ||
| 624 | + U - the `BV` factor | ||
| 625 | - S - the `Mat` factor | ||
| 626 | |||
| 627 | Notes: | ||
| 628 | The returned factors are references (not copies) of the internal factors, | ||
| 629 | so modifying them will change the tensor `BV` as well. Some operations of the | ||
| 630 | tensor `BV` assume that `U` has orthonormal columns, so if the user modifies `U` | ||
| 631 | this restriction must be taken into account. | ||
| 632 | |||
| 633 | The returned factors must not be destroyed. `BVTensorRestoreFactors()` must | ||
| 634 | be called when they are no longer needed. | ||
| 635 | |||
| 636 | Pass `NULL` for any of the arguments that is not needed. | ||
| 637 | |||
| 638 | Level: advanced | ||
| 639 | |||
| 640 | .seealso: [](sec:bv), `BVTENSOR`, `BVCreateTensor()`, `BVTensorRestoreFactors()` | ||
| 641 | @*/ | ||
| 642 | 60351 | PetscErrorCode BVTensorGetFactors(BV V,BV *U,Mat *S) | |
| 643 | { | ||
| 644 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
60351 | PetscFunctionBegin; |
| 645 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
|
60351 | PetscValidHeaderSpecific(V,BV_CLASSID,1); |
| 646 |
9/16✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 10 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
60351 | PetscUseMethod(V,"BVTensorGetFactors_C",(BV,BV*,Mat*),(V,U,S)); |
| 647 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
60351 | PetscFunctionReturn(PETSC_SUCCESS); |
| 648 | } | ||
| 649 | |||
| 650 | 60351 | static PetscErrorCode BVTensorRestoreFactors_Tensor(BV V,BV *U,Mat *S) | |
| 651 | { | ||
| 652 | 60351 | BV_TENSOR *ctx = (BV_TENSOR*)V->data; | |
| 653 | |||
| 654 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
60351 | PetscFunctionBegin; |
| 655 |
2/4✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
60351 | PetscCall(PetscObjectStateIncrease((PetscObject)V)); |
| 656 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
60351 | if (U) *U = NULL; |
| 657 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
60351 | if (S) *S = NULL; |
| 658 |
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.
|
60351 | PetscCall(BVTensorUpdateMatrix(V,ctx->puk,ctx->U->k)); |
| 659 | 60351 | ctx->puk = -1; | |
| 660 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
60351 | PetscFunctionReturn(PETSC_SUCCESS); |
| 661 | } | ||
| 662 | |||
| 663 | /*@ | ||
| 664 | BVTensorRestoreFactors - Restore the two factors that were obtained with | ||
| 665 | `BVTensorGetFactors()`. | ||
| 666 | |||
| 667 | Logically Collective | ||
| 668 | |||
| 669 | Input Parameters: | ||
| 670 | + V - the basis vectors context | ||
| 671 | . U - the `BV` factor (or `NULL`) | ||
| 672 | - S - the `Mat` factor (or `NULL`) | ||
| 673 | |||
| 674 | Notes: | ||
| 675 | The arguments must match the corresponding call to `BVTensorGetFactors()`. | ||
| 676 | |||
| 677 | Level: advanced | ||
| 678 | |||
| 679 | .seealso: [](sec:bv), `BVTENSOR`, `BVTensorGetFactors()` | ||
| 680 | @*/ | ||
| 681 | 60351 | PetscErrorCode BVTensorRestoreFactors(BV V,BV *U,Mat *S) | |
| 682 | { | ||
| 683 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
60351 | PetscFunctionBegin; |
| 684 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
|
60351 | PetscValidHeaderSpecific(V,BV_CLASSID,1); |
| 685 |
5/18✓ Branch 0 taken 2 times.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
|
60351 | if (U) PetscValidHeaderSpecific(*U,BV_CLASSID,2); |
| 686 |
4/18✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
|
60351 | if (S) PetscValidHeaderSpecific(*S,MAT_CLASSID,3); |
| 687 |
9/16✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 10 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
60351 | PetscUseMethod(V,"BVTensorRestoreFactors_C",(BV,BV*,Mat*),(V,U,S)); |
| 688 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
60351 | PetscFunctionReturn(PETSC_SUCCESS); |
| 689 | } | ||
| 690 | |||
| 691 | 1813 | static PetscErrorCode BVDestroy_Tensor(BV bv) | |
| 692 | { | ||
| 693 | 1813 | BV_TENSOR *ctx = (BV_TENSOR*)bv->data; | |
| 694 | |||
| 695 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
1813 | PetscFunctionBegin; |
| 696 |
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.
|
1813 | PetscCall(BVDestroy(&ctx->U)); |
| 697 |
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.
|
1813 | PetscCall(MatDestroy(&ctx->S)); |
| 698 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1813 | if (ctx->u) { |
| 699 |
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.
|
367 | PetscCall(PetscFree2(ctx->qB,ctx->sw)); |
| 700 |
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.
|
367 | PetscCall(VecDestroy(&ctx->u)); |
| 701 | } | ||
| 702 |
6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 2 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.
|
1813 | PetscCall(PetscFree(bv->data)); |
| 703 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorBuildFirstColumn_C",NULL)); |
| 704 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorCompress_C",NULL)); |
| 705 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorGetDegree_C",NULL)); |
| 706 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorGetFactors_C",NULL)); |
| 707 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorRestoreFactors_C",NULL)); |
| 708 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
343 | PetscFunctionReturn(PETSC_SUCCESS); |
| 709 | } | ||
| 710 | |||
| 711 | /*MC | ||
| 712 | BVTENSOR - BVTENSOR = "tensor" - A special basis vectors object represented | ||
| 713 | in compact form as $V = (I \otimes U) S$. | ||
| 714 | |||
| 715 | Level: beginner | ||
| 716 | |||
| 717 | Note: | ||
| 718 | This is not a general-purpose `BV` and cannot be used interchangeably with other | ||
| 719 | `BV` types. It is intended to store structured Krylov bases for eigensolvers | ||
| 720 | based on linearization. | ||
| 721 | |||
| 722 | .seealso: [](sec:bv), `BV`, `BVType`, `BVSetType()`, `BVCreateTensor()` | ||
| 723 | M*/ | ||
| 724 | |||
| 725 | 1813 | SLEPC_EXTERN PetscErrorCode BVCreate_Tensor(BV bv) | |
| 726 | { | ||
| 727 | 1813 | BV_TENSOR *ctx; | |
| 728 | |||
| 729 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
1813 | PetscFunctionBegin; |
| 730 |
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.
|
1813 | PetscCall(PetscNew(&ctx)); |
| 731 | 1813 | bv->data = (void*)ctx; | |
| 732 | 1813 | ctx->puk = -1; | |
| 733 | |||
| 734 | 1813 | bv->ops->multinplace = BVMultInPlace_Tensor; | |
| 735 | 1813 | bv->ops->multinplacetrans = BVMultInPlaceHermitianTranspose_Tensor; | |
| 736 | 1813 | bv->ops->dot = BVDot_Tensor; | |
| 737 | 1813 | bv->ops->scale = BVScale_Tensor; | |
| 738 | 1813 | bv->ops->norm = BVNorm_Tensor; | |
| 739 | 1813 | bv->ops->copycolumn = BVCopyColumn_Tensor; | |
| 740 | 1813 | bv->ops->gramschmidt = BVOrthogonalizeGS1_Tensor; | |
| 741 | 1813 | bv->ops->destroy = BVDestroy_Tensor; | |
| 742 | 1813 | bv->ops->view = BVView_Tensor; | |
| 743 | |||
| 744 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorBuildFirstColumn_C",BVTensorBuildFirstColumn_Tensor)); |
| 745 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorCompress_C",BVTensorCompress_Tensor)); |
| 746 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorGetDegree_C",BVTensorGetDegree_Tensor)); |
| 747 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorGetFactors_C",BVTensorGetFactors_Tensor)); |
| 748 |
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.
|
1813 | PetscCall(PetscObjectComposeFunction((PetscObject)bv,"BVTensorRestoreFactors_C",BVTensorRestoreFactors_Tensor)); |
| 749 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
343 | PetscFunctionReturn(PETSC_SUCCESS); |
| 750 | } | ||
| 751 | |||
| 752 | /*@ | ||
| 753 | BVCreateTensor - Creates a tensor `BV` that is represented in compact form | ||
| 754 | as $V = (I \otimes U) S$, where $U$ has orthonormal columns. | ||
| 755 | |||
| 756 | Collective | ||
| 757 | |||
| 758 | Input Parameters: | ||
| 759 | + U - a basis vectors object | ||
| 760 | - d - the number of blocks (degree) of the tensor `BV` | ||
| 761 | |||
| 762 | Output Parameter: | ||
| 763 | . V - the new basis vectors context | ||
| 764 | |||
| 765 | Notes: | ||
| 766 | The new basis vectors object is $V = (I \otimes U) S$, where $\otimes$ denotes | ||
| 767 | the Kronecker product, $I$ is the identity matrix of order $d$, and $S$ is a | ||
| 768 | sequential matrix allocated internally. This compact representation is | ||
| 769 | used, e.g., to represent the Krylov basis generated with the linearization | ||
| 770 | of a matrix polynomial of degree $d$. | ||
| 771 | |||
| 772 | The size of `V` (number of rows) is equal to `d` times `n`, where `n` is the size | ||
| 773 | of `U`. Matrix $S$ has $d \times m$ rows and $m-d+1$ columns, where `m` is | ||
| 774 | the number of columns of `U`, so `m` should be at least `d`. | ||
| 775 | |||
| 776 | The communicator of `V` will be the same as `U`. | ||
| 777 | |||
| 778 | On input, the content of `U` is irrelevant. Alternatively, it may contain | ||
| 779 | some nonzero columns that will be used by `BVTensorBuildFirstColumn()`. | ||
| 780 | |||
| 781 | Level: advanced | ||
| 782 | |||
| 783 | .seealso: [](sec:bv), `BVTensorGetDegree()`, `BVTensorGetFactors()`, `BVTensorBuildFirstColumn()` | ||
| 784 | @*/ | ||
| 785 | 1813 | PetscErrorCode BVCreateTensor(BV U,PetscInt d,BV *V) | |
| 786 | { | ||
| 787 | 1813 | PetscBool match; | |
| 788 | 1813 | PetscInt n,N,m; | |
| 789 | 1813 | VecType vtype; | |
| 790 | 1813 | BV_TENSOR *ctx; | |
| 791 | |||
| 792 |
2/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
|
1813 | PetscFunctionBegin; |
| 793 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
|
1813 | PetscValidHeaderSpecific(U,BV_CLASSID,1); |
| 794 |
31/70✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✗ Branch 16 not taken.
✓ Branch 17 taken 2 times.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✓ Branch 35 taken 2 times.
✓ Branch 36 taken 2 times.
✗ Branch 37 not taken.
✓ Branch 38 taken 2 times.
✗ Branch 39 not taken.
✓ Branch 40 taken 2 times.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✓ Branch 47 taken 2 times.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✓ Branch 51 taken 2 times.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✓ Branch 56 taken 2 times.
✗ Branch 57 not taken.
✓ Branch 58 taken 2 times.
✗ Branch 59 not taken.
✓ Branch 60 taken 2 times.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✓ Branch 63 taken 2 times.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✓ Branch 67 taken 2 times.
✗ Branch 68 not taken.
✗ Branch 69 not taken.
|
1813 | PetscValidLogicalCollectiveInt(U,d,2); |
| 795 |
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.
|
1813 | PetscCall(PetscObjectTypeCompare((PetscObject)U,BVTENSOR,&match)); |
| 796 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1813 | PetscCheck(!match,PetscObjectComm((PetscObject)U),PETSC_ERR_SUP,"U cannot be of type tensor"); |
| 797 | |||
| 798 |
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.
|
1813 | PetscCall(BVCreate(PetscObjectComm((PetscObject)U),V)); |
| 799 |
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.
|
1813 | PetscCall(BVGetSizes(U,&n,&N,&m)); |
| 800 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1813 | PetscCheck(m>=d,PetscObjectComm((PetscObject)U),PETSC_ERR_ARG_SIZ,"U has %" PetscInt_FMT " columns, it should have at least d=%" PetscInt_FMT,m,d); |
| 801 |
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.
|
1813 | PetscCall(BVSetSizes(*V,d*n,d*N,m-d+1)); |
| 802 |
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.
|
1813 | PetscCall(BVGetVecType(U,&vtype)); |
| 803 |
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.
|
1813 | PetscCall(BVSetVecType(*V,vtype)); |
| 804 |
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.
|
1813 | PetscCall(PetscObjectChangeTypeName((PetscObject)*V,BVTENSOR)); |
| 805 |
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.
|
1813 | PetscCall(PetscLogEventBegin(BV_Create,*V,0,0,0)); |
| 806 |
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.
|
1813 | PetscCall(BVCreate_Tensor(*V)); |
| 807 |
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.
|
1813 | PetscCall(PetscLogEventEnd(BV_Create,*V,0,0,0)); |
| 808 | |||
| 809 | 1813 | ctx = (BV_TENSOR*)(*V)->data; | |
| 810 | 1813 | ctx->U = U; | |
| 811 | 1813 | ctx->d = d; | |
| 812 | 1813 | ctx->ld = m; | |
| 813 |
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.
|
1813 | PetscCall(PetscObjectReference((PetscObject)U)); |
| 814 |
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.
|
1813 | PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,d*m,m-d+1,NULL,&ctx->S)); |
| 815 |
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.
|
1813 | PetscCall(PetscObjectSetName((PetscObject)ctx->S,"S")); |
| 816 | |||
| 817 | /* Copy user-provided attributes of U */ | ||
| 818 | 1813 | (*V)->orthog_type = U->orthog_type; | |
| 819 | 1813 | (*V)->orthog_ref = U->orthog_ref; | |
| 820 | 1813 | (*V)->orthog_eta = U->orthog_eta; | |
| 821 | 1813 | (*V)->orthog_block = U->orthog_block; | |
| 822 | 1813 | (*V)->vmm = U->vmm; | |
| 823 | 1813 | (*V)->rrandom = U->rrandom; | |
| 824 |
7/14✓ 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 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
1813 | PetscFunctionReturn(PETSC_SUCCESS); |
| 825 | } | ||
| 826 |