| 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 | Common subroutines for all Krylov-type solvers | ||
| 12 | */ | ||
| 13 | |||
| 14 | #include <slepc/private/epsimpl.h> | ||
| 15 | #include <slepc/private/slepcimpl.h> | ||
| 16 | #include <slepcblaslapack.h> | ||
| 17 | |||
| 18 | /* | ||
| 19 | EPSDelayedArnoldi - This function is equivalent to BVMatArnoldi but | ||
| 20 | performs the computation in a different way. The main idea is that | ||
| 21 | reorthogonalization is delayed to the next Arnoldi step. This version is | ||
| 22 | more scalable but in some cases convergence may stagnate. | ||
| 23 | */ | ||
| 24 | 430 | PetscErrorCode EPSDelayedArnoldi(EPS eps,PetscScalar *H,PetscInt ldh,PetscInt k,PetscInt *M,PetscReal *beta,PetscBool *breakdown) | |
| 25 | { | ||
| 26 | 430 | PetscInt i,j,m=*M; | |
| 27 | 430 | Vec u,t; | |
| 28 | 430 | PetscScalar shh[100],*lhh,dot,dot2; | |
| 29 | 430 | PetscReal norm1=0.0,norm2=1.0; | |
| 30 | 430 | Vec vj,vj1,vj2=NULL; | |
| 31 | |||
| 32 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
430 | PetscFunctionBegin; |
| 33 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
430 | if (m<=100) lhh = shh; |
| 34 | ✗ | else PetscCall(PetscMalloc1(m,&lhh)); | |
| 35 |
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.
|
430 | PetscCall(BVCreateVec(eps->V,&u)); |
| 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.
|
430 | PetscCall(BVCreateVec(eps->V,&t)); |
| 37 | |||
| 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.
|
430 | PetscCall(BVSetActiveColumns(eps->V,0,m)); |
| 39 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7885 | for (j=k;j<m;j++) { |
| 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.
|
7455 | PetscCall(BVGetColumn(eps->V,j,&vj)); |
| 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.
|
7455 | PetscCall(BVGetColumn(eps->V,j+1,&vj1)); |
| 42 |
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.
|
7455 | PetscCall(STApply(eps->st,vj,vj1)); |
| 43 |
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.
|
7455 | PetscCall(BVRestoreColumn(eps->V,j,&vj)); |
| 44 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7455 | PetscCall(BVRestoreColumn(eps->V,j+1,&vj1)); |
| 45 | |||
| 46 |
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.
|
7455 | PetscCall(BVDotColumnBegin(eps->V,j+1,H+ldh*j)); |
| 47 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7455 | if (j>k) { |
| 48 |
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.
|
7025 | PetscCall(BVDotColumnBegin(eps->V,j,lhh)); |
| 49 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7025 | PetscCall(BVGetColumn(eps->V,j,&vj)); |
| 50 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7025 | PetscCall(VecDotBegin(vj,vj,&dot)); |
| 51 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7025 | if (j>k+1) { |
| 52 |
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.
|
6595 | PetscCall(BVNormVecBegin(eps->V,u,NORM_2,&norm2)); |
| 53 |
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.
|
6595 | PetscCall(BVGetColumn(eps->V,j-2,&vj2)); |
| 54 |
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.
|
6595 | PetscCall(VecDotBegin(u,vj2,&dot2)); |
| 55 | } | ||
| 56 |
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.
|
7025 | PetscCall(BVDotColumnEnd(eps->V,j+1,H+ldh*j)); |
| 57 |
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.
|
7025 | PetscCall(BVDotColumnEnd(eps->V,j,lhh)); |
| 58 |
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.
|
7025 | PetscCall(VecDotEnd(vj,vj,&dot)); |
| 59 |
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.
|
7025 | PetscCall(BVRestoreColumn(eps->V,j,&vj)); |
| 60 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7025 | if (j>k+1) { |
| 61 |
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.
|
6595 | PetscCall(BVNormVecEnd(eps->V,u,NORM_2,&norm2)); |
| 62 |
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.
|
6595 | PetscCall(VecDotEnd(u,vj2,&dot2)); |
| 63 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6595 | PetscCall(BVRestoreColumn(eps->V,j-2,&vj2)); |
| 64 | } | ||
| 65 | 7025 | norm1 = PetscSqrtReal(PetscRealPart(dot)); | |
| 66 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
76470 | for (i=0;i<j;i++) H[ldh*j+i] = H[ldh*j+i]/norm1; |
| 67 | 7025 | H[ldh*j+j] = H[ldh*j+j]/dot; | |
| 68 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7025 | PetscCall(BVCopyVec(eps->V,j,t)); |
| 69 |
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.
|
7025 | PetscCall(BVScaleColumn(eps->V,j,1.0/norm1)); |
| 70 |
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.
|
7025 | PetscCall(BVScaleColumn(eps->V,j+1,1.0/norm1)); |
| 71 |
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.
|
430 | } else PetscCall(BVDotColumnEnd(eps->V,j+1,H+ldh*j)); /* j==k */ |
| 72 | |||
| 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.
|
7455 | PetscCall(BVMultColumn(eps->V,-1.0,1.0,j+1,H+ldh*j)); |
| 74 | |||
| 75 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7455 | if (j>k) { |
| 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.
|
7025 | PetscCall(BVSetActiveColumns(eps->V,0,j)); |
| 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.
|
7025 | PetscCall(BVMultVec(eps->V,-1.0,1.0,t,lhh)); |
| 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.
|
7025 | PetscCall(BVSetActiveColumns(eps->V,0,m)); |
| 79 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
76470 | for (i=0;i<j;i++) H[ldh*(j-1)+i] += lhh[i]; |
| 80 | } | ||
| 81 | |||
| 82 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7455 | if (j>k+1) { |
| 83 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6595 | PetscCall(BVGetColumn(eps->V,j-1,&vj1)); |
| 84 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6595 | PetscCall(VecCopy(u,vj1)); |
| 85 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6595 | PetscCall(BVRestoreColumn(eps->V,j-1,&vj1)); |
| 86 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6595 | PetscCall(BVScaleColumn(eps->V,j-1,1.0/norm2)); |
| 87 | 6595 | H[ldh*(j-2)+j-1] = norm2; | |
| 88 | } | ||
| 89 | |||
| 90 |
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.
|
7455 | if (j<m-1) PetscCall(VecCopy(t,u)); |
| 91 | } | ||
| 92 | |||
| 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.
|
430 | PetscCall(BVNormVec(eps->V,t,NORM_2,&norm2)); |
| 94 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
430 | PetscCall(VecScale(t,1.0/norm2)); |
| 95 |
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.
|
430 | PetscCall(BVGetColumn(eps->V,m-1,&vj1)); |
| 96 |
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.
|
430 | PetscCall(VecCopy(t,vj1)); |
| 97 |
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.
|
430 | PetscCall(BVRestoreColumn(eps->V,m-1,&vj1)); |
| 98 | 430 | H[ldh*(m-2)+m-1] = norm2; | |
| 99 | |||
| 100 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
430 | PetscCall(BVDotColumn(eps->V,m,lhh)); |
| 101 | |||
| 102 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
430 | PetscCall(BVMultColumn(eps->V,-1.0,1.0,m,lhh)); |
| 103 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
8420 | for (i=0;i<m;i++) |
| 104 | 7990 | H[ldh*(m-1)+i] += lhh[i]; | |
| 105 | |||
| 106 |
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.
|
430 | PetscCall(BVNormColumn(eps->V,m,NORM_2,beta)); |
| 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.
|
430 | PetscCall(BVScaleColumn(eps->V,m,1.0 / *beta)); |
| 108 | 430 | *breakdown = PETSC_FALSE; | |
| 109 | |||
| 110 |
1/10✗ 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.
|
430 | if (m>100) PetscCall(PetscFree(lhh)); |
| 111 |
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.
|
430 | PetscCall(VecDestroy(&u)); |
| 112 |
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.
|
430 | PetscCall(VecDestroy(&t)); |
| 113 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
86 | PetscFunctionReturn(PETSC_SUCCESS); |
| 114 | } | ||
| 115 | |||
| 116 | /* | ||
| 117 | EPSDelayedArnoldi1 - This function is similar to EPSDelayedArnoldi, | ||
| 118 | but without reorthogonalization (only delayed normalization). | ||
| 119 | */ | ||
| 120 | 180 | PetscErrorCode EPSDelayedArnoldi1(EPS eps,PetscScalar *H,PetscInt ldh,PetscInt k,PetscInt *M,PetscReal *beta,PetscBool *breakdown) | |
| 121 | { | ||
| 122 | 180 | PetscInt i,j,m=*M; | |
| 123 | 180 | PetscScalar dot; | |
| 124 | 180 | PetscReal norm=0.0; | |
| 125 | 180 | Vec vj,vj1; | |
| 126 | |||
| 127 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
180 | PetscFunctionBegin; |
| 128 |
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.
|
180 | PetscCall(BVSetActiveColumns(eps->V,0,m)); |
| 129 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
3260 | for (j=k;j<m;j++) { |
| 130 |
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.
|
3080 | PetscCall(BVGetColumn(eps->V,j,&vj)); |
| 131 |
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.
|
3080 | PetscCall(BVGetColumn(eps->V,j+1,&vj1)); |
| 132 |
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.
|
3080 | PetscCall(STApply(eps->st,vj,vj1)); |
| 133 |
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.
|
3080 | PetscCall(BVRestoreColumn(eps->V,j+1,&vj1)); |
| 134 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
3080 | if (j>k) { |
| 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.
|
2900 | PetscCall(BVDotColumnBegin(eps->V,j+1,H+ldh*j)); |
| 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.
|
2900 | PetscCall(VecDotBegin(vj,vj,&dot)); |
| 137 |
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.
|
2900 | PetscCall(BVDotColumnEnd(eps->V,j+1,H+ldh*j)); |
| 138 |
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.
|
2900 | PetscCall(VecDotEnd(vj,vj,&dot)); |
| 139 | 2900 | norm = PetscSqrtReal(PetscRealPart(dot)); | |
| 140 |
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.
|
2900 | PetscCall(BVScaleColumn(eps->V,j,1.0/norm)); |
| 141 | 2900 | H[ldh*(j-1)+j] = norm; | |
| 142 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
30220 | for (i=0;i<j;i++) H[ldh*j+i] = H[ldh*j+i]/norm; |
| 143 | 2900 | H[ldh*j+j] = H[ldh*j+j]/dot; | |
| 144 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2900 | PetscCall(BVScaleColumn(eps->V,j+1,1.0/norm)); |
| 145 | 2900 | *beta = norm; | |
| 146 | } else { /* j==k */ | ||
| 147 |
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.
|
180 | PetscCall(BVDotColumn(eps->V,j+1,H+ldh*j)); |
| 148 | } | ||
| 149 |
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.
|
3080 | PetscCall(BVRestoreColumn(eps->V,j,&vj)); |
| 150 |
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.
|
3080 | PetscCall(BVMultColumn(eps->V,-1.0,1.0,j+1,H+ldh*j)); |
| 151 | } | ||
| 152 | |||
| 153 | 180 | *breakdown = PETSC_FALSE; | |
| 154 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
180 | PetscFunctionReturn(PETSC_SUCCESS); |
| 155 | } | ||
| 156 | |||
| 157 | /* | ||
| 158 | EPSKrylovConvergence - Implements the loop that checks for convergence | ||
| 159 | in Krylov methods. | ||
| 160 | |||
| 161 | Input Parameters: | ||
| 162 | eps - the eigensolver; some error estimates are updated in eps->errest | ||
| 163 | getall - whether all residuals must be computed | ||
| 164 | kini - initial value of k (the loop variable) | ||
| 165 | nits - number of iterations of the loop | ||
| 166 | V - set of basis vectors (used only if trueresidual is activated) | ||
| 167 | nv - number of vectors to process (dimension of Q, columns of V) | ||
| 168 | beta - norm of f (the residual vector of the Arnoldi/Lanczos factorization) | ||
| 169 | corrf - correction factor for residual estimates (only in harmonic KS) | ||
| 170 | |||
| 171 | Output Parameters: | ||
| 172 | kout - the first index where the convergence test failed | ||
| 173 | */ | ||
| 174 | 56602 | PetscErrorCode EPSKrylovConvergence(EPS eps,PetscBool getall,PetscInt kini,PetscInt nits,PetscReal beta,PetscReal betat,PetscReal corrf,PetscInt *kout) | |
| 175 | { | ||
| 176 | 56602 | PetscInt k,newk,newk2,marker,ld,inside; | |
| 177 | 56602 | PetscScalar re,im,*Zr,*Zi,*X; | |
| 178 | 56602 | PetscReal resnorm,lerrest; | |
| 179 | 56602 | PetscBool isshift,refined,istrivial; | |
| 180 | 56602 | Vec x=NULL,y=NULL,w[3]; | |
| 181 | |||
| 182 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
56602 | PetscFunctionBegin; |
| 183 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
56602 | PetscCall(RGIsTrivial(eps->rg,&istrivial)); |
| 184 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
56602 | if (PetscUnlikely(eps->trueres)) { |
| 185 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
420 | PetscCall(BVCreateVec(eps->V,&x)); |
| 186 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
420 | PetscCall(BVCreateVec(eps->V,&y)); |
| 187 |
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.
|
420 | PetscCall(BVCreateVec(eps->V,&w[0])); |
| 188 |
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.
|
420 | PetscCall(BVCreateVec(eps->V,&w[2])); |
| 189 | #if !defined(PETSC_USE_COMPLEX) | ||
| 190 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
240 | PetscCall(BVCreateVec(eps->V,&w[1])); |
| 191 | #else | ||
| 192 | 180 | w[1] = NULL; | |
| 193 | #endif | ||
| 194 | } | ||
| 195 |
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.
|
56602 | PetscCall(DSGetLeadingDimension(eps->ds,&ld)); |
| 196 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
56602 | PetscCall(DSGetRefined(eps->ds,&refined)); |
| 197 |
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.
|
56602 | PetscCall(PetscObjectTypeCompare((PetscObject)eps->st,STSHIFT,&isshift)); |
| 198 | 56602 | marker = -1; | |
| 199 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
56602 | if (eps->trackall) getall = PETSC_TRUE; |
| 200 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
151841 | for (k=kini;k<kini+nits;k++) { |
| 201 | /* eigenvalue */ | ||
| 202 | 150268 | re = eps->eigr[k]; | |
| 203 | 150268 | im = eps->eigi[k]; | |
| 204 |
12/14✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 10 times.
✓ Branch 5 taken 10 times.
✓ Branch 6 taken 10 times.
✓ Branch 7 taken 10 times.
✓ Branch 8 taken 2 times.
✓ Branch 9 taken 8 times.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
150268 | if (!istrivial || eps->trueres || isshift || eps->conv==EPS_CONV_NORM) PetscCall(STBackTransform(eps->st,1,&re,&im)); |
| 205 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
150268 | if (PetscUnlikely(!istrivial)) { |
| 206 |
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.
|
1150 | PetscCall(RGCheckInside(eps->rg,1,&re,&im,&inside)); |
| 207 |
3/4✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
1150 | if (marker==-1 && inside<0) marker = k; |
| 208 |
4/6✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 10 times.
✗ Branch 5 not taken.
|
1150 | if (!(eps->trueres || isshift || eps->conv==EPS_CONV_NORM)) { /* make sure eps->converged below uses the right value */ |
| 209 | 425 | re = eps->eigr[k]; | |
| 210 | 425 | im = eps->eigi[k]; | |
| 211 | } | ||
| 212 | } | ||
| 213 | 150268 | newk = k; | |
| 214 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
150268 | PetscCall(DSVectors(eps->ds,DS_MAT_X,&newk,&resnorm)); |
| 215 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
150268 | if (PetscUnlikely(eps->trueres)) { |
| 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.
|
655 | PetscCall(DSGetArray(eps->ds,DS_MAT_X,&X)); |
| 217 | 655 | Zr = X+k*ld; | |
| 218 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 10 times.
|
655 | if (newk==k+1) Zi = X+newk*ld; |
| 219 | else Zi = NULL; | ||
| 220 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
655 | PetscCall(EPSComputeRitzVector(eps,Zr,Zi,eps->V,x,y)); |
| 221 |
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.
|
655 | PetscCall(DSRestoreArray(eps->ds,DS_MAT_X,&X)); |
| 222 |
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.
|
655 | PetscCall(EPSComputeResidualNorm_Private(eps,PETSC_FALSE,re,im,x,y,w,&resnorm)); |
| 223 | } | ||
| 224 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
149613 | else if (!refined) resnorm *= beta*corrf; |
| 225 | /* error estimate */ | ||
| 226 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
150268 | PetscCall((*eps->converged)(eps,re,im,resnorm,&eps->errest[k],eps->convergedctx)); |
| 227 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
150268 | if (marker==-1 && eps->errest[k] >= eps->tol) marker = k; |
| 228 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
150268 | if (PetscUnlikely(eps->twosided)) { |
| 229 | 1942 | newk2 = k; | |
| 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.
|
1942 | PetscCall(DSVectors(eps->ds,DS_MAT_Y,&newk2,&resnorm)); |
| 231 | 1942 | resnorm *= betat; | |
| 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.
|
1942 | PetscCall((*eps->converged)(eps,re,im,resnorm,&lerrest,eps->convergedctx)); |
| 233 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1942 | eps->errest[k] = PetscMax(eps->errest[k],lerrest); |
| 234 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
1942 | if (marker==-1 && lerrest >= eps->tol) marker = k; |
| 235 | } | ||
| 236 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
150268 | if (PetscUnlikely(newk==k+1)) { |
| 237 | 2981 | eps->errest[k+1] = eps->errest[k]; | |
| 238 | 2981 | k++; | |
| 239 | } | ||
| 240 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
150268 | if (marker!=-1 && !getall) break; |
| 241 | } | ||
| 242 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
56602 | if (marker!=-1) k = marker; |
| 243 | 56602 | *kout = k; | |
| 244 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
56602 | if (PetscUnlikely(eps->trueres)) { |
| 245 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
420 | PetscCall(VecDestroy(&x)); |
| 246 |
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.
|
420 | PetscCall(VecDestroy(&y)); |
| 247 |
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.
|
420 | PetscCall(VecDestroy(&w[0])); |
| 248 |
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.
|
420 | PetscCall(VecDestroy(&w[2])); |
| 249 | #if !defined(PETSC_USE_COMPLEX) | ||
| 250 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
240 | PetscCall(VecDestroy(&w[1])); |
| 251 | #endif | ||
| 252 | } | ||
| 253 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
11090 | PetscFunctionReturn(PETSC_SUCCESS); |
| 254 | } | ||
| 255 | |||
| 256 | 7064 | PetscErrorCode EPSPseudoLanczos(EPS eps,PetscReal *alpha,PetscReal *beta,PetscReal *omega,PetscInt k,PetscInt *M,PetscBool *breakdown,PetscBool *symmlost,PetscReal *cos,Vec w) | |
| 257 | { | ||
| 258 | 7064 | PetscInt j,m = *M,i,ld,l; | |
| 259 | 7064 | Vec vj,vj1; | |
| 260 | 7064 | PetscScalar *hwork,lhwork[100]; | |
| 261 | 7064 | PetscReal norm,norm1,norm2,t,sym=0.0,fro=0.0; | |
| 262 | 7064 | PetscBLASInt j_,one=1; | |
| 263 | |||
| 264 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
7064 | PetscFunctionBegin; |
| 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.
|
7064 | PetscCall(DSGetLeadingDimension(eps->ds,&ld)); |
| 266 |
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.
|
7064 | PetscCall(DSGetDimensions(eps->ds,NULL,&l,NULL,NULL)); |
| 267 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
7064 | if (cos) *cos = 1.0; |
| 268 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
7064 | if (m > 100) PetscCall(PetscMalloc1(m,&hwork)); |
| 269 | 7016 | else hwork = lhwork; | |
| 270 | |||
| 271 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7064 | PetscCall(BVSetActiveColumns(eps->V,0,m)); |
| 272 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
87676 | for (j=k;j<m;j++) { |
| 273 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
80612 | PetscCall(BVGetColumn(eps->V,j,&vj)); |
| 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.
|
80612 | PetscCall(BVGetColumn(eps->V,j+1,&vj1)); |
| 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.
|
80612 | PetscCall(STApply(eps->st,vj,vj1)); |
| 276 |
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.
|
80612 | PetscCall(BVRestoreColumn(eps->V,j,&vj)); |
| 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.
|
80612 | PetscCall(BVRestoreColumn(eps->V,j+1,&vj1)); |
| 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.
|
80612 | PetscCall(BVOrthogonalizeColumn(eps->V,j+1,hwork,&norm,breakdown)); |
| 279 | 80612 | alpha[j] = PetscRealPart(hwork[j]); | |
| 280 | 80612 | beta[j] = PetscAbsReal(norm); | |
| 281 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
80612 | if (j==k) { |
| 282 | 7064 | PetscReal *f; | |
| 283 | |||
| 284 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7064 | PetscCall(DSGetArrayReal(eps->ds,DS_MAT_T,&f)); |
| 285 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11109 | for (i=0;i<l;i++) hwork[i] = 0.0; |
| 286 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
176081 | for (;i<j-1;i++) hwork[i] -= f[2*ld+i]; |
| 287 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7064 | PetscCall(DSRestoreArrayReal(eps->ds,DS_MAT_T,&f)); |
| 288 | } | ||
| 289 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
80612 | if (j>0) { |
| 290 | 80491 | hwork[j-1] -= beta[j-1]; | |
| 291 |
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.
|
80491 | PetscCall(PetscBLASIntCast(j,&j_)); |
| 292 | 80491 | sym = SlepcAbs(BLASnrm2_(&j_,hwork,&one),sym); | |
| 293 | } | ||
| 294 | 80612 | fro = SlepcAbs(fro,SlepcAbs(alpha[j],beta[j])); | |
| 295 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
80612 | if (j>0) fro = SlepcAbs(fro,beta[j-1]); |
| 296 |
3/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
|
81507 | if (sym/fro>PetscMax(PETSC_SQRT_MACHINE_EPSILON,10*eps->tol)) { *symmlost = PETSC_TRUE; *M=j; break; } |
| 297 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
80612 | omega[j+1] = (norm<0.0)? -1.0: 1.0; |
| 298 |
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.
|
80612 | PetscCall(BVScaleColumn(eps->V,j+1,1.0/norm)); |
| 299 | /* */ | ||
| 300 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
80612 | if (cos) { |
| 301 | ✗ | PetscCall(BVGetColumn(eps->V,j+1,&vj1)); | |
| 302 | ✗ | PetscCall(VecNorm(vj1,NORM_2,&norm1)); | |
| 303 | ✗ | PetscCall(BVApplyMatrix(eps->V,vj1,w)); | |
| 304 | ✗ | PetscCall(BVRestoreColumn(eps->V,j+1,&vj1)); | |
| 305 | ✗ | PetscCall(VecNorm(w,NORM_2,&norm2)); | |
| 306 | ✗ | t = 1.0/(norm1*norm2); | |
| 307 | ✗ | if (*cos>t) *cos = t; | |
| 308 | } | ||
| 309 | } | ||
| 310 |
7/10✓ Branch 0 taken 6 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 4 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
7064 | if (m > 100) PetscCall(PetscFree(hwork)); |
| 311 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
1450 | PetscFunctionReturn(PETSC_SUCCESS); |
| 312 | } | ||
| 313 |