| 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 | EPS routines related to the solution process | ||
| 12 | */ | ||
| 13 | |||
| 14 | #include <slepc/private/epsimpl.h> /*I "slepceps.h" I*/ | ||
| 15 | #include <slepc/private/bvimpl.h> | ||
| 16 | #include <petscdraw.h> | ||
| 17 | |||
| 18 | 101238 | PetscErrorCode EPSComputeVectors(EPS eps) | |
| 19 | { | ||
| 20 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
101238 | PetscFunctionBegin; |
| 21 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
101238 | EPSCheckSolved(eps,1); |
| 22 |
8/10✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
101238 | if (eps->state==EPS_STATE_SOLVED) PetscTryTypeMethod(eps,computevectors); |
| 23 | 101238 | eps->state = EPS_STATE_EIGENVECTORS; | |
| 24 |
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.
|
101238 | PetscFunctionReturn(PETSC_SUCCESS); |
| 25 | } | ||
| 26 | |||
| 27 | 11021 | static PetscErrorCode EPSComputeValues(EPS eps) | |
| 28 | { | ||
| 29 | 11021 | PetscBool injective,iscomp,isfilter; | |
| 30 | 11021 | PetscInt i,n,aux,nconv0; | |
| 31 | 11021 | Mat A,B=NULL,G,Z; | |
| 32 | |||
| 33 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
11021 | PetscFunctionBegin; |
| 34 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11021 | switch (eps->categ) { |
| 35 | 8901 | case EPS_CATEGORY_KRYLOV: | |
| 36 | case EPS_CATEGORY_OTHER: | ||
| 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.
|
8901 | PetscCall(STIsInjective(eps->st,&injective)); |
| 38 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
8901 | if (injective) { |
| 39 | /* one-to-one mapping: backtransform eigenvalues */ | ||
| 40 |
5/10✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
8685 | PetscUseTypeMethod(eps,backtransform); |
| 41 | } else { | ||
| 42 | /* compute eigenvalues from Rayleigh quotient */ | ||
| 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.
|
216 | PetscCall(DSGetDimensions(eps->ds,&n,NULL,NULL,NULL)); |
| 44 |
2/2✓ Branch 0 taken 9 times.
✓ Branch 1 taken 1 times.
|
216 | if (!n) break; |
| 45 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
216 | PetscCall(EPSGetOperators(eps,&A,&B)); |
| 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.
|
216 | PetscCall(BVSetActiveColumns(eps->V,0,n)); |
| 47 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
216 | PetscCall(DSGetCompact(eps->ds,&iscomp)); |
| 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.
|
216 | PetscCall(DSSetCompact(eps->ds,PETSC_FALSE)); |
| 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.
|
216 | PetscCall(DSGetMat(eps->ds,DS_MAT_A,&G)); |
| 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.
|
216 | PetscCall(BVMatProject(eps->V,A,eps->V,G)); |
| 51 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
216 | PetscCall(DSRestoreMat(eps->ds,DS_MAT_A,&G)); |
| 52 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
216 | if (B) { |
| 53 | ✗ | PetscCall(DSGetMat(eps->ds,DS_MAT_B,&G)); | |
| 54 | ✗ | PetscCall(BVMatProject(eps->V,B,eps->V,G)); | |
| 55 | ✗ | PetscCall(DSRestoreMat(eps->ds,DS_MAT_B,&G)); | |
| 56 | } | ||
| 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.
|
216 | PetscCall(DSSolve(eps->ds,eps->eigr,eps->eigi)); |
| 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.
|
216 | PetscCall(DSSort(eps->ds,eps->eigr,eps->eigi,NULL,NULL,NULL)); |
| 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.
|
216 | PetscCall(DSSynchronize(eps->ds,eps->eigr,eps->eigi)); |
| 60 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
216 | PetscCall(DSSetCompact(eps->ds,iscomp)); |
| 61 |
2/6✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
|
216 | if (eps->ishermitian && (!eps->isgeneralized || eps->ispositive)) { /* V = V * Z */ |
| 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.
|
216 | PetscCall(DSVectors(eps->ds,DS_MAT_X,NULL,NULL)); |
| 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.
|
216 | PetscCall(DSGetMat(eps->ds,DS_MAT_X,&Z)); |
| 64 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
216 | PetscCall(BVMultInPlace(eps->V,Z,0,n)); |
| 65 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
216 | PetscCall(DSRestoreMat(eps->ds,DS_MAT_X,&Z)); |
| 66 | } | ||
| 67 | /* in case of STFILTER discard computed eigenvalues that lie outside the wanted interval */ | ||
| 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.
|
216 | PetscCall(PetscObjectTypeCompare((PetscObject)eps->st,STFILTER,&isfilter)); |
| 69 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
216 | if (isfilter) { |
| 70 | 206 | nconv0 = eps->nconv; | |
| 71 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
2308 | for (i=0;i<eps->nconv;i++) { |
| 72 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
2102 | if (PetscRealPart(eps->eigr[eps->perm[i]])<eps->inta || PetscRealPart(eps->eigr[eps->perm[i]])>eps->intb) { |
| 73 | 412 | eps->nconv--; | |
| 74 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
412 | if (i<eps->nconv) { SlepcSwap(eps->perm[i],eps->perm[eps->nconv],aux); i--; } |
| 75 | } | ||
| 76 | } | ||
| 77 |
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.
|
206 | if (nconv0>eps->nconv) PetscCall(PetscInfo(eps,"Discarded %" PetscInt_FMT " computed eigenvalues lying outside the interval\n",nconv0-eps->nconv)); |
| 78 | } | ||
| 79 | } | ||
| 80 | break; | ||
| 81 | case EPS_CATEGORY_PRECOND: | ||
| 82 | case EPS_CATEGORY_CONTOUR: | ||
| 83 | /* eigenvalues already available as an output of the solver */ | ||
| 84 | break; | ||
| 85 | } | ||
| 86 |
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.
|
2125 | PetscFunctionReturn(PETSC_SUCCESS); |
| 87 | } | ||
| 88 | |||
| 89 | /*@ | ||
| 90 | EPSSolve - Solves the eigensystem. | ||
| 91 | |||
| 92 | Collective | ||
| 93 | |||
| 94 | Input Parameter: | ||
| 95 | . eps - the linear eigensolver context | ||
| 96 | |||
| 97 | Options Database Keys: | ||
| 98 | + -eps_view - print information about the solver once the solve is complete | ||
| 99 | . -eps_view_pre - print information about the solver before the solve starts | ||
| 100 | . -eps_view_mat0 - view the first matrix ($A$) | ||
| 101 | . -eps_view_mat1 - view the second matrix ($B$) | ||
| 102 | . -eps_view_vectors - view the computed eigenvectors | ||
| 103 | . -eps_view_values - view the computed eigenvalues | ||
| 104 | . -eps_converged_reason - print reason for convergence/divergence, and number of iterations | ||
| 105 | . -eps_error_absolute - print absolute errors of each eigenpair | ||
| 106 | . -eps_error_relative - print relative errors of each eigenpair | ||
| 107 | - -eps_error_backward - print backward errors of each eigenpair | ||
| 108 | |||
| 109 | Notes: | ||
| 110 | The problem matrices are specified with `EPSSetOperators()`. | ||
| 111 | |||
| 112 | `EPSSolve()` will return without generating an error regardless of whether | ||
| 113 | all requested solutions were computed or not. Call `EPSGetConverged()` to get the | ||
| 114 | actual number of computed solutions, and `EPSGetConvergedReason()` to determine if | ||
| 115 | the solver converged or failed and why. | ||
| 116 | |||
| 117 | All the command-line options listed above admit an optional argument specifying | ||
| 118 | the viewer type and options. For instance, use `-eps_view_mat0 binary:amatrix.bin` | ||
| 119 | to save the $A$ matrix to a binary file, `-eps_view_values draw` to draw the computed | ||
| 120 | eigenvalues graphically, or `-eps_error_relative :myerr.m:ascii_matlab` to save | ||
| 121 | the errors in a file that can be executed in Matlab. | ||
| 122 | See `PetscObjectViewFromOptions()` for more details. | ||
| 123 | |||
| 124 | Level: beginner | ||
| 125 | |||
| 126 | .seealso: [](ch:eps), `EPSCreate()`, `EPSSetUp()`, `EPSDestroy()`, `EPSSetOperators()`, `EPSGetConverged()`, `EPSGetConvergedReason()` | ||
| 127 | @*/ | ||
| 128 | 11021 | PetscErrorCode EPSSolve(EPS eps) | |
| 129 | { | ||
| 130 | 11021 | PetscInt i; | |
| 131 | 11021 | PetscBool hasname; | |
| 132 | 11021 | STMatMode matmode; | |
| 133 | 11021 | Mat A,B; | |
| 134 | |||
| 135 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
11021 | PetscFunctionBegin; |
| 136 |
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.
|
11021 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 137 |
2/14✓ Branch 0 taken 8 times.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
|
11021 | if (eps->state>=EPS_STATE_SOLVED) PetscFunctionReturn(PETSC_SUCCESS); |
| 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.
|
11021 | PetscCall(PetscLogEventBegin(EPS_Solve,eps,0,0,0)); |
| 139 | |||
| 140 | /* Call setup */ | ||
| 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.
|
11021 | PetscCall(EPSSetUp(eps)); |
| 142 | |||
| 143 | /* Safeguard for matrices of size 0 */ | ||
| 144 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
11021 | if (eps->n == 0) { |
| 145 | ✗ | eps->nconv = 0; | |
| 146 | ✗ | eps->reason = EPS_CONVERGED_TOL; | |
| 147 | ✗ | eps->state = EPS_STATE_SOLVED; | |
| 148 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 149 | } | ||
| 150 | |||
| 151 | 11021 | eps->nconv = 0; | |
| 152 | 11021 | eps->its = 0; | |
| 153 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
311575 | for (i=0;i<eps->ncv;i++) { |
| 154 | 300554 | eps->eigr[i] = 0.0; | |
| 155 | 300554 | eps->eigi[i] = 0.0; | |
| 156 | 300554 | eps->errest[i] = 0.0; | |
| 157 | 300554 | eps->perm[i] = i; | |
| 158 | } | ||
| 159 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(EPSViewFromOptions(eps,NULL,"-eps_view_pre")); |
| 160 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(RGViewFromOptions(eps->rg,NULL,"-rg_view")); |
| 161 | |||
| 162 | /* Call solver */ | ||
| 163 |
5/10✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
11021 | PetscUseTypeMethod(eps,solve); |
| 164 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
11021 | PetscCheck(eps->reason,PetscObjectComm((PetscObject)eps),PETSC_ERR_PLIB,"Internal error, solver returned without setting converged reason"); |
| 165 | 11021 | eps->state = EPS_STATE_SOLVED; | |
| 166 | |||
| 167 | /* Only the first nconv columns contain useful information (except in CISS) */ | ||
| 168 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(BVSetActiveColumns(eps->V,0,eps->nconv)); |
| 169 |
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.
|
11021 | if (eps->twosided) PetscCall(BVSetActiveColumns(eps->W,0,eps->nconv)); |
| 170 | |||
| 171 | /* If inplace, purify eigenvectors before reverting operator */ | ||
| 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.
|
11021 | PetscCall(STGetMatMode(eps->st,&matmode)); |
| 173 |
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.
|
11021 | if (matmode == ST_MATMODE_INPLACE && eps->ispositive) PetscCall(EPSComputeVectors(eps)); |
| 174 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(STPostSolve(eps->st)); |
| 175 | |||
| 176 | /* Map eigenvalues back to the original problem if appropriate */ | ||
| 177 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(EPSComputeValues(eps)); |
| 178 | |||
| 179 | #if !PetscDefined(USE_COMPLEX) | ||
| 180 | /* Reorder conjugate eigenvalues (positive imaginary first) */ | ||
| 181 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
59332 | for (i=0;i<eps->nconv-1;i++) { |
| 182 |
6/6✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 5 times.
✓ Branch 3 taken 5 times.
✓ Branch 4 taken 5 times.
✓ Branch 5 taken 5 times.
|
53701 | if (eps->eigi[i] != 0 && (eps->problem_type!=EPS_HAMILT || eps->eigr[i]!=0)) { |
| 183 | /* conjugate eigenvalues */ | ||
| 184 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
925 | if (eps->eigi[i] < 0) { |
| 185 | 116 | eps->eigi[i] = -eps->eigi[i]; | |
| 186 | 116 | eps->eigi[i+1] = -eps->eigi[i+1]; | |
| 187 | /* the next correction only works with eigenvectors */ | ||
| 188 |
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.
|
116 | PetscCall(EPSComputeVectors(eps)); |
| 189 |
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.
|
116 | PetscCall(BVScaleColumn(eps->V,i+1,-1.0)); |
| 190 |
6/8✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
116 | if (eps->W) PetscCall(BVScaleColumn(eps->W,i+1,-1.0)); |
| 191 | } | ||
| 192 | 925 | i++; | |
| 193 | } | ||
| 194 | } | ||
| 195 | #endif | ||
| 196 | |||
| 197 | /* Sort eigenvalues according to eps->which parameter */ | ||
| 198 |
6/8✓ Branch 0 taken 5 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
11021 | if (eps->problem_type==EPS_HAMILT) PetscCall(SlepcSortEigenvaluesSpecial(eps->sc,eps->nconv,eps->eigr,eps->eigi,eps->perm)); |
| 199 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10986 | else PetscCall(SlepcSortEigenvalues(eps->sc,eps->nconv,eps->eigr,eps->eigi,eps->perm)); |
| 200 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(PetscLogEventEnd(EPS_Solve,eps,0,0,0)); |
| 201 | |||
| 202 | /* Various viewers */ | ||
| 203 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(EPSViewFromOptions(eps,NULL,"-eps_view")); |
| 204 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(EPSConvergedReasonViewFromOptions(eps)); |
| 205 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(EPSErrorViewFromOptions(eps)); |
| 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.
|
11021 | PetscCall(EPSValuesViewFromOptions(eps)); |
| 207 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(EPSVectorsViewFromOptions(eps)); |
| 208 | |||
| 209 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11021 | PetscCall(PetscOptionsHasName(((PetscObject)eps)->options,((PetscObject)eps)->prefix,"-eps_view_mat0",&hasname)); |
| 210 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
11021 | if (hasname) { |
| 211 | ✗ | PetscCall(EPSGetOperators(eps,&A,NULL)); | |
| 212 | ✗ | PetscCall(MatViewFromOptions(A,(PetscObject)eps,"-eps_view_mat0")); | |
| 213 | } | ||
| 214 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11021 | if (eps->isgeneralized) { |
| 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.
|
2080 | PetscCall(PetscOptionsHasName(((PetscObject)eps)->options,((PetscObject)eps)->prefix,"-eps_view_mat1",&hasname)); |
| 216 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
2080 | if (hasname) { |
| 217 | ✗ | PetscCall(EPSGetOperators(eps,NULL,&B)); | |
| 218 | ✗ | PetscCall(MatViewFromOptions(B,(PetscObject)eps,"-eps_view_mat1")); | |
| 219 | } | ||
| 220 | } | ||
| 221 | |||
| 222 | /* Remove deflation and initial subspaces */ | ||
| 223 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11021 | if (eps->nds) { |
| 224 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
153 | PetscCall(BVSetNumConstraints(eps->V,0)); |
| 225 | 153 | eps->nds = 0; | |
| 226 | } | ||
| 227 | 11021 | eps->nini = 0; | |
| 228 |
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.
|
11021 | PetscFunctionReturn(PETSC_SUCCESS); |
| 229 | } | ||
| 230 | |||
| 231 | /*@ | ||
| 232 | EPSGetIterationNumber - Gets the current iteration number. If the | ||
| 233 | call to `EPSSolve()` is complete, then it returns the number of iterations | ||
| 234 | carried out by the solution method. | ||
| 235 | |||
| 236 | Not Collective | ||
| 237 | |||
| 238 | Input Parameter: | ||
| 239 | . eps - the linear eigensolver context | ||
| 240 | |||
| 241 | Output Parameter: | ||
| 242 | . its - number of iterations | ||
| 243 | |||
| 244 | Note: | ||
| 245 | During the $i$-th iteration this call returns $i-1$. If `EPSSolve()` is | ||
| 246 | complete, then parameter `its` contains either the iteration number at | ||
| 247 | which convergence was successfully reached, or failure was detected. | ||
| 248 | Call `EPSGetConvergedReason()` to determine if the solver converged or | ||
| 249 | failed and why. | ||
| 250 | |||
| 251 | Level: intermediate | ||
| 252 | |||
| 253 | .seealso: [](ch:eps), `EPSGetConvergedReason()`, `EPSSetTolerances()` | ||
| 254 | @*/ | ||
| 255 | 1564 | PetscErrorCode EPSGetIterationNumber(EPS eps,PetscInt *its) | |
| 256 | { | ||
| 257 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1564 | PetscFunctionBegin; |
| 258 |
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.
|
1564 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 259 |
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.
|
1564 | PetscAssertPointer(its,2); |
| 260 | 1564 | *its = eps->its; | |
| 261 |
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.
|
1564 | PetscFunctionReturn(PETSC_SUCCESS); |
| 262 | } | ||
| 263 | |||
| 264 | /*@ | ||
| 265 | EPSGetConverged - Gets the number of converged eigenpairs. | ||
| 266 | |||
| 267 | Not Collective | ||
| 268 | |||
| 269 | Input Parameter: | ||
| 270 | . eps - the linear eigensolver context | ||
| 271 | |||
| 272 | Output Parameter: | ||
| 273 | . nconv - number of converged eigenpairs | ||
| 274 | |||
| 275 | Notes: | ||
| 276 | This function should be called after `EPSSolve()` has finished. | ||
| 277 | |||
| 278 | The value `nconv` may be different from the number of requested solutions | ||
| 279 | `nev`, but not larger than `ncv`, see `EPSSetDimensions()`. | ||
| 280 | |||
| 281 | Level: beginner | ||
| 282 | |||
| 283 | .seealso: [](ch:eps), `EPSSetDimensions()`, `EPSSolve()`, `EPSGetEigenpair()` | ||
| 284 | @*/ | ||
| 285 | 7180 | PetscErrorCode EPSGetConverged(EPS eps,PetscInt *nconv) | |
| 286 | { | ||
| 287 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
7180 | PetscFunctionBegin; |
| 288 |
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.
|
7180 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 289 |
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.
|
7180 | PetscAssertPointer(nconv,2); |
| 290 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
7180 | EPSCheckSolved(eps,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.
|
7180 | PetscCall(EPS_GetActualConverged(eps,nconv)); |
| 292 |
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.
|
1383 | PetscFunctionReturn(PETSC_SUCCESS); |
| 293 | } | ||
| 294 | |||
| 295 | /*@ | ||
| 296 | EPSGetConvergedReason - Gets the reason why the `EPSSolve()` iteration was | ||
| 297 | stopped. | ||
| 298 | |||
| 299 | Not Collective | ||
| 300 | |||
| 301 | Input Parameter: | ||
| 302 | . eps - the linear eigensolver context | ||
| 303 | |||
| 304 | Output Parameter: | ||
| 305 | . reason - negative value indicates diverged, positive value converged, see | ||
| 306 | `EPSConvergedReason` for the possible values | ||
| 307 | |||
| 308 | Options Database Key: | ||
| 309 | . -eps_converged_reason - print reason for convergence/divergence, and number of iterations | ||
| 310 | |||
| 311 | Note: | ||
| 312 | If this routine is called before or doing the `EPSSolve()` the value of | ||
| 313 | `EPS_CONVERGED_ITERATING` is returned. | ||
| 314 | |||
| 315 | Level: intermediate | ||
| 316 | |||
| 317 | .seealso: [](ch:eps), `EPSSetTolerances()`, `EPSSolve()`, `EPSConvergedReason` | ||
| 318 | @*/ | ||
| 319 | 1513 | PetscErrorCode EPSGetConvergedReason(EPS eps,EPSConvergedReason *reason) | |
| 320 | { | ||
| 321 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1513 | PetscFunctionBegin; |
| 322 |
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.
|
1513 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 323 |
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.
|
1513 | PetscAssertPointer(reason,2); |
| 324 | 1513 | *reason = eps->reason; | |
| 325 |
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.
|
1513 | PetscFunctionReturn(PETSC_SUCCESS); |
| 326 | } | ||
| 327 | |||
| 328 | /*@ | ||
| 329 | EPSGetInvariantSubspace - Gets an orthonormal basis of the computed invariant | ||
| 330 | subspace. | ||
| 331 | |||
| 332 | Collective | ||
| 333 | |||
| 334 | Input Parameter: | ||
| 335 | . eps - the linear eigensolver context | ||
| 336 | |||
| 337 | Output Parameter: | ||
| 338 | . v - an array of vectors | ||
| 339 | |||
| 340 | Notes: | ||
| 341 | This function should be called after `EPSSolve()` has finished. | ||
| 342 | |||
| 343 | The user should provide in `v` an array of `nconv` vectors, where `nconv` | ||
| 344 | is the value returned by `EPSGetConverged()`. | ||
| 345 | |||
| 346 | The first $k$ vectors returned in `v` span an invariant subspace associated | ||
| 347 | with the first $k$ computed eigenvalues (note that this is not true if the | ||
| 348 | $k$-th eigenvalue is complex and matrix $A$ is real; in this case the first | ||
| 349 | $k+1$ vectors should be used). An invariant subspace $X$ of $A$ satisfies | ||
| 350 | $Ax\in X, \forall x\in X$ (a similar definition applies for generalized | ||
| 351 | eigenproblems). | ||
| 352 | |||
| 353 | Level: intermediate | ||
| 354 | |||
| 355 | .seealso: [](ch:eps), `EPSGetEigenpair()`, `EPSGetConverged()`, `EPSSolve()` | ||
| 356 | @*/ | ||
| 357 | 80 | PetscErrorCode EPSGetInvariantSubspace(EPS eps,Vec v[]) | |
| 358 | { | ||
| 359 | 80 | PetscInt i; | |
| 360 | 80 | BV V=eps->V; | |
| 361 | 80 | Vec w; | |
| 362 | |||
| 363 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
80 | PetscFunctionBegin; |
| 364 |
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.
|
80 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 365 |
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.
|
80 | PetscAssertPointer(v,2); |
| 366 |
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.
|
80 | PetscValidHeaderSpecific(*v,VEC_CLASSID,2); |
| 367 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
80 | EPSCheckSolved(eps,1); |
| 368 |
2/6✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
|
80 | PetscCheck(eps->ishermitian || eps->state!=EPS_STATE_EIGENVECTORS,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONGSTATE,"EPSGetInvariantSubspace must be called before EPSGetEigenpair,EPSGetEigenvector or EPSComputeError"); |
| 369 |
3/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
|
80 | if (eps->balance!=EPS_BALANCE_NONE && eps->D) { |
| 370 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10 | PetscCall(BVDuplicateResize(eps->V,eps->nconv,&V)); |
| 371 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10 | PetscCall(BVSetActiveColumns(eps->V,0,eps->nconv)); |
| 372 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10 | PetscCall(BVCopy(eps->V,V)); |
| 373 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
55 | for (i=0;i<eps->nconv;i++) { |
| 374 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
45 | PetscCall(BVGetColumn(V,i,&w)); |
| 375 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
45 | PetscCall(VecPointwiseDivide(w,w,eps->D)); |
| 376 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
45 | PetscCall(BVRestoreColumn(V,i,&w)); |
| 377 | } | ||
| 378 |
3/6✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
|
10 | PetscCall(BVOrthogonalize(V,NULL)); |
| 379 | } | ||
| 380 |
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.
|
625 | for (i=0;i<eps->nconv;i++) PetscCall(BVCopyVec(V,i,v[i])); |
| 381 |
7/10✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
80 | if (eps->balance!=EPS_BALANCE_NONE && eps->D) PetscCall(BVDestroy(&V)); |
| 382 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
16 | PetscFunctionReturn(PETSC_SUCCESS); |
| 383 | } | ||
| 384 | |||
| 385 | /*@ | ||
| 386 | EPSGetEigenpair - Gets the `i`-th solution of the eigenproblem as computed by | ||
| 387 | `EPSSolve()`. The solution consists in both the eigenvalue and the eigenvector. | ||
| 388 | |||
| 389 | Collective | ||
| 390 | |||
| 391 | Input Parameters: | ||
| 392 | + eps - the linear eigensolver context | ||
| 393 | - i - index of the solution | ||
| 394 | |||
| 395 | Output Parameters: | ||
| 396 | + eigr - real part of eigenvalue | ||
| 397 | . eigi - imaginary part of eigenvalue | ||
| 398 | . Vr - real part of eigenvector | ||
| 399 | - Vi - imaginary part of eigenvector | ||
| 400 | |||
| 401 | Notes: | ||
| 402 | It is allowed to pass `NULL` for `Vr` and `Vi`, if the eigenvector is not | ||
| 403 | required. Otherwise, the caller must provide valid `Vec` objects, i.e., | ||
| 404 | they must be created by the calling program with e.g. `MatCreateVecs()`. | ||
| 405 | |||
| 406 | If the eigenvalue is real, then `eigi` and `Vi` are set to zero. If PETSc is | ||
| 407 | configured with complex scalars the eigenvalue is stored | ||
| 408 | directly in `eigr` (`eigi` is set to zero) and the eigenvector in `Vr` (`Vi` is | ||
| 409 | set to zero). In both cases, the user can pass `NULL` in `eigi` and `Vi`. | ||
| 410 | |||
| 411 | The index `i` should be a value between 0 and `nconv`-1 (see `EPSGetConverged()`). | ||
| 412 | Eigenpairs are indexed according to the ordering criterion established | ||
| 413 | with `EPSSetWhichEigenpairs()`. | ||
| 414 | |||
| 415 | The 2-norm of the eigenvector is one unless the problem is generalized | ||
| 416 | Hermitian. In this case the eigenvector is normalized with respect to the | ||
| 417 | norm defined by the $B$ matrix. | ||
| 418 | |||
| 419 | In case of structured eigenproblems such as `EPS_BSE`, see the discussion about | ||
| 420 | [](#sec:structured-vectors). | ||
| 421 | |||
| 422 | Level: beginner | ||
| 423 | |||
| 424 | .seealso: [](ch:eps), `EPSGetEigenvalue()`, `EPSGetEigenvector()`, `EPSGetLeftEigenvector()`, `EPSSolve()`, `EPSGetConverged()`, `EPSSetWhichEigenpairs()`, `EPSGetInvariantSubspace()` | ||
| 425 | @*/ | ||
| 426 | 56993 | PetscErrorCode EPSGetEigenpair(EPS eps,PetscInt i,PetscScalar *eigr,PetscScalar *eigi,Vec Vr,Vec Vi) | |
| 427 | { | ||
| 428 | 56993 | PetscInt nconv; | |
| 429 | |||
| 430 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
56993 | PetscFunctionBegin; |
| 431 |
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.
|
56993 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 432 |
27/62✓ 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 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✓ Branch 59 taken 2 times.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
|
56993 | PetscValidLogicalCollectiveInt(eps,i,2); |
| 433 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
56993 | EPSCheckSolved(eps,1); |
| 434 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
56993 | PetscCheck(i>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index cannot be negative"); |
| 435 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
56993 | PetscCall(EPS_GetActualConverged(eps,&nconv)); |
| 436 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
56993 | PetscCheck(i<nconv,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index can be nconv-1 at most, see EPSGetConverged()"); |
| 437 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
56993 | PetscCall(EPSGetEigenvalue(eps,i,eigr,eigi)); |
| 438 |
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.
|
56993 | if (Vr || Vi) PetscCall(EPSGetEigenvector(eps,i,Vr,Vi)); |
| 439 |
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.
|
11399 | PetscFunctionReturn(PETSC_SUCCESS); |
| 440 | } | ||
| 441 | |||
| 442 | /*@ | ||
| 443 | EPSGetEigenvalue - Gets the `i`-th eigenvalue as computed by `EPSSolve()`. | ||
| 444 | |||
| 445 | Not Collective | ||
| 446 | |||
| 447 | Input Parameters: | ||
| 448 | + eps - the linear eigensolver context | ||
| 449 | - i - index of the solution | ||
| 450 | |||
| 451 | Output Parameters: | ||
| 452 | + eigr - real part of eigenvalue | ||
| 453 | - eigi - imaginary part of eigenvalue | ||
| 454 | |||
| 455 | Notes: | ||
| 456 | If the eigenvalue is real, then `eigi` is set to zero. If PETSc is | ||
| 457 | configured with complex scalars the eigenvalue is stored | ||
| 458 | directly in `eigr` (`eigi` is set to zero). | ||
| 459 | |||
| 460 | The index `i` should be a value between 0 and `nconv`-1 (see `EPSGetConverged()`). | ||
| 461 | Eigenpairs are indexed according to the ordering criterion established | ||
| 462 | with `EPSSetWhichEigenpairs()`. | ||
| 463 | |||
| 464 | Level: beginner | ||
| 465 | |||
| 466 | .seealso: [](ch:eps), `EPSSolve()`, `EPSGetConverged()`, `EPSSetWhichEigenpairs()`, `EPSGetEigenpair()` | ||
| 467 | @*/ | ||
| 468 | 112058 | PetscErrorCode EPSGetEigenvalue(EPS eps,PetscInt i,PetscScalar *eigr,PetscScalar *eigi) | |
| 469 | { | ||
| 470 | 112058 | PetscInt k,nconv; | |
| 471 | #if !PetscDefined(USE_COMPLEX) | ||
| 472 | 59101 | PetscInt k2, iquad; | |
| 473 | #endif | ||
| 474 | |||
| 475 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
112058 | PetscFunctionBegin; |
| 476 |
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.
|
112058 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 477 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
112058 | EPSCheckSolved(eps,1); |
| 478 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
112058 | PetscCheck(i>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index cannot be negative"); |
| 479 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
112058 | PetscCall(EPS_GetActualConverged(eps,&nconv)); |
| 480 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
112058 | PetscCheck(i<nconv,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index can be nconv-1 at most, see EPSGetConverged()"); |
| 481 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
112058 | if (nconv==eps->nconv) { |
| 482 | 101850 | k = eps->perm[i]; | |
| 483 | #if PetscDefined(USE_COMPLEX) | ||
| 484 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
49229 | if (eigr) *eigr = eps->eigr[k]; |
| 485 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
49229 | if (eigi) *eigi = 0; |
| 486 | #else | ||
| 487 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
52621 | if (eigr) *eigr = eps->eigr[k]; |
| 488 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
52621 | if (eigi) *eigi = eps->eigi[k]; |
| 489 | #endif | ||
| 490 | } else { | ||
| 491 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
10208 | PetscCheck(eps->problem_type==EPS_BSE || eps->problem_type==EPS_HAMILT || eps->problem_type==EPS_LREP,PetscObjectComm((PetscObject)eps),PETSC_ERR_PLIB,"Problem type should be BSE, Hamiltonian, or LREP"); |
| 492 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 5 times.
|
10208 | if (eps->problem_type==EPS_BSE || eps->problem_type==EPS_LREP) { |
| 493 | /* BSE problem, even index is +lambda, odd index is -lambda */ | ||
| 494 | 9688 | k = eps->perm[i/2]; | |
| 495 | #if PetscDefined(USE_COMPLEX) | ||
| 496 |
3/4✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 5 times.
✓ Branch 3 taken 5 times.
|
3728 | if (eigr) *eigr = (i%2)? -eps->eigr[k]: eps->eigr[k]; |
| 497 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
3728 | if (eigi) *eigi = 0; |
| 498 | #else | ||
| 499 |
3/4✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 5 times.
✓ Branch 3 taken 5 times.
|
5960 | if (eigr) *eigr = (i%2)? -eps->eigr[k]: eps->eigr[k]; |
| 500 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
5960 | if (eigi) *eigi = eps->eigi[k]; |
| 501 | #endif | ||
| 502 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
520 | } else if (eps->problem_type==EPS_HAMILT) { |
| 503 | /* Hamiltonian eigenproblem */ | ||
| 504 | 520 | k = eps->perm[i/2]; | |
| 505 | #if PetscDefined(USE_COMPLEX) | ||
| 506 | ✗ | if (eigr) *eigr = (i%2)? -eps->eigr[k]: eps->eigr[k]; | |
| 507 | ✗ | if (eigi) *eigi = 0; | |
| 508 | #else | ||
| 509 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 5 times.
|
520 | if (eps->eigi[k]==0.0) { /* real eigenvalue */ |
| 510 | ✗ | if (eigr) *eigr = (i%2)? -eps->eigr[k]: eps->eigr[k]; | |
| 511 | ✗ | if (eigi) *eigi = 0.0; | |
| 512 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
520 | } else if (eps->eigr[k]==0.0) { /* purely imaginary eigenvalue */ |
| 513 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
120 | if (eigr) *eigr = 0.0; |
| 514 |
3/4✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 5 times.
✓ Branch 3 taken 5 times.
|
120 | if (eigi) *eigi = (i%2)? -eps->eigi[k]: eps->eigi[k]; |
| 515 | } else { /* quadruple eigenvalue (-conj(lambda),-lambda,lambda,conj(lambda)) */ | ||
| 516 | 400 | iquad = i%2; /* index within the 4 values */ | |
| 517 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
400 | if (i>1) { |
| 518 | 320 | k2 = eps->perm[(i-2)/2]; | |
| 519 |
3/4✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 5 times.
✗ Branch 3 not taken.
|
320 | if (eps->eigr[k]==eps->eigr[k2] && eps->eigi[k]==-eps->eigi[k2]) iquad += 2; |
| 520 | } | ||
| 521 |
3/4✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 5 times.
✓ Branch 3 taken 5 times.
|
400 | if (eigr) *eigr = (iquad<2)? -eps->eigr[k]: eps->eigr[k]; |
| 522 |
3/4✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 5 times.
✓ Branch 3 taken 5 times.
|
400 | if (eigi) *eigi = (iquad%3)? -eps->eigi[k]: eps->eigi[k]; |
| 523 | } | ||
| 524 | #endif | ||
| 525 | } | ||
| 526 | } | ||
| 527 |
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.
|
22777 | PetscFunctionReturn(PETSC_SUCCESS); |
| 528 | } | ||
| 529 | |||
| 530 | /*@ | ||
| 531 | EPSGetEigenvector - Gets the `i`-th right eigenvector as computed by `EPSSolve()`. | ||
| 532 | |||
| 533 | Collective | ||
| 534 | |||
| 535 | Input Parameters: | ||
| 536 | + eps - the linear eigensolver context | ||
| 537 | - i - index of the solution | ||
| 538 | |||
| 539 | Output Parameters: | ||
| 540 | + Vr - real part of eigenvector | ||
| 541 | - Vi - imaginary part of eigenvector | ||
| 542 | |||
| 543 | Notes: | ||
| 544 | The caller must provide valid `Vec` objects, i.e., they must be created | ||
| 545 | by the calling program with e.g. `MatCreateVecs()`. | ||
| 546 | |||
| 547 | If the corresponding eigenvalue is real, then `Vi` is set to zero. If PETSc is | ||
| 548 | configured with complex scalars the eigenvector is stored | ||
| 549 | directly in `Vr` (`Vi` is set to zero). In any case, the user can pass `NULL` in `Vr` | ||
| 550 | or `Vi` if one of them is not required. | ||
| 551 | |||
| 552 | The index `i` should be a value between 0 and `nconv`-1 (see `EPSGetConverged()`). | ||
| 553 | Eigenpairs are indexed according to the ordering criterion established | ||
| 554 | with `EPSSetWhichEigenpairs()`. | ||
| 555 | |||
| 556 | The 2-norm of the eigenvector is one unless the problem is generalized | ||
| 557 | Hermitian. In this case the eigenvector is normalized with respect to the | ||
| 558 | norm defined by the $B$ matrix. | ||
| 559 | |||
| 560 | In case of structured eigenproblems such as `EPS_BSE`, see the discussion about | ||
| 561 | [](#sec:structured-vectors). | ||
| 562 | |||
| 563 | Level: beginner | ||
| 564 | |||
| 565 | .seealso: [](ch:eps), `EPSSolve()`, `EPSGetConverged()`, `EPSSetWhichEigenpairs()`, `EPSGetEigenpair()`, `EPSGetLeftEigenvector()` | ||
| 566 | @*/ | ||
| 567 | 89945 | PetscErrorCode EPSGetEigenvector(EPS eps,PetscInt i,Vec Vr,Vec Vi) | |
| 568 | { | ||
| 569 | 89945 | PetscInt nconv; | |
| 570 | |||
| 571 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
89945 | PetscFunctionBegin; |
| 572 |
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.
|
89945 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 573 |
27/62✓ 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 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✓ Branch 59 taken 2 times.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
|
89945 | PetscValidLogicalCollectiveInt(eps,i,2); |
| 574 |
16/46✓ 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 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✓ 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 taken 2 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✓ Branch 26 taken 2 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✓ Branch 30 taken 2 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✓ Branch 35 taken 2 times.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✓ Branch 39 taken 2 times.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
|
89945 | if (Vr) { PetscValidHeaderSpecific(Vr,VEC_CLASSID,3); PetscCheckSameComm(eps,1,Vr,3); } |
| 575 |
17/46✓ 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 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✓ 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 taken 2 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✓ Branch 26 taken 2 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✓ Branch 30 taken 2 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✓ Branch 35 taken 2 times.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✓ Branch 39 taken 2 times.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
|
89945 | if (Vi) { PetscValidHeaderSpecific(Vi,VEC_CLASSID,4); PetscCheckSameComm(eps,1,Vi,4); } |
| 576 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
89945 | EPSCheckSolved(eps,1); |
| 577 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
89945 | PetscCheck(i>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index cannot be negative"); |
| 578 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
89945 | PetscCall(EPS_GetActualConverged(eps,&nconv)); |
| 579 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
89945 | PetscCheck(i<nconv,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index can be nconv-1 at most, see EPSGetConverged()"); |
| 580 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
89945 | PetscCall(EPSComputeVectors(eps)); |
| 581 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
89945 | PetscCall(EPS_GetEigenvector(eps,eps->V,i,Vr,Vi)); |
| 582 |
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.
|
16053 | PetscFunctionReturn(PETSC_SUCCESS); |
| 583 | } | ||
| 584 | |||
| 585 | /*@ | ||
| 586 | EPSGetLeftEigenvector - Gets the `i`-th left eigenvector as computed by `EPSSolve()`. | ||
| 587 | |||
| 588 | Collective | ||
| 589 | |||
| 590 | Input Parameters: | ||
| 591 | + eps - the linear eigensolver context | ||
| 592 | - i - index of the solution | ||
| 593 | |||
| 594 | Output Parameters: | ||
| 595 | + Wr - real part of left eigenvector | ||
| 596 | - Wi - imaginary part of left eigenvector | ||
| 597 | |||
| 598 | Notes: | ||
| 599 | The caller must provide valid `Vec` objects, i.e., they must be created | ||
| 600 | by the calling program with e.g. `MatCreateVecs()`. | ||
| 601 | |||
| 602 | If the corresponding eigenvalue is real, then `Wi` is set to zero. If PETSc is | ||
| 603 | configured with complex scalars the eigenvector is stored directly in `Wr` | ||
| 604 | (`Wi` is set to zero). In any case, the user can pass `NULL` in `Wr` or `Wi` if | ||
| 605 | one of them is not required. | ||
| 606 | |||
| 607 | The index `i` should be a value between 0 and `nconv`-1 (see `EPSGetConverged()`). | ||
| 608 | Eigensolutions are indexed according to the ordering criterion established | ||
| 609 | with `EPSSetWhichEigenpairs()`. | ||
| 610 | |||
| 611 | Left eigenvectors are available only if the `twosided` flag was set, see | ||
| 612 | `EPSSetTwoSided()`. | ||
| 613 | |||
| 614 | In case of structured eigenproblems such as `EPS_BSE`, see the discussion about | ||
| 615 | [](#sec:structured-vectors). | ||
| 616 | |||
| 617 | Level: intermediate | ||
| 618 | |||
| 619 | .seealso: [](ch:eps), `EPSGetEigenvector()`, `EPSGetConverged()`, `EPSSetWhichEigenpairs()`, `EPSSetTwoSided()` | ||
| 620 | @*/ | ||
| 621 | 10665 | PetscErrorCode EPSGetLeftEigenvector(EPS eps,PetscInt i,Vec Wr,Vec Wi) | |
| 622 | { | ||
| 623 | 10665 | PetscInt nconv; | |
| 624 | 10665 | PetscBool trivial,lrepreduced=PETSC_FALSE; | |
| 625 | 10665 | Mat H; | |
| 626 | 10665 | IS is[2]; | |
| 627 | 10665 | Vec v0,v1,z,z0,z1; | |
| 628 | |||
| 629 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
10665 | PetscFunctionBegin; |
| 630 |
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.
|
10665 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 631 |
27/62✓ 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 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✓ Branch 59 taken 2 times.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
|
10665 | PetscValidLogicalCollectiveInt(eps,i,2); |
| 632 |
16/46✓ 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 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✓ 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 taken 2 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✓ Branch 26 taken 2 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✓ Branch 30 taken 2 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✓ Branch 35 taken 2 times.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✓ Branch 39 taken 2 times.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
|
10665 | if (Wr) { PetscValidHeaderSpecific(Wr,VEC_CLASSID,3); PetscCheckSameComm(eps,1,Wr,3); } |
| 633 |
17/46✓ 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 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✓ 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 taken 2 times.
✗ Branch 21 not taken.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✓ Branch 26 taken 2 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✓ Branch 30 taken 2 times.
✗ Branch 31 not taken.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✓ Branch 35 taken 2 times.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✓ Branch 39 taken 2 times.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
|
10665 | if (Wi) { PetscValidHeaderSpecific(Wi,VEC_CLASSID,4); PetscCheckSameComm(eps,1,Wi,4); } |
| 634 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
10665 | EPSCheckSolved(eps,1); |
| 635 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
10665 | PetscCheck(i>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index cannot be negative"); |
| 636 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10665 | PetscCall(EPS_GetActualConverged(eps,&nconv)); |
| 637 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
10665 | PetscCheck(i<nconv,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index can be nconv-1 at most, see EPSGetConverged()"); |
| 638 | |||
| 639 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 1 times.
|
10665 | trivial = (eps->problem_type==EPS_HEP || eps->problem_type==EPS_GHEP || eps->problem_type==EPS_BSE || eps->problem_type==EPS_LREP)? PETSC_TRUE: PETSC_FALSE; |
| 640 |
3/6✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
|
10665 | if (!trivial) PetscCheck(eps->twosided,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONGSTATE,"Must request left vectors with EPSSetTwoSided"); |
| 641 | |||
| 642 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10665 | PetscCall(EPSComputeVectors(eps)); |
| 643 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10665 | if (trivial) { |
| 644 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
9772 | PetscCall(EPS_GetEigenvector(eps,eps->V,i,Wr,Wi)); |
| 645 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
9772 | PetscCall(STGetMatrix(eps->st,0,&H)); |
| 646 |
6/8✓ Branch 0 taken 5 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
9772 | if (eps->problem_type==EPS_LREP) PetscCall(SlepcCheckMatLREPReduced(H,&lrepreduced)); |
| 647 |
4/4✓ Branch 0 taken 5 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 5 times.
✓ Branch 3 taken 5 times.
|
9772 | if (eps->problem_type==EPS_BSE || !lrepreduced) { /* change sign of bottom part of the vector */ |
| 648 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7524 | PetscCall(MatNestGetISs(H,is,NULL)); |
| 649 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
7524 | if (Wr) { |
| 650 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7524 | PetscCall(VecGetSubVector(Wr,is[1],&v1)); |
| 651 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7524 | PetscCall(VecScale(v1,-1.0)); |
| 652 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7524 | PetscCall(VecRestoreSubVector(Wr,is[1],&v1)); |
| 653 | } | ||
| 654 | #if !PetscDefined(USE_COMPLEX) | ||
| 655 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 5 times.
|
4864 | if (Wi) { |
| 656 | ✗ | PetscCall(VecGetSubVector(Wi,is[1],&v1)); | |
| 657 | ✗ | PetscCall(VecScale(v1,-1.0)); | |
| 658 | ✗ | PetscCall(VecRestoreSubVector(Wi,is[1],&v1)); | |
| 659 | } | ||
| 660 | #endif | ||
| 661 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
2248 | } else if (eps->problem_type==EPS_LREP) { /* swap the two parts of the vector */ |
| 662 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
2248 | if (Wr) { |
| 663 |
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.
|
2248 | PetscCall(VecDuplicate(Wr,&z)); |
| 664 |
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.
|
2248 | PetscCall(VecCopy(Wr,z)); |
| 665 |
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.
|
2248 | PetscCall(MatNestGetISs(H,is,NULL)); |
| 666 |
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.
|
2248 | PetscCall(VecGetSubVector(Wr,is[0],&v0)); |
| 667 |
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.
|
2248 | PetscCall(VecGetSubVector(Wr,is[1],&v1)); |
| 668 |
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.
|
2248 | PetscCall(VecGetSubVector(z,is[0],&z0)); |
| 669 |
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.
|
2248 | PetscCall(VecGetSubVector(z,is[1],&z1)); |
| 670 |
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.
|
2248 | PetscCall(VecCopy(z0,v1)); |
| 671 |
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.
|
2248 | PetscCall(VecCopy(z1,v0)); |
| 672 |
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.
|
2248 | PetscCall(VecRestoreSubVector(Wr,is[0],&v0)); |
| 673 |
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.
|
2248 | PetscCall(VecRestoreSubVector(Wr,is[1],&v1)); |
| 674 |
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.
|
2248 | PetscCall(VecRestoreSubVector(z,is[0],&z0)); |
| 675 |
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.
|
2248 | PetscCall(VecRestoreSubVector(z,is[1],&z1)); |
| 676 |
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.
|
2248 | PetscCall(VecDestroy(&z)); |
| 677 | } | ||
| 678 | } | ||
| 679 | } else { | ||
| 680 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
893 | PetscCall(EPS_GetEigenvector(eps,eps->W,i,Wr,Wi)); |
| 681 | } | ||
| 682 |
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.
|
1831 | PetscFunctionReturn(PETSC_SUCCESS); |
| 683 | } | ||
| 684 | |||
| 685 | /*@ | ||
| 686 | EPSGetErrorEstimate - Returns the error estimate associated to the `i`-th | ||
| 687 | computed eigenpair. | ||
| 688 | |||
| 689 | Not Collective | ||
| 690 | |||
| 691 | Input Parameters: | ||
| 692 | + eps - the linear eigensolver context | ||
| 693 | - i - index of eigenpair | ||
| 694 | |||
| 695 | Output Parameter: | ||
| 696 | . errest - the error estimate | ||
| 697 | |||
| 698 | Note: | ||
| 699 | This is the error estimate used internally by the eigensolver. The actual | ||
| 700 | error bound can be computed with `EPSComputeError()`. See discussion at | ||
| 701 | section [](#sec:errbnd). | ||
| 702 | |||
| 703 | Level: advanced | ||
| 704 | |||
| 705 | .seealso: [](ch:eps), [](#sec:errbnd), `EPSComputeError()` | ||
| 706 | @*/ | ||
| 707 | 612 | PetscErrorCode EPSGetErrorEstimate(EPS eps,PetscInt i,PetscReal *errest) | |
| 708 | { | ||
| 709 | 612 | PetscInt nconv; | |
| 710 | |||
| 711 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
612 | PetscFunctionBegin; |
| 712 |
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.
|
612 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 713 |
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.
|
612 | PetscAssertPointer(errest,3); |
| 714 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
612 | EPSCheckSolved(eps,1); |
| 715 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
612 | PetscCheck(i>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index cannot be negative"); |
| 716 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
612 | PetscCall(EPS_GetActualConverged(eps,&nconv)); |
| 717 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
612 | PetscCheck(i<nconv,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"The index can be nconv-1 at most, see EPSGetConverged()"); |
| 718 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
612 | if (nconv==eps->nconv) { |
| 719 | 612 | *errest = eps->errest[eps->perm[i]]; | |
| 720 | } else { | ||
| 721 | ✗ | PetscCheck(eps->problem_type==EPS_BSE || eps->problem_type==EPS_LREP || eps->problem_type==EPS_HAMILT,PetscObjectComm((PetscObject)eps),PETSC_ERR_PLIB,"Wrong problem type"); | |
| 722 | /* even index is +lambda, odd index is -lambda, assume both have same error */ | ||
| 723 | ✗ | *errest = eps->errest[eps->perm[i/2]]; | |
| 724 | } | ||
| 725 |
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.
|
116 | PetscFunctionReturn(PETSC_SUCCESS); |
| 726 | } | ||
| 727 | |||
| 728 | /* | ||
| 729 | EPSComputeResidualNorm_Private - Computes the norm of the residual vector | ||
| 730 | associated with an eigenpair. | ||
| 731 | |||
| 732 | Input Parameters: | ||
| 733 | trans - whether A' must be used instead of A | ||
| 734 | kr,ki - eigenvalue | ||
| 735 | xr,xi - eigenvector | ||
| 736 | z - three work vectors (the second one not referenced in complex scalars) | ||
| 737 | */ | ||
| 738 | 46347 | PetscErrorCode EPSComputeResidualNorm_Private(EPS eps,PetscBool trans,PetscScalar kr,PetscScalar ki,Vec xr,Vec xi,Vec *z,PetscReal *norm) | |
| 739 | { | ||
| 740 | 46347 | PetscInt nmat; | |
| 741 | 46347 | Mat A,B; | |
| 742 | 46347 | Vec u,w; | |
| 743 | 46347 | PetscScalar alpha; | |
| 744 | #if !PetscDefined(USE_COMPLEX) | ||
| 745 | 24152 | Vec v; | |
| 746 | 24152 | PetscReal ni,nr; | |
| 747 | #endif | ||
| 748 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
46347 | PetscErrorCode (*matmult)(Mat,Vec,Vec) = trans? MatMultHermitianTranspose: MatMult; |
| 749 | |||
| 750 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
46347 | PetscFunctionBegin; |
| 751 | 46347 | u = z[0]; w = z[2]; | |
| 752 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
46347 | PetscCall(STGetNumMatrices(eps->st,&nmat)); |
| 753 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
46347 | PetscCall(STGetMatrix(eps->st,0,&A)); |
| 754 |
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.
|
46347 | if (nmat>1) PetscCall(STGetMatrix(eps->st,1,&B)); |
| 755 | |||
| 756 | #if !PetscDefined(USE_COMPLEX) | ||
| 757 | 24152 | v = z[1]; | |
| 758 |
3/4✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 5 times.
|
24152 | if (ki == 0 || PetscAbsScalar(ki) < PetscAbsScalar(kr*PETSC_MACHINE_EPSILON)) { |
| 759 | #endif | ||
| 760 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
45177 | PetscCall((*matmult)(A,xr,u)); /* u=A*x */ |
| 761 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
45177 | if (PetscAbsScalar(kr) > PETSC_MACHINE_EPSILON) { |
| 762 |
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.
|
45177 | if (nmat>1) PetscCall((*matmult)(B,xr,w)); |
| 763 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
23406 | else PetscCall(VecCopy(xr,w)); /* w=B*x */ |
| 764 |
2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
|
45177 | alpha = trans? -PetscConj(kr): -kr; |
| 765 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
45177 | PetscCall(VecAXPY(u,alpha,w)); /* u=A*x-k*B*x */ |
| 766 | } | ||
| 767 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
45177 | PetscCall(VecNorm(u,NORM_2,norm)); |
| 768 | #if !PetscDefined(USE_COMPLEX) | ||
| 769 | } else { | ||
| 770 |
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.
|
1170 | PetscCall((*matmult)(A,xr,u)); /* u=A*xr */ |
| 771 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
1170 | if (SlepcAbsEigenvalue(kr,ki) > PETSC_MACHINE_EPSILON) { |
| 772 |
6/8✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
1170 | if (nmat>1) PetscCall((*matmult)(B,xr,v)); |
| 773 |
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.
|
1156 | else PetscCall(VecCopy(xr,v)); /* v=B*xr */ |
| 774 |
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.
|
1170 | PetscCall(VecAXPY(u,-kr,v)); /* u=A*xr-kr*B*xr */ |
| 775 |
6/8✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
1170 | if (nmat>1) PetscCall((*matmult)(B,xi,w)); |
| 776 |
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.
|
1156 | else PetscCall(VecCopy(xi,w)); /* w=B*xi */ |
| 777 |
7/8✓ Branch 0 taken 5 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 1 times.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
1170 | PetscCall(VecAXPY(u,trans?-ki:ki,w)); /* u=A*xr-kr*B*xr+ki*B*xi */ |
| 778 | } | ||
| 779 |
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.
|
1170 | PetscCall(VecNorm(u,NORM_2,&nr)); |
| 780 |
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.
|
1170 | PetscCall((*matmult)(A,xi,u)); /* u=A*xi */ |
| 781 |
1/2✓ Branch 0 taken 5 times.
✗ Branch 1 not taken.
|
1170 | if (SlepcAbsEigenvalue(kr,ki) > PETSC_MACHINE_EPSILON) { |
| 782 |
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.
|
1170 | PetscCall(VecAXPY(u,-kr,w)); /* u=A*xi-kr*B*xi */ |
| 783 |
7/8✓ Branch 0 taken 5 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 1 times.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
1170 | PetscCall(VecAXPY(u,trans?ki:-ki,v)); /* u=A*xi-kr*B*xi-ki*B*xr */ |
| 784 | } | ||
| 785 |
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.
|
1170 | PetscCall(VecNorm(u,NORM_2,&ni)); |
| 786 | 1170 | *norm = SlepcAbsEigenvalue(nr,ni); | |
| 787 | } | ||
| 788 | #endif | ||
| 789 |
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.
|
9139 | PetscFunctionReturn(PETSC_SUCCESS); |
| 790 | } | ||
| 791 | |||
| 792 | /*@ | ||
| 793 | EPSComputeError - Computes the error (based on the residual norm) associated | ||
| 794 | with the `i`-th computed eigenpair. | ||
| 795 | |||
| 796 | Collective | ||
| 797 | |||
| 798 | Input Parameters: | ||
| 799 | + eps - the linear eigensolver context | ||
| 800 | . i - the solution index | ||
| 801 | - type - the type of error to compute, see `EPSErrorType` | ||
| 802 | |||
| 803 | Output Parameter: | ||
| 804 | . error - the error | ||
| 805 | |||
| 806 | Notes: | ||
| 807 | The error can be computed in various ways, all of them based on the residual | ||
| 808 | norm $\|Ax-\lambda Bx\|_2$ where $(\lambda,x)$ is the approximate eigenpair. | ||
| 809 | |||
| 810 | If the computation of left eigenvectors was enabled with `EPSSetTwoSided()`, | ||
| 811 | then the error will be computed using the maximum of the value above and | ||
| 812 | the left residual norm $\|y^*A-\lambda y^*B\|_2$, where $y$ is the approximate left | ||
| 813 | eigenvector. | ||
| 814 | |||
| 815 | Level: beginner | ||
| 816 | |||
| 817 | .seealso: [](ch:eps), `EPSErrorType`, `EPSSolve()`, `EPSGetErrorEstimate()`, `EPSSetTwoSided()` | ||
| 818 | @*/ | ||
| 819 | 42962 | PetscErrorCode EPSComputeError(EPS eps,PetscInt i,EPSErrorType type,PetscReal *error) | |
| 820 | { | ||
| 821 | 42962 | Mat A,B; | |
| 822 | 42962 | Vec xr,xi,w[3]; | |
| 823 | 42962 | PetscReal t,vecnorm=1.0,errorl; | |
| 824 | 42962 | PetscScalar kr,ki; | |
| 825 | 42962 | PetscBool flg; | |
| 826 | |||
| 827 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
42962 | PetscFunctionBegin; |
| 828 |
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.
|
42962 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 829 |
27/62✓ 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 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✓ Branch 59 taken 2 times.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
|
42962 | PetscValidLogicalCollectiveInt(eps,i,2); |
| 830 |
27/62✓ 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 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✓ Branch 59 taken 2 times.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
|
42962 | PetscValidLogicalCollectiveEnum(eps,type,3); |
| 831 |
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.
|
42962 | PetscAssertPointer(error,4); |
| 832 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
42962 | EPSCheckSolved(eps,1); |
| 833 | |||
| 834 | /* allocate work vectors */ | ||
| 835 | #if PetscDefined(USE_COMPLEX) | ||
| 836 |
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.
|
20440 | PetscCall(EPSSetWorkVecs(eps,3)); |
| 837 | 20440 | xi = NULL; | |
| 838 | 20440 | w[1] = NULL; | |
| 839 | #else | ||
| 840 |
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.
|
22522 | PetscCall(EPSSetWorkVecs(eps,5)); |
| 841 | 22522 | xi = eps->work[3]; | |
| 842 | 22522 | w[1] = eps->work[4]; | |
| 843 | #endif | ||
| 844 | 42962 | xr = eps->work[0]; | |
| 845 | 42962 | w[0] = eps->work[1]; | |
| 846 | 42962 | w[2] = eps->work[2]; | |
| 847 | |||
| 848 | /* compute residual norm */ | ||
| 849 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42962 | PetscCall(EPSGetEigenpair(eps,i,&kr,&ki,xr,xi)); |
| 850 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42962 | PetscCall(EPSComputeResidualNorm_Private(eps,PETSC_FALSE,kr,ki,xr,xi,w,error)); |
| 851 | |||
| 852 | /* compute 2-norm of eigenvector */ | ||
| 853 |
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.
|
42962 | if (eps->problem_type==EPS_GHEP) PetscCall(VecNorm(xr,NORM_2,&vecnorm)); |
| 854 | |||
| 855 | /* if two-sided, compute left residual norm and take the maximum */ | ||
| 856 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
42962 | if (eps->twosided) { |
| 857 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
450 | PetscCall(EPSGetLeftEigenvector(eps,i,xr,xi)); |
| 858 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
450 | PetscCall(EPSComputeResidualNorm_Private(eps,PETSC_TRUE,kr,ki,xr,xi,w,&errorl)); |
| 859 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
667 | *error = PetscMax(*error,errorl); |
| 860 | } | ||
| 861 | |||
| 862 | /* compute error */ | ||
| 863 |
3/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
|
42962 | switch (type) { |
| 864 | case EPS_ERROR_ABSOLUTE: | ||
| 865 | break; | ||
| 866 | 25512 | case EPS_ERROR_RELATIVE: | |
| 867 | 25512 | *error /= SlepcAbsEigenvalue(kr,ki)*vecnorm; | |
| 868 | 25512 | break; | |
| 869 | 17198 | case EPS_ERROR_BACKWARD: | |
| 870 | /* initialization of matrix norms */ | ||
| 871 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
17198 | if (!eps->nrma) { |
| 872 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
395 | PetscCall(STGetMatrix(eps->st,0,&A)); |
| 873 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
395 | PetscCall(MatHasOperation(A,MATOP_NORM,&flg)); |
| 874 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
395 | PetscCheck(flg,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONG,"The computation of backward errors requires a matrix norm operation"); |
| 875 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
395 | PetscCall(MatNorm(A,NORM_INFINITY,&eps->nrma)); |
| 876 | } | ||
| 877 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
17198 | if (eps->isgeneralized) { |
| 878 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
16945 | if (!eps->nrmb) { |
| 879 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
336 | PetscCall(STGetMatrix(eps->st,1,&B)); |
| 880 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
336 | PetscCall(MatHasOperation(B,MATOP_NORM,&flg)); |
| 881 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
336 | PetscCheck(flg,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONG,"The computation of backward errors requires a matrix norm operation"); |
| 882 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
336 | PetscCall(MatNorm(B,NORM_INFINITY,&eps->nrmb)); |
| 883 | } | ||
| 884 | 253 | } else eps->nrmb = 1.0; | |
| 885 | 17198 | t = SlepcAbsEigenvalue(kr,ki); | |
| 886 | 17198 | *error /= (eps->nrma+t*eps->nrmb)*vecnorm; | |
| 887 | 17198 | break; | |
| 888 | ✗ | default: | |
| 889 | − | SETERRQ(PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"Invalid error type"); | |
| 890 | } | ||
| 891 |
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.
|
8423 | PetscFunctionReturn(PETSC_SUCCESS); |
| 892 | } | ||
| 893 | |||
| 894 | /* | ||
| 895 | EPSGetStartVector - Generate a suitable vector to be used as the starting vector | ||
| 896 | for the recurrence that builds the right subspace. | ||
| 897 | |||
| 898 | Collective | ||
| 899 | |||
| 900 | Input Parameters: | ||
| 901 | + eps - the linear eigensolver context | ||
| 902 | - i - iteration number | ||
| 903 | |||
| 904 | Output Parameter: | ||
| 905 | . breakdown - flag indicating that a breakdown has occurred | ||
| 906 | |||
| 907 | Notes: | ||
| 908 | The start vector is computed from another vector: for the first step (i=0), | ||
| 909 | the first initial vector is used (see EPSSetInitialSpace()); otherwise a random | ||
| 910 | vector is created. Then this vector is forced to be in the range of OP (only | ||
| 911 | for generalized definite problems) and orthonormalized with respect to all | ||
| 912 | V-vectors up to i-1. The resulting vector is placed in V[i]. | ||
| 913 | |||
| 914 | The flag breakdown is set to true if either i=0 and the vector belongs to the | ||
| 915 | deflation space, or i>0 and the vector is linearly dependent with respect | ||
| 916 | to the V-vectors. | ||
| 917 | */ | ||
| 918 | 7984 | PetscErrorCode EPSGetStartVector(EPS eps,PetscInt i,PetscBool *breakdown) | |
| 919 | { | ||
| 920 | 7984 | PetscReal norm; | |
| 921 | 7984 | PetscBool lindep; | |
| 922 | 7984 | Vec w,z; | |
| 923 | |||
| 924 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
7984 | PetscFunctionBegin; |
| 925 |
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.
|
7984 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 926 |
27/62✓ 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 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✓ Branch 59 taken 2 times.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
|
7984 | PetscValidLogicalCollectiveInt(eps,i,2); |
| 927 | |||
| 928 | /* For the first step, use the first initial vector, otherwise a random one */ | ||
| 929 |
8/10✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
7984 | if (i>0 || eps->nini==0) PetscCall(BVSetRandomColumn(eps->V,i)); |
| 930 | |||
| 931 | /* Force the vector to be in the range of OP for generalized problems with B-inner product */ | ||
| 932 |
6/6✓ 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.
|
7984 | if (eps->ispositive || (eps->isgeneralized && eps->ishermitian)) { |
| 933 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1033 | PetscCall(BVCreateVec(eps->V,&w)); |
| 934 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1033 | PetscCall(BVCopyVec(eps->V,i,w)); |
| 935 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1033 | PetscCall(BVGetColumn(eps->V,i,&z)); |
| 936 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1033 | PetscCall(STApply(eps->st,w,z)); |
| 937 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1033 | PetscCall(BVRestoreColumn(eps->V,i,&z)); |
| 938 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1033 | PetscCall(VecDestroy(&w)); |
| 939 | } | ||
| 940 | |||
| 941 | /* Orthonormalize the vector with respect to previous vectors */ | ||
| 942 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7984 | PetscCall(BVOrthogonalizeColumn(eps->V,i,NULL,&norm,&lindep)); |
| 943 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7984 | if (breakdown) *breakdown = lindep; |
| 944 |
2/4✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
|
7627 | else if (lindep || norm == 0.0) { |
| 945 | ✗ | PetscCheck(i,PetscObjectComm((PetscObject)eps),PETSC_ERR_PLIB,"Initial vector is zero or belongs to the deflation space"); | |
| 946 | ✗ | PetscCheck(!i,PetscObjectComm((PetscObject)eps),PETSC_ERR_CONV_FAILED,"Unable to generate more start vectors"); | |
| 947 | } | ||
| 948 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7984 | PetscCall(BVScaleColumn(eps->V,i,1.0/norm)); |
| 949 |
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.
|
1549 | PetscFunctionReturn(PETSC_SUCCESS); |
| 950 | } | ||
| 951 | |||
| 952 | /* | ||
| 953 | EPSGetLeftStartVector - Generate a suitable vector to be used as the left starting | ||
| 954 | vector for the recurrence that builds the left subspace. See EPSGetStartVector(). | ||
| 955 | */ | ||
| 956 | 254 | PetscErrorCode EPSGetLeftStartVector(EPS eps,PetscInt i,PetscBool *breakdown) | |
| 957 | { | ||
| 958 | 254 | PetscReal norm; | |
| 959 | 254 | PetscBool lindep; | |
| 960 | 254 | Vec w,z; | |
| 961 | |||
| 962 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
254 | PetscFunctionBegin; |
| 963 |
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.
|
254 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 964 |
27/62✓ 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 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 2 times.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✓ Branch 32 taken 2 times.
✗ Branch 33 not taken.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✓ Branch 37 taken 2 times.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✓ Branch 41 taken 2 times.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✓ Branch 45 taken 2 times.
✓ Branch 46 taken 2 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2 times.
✓ Branch 50 taken 2 times.
✗ Branch 51 not taken.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✓ Branch 55 taken 2 times.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✓ Branch 59 taken 2 times.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
|
254 | PetscValidLogicalCollectiveInt(eps,i,2); |
| 965 | |||
| 966 | /* For the first step, use the first initial vector, otherwise a random one */ | ||
| 967 |
8/10✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
254 | if (i>0 || eps->ninil==0) PetscCall(BVSetRandomColumn(eps->W,i)); |
| 968 | |||
| 969 | /* Force the vector to be in the range of OP' for generalized problems with B-inner product */ | ||
| 970 |
5/6✓ Branch 0 taken 10 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 8 times.
|
254 | if (eps->ispositive || (eps->isgeneralized && eps->ishermitian)) { |
| 971 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(BVCreateVec(eps->W,&w)); |
| 972 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(BVCopyVec(eps->W,i,w)); |
| 973 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(BVGetColumn(eps->W,i,&z)); |
| 974 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(STApplyHermitianTranspose(eps->st,w,z)); |
| 975 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(BVRestoreColumn(eps->W,i,&z)); |
| 976 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(VecDestroy(&w)); |
| 977 | } | ||
| 978 | |||
| 979 | /* Orthonormalize the vector with respect to previous vectors */ | ||
| 980 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
254 | PetscCall(BVOrthogonalizeColumn(eps->W,i,NULL,&norm,&lindep)); |
| 981 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
254 | if (breakdown) *breakdown = lindep; |
| 982 |
2/4✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
|
199 | else if (lindep || norm == 0.0) { |
| 983 | ✗ | PetscCheck(i,PetscObjectComm((PetscObject)eps),PETSC_ERR_PLIB,"Left initial vector is zero"); | |
| 984 | ✗ | PetscCheck(!i,PetscObjectComm((PetscObject)eps),PETSC_ERR_CONV_FAILED,"Unable to generate more left start vectors"); | |
| 985 | } | ||
| 986 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
254 | PetscCall(BVScaleColumn(eps->W,i,1.0/norm)); |
| 987 |
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.
|
50 | PetscFunctionReturn(PETSC_SUCCESS); |
| 988 | } | ||
| 989 |