| 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 problem setup | ||
| 12 | */ | ||
| 13 | |||
| 14 | #include <slepc/private/epsimpl.h> /*I "slepceps.h" I*/ | ||
| 15 | |||
| 16 | /* | ||
| 17 | Let the solver choose the ST type that should be used by default, | ||
| 18 | otherwise set it to SHIFT. | ||
| 19 | This is called at EPSSetFromOptions (before STSetFromOptions) | ||
| 20 | and also at EPSSetUp (in case EPSSetFromOptions was not called). | ||
| 21 | */ | ||
| 22 | 16823 | PetscErrorCode EPSSetDefaultST(EPS eps) | |
| 23 | { | ||
| 24 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
16823 | PetscFunctionBegin; |
| 25 |
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.
|
16823 | PetscTryTypeMethod(eps,setdefaultst); |
| 26 |
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.
|
16823 | if (!((PetscObject)eps->st)->type_name) PetscCall(STSetType(eps->st,STSHIFT)); |
| 27 |
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.
|
3252 | PetscFunctionReturn(PETSC_SUCCESS); |
| 28 | } | ||
| 29 | |||
| 30 | /* | ||
| 31 | This is done by preconditioned eigensolvers that use the PC only. | ||
| 32 | It sets STPRECOND with KSPPREONLY. | ||
| 33 | */ | ||
| 34 | 1048 | PetscErrorCode EPSSetDefaultST_Precond(EPS eps) | |
| 35 | { | ||
| 36 | 1048 | KSP ksp; | |
| 37 | |||
| 38 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1048 | PetscFunctionBegin; |
| 39 |
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.
|
1048 | if (!((PetscObject)eps->st)->type_name) PetscCall(STSetType(eps->st,STPRECOND)); |
| 40 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1048 | PetscCall(STGetKSP(eps->st,&ksp)); |
| 41 |
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.
|
1048 | if (!((PetscObject)ksp)->type_name) PetscCall(KSPSetType(ksp,KSPPREONLY)); |
| 42 |
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.
|
208 | PetscFunctionReturn(PETSC_SUCCESS); |
| 43 | } | ||
| 44 | |||
| 45 | /* | ||
| 46 | This is done by preconditioned eigensolvers that can also use the KSP. | ||
| 47 | It sets STPRECOND with the default KSP (GMRES) and maxit=5. | ||
| 48 | */ | ||
| 49 | 952 | PetscErrorCode EPSSetDefaultST_GMRES(EPS eps) | |
| 50 | { | ||
| 51 | 952 | KSP ksp; | |
| 52 | |||
| 53 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
952 | PetscFunctionBegin; |
| 54 |
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.
|
952 | if (!((PetscObject)eps->st)->type_name) PetscCall(STSetType(eps->st,STPRECOND)); |
| 55 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
952 | PetscCall(STPrecondSetKSPHasMat(eps->st,PETSC_TRUE)); |
| 56 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
952 | PetscCall(STGetKSP(eps->st,&ksp)); |
| 57 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
952 | if (!((PetscObject)ksp)->type_name) { |
| 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.
|
391 | PetscCall(KSPSetType(ksp,KSPGMRES)); |
| 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.
|
391 | PetscCall(KSPSetTolerances(ksp,PETSC_CURRENT,PETSC_CURRENT,PETSC_CURRENT,5)); |
| 60 | } | ||
| 61 |
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.
|
216 | PetscFunctionReturn(PETSC_SUCCESS); |
| 62 | } | ||
| 63 | |||
| 64 | #if defined(SLEPC_HAVE_SCALAPACK) || defined(SLEPC_HAVE_ELPA) || defined(SLEPC_HAVE_ELEMENTAL) || defined(SLEPC_HAVE_EVSL) | ||
| 65 | /* | ||
| 66 | This is for direct eigensolvers that work with A and B directly, so | ||
| 67 | no need to factorize B. | ||
| 68 | */ | ||
| 69 | 192 | PetscErrorCode EPSSetDefaultST_NoFactor(EPS eps) | |
| 70 | { | ||
| 71 | 192 | KSP ksp; | |
| 72 | 192 | PC pc; | |
| 73 | |||
| 74 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
192 | PetscFunctionBegin; |
| 75 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 5 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
192 | if (!((PetscObject)eps->st)->type_name) PetscCall(STSetType(eps->st,STSHIFT)); |
| 76 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
192 | PetscCall(STGetKSP(eps->st,&ksp)); |
| 77 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 5 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
192 | if (!((PetscObject)ksp)->type_name) PetscCall(KSPSetType(ksp,KSPPREONLY)); |
| 78 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
192 | PetscCall(KSPGetPC(ksp,&pc)); |
| 79 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 5 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
|
192 | if (!((PetscObject)pc)->type_name) PetscCall(PCSetType(pc,PCNONE)); |
| 80 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 1 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 1 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 1 times.
|
24 | PetscFunctionReturn(PETSC_SUCCESS); |
| 81 | } | ||
| 82 | #endif | ||
| 83 | |||
| 84 | /* | ||
| 85 | Check that the ST selected by the user is compatible with the EPS solver and options | ||
| 86 | */ | ||
| 87 | 10005 | static PetscErrorCode EPSCheckCompatibleST(EPS eps) | |
| 88 | { | ||
| 89 | 10005 | PetscBool precond,shift,sinvert,cayley,lyapii; | |
| 90 | #if defined(PETSC_USE_COMPLEX) | ||
| 91 | 4828 | PetscScalar sigma; | |
| 92 | #endif | ||
| 93 | |||
| 94 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
10005 | PetscFunctionBegin; |
| 95 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(PetscObjectTypeCompare((PetscObject)eps->st,STPRECOND,&precond)); |
| 96 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(PetscObjectTypeCompare((PetscObject)eps->st,STSHIFT,&shift)); |
| 97 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(PetscObjectTypeCompare((PetscObject)eps->st,STSINVERT,&sinvert)); |
| 98 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(PetscObjectTypeCompare((PetscObject)eps->st,STCAYLEY,&cayley)); |
| 99 | |||
| 100 | /* preconditioned eigensolvers */ | ||
| 101 |
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.
|
10005 | PetscCheck(eps->categ!=EPS_CATEGORY_PRECOND || precond,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"This solver requires ST=PRECOND"); |
| 102 |
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.
|
10005 | PetscCheck(eps->categ==EPS_CATEGORY_PRECOND || !precond,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"STPRECOND is intended for preconditioned eigensolvers only"); |
| 103 | |||
| 104 | /* harmonic extraction */ | ||
| 105 |
5/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 10 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
10005 | PetscCheck(precond || shift || !eps->extraction || eps->extraction==EPS_RITZ,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Cannot use a spectral transformation combined with harmonic extraction"); |
| 106 | |||
| 107 | /* real shifts in Hermitian problems */ | ||
| 108 | #if defined(PETSC_USE_COMPLEX) | ||
| 109 |
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.
|
4828 | PetscCall(STGetShift(eps->st,&sigma)); |
| 110 |
3/6✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 5 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
|
4828 | PetscCheck(!eps->ishermitian || PetscImaginaryPart(sigma)==0.0,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Hermitian problems are not compatible with complex shifts"); |
| 111 | #endif | ||
| 112 | |||
| 113 | /* Cayley with PGNHEP */ | ||
| 114 |
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.
|
10005 | PetscCheck(!cayley || eps->problem_type!=EPS_PGNHEP,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Cayley spectral transformation is not compatible with PGNHEP"); |
| 115 | |||
| 116 | /* make sure that the user does not specify smallest magnitude with shift-and-invert */ | ||
| 117 |
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.
|
10005 | if ((cayley || sinvert) && (eps->categ==EPS_CATEGORY_KRYLOV || eps->categ==EPS_CATEGORY_OTHER)) { |
| 118 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1292 | PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSLYAPII,&lyapii)); |
| 119 |
10/14✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 10 times.
✓ Branch 5 taken 10 times.
✓ Branch 6 taken 10 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 10 times.
✓ Branch 9 taken 10 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 10 times.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
|
1292 | PetscCheck(lyapii || eps->which==EPS_TARGET_MAGNITUDE || eps->which==EPS_TARGET_REAL || eps->which==EPS_TARGET_IMAGINARY || eps->which==EPS_ALL || eps->which==EPS_WHICH_USER,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"Shift-and-invert requires a target 'which' (see EPSSetWhichEigenpairs), for instance -st_type sinvert -eps_target 0 -eps_target_magnitude"); |
| 120 | } | ||
| 121 |
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.
|
1926 | PetscFunctionReturn(PETSC_SUCCESS); |
| 122 | } | ||
| 123 | |||
| 124 | /* | ||
| 125 | MatEstimateSpectralRange_EPS: estimate the spectral range [left,right] of a | ||
| 126 | symmetric/Hermitian matrix A using an auxiliary EPS object | ||
| 127 | */ | ||
| 128 | 146 | PetscErrorCode MatEstimateSpectralRange_EPS(Mat A,PetscReal *left,PetscReal *right) | |
| 129 | { | ||
| 130 | 146 | PetscInt nconv; | |
| 131 | 146 | PetscScalar eig0; | |
| 132 | 146 | PetscReal tol=1e-3,errest=tol; | |
| 133 | 146 | EPS eps; | |
| 134 | |||
| 135 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
146 | PetscFunctionBegin; |
| 136 | 146 | *left = 0.0; *right = 0.0; | |
| 137 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSCreate(PetscObjectComm((PetscObject)A),&eps)); |
| 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.
|
146 | PetscCall(EPSSetOptionsPrefix(eps,"eps_filter_")); |
| 139 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSSetOperators(eps,A,NULL)); |
| 140 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSSetProblemType(eps,EPS_HEP)); |
| 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.
|
146 | PetscCall(EPSSetTolerances(eps,tol,50)); |
| 142 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSSetConvergenceTest(eps,EPS_CONV_ABS)); |
| 143 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSSetWhichEigenpairs(eps,EPS_SMALLEST_REAL)); |
| 144 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSSolve(eps)); |
| 145 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSGetConverged(eps,&nconv)); |
| 146 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
146 | if (nconv>0) { |
| 147 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSGetEigenvalue(eps,0,&eig0,NULL)); |
| 148 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSGetErrorEstimate(eps,0,&errest)); |
| 149 | ✗ | } else eig0 = eps->eigr[0]; | |
| 150 | 146 | *left = PetscRealPart(eig0)-errest; | |
| 151 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSSetWhichEigenpairs(eps,EPS_LARGEST_REAL)); |
| 152 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSSolve(eps)); |
| 153 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSGetConverged(eps,&nconv)); |
| 154 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
146 | if (nconv>0) { |
| 155 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSGetEigenvalue(eps,0,&eig0,NULL)); |
| 156 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
146 | PetscCall(EPSGetErrorEstimate(eps,0,&errest)); |
| 157 | ✗ | } else eig0 = eps->eigr[0]; | |
| 158 | 146 | *right = PetscRealPart(eig0)+errest; | |
| 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.
|
146 | PetscCall(EPSDestroy(&eps)); |
| 160 |
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.
|
26 | PetscFunctionReturn(PETSC_SUCCESS); |
| 161 | } | ||
| 162 | |||
| 163 | /* | ||
| 164 | EPSSetUpSort_Basic: configure the EPS sorting criterion according to 'which' | ||
| 165 | */ | ||
| 166 | 9400 | PetscErrorCode EPSSetUpSort_Basic(EPS eps) | |
| 167 | { | ||
| 168 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
9400 | PetscFunctionBegin; |
| 169 |
11/11✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 10 times.
✓ Branch 5 taken 10 times.
✓ Branch 6 taken 10 times.
✓ Branch 7 taken 10 times.
✓ Branch 8 taken 5 times.
✓ Branch 9 taken 10 times.
✓ Branch 10 taken 10 times.
|
9400 | switch (eps->which) { |
| 170 | 3836 | case EPS_LARGEST_MAGNITUDE: | |
| 171 | 3836 | eps->sc->comparison = SlepcCompareLargestMagnitude; | |
| 172 | 3836 | eps->sc->comparisonctx = NULL; | |
| 173 | 3836 | break; | |
| 174 | 493 | case EPS_SMALLEST_MAGNITUDE: | |
| 175 | 493 | eps->sc->comparison = SlepcCompareSmallestMagnitude; | |
| 176 | 493 | eps->sc->comparisonctx = NULL; | |
| 177 | 493 | break; | |
| 178 | 1199 | case EPS_LARGEST_REAL: | |
| 179 | 1199 | eps->sc->comparison = SlepcCompareLargestReal; | |
| 180 | 1199 | eps->sc->comparisonctx = NULL; | |
| 181 | 1199 | break; | |
| 182 | 1634 | case EPS_SMALLEST_REAL: | |
| 183 | 1634 | eps->sc->comparison = SlepcCompareSmallestReal; | |
| 184 | 1634 | eps->sc->comparisonctx = NULL; | |
| 185 | 1634 | break; | |
| 186 | 14 | case EPS_LARGEST_IMAGINARY: | |
| 187 | 14 | eps->sc->comparison = SlepcCompareLargestImaginary; | |
| 188 | 14 | eps->sc->comparisonctx = NULL; | |
| 189 | 14 | break; | |
| 190 | 10 | case EPS_SMALLEST_IMAGINARY: | |
| 191 | 10 | eps->sc->comparison = SlepcCompareSmallestImaginary; | |
| 192 | 10 | eps->sc->comparisonctx = NULL; | |
| 193 | 10 | break; | |
| 194 | 1257 | case EPS_TARGET_MAGNITUDE: | |
| 195 | 1257 | eps->sc->comparison = SlepcCompareTargetMagnitude; | |
| 196 | 1257 | eps->sc->comparisonctx = &eps->target; | |
| 197 | 1257 | break; | |
| 198 | 52 | case EPS_TARGET_REAL: | |
| 199 | 52 | eps->sc->comparison = SlepcCompareTargetReal; | |
| 200 | 52 | eps->sc->comparisonctx = &eps->target; | |
| 201 | 52 | break; | |
| 202 | 5 | case EPS_TARGET_IMAGINARY: | |
| 203 | 5 | eps->sc->comparison = SlepcCompareTargetImaginary; | |
| 204 | 5 | eps->sc->comparisonctx = &eps->target; | |
| 205 | 5 | break; | |
| 206 | 645 | case EPS_ALL: | |
| 207 | 645 | eps->sc->comparison = SlepcCompareSmallestReal; | |
| 208 | 645 | eps->sc->comparisonctx = NULL; | |
| 209 | 645 | break; | |
| 210 | case EPS_WHICH_USER: | ||
| 211 | break; | ||
| 212 | } | ||
| 213 | 9400 | eps->sc->map = NULL; | |
| 214 | 9400 | eps->sc->mapobj = NULL; | |
| 215 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
9400 | PetscFunctionReturn(PETSC_SUCCESS); |
| 216 | } | ||
| 217 | |||
| 218 | /* | ||
| 219 | EPSSetUpSort_Default: configure both EPS and DS sorting criterion | ||
| 220 | */ | ||
| 221 | 9146 | PetscErrorCode EPSSetUpSort_Default(EPS eps) | |
| 222 | { | ||
| 223 | 9146 | SlepcSC sc; | |
| 224 | 9146 | PetscBool istrivial; | |
| 225 | |||
| 226 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
9146 | PetscFunctionBegin; |
| 227 | /* fill sorting criterion context */ | ||
| 228 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
9146 | PetscCall(EPSSetUpSort_Basic(eps)); |
| 229 | /* fill sorting criterion for DS */ | ||
| 230 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
9146 | PetscCall(DSGetSlepcSC(eps->ds,&sc)); |
| 231 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
9146 | PetscCall(RGIsTrivial(eps->rg,&istrivial)); |
| 232 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
9146 | sc->rg = istrivial? NULL: eps->rg; |
| 233 | 9146 | sc->comparison = eps->sc->comparison; | |
| 234 | 9146 | sc->comparisonctx = eps->sc->comparisonctx; | |
| 235 | 9146 | sc->map = SlepcMap_ST; | |
| 236 | 9146 | sc->mapobj = (PetscObject)eps->st; | |
| 237 |
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.
|
9146 | PetscFunctionReturn(PETSC_SUCCESS); |
| 238 | } | ||
| 239 | |||
| 240 | /*@ | ||
| 241 | EPSSetDSType - Sets the type of the internal `DS` object based on the current | ||
| 242 | settings of the eigenvalue solver. | ||
| 243 | |||
| 244 | Collective | ||
| 245 | |||
| 246 | Input Parameter: | ||
| 247 | . eps - the linear eigensolver context | ||
| 248 | |||
| 249 | Note: | ||
| 250 | This function need not be called explicitly, since it will be called at | ||
| 251 | both `EPSSetFromOptions()` and `EPSSetUp()`. | ||
| 252 | |||
| 253 | Level: developer | ||
| 254 | |||
| 255 | .seealso: [](ch:eps), `EPSSetFromOptions()`, `EPSSetUp()` | ||
| 256 | @*/ | ||
| 257 | 15464 | PetscErrorCode EPSSetDSType(EPS eps) | |
| 258 | { | ||
| 259 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
15464 | PetscFunctionBegin; |
| 260 |
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.
|
15464 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 261 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
15464 | PetscTryTypeMethod(eps,setdstype); |
| 262 |
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.
|
2964 | PetscFunctionReturn(PETSC_SUCCESS); |
| 263 | } | ||
| 264 | |||
| 265 | /*@ | ||
| 266 | EPSSetUp - Sets up all the internal data structures necessary for the | ||
| 267 | execution of the eigensolver. Then calls `STSetUp()` for any set-up | ||
| 268 | operations associated to the internal `ST` object. | ||
| 269 | |||
| 270 | Collective | ||
| 271 | |||
| 272 | Input Parameter: | ||
| 273 | . eps - the linear eigensolver context | ||
| 274 | |||
| 275 | Notes: | ||
| 276 | This function need not be called explicitly in most cases, since `EPSSolve()` | ||
| 277 | calls it. It can be useful when one wants to measure the set-up time | ||
| 278 | separately from the solve time. | ||
| 279 | |||
| 280 | Level: developer | ||
| 281 | |||
| 282 | .seealso: [](ch:eps), `EPSCreate()`, `EPSSolve()`, `EPSDestroy()`, `STSetUp()`, `EPSSetInitialSpace()`, `EPSSetDeflationSpace()` | ||
| 283 | @*/ | ||
| 284 | 11849 | PetscErrorCode EPSSetUp(EPS eps) | |
| 285 | { | ||
| 286 | 11849 | Mat A,B; | |
| 287 | 11849 | PetscInt k,nmat; | |
| 288 | 11849 | PetscBool flg; | |
| 289 | 11849 | EPSStoppingCtx ctx; | |
| 290 | |||
| 291 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
11849 | PetscFunctionBegin; |
| 292 |
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.
|
11849 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 293 |
8/14✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
|
11849 | if (eps->state) PetscFunctionReturn(PETSC_SUCCESS); |
| 294 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(PetscLogEventBegin(EPS_SetUp,eps,0,0,0)); |
| 295 | |||
| 296 | /* reset the convergence flag from the previous solves */ | ||
| 297 | 10005 | eps->reason = EPS_CONVERGED_ITERATING; | |
| 298 | |||
| 299 | /* Set default solver type (EPSSetFromOptions was not called) */ | ||
| 300 |
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.
|
10005 | if (!((PetscObject)eps)->type_name) PetscCall(EPSSetType(eps,EPSKRYLOVSCHUR)); |
| 301 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
10005 | if (!eps->st) PetscCall(EPSGetST(eps,&eps->st)); |
| 302 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(EPSSetDefaultST(eps)); |
| 303 | |||
| 304 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(STSetTransform(eps->st,PETSC_TRUE)); |
| 305 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(STSetStructured(eps->st,PETSC_FALSE)); |
| 306 |
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.
|
10005 | if (eps->useds && !eps->ds) PetscCall(EPSGetDS(eps,&eps->ds)); |
| 307 |
6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
10005 | if (eps->useds) PetscCall(EPSSetDSType(eps)); |
| 308 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->twosided) { |
| 309 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
209 | PetscCheck(!eps->ishermitian || (eps->isgeneralized && !eps->ispositive),PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Two-sided methods are not intended for %s problems",SLEPC_STRING_HERMITIAN); |
| 310 | } | ||
| 311 |
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.
|
10005 | if (!eps->rg) PetscCall(EPSGetRG(eps,&eps->rg)); |
| 312 |
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.
|
10005 | if (!((PetscObject)eps->rg)->type_name) PetscCall(RGSetType(eps->rg,RGINTERVAL)); |
| 313 | |||
| 314 | /* Set problem dimensions */ | ||
| 315 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(STGetNumMatrices(eps->st,&nmat)); |
| 316 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
10005 | PetscCheck(nmat,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONGSTATE,"EPSSetOperators must be called first"); |
| 317 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(STMatGetSize(eps->st,&eps->n,NULL)); |
| 318 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(STMatGetLocalSize(eps->st,&eps->nloc,NULL)); |
| 319 | |||
| 320 | /* Set default problem type */ | ||
| 321 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (!eps->problem_type) { |
| 322 |
6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 8 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.
|
457 | if (nmat==1) PetscCall(EPSSetProblemType(eps,EPS_NHEP)); |
| 323 |
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.
|
41 | else PetscCall(EPSSetProblemType(eps,EPS_GNHEP)); |
| 324 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 3 times.
✓ Branch 3 taken 10 times.
|
9548 | } else if (nmat==1 && eps->isgeneralized) { |
| 325 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
3 | PetscCall(PetscInfo(eps,"Eigenproblem set as generalized but no matrix B was provided; reverting to a standard eigenproblem\n")); |
| 326 | 3 | eps->isgeneralized = PETSC_FALSE; | |
| 327 |
1/2✓ Branch 0 taken 3 times.
✗ Branch 1 not taken.
|
6 | eps->problem_type = eps->ishermitian? EPS_HEP: EPS_NHEP; |
| 328 |
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.
|
9545 | } else PetscCheck(nmat==1 || eps->isgeneralized,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_INCOMP,"Inconsistent EPS state: the problem type does not match the number of matrices"); |
| 329 | |||
| 330 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->isstructured) { |
| 331 | /* make sure the user has set the appropriate matrix */ | ||
| 332 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
540 | PetscCall(STGetMatrix(eps->st,0,&A)); |
| 333 |
6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 5 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.
|
540 | if (eps->problem_type==EPS_BSE) PetscCall(SlepcCheckMatStruct(A,SLEPC_MAT_STRUCT_BSE,NULL)); |
| 334 |
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.
|
540 | if (eps->problem_type==EPS_HAMILT) PetscCall(SlepcCheckMatStruct(A,SLEPC_MAT_STRUCT_HAMILT,NULL)); |
| 335 | } | ||
| 336 | |||
| 337 | /* safeguard for small problems */ | ||
| 338 |
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.
|
10005 | if (eps->n == 0) PetscFunctionReturn(PETSC_SUCCESS); |
| 339 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->nev > eps->n) eps->nev = eps->n; |
| 340 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->problem_type == EPS_BSE) { |
| 341 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
505 | if (2*eps->ncv > eps->n) eps->ncv = eps->n/2; |
| 342 | } else { | ||
| 343 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
9500 | if (eps->ncv > eps->n) eps->ncv = eps->n; |
| 344 | } | ||
| 345 | |||
| 346 | /* check some combinations of eps->which */ | ||
| 347 |
7/10✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 10 times.
✓ Branch 5 taken 10 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 10 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
|
10005 | PetscCheck(!eps->ishermitian || (eps->isgeneralized && !eps->ispositive) || (eps->which!=EPS_LARGEST_IMAGINARY && eps->which!=EPS_SMALLEST_IMAGINARY && eps->which!=EPS_TARGET_IMAGINARY),PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Sorting the eigenvalues along the imaginary axis is not allowed when all eigenvalues are real"); |
| 348 | |||
| 349 | /* initialization of matrix norms */ | ||
| 350 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->conv==EPS_CONV_NORM) { |
| 351 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
496 | if (!eps->nrma) { |
| 352 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
496 | PetscCall(STGetMatrix(eps->st,0,&A)); |
| 353 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
496 | PetscCall(MatNorm(A,NORM_INFINITY,&eps->nrma)); |
| 354 | } | ||
| 355 |
3/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✗ Branch 3 not taken.
|
496 | if (nmat>1 && !eps->nrmb) { |
| 356 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
486 | PetscCall(STGetMatrix(eps->st,1,&B)); |
| 357 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
486 | PetscCall(MatNorm(B,NORM_INFINITY,&eps->nrmb)); |
| 358 | } | ||
| 359 | } | ||
| 360 | |||
| 361 | /* call specific solver setup */ | ||
| 362 |
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.
|
10005 | PetscUseTypeMethod(eps,setup); |
| 363 | |||
| 364 | /* threshold stopping test */ | ||
| 365 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->stop==EPS_STOP_THRESHOLD) { |
| 366 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
146 | PetscCheck(eps->thres!=PETSC_MIN_REAL,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"Must give a threshold value with EPSSetThreshold()"); |
| 367 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
146 | PetscCheck(eps->which==EPS_LARGEST_MAGNITUDE || eps->which==EPS_SMALLEST_MAGNITUDE || eps->which==EPS_LARGEST_REAL || eps->which==EPS_SMALLEST_REAL || eps->which==EPS_TARGET_MAGNITUDE,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Threshold stopping test can only be used with largest/smallest/target magnitude or largest/smallest real selection of eigenvalues"); |
| 368 |
1/6✗ 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.
|
146 | if (eps->which==EPS_LARGEST_REAL || eps->which==EPS_SMALLEST_REAL) PetscCheck(eps->problem_type==EPS_HEP || eps->problem_type==EPS_GHEP || eps->problem_type==EPS_BSE,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Threshold stopping test with largest/smallest real can only be used in problems that have all eigenvaues real"); |
| 369 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
146 | else PetscCheck(eps->thres>0.0,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"In case of largest/smallest/target magnitude the threshold value must be positive"); |
| 370 |
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.
|
146 | PetscCheck(eps->which==EPS_LARGEST_MAGNITUDE || eps->which==EPS_TARGET_MAGNITUDE || !eps->threlative,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"Can only use a relative threshold with largest/target magnitude selection of eigenvalues"); |
| 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.
|
146 | PetscCall(PetscNew(&ctx)); |
| 372 | 146 | ctx->thres = eps->thres; | |
| 373 | 146 | ctx->threlative = eps->threlative; | |
| 374 | 146 | ctx->which = eps->which; | |
| 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.
|
146 | PetscCall(EPSSetStoppingTestFunction(eps,EPSStoppingThreshold,ctx,PetscCtxDestroyDefault)); |
| 376 | } | ||
| 377 | |||
| 378 | /* if purification is set, check that it really makes sense */ | ||
| 379 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->purify) { |
| 380 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7126 | if (eps->categ==EPS_CATEGORY_PRECOND || eps->categ==EPS_CATEGORY_CONTOUR) eps->purify = PETSC_FALSE; |
| 381 | else { | ||
| 382 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
5464 | if (!eps->isgeneralized) eps->purify = PETSC_FALSE; |
| 383 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 8 times.
|
1490 | else if (!eps->ishermitian && !eps->ispositive) eps->purify = PETSC_FALSE; |
| 384 | else { | ||
| 385 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1164 | PetscCall(PetscObjectTypeCompare((PetscObject)eps->st,STCAYLEY,&flg)); |
| 386 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1164 | if (flg) eps->purify = PETSC_FALSE; |
| 387 | } | ||
| 388 | } | ||
| 389 | } | ||
| 390 | |||
| 391 | /* set tolerance if not yet set */ | ||
| 392 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->tol==(PetscReal)PETSC_DETERMINE) eps->tol = SLEPC_DEFAULT_TOL; |
| 393 | |||
| 394 | /* set up sorting criterion */ | ||
| 395 |
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.
|
10005 | PetscTryTypeMethod(eps,setupsort); |
| 396 | |||
| 397 | /* Build balancing matrix if required */ | ||
| 398 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
10005 | if (eps->balance!=EPS_BALANCE_USER) { |
| 399 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(STSetBalanceMatrix(eps->st,NULL)); |
| 400 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
10005 | if (!eps->ishermitian && (eps->balance==EPS_BALANCE_ONESIDE || eps->balance==EPS_BALANCE_TWOSIDE)) { |
| 401 |
5/8✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
130 | if (!eps->D) PetscCall(BVCreateVec(eps->V,&eps->D)); |
| 402 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
130 | PetscCall(EPSBuildBalance_Krylov(eps)); |
| 403 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
130 | PetscCall(STSetBalanceMatrix(eps->st,eps->D)); |
| 404 | } | ||
| 405 | } | ||
| 406 | |||
| 407 | /* Setup ST */ | ||
| 408 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(STSetUp(eps->st)); |
| 409 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(EPSCheckCompatibleST(eps)); |
| 410 | |||
| 411 | /* process deflation and initial vectors */ | ||
| 412 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->nds<0) { |
| 413 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
153 | PetscCheck(!eps->isstructured,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Deflation space is not supported in structured eigenproblems"); |
| 414 | 153 | k = -eps->nds; | |
| 415 |
4/6✓ Branch 0 taken 2 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(BVInsertConstraints(eps->V,&k,eps->defl)); |
| 416 |
4/6✓ Branch 0 taken 2 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(SlepcBasisDestroy_Private(&eps->nds,&eps->defl)); |
| 417 | 153 | eps->nds = k; | |
| 418 |
4/6✓ Branch 0 taken 2 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(STCheckNullSpace(eps->st,eps->V)); |
| 419 | } | ||
| 420 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->nini<0) { |
| 421 | 2624 | k = -eps->nini; | |
| 422 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
2624 | PetscCheck(k<=eps->ncv,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"The number of initial vectors is larger than ncv"); |
| 423 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2624 | PetscCall(BVInsertVecs(eps->V,0,&k,eps->IS,PETSC_TRUE)); |
| 424 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2624 | PetscCall(SlepcBasisDestroy_Private(&eps->nini,&eps->IS)); |
| 425 | 2624 | eps->nini = k; | |
| 426 | } | ||
| 427 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
10005 | if (eps->twosided && eps->ninil<0) { |
| 428 | 40 | k = -eps->ninil; | |
| 429 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
40 | PetscCheck(k<=eps->ncv,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"The number of left initial vectors is larger than ncv"); |
| 430 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
40 | PetscCall(BVInsertVecs(eps->W,0,&k,eps->ISL,PETSC_TRUE)); |
| 431 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
40 | PetscCall(SlepcBasisDestroy_Private(&eps->ninil,&eps->ISL)); |
| 432 | 40 | eps->ninil = k; | |
| 433 | } | ||
| 434 | |||
| 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.
|
10005 | PetscCall(PetscLogEventEnd(EPS_SetUp,eps,0,0,0)); |
| 436 | 10005 | eps->state = EPS_STATE_SETUP; | |
| 437 |
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.
|
10005 | PetscFunctionReturn(PETSC_SUCCESS); |
| 438 | } | ||
| 439 | |||
| 440 | /*@ | ||
| 441 | EPSSetOperators - Sets the matrices associated with the eigenvalue problem. | ||
| 442 | |||
| 443 | Collective | ||
| 444 | |||
| 445 | Input Parameters: | ||
| 446 | + eps - the linear eigensolver context | ||
| 447 | . A - the matrix associated with the eigensystem | ||
| 448 | - B - the second matrix in the case of generalized eigenproblems | ||
| 449 | |||
| 450 | Notes: | ||
| 451 | To specify a standard eigenproblem, use `NULL` for parameter `B`. | ||
| 452 | |||
| 453 | It must be called before `EPSSetUp()`. If it is called again after `EPSSetUp()` and | ||
| 454 | the matrix sizes have changed then the `EPS` object is reset. | ||
| 455 | |||
| 456 | For structured eigenproblem types such as `EPS_BSE` (see `EPSSetProblemType()`), the | ||
| 457 | provided matrices must have been created with the corresponding helper function, | ||
| 458 | i.e., `MatCreateBSE()`. | ||
| 459 | |||
| 460 | Level: beginner | ||
| 461 | |||
| 462 | .seealso: [](ch:eps), `EPSSolve()`, `EPSSetUp()`, `EPSReset()`, `EPSSetProblemType()` | ||
| 463 | @*/ | ||
| 464 | 8382 | PetscErrorCode EPSSetOperators(EPS eps,Mat A,Mat B) | |
| 465 | { | ||
| 466 | 8382 | PetscInt m,n,m0,mloc,nloc,mloc0,nmat; | |
| 467 | 8382 | Mat mat[2]; | |
| 468 | 8382 | VecType ta,tb; | |
| 469 | 8382 | PetscBool same; | |
| 470 | |||
| 471 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
8382 | PetscFunctionBegin; |
| 472 |
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.
|
8382 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 473 |
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.
|
8382 | PetscValidHeaderSpecific(A,MAT_CLASSID,2); |
| 474 |
4/14✓ 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.
|
8382 | if (B) PetscValidHeaderSpecific(B,MAT_CLASSID,3); |
| 475 |
13/32✓ 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 not taken.
✓ Branch 11 taken 2 times.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
|
8382 | PetscCheckSameComm(eps,1,A,2); |
| 476 |
15/34✓ Branch 0 taken 2 times.
✓ Branch 1 taken 2 times.
✓ 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 taken 2 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✓ Branch 17 taken 2 times.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 2 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 2 times.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✓ Branch 27 taken 2 times.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✓ Branch 31 taken 2 times.
✗ Branch 32 not taken.
✗ Branch 33 not taken.
|
8382 | if (B) PetscCheckSameComm(eps,1,B,3); |
| 477 | |||
| 478 | /* Check matrix sizes */ | ||
| 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.
|
8382 | PetscCall(MatGetSize(A,&m,&n)); |
| 480 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
8382 | PetscCall(MatGetLocalSize(A,&mloc,&nloc)); |
| 481 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
8382 | PetscCheck(m==n,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONG,"A is a non-square matrix (%" PetscInt_FMT " rows, %" PetscInt_FMT " cols)",m,n); |
| 482 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
8382 | PetscCheck(mloc==nloc,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONG,"A does not have equal row and column sizes (%" PetscInt_FMT ", %" PetscInt_FMT ")",mloc,nloc); |
| 483 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
8382 | if (B) { |
| 484 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2311 | PetscCall(MatGetSize(B,&m0,&n)); |
| 485 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2311 | PetscCall(MatGetLocalSize(B,&mloc0,&nloc)); |
| 486 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
2311 | PetscCheck(m0==n,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONG,"B is a non-square matrix (%" PetscInt_FMT " rows, %" PetscInt_FMT " cols)",m0,n); |
| 487 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
2311 | PetscCheck(mloc0==nloc,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_WRONG,"B does not have equal row and column local sizes (%" PetscInt_FMT ", %" PetscInt_FMT ")",mloc0,nloc); |
| 488 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
2311 | PetscCheck(m==m0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_INCOMP,"Dimensions of A and B do not match (%" PetscInt_FMT ", %" PetscInt_FMT ")",m,m0); |
| 489 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
2311 | PetscCheck(mloc==mloc0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_INCOMP,"Local dimensions of A and B do not match (%" PetscInt_FMT ", %" PetscInt_FMT ")",mloc,mloc0); |
| 490 | /* make sure both matrices have compatible VecType */ | ||
| 491 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2311 | PetscCall(MatGetVecType(A,&ta)); |
| 492 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2311 | PetscCall(MatGetVecType(B,&tb)); |
| 493 |
3/6✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2311 | PetscCall(PetscStrcmp(ta,tb,&same)); |
| 494 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
2311 | PetscCheck(same,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_INCOMP,"The provided matrices have different vector types (%s vs %s), consider calling MatSetVecType() in the one that is not on GPU",ta,tb); |
| 495 | } | ||
| 496 |
10/12✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
✓ Branch 4 taken 10 times.
✓ Branch 5 taken 10 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 8 times.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
8382 | if (eps->state && (n!=eps->n || nloc!=eps->nloc)) PetscCall(EPSReset(eps)); |
| 497 | 8382 | eps->nrma = 0.0; | |
| 498 | 8382 | eps->nrmb = 0.0; | |
| 499 |
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.
|
8382 | if (!eps->st) PetscCall(EPSGetST(eps,&eps->st)); |
| 500 | 8382 | mat[0] = A; | |
| 501 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
8382 | if (B) { |
| 502 | 2311 | mat[1] = B; | |
| 503 | 2311 | nmat = 2; | |
| 504 | } else nmat = 1; | ||
| 505 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
8382 | PetscCall(STSetMatrices(eps->st,nmat,mat)); |
| 506 | 8382 | eps->state = EPS_STATE_INITIAL; | |
| 507 |
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.
|
8382 | PetscFunctionReturn(PETSC_SUCCESS); |
| 508 | } | ||
| 509 | |||
| 510 | /*@ | ||
| 511 | EPSGetOperators - Gets the matrices associated with the eigensystem. | ||
| 512 | |||
| 513 | Collective | ||
| 514 | |||
| 515 | Input Parameter: | ||
| 516 | . eps - the linear eigensolver context | ||
| 517 | |||
| 518 | Output Parameters: | ||
| 519 | + A - the matrix associated with the eigensystem | ||
| 520 | - B - the second matrix in the case of generalized eigenproblems | ||
| 521 | |||
| 522 | Note: | ||
| 523 | Does not increase the reference count of the matrices, so you should not destroy them. | ||
| 524 | |||
| 525 | Level: intermediate | ||
| 526 | |||
| 527 | .seealso: [](ch:eps), `EPSSolve()`, `EPSGetST()`, `STGetMatrix()`, `STSetMatrices()` | ||
| 528 | @*/ | ||
| 529 | 11369 | PetscErrorCode EPSGetOperators(EPS eps,Mat *A,Mat *B) | |
| 530 | { | ||
| 531 | 11369 | ST st; | |
| 532 | 11369 | PetscInt k; | |
| 533 | |||
| 534 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
11369 | PetscFunctionBegin; |
| 535 |
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.
|
11369 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 536 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11369 | PetscCall(EPSGetST(eps,&st)); |
| 537 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
11369 | PetscCall(STGetNumMatrices(st,&k)); |
| 538 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11369 | if (A) { |
| 539 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
|
5881 | if (k<1) *A = NULL; |
| 540 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
5881 | else PetscCall(STGetMatrix(st,0,A)); |
| 541 | } | ||
| 542 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
11369 | if (B) { |
| 543 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7273 | if (k<2) *B = NULL; |
| 544 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6861 | else PetscCall(STGetMatrix(st,1,B)); |
| 545 | } | ||
| 546 |
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.
|
2239 | PetscFunctionReturn(PETSC_SUCCESS); |
| 547 | } | ||
| 548 | |||
| 549 | /*@ | ||
| 550 | EPSSetDeflationSpace - Specify a basis of vectors that constitute the deflation | ||
| 551 | space. | ||
| 552 | |||
| 553 | Collective | ||
| 554 | |||
| 555 | Input Parameters: | ||
| 556 | + eps - the linear eigensolver context | ||
| 557 | . n - number of vectors | ||
| 558 | - v - set of basis vectors of the deflation space | ||
| 559 | |||
| 560 | Notes: | ||
| 561 | When a deflation space is given, the eigensolver seeks the eigensolution | ||
| 562 | in the restriction of the problem to the orthogonal complement of this | ||
| 563 | space. This can be used for instance in the case that an invariant | ||
| 564 | subspace is known beforehand (such as the nullspace of the matrix). | ||
| 565 | |||
| 566 | These vectors do not persist from one `EPSSolve()` call to the other, so the | ||
| 567 | deflation space should be set every time. | ||
| 568 | |||
| 569 | The vectors do not need to be mutually orthonormal, since they are explicitly | ||
| 570 | orthonormalized internally. | ||
| 571 | |||
| 572 | Level: intermediate | ||
| 573 | |||
| 574 | .seealso: [](ch:eps), `EPSSetInitialSpace()` | ||
| 575 | @*/ | ||
| 576 | 171 | PetscErrorCode EPSSetDeflationSpace(EPS eps,PetscInt n,Vec v[]) | |
| 577 | { | ||
| 578 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
171 | PetscFunctionBegin; |
| 579 |
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.
|
171 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 580 |
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.
|
171 | PetscValidLogicalCollectiveInt(eps,n,2); |
| 581 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
171 | PetscCheck(n>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"Argument n cannot be negative"); |
| 582 |
2/2✓ Branch 0 taken 2 times.
✓ Branch 1 taken 2 times.
|
171 | if (n>0) { |
| 583 |
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.
|
167 | PetscAssertPointer(v,3); |
| 584 |
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.
|
167 | PetscValidHeaderSpecific(*v,VEC_CLASSID,3); |
| 585 | } | ||
| 586 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
171 | PetscCall(SlepcBasisReference_Private(n,v,&eps->nds,&eps->defl)); |
| 587 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 8 times.
|
171 | if (n>0) eps->state = EPS_STATE_INITIAL; |
| 588 |
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.
|
33 | PetscFunctionReturn(PETSC_SUCCESS); |
| 589 | } | ||
| 590 | |||
| 591 | /*@ | ||
| 592 | EPSSetInitialSpace - Specify a basis of vectors that constitute the initial | ||
| 593 | space, that is, the subspace from which the solver starts to iterate. | ||
| 594 | |||
| 595 | Collective | ||
| 596 | |||
| 597 | Input Parameters: | ||
| 598 | + eps - the linear eigensolver context | ||
| 599 | . n - number of vectors | ||
| 600 | - is - set of basis vectors of the initial space | ||
| 601 | |||
| 602 | Notes: | ||
| 603 | Some solvers such as `EPSKRYLOVSCHUR` start to iterate on a single vector | ||
| 604 | (initial vector). In that case, only `is[0]` is taken into account and the | ||
| 605 | other vectors are ignored. But other solvers such as `EPSSUBSPACE` are | ||
| 606 | able to make use of the whole initial subspace as an initial guess. | ||
| 607 | |||
| 608 | These vectors do not persist from one `EPSSolve()` call to the other, so the | ||
| 609 | initial space should be set every time. | ||
| 610 | |||
| 611 | The vectors do not need to be mutually orthonormal, since they are explicitly | ||
| 612 | orthonormalized internally. | ||
| 613 | |||
| 614 | Common usage of this function is when the user can provide a rough approximation | ||
| 615 | of the wanted eigenspace. Then, convergence may be faster. | ||
| 616 | |||
| 617 | Level: intermediate | ||
| 618 | |||
| 619 | .seealso: [](ch:eps), `EPSSetLeftInitialSpace()`, `EPSSetDeflationSpace()` | ||
| 620 | @*/ | ||
| 621 | 2646 | PetscErrorCode EPSSetInitialSpace(EPS eps,PetscInt n,Vec is[]) | |
| 622 | { | ||
| 623 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
2646 | PetscFunctionBegin; |
| 624 |
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.
|
2646 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 625 |
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.
|
2646 | PetscValidLogicalCollectiveInt(eps,n,2); |
| 626 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
2646 | PetscCheck(n>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"Argument n cannot be negative"); |
| 627 |
2/2✓ Branch 0 taken 2 times.
✓ Branch 1 taken 2 times.
|
2646 | if (n>0) { |
| 628 |
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.
|
2638 | PetscAssertPointer(is,3); |
| 629 |
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.
|
2638 | PetscValidHeaderSpecific(*is,VEC_CLASSID,3); |
| 630 | } | ||
| 631 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2646 | PetscCall(SlepcBasisReference_Private(n,is,&eps->nini,&eps->IS)); |
| 632 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 8 times.
|
2646 | if (n>0) eps->state = EPS_STATE_INITIAL; |
| 633 |
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.
|
501 | PetscFunctionReturn(PETSC_SUCCESS); |
| 634 | } | ||
| 635 | |||
| 636 | /*@ | ||
| 637 | EPSSetLeftInitialSpace - Specify a basis of vectors that constitute the left | ||
| 638 | initial space, used by two-sided solvers to start the left subspace. | ||
| 639 | |||
| 640 | Collective | ||
| 641 | |||
| 642 | Input Parameters: | ||
| 643 | + eps - the linear eigensolver context | ||
| 644 | . n - number of vectors | ||
| 645 | - isl - set of basis vectors of the left initial space | ||
| 646 | |||
| 647 | Notes: | ||
| 648 | Left initial vectors are used to initiate the left search space in two-sided | ||
| 649 | eigensolvers. Users should pass here an approximation of the left eigenspace, | ||
| 650 | if available. | ||
| 651 | |||
| 652 | The same comments in `EPSSetInitialSpace()` are applicable here. | ||
| 653 | |||
| 654 | Level: intermediate | ||
| 655 | |||
| 656 | .seealso: [](ch:eps), `EPSSetInitialSpace()`, `EPSSetTwoSided()` | ||
| 657 | @*/ | ||
| 658 | 40 | PetscErrorCode EPSSetLeftInitialSpace(EPS eps,PetscInt n,Vec isl[]) | |
| 659 | { | ||
| 660 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
40 | PetscFunctionBegin; |
| 661 |
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.
|
40 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 662 |
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.
|
40 | PetscValidLogicalCollectiveInt(eps,n,2); |
| 663 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
40 | PetscCheck(n>=0,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"Argument n cannot be negative"); |
| 664 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
40 | if (n>0) { |
| 665 |
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.
|
40 | PetscAssertPointer(isl,3); |
| 666 |
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.
|
40 | PetscValidHeaderSpecific(*isl,VEC_CLASSID,3); |
| 667 | } | ||
| 668 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
40 | PetscCall(SlepcBasisReference_Private(n,isl,&eps->ninil,&eps->ISL)); |
| 669 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
40 | if (n>0) eps->state = EPS_STATE_INITIAL; |
| 670 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
8 | PetscFunctionReturn(PETSC_SUCCESS); |
| 671 | } | ||
| 672 | |||
| 673 | /* | ||
| 674 | EPSSetDimensions_Default - Set reasonable values for ncv, mpd if not set | ||
| 675 | by the user. This is called at setup. | ||
| 676 | */ | ||
| 677 | 7224 | PetscErrorCode EPSSetDimensions_Default(EPS eps,PetscInt *nev,PetscInt *ncv,PetscInt *mpd) | |
| 678 | { | ||
| 679 | 7224 | PetscBool krylov; | |
| 680 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7224 | PetscInt nev2, n = eps->problem_type==EPS_BSE? eps->n/2: eps->n; |
| 681 | |||
| 682 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
7224 | PetscFunctionBegin; |
| 683 |
4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
7224 | if (*nev==0 && eps->stop!=EPS_STOP_THRESHOLD) *nev = 1; |
| 684 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7224 | nev2 = eps->problem_type==EPS_BSE? (*nev+1)/2: *nev; |
| 685 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7224 | if (*ncv!=PETSC_DETERMINE) { /* ncv set */ |
| 686 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
3952 | PetscCall(PetscObjectTypeCompareAny((PetscObject)eps,&krylov,EPSKRYLOVSCHUR,EPSARNOLDI,EPSLANCZOS,"")); |
| 687 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
3952 | if (krylov) { |
| 688 |
4/8✓ Branch 0 taken 1 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
3774 | PetscCheck(*ncv>=nev2+1 || (*ncv==nev2 && *ncv==n),PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"The value of ncv must be at least nev+1"); |
| 689 | } else { | ||
| 690 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
178 | PetscCheck(*ncv>=nev2,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"The value of ncv must be at least nev"); |
| 691 | } | ||
| 692 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
3272 | } else if (*mpd!=PETSC_DETERMINE) { /* mpd set */ |
| 693 | 60 | *ncv = PetscMin(n,nev2+(*mpd)); | |
| 694 | } else { /* neither set: defaults depend on nev being small or large */ | ||
| 695 |
3/4✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.
|
3212 | if (nev2<500) *ncv = PetscMin(n,PetscMax(2*(nev2),nev2+15)); |
| 696 | else { | ||
| 697 | ✗ | *mpd = 500; | |
| 698 | ✗ | *ncv = PetscMin(n,nev2+(*mpd)); | |
| 699 | } | ||
| 700 | } | ||
| 701 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7224 | if (*mpd==PETSC_DETERMINE) *mpd = *ncv; |
| 702 |
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.
|
1360 | PetscFunctionReturn(PETSC_SUCCESS); |
| 703 | } | ||
| 704 | |||
| 705 | /*@ | ||
| 706 | EPSAllocateSolution - Allocate memory storage for common variables such | ||
| 707 | as eigenvalues and eigenvectors. | ||
| 708 | |||
| 709 | Collective | ||
| 710 | |||
| 711 | Input Parameters: | ||
| 712 | + eps - the linear eigensolver context | ||
| 713 | - extra - number of additional positions, used for methods that require a | ||
| 714 | working basis slightly larger than `ncv` | ||
| 715 | |||
| 716 | Developer Note: | ||
| 717 | This is `SLEPC_EXTERN` because it may be required by user plugin `EPS` | ||
| 718 | implementations. | ||
| 719 | |||
| 720 | Level: developer | ||
| 721 | |||
| 722 | .seealso: [](ch:eps), `EPSSetUp()`, `EPSSetDimensions()` | ||
| 723 | @*/ | ||
| 724 | 10005 | PetscErrorCode EPSAllocateSolution(EPS eps,PetscInt extra) | |
| 725 | { | ||
| 726 | 10005 | PetscInt oldsize,requested; | |
| 727 | 10005 | PetscRandom rand; | |
| 728 | 10005 | Vec t; | |
| 729 | |||
| 730 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
10005 | PetscFunctionBegin; |
| 731 | 10005 | requested = eps->ncv + extra; | |
| 732 | |||
| 733 | /* oldsize is zero if this is the first time setup is called */ | ||
| 734 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10005 | PetscCall(BVGetSizes(eps->V,NULL,NULL,&oldsize)); |
| 735 | |||
| 736 | /* allocate space for eigenvalues and friends */ | ||
| 737 |
3/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
|
10005 | if (requested != oldsize || !eps->eigr) { |
| 738 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7499 | PetscCall(PetscFree4(eps->eigr,eps->eigi,eps->errest,eps->perm)); |
| 739 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7499 | PetscCall(PetscMalloc4(requested,&eps->eigr,requested,&eps->eigi,requested,&eps->errest,requested,&eps->perm)); |
| 740 | } | ||
| 741 | |||
| 742 | /* workspace for the case of arbitrary selection */ | ||
| 743 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->arbitrary) { |
| 744 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
82 | if (eps->rr) PetscCall(PetscFree2(eps->rr,eps->ri)); |
| 745 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
82 | PetscCall(PetscMalloc2(requested,&eps->rr,requested,&eps->ri)); |
| 746 | } | ||
| 747 | |||
| 748 | /* allocate V */ | ||
| 749 |
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.
|
10005 | if (!eps->V) PetscCall(EPSGetBV(eps,&eps->V)); |
| 750 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (!oldsize) { |
| 751 |
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.
|
7358 | if (!((PetscObject)eps->V)->type_name) PetscCall(BVSetType(eps->V,BVMAT)); |
| 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.
|
7358 | PetscCall(STMatCreateVecsEmpty(eps->st,&t,NULL)); |
| 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.
|
7358 | PetscCall(BVSetSizesFromVec(eps->V,t,requested)); |
| 754 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
7358 | PetscCall(VecDestroy(&t)); |
| 755 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2647 | } else PetscCall(BVResize(eps->V,requested,PETSC_FALSE)); |
| 756 | |||
| 757 | /* allocate W */ | ||
| 758 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
10005 | if (eps->twosided) { |
| 759 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
209 | PetscCall(BVGetRandomContext(eps->V,&rand)); /* make sure the random context is available when duplicating */ |
| 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.
|
209 | PetscCall(BVDestroy(&eps->W)); |
| 761 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
209 | PetscCall(BVDuplicate(eps->V,&eps->W)); |
| 762 | } | ||
| 763 |
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.
|
1926 | PetscFunctionReturn(PETSC_SUCCESS); |
| 764 | } | ||
| 765 | |||
| 766 | /*@ | ||
| 767 | EPSReallocateSolution - Reallocate memory storage for common variables such | ||
| 768 | as the eigenvalues and the basis vectors. | ||
| 769 | |||
| 770 | Collective | ||
| 771 | |||
| 772 | Input Parameters: | ||
| 773 | + eps - the linear eigensolver context | ||
| 774 | - newsize - new size | ||
| 775 | |||
| 776 | Developer Notes: | ||
| 777 | This is `SLEPC_EXTERN` because it may be required by user plugin `EPS` | ||
| 778 | implementations. | ||
| 779 | |||
| 780 | This is called during the iteration in case the threshold stopping test has | ||
| 781 | been selected. | ||
| 782 | |||
| 783 | Level: developer | ||
| 784 | |||
| 785 | .seealso: [](ch:eps), `EPSAllocateSolution()`, `EPSSetThreshold()` | ||
| 786 | @*/ | ||
| 787 | 73 | PetscErrorCode EPSReallocateSolution(EPS eps,PetscInt newsize) | |
| 788 | { | ||
| 789 | 73 | PetscInt oldsize,*nperm; | |
| 790 | 73 | PetscReal *nerrest; | |
| 791 | 73 | PetscScalar *neigr,*neigi; | |
| 792 | |||
| 793 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
73 | PetscFunctionBegin; |
| 794 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(BVGetSizes(eps->V,NULL,NULL,&oldsize)); |
| 795 |
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.
|
73 | if (oldsize>=newsize) PetscFunctionReturn(PETSC_SUCCESS); |
| 796 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(PetscInfo(eps,"Reallocating basis vectors to %" PetscInt_FMT " columns\n",newsize)); |
| 797 | |||
| 798 | /* reallocate eigenvalues */ | ||
| 799 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(PetscMalloc4(newsize,&neigr,newsize,&neigi,newsize,&nerrest,newsize,&nperm)); |
| 800 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(PetscArraycpy(neigr,eps->eigr,oldsize)); |
| 801 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(PetscArraycpy(neigi,eps->eigi,oldsize)); |
| 802 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(PetscArraycpy(nerrest,eps->errest,oldsize)); |
| 803 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(PetscArraycpy(nperm,eps->perm,oldsize)); |
| 804 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(PetscFree4(eps->eigr,eps->eigi,eps->errest,eps->perm)); |
| 805 | 73 | eps->eigr = neigr; | |
| 806 | 73 | eps->eigi = neigi; | |
| 807 | 73 | eps->errest = nerrest; | |
| 808 | 73 | eps->perm = nperm; | |
| 809 | /* reallocate V,W */ | ||
| 810 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
73 | PetscCall(BVResize(eps->V,newsize,PETSC_TRUE)); |
| 811 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
73 | if (eps->twosided) PetscCall(BVResize(eps->W,newsize,PETSC_TRUE)); |
| 812 |
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.
|
15 | PetscFunctionReturn(PETSC_SUCCESS); |
| 813 | } | ||
| 814 |