| 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 | This file implements a wrapper to the BLOPEX package | ||
| 12 | */ | ||
| 13 | |||
| 14 | #include <slepc/private/epsimpl.h> /*I "slepceps.h" I*/ | ||
| 15 | #include "blopex.h" | ||
| 16 | #include <lobpcg.h> | ||
| 17 | #include <interpreter.h> | ||
| 18 | #include <multivector.h> | ||
| 19 | #include <temp_multivector.h> | ||
| 20 | |||
| 21 | PetscInt slepc_blopex_useconstr = -1; | ||
| 22 | |||
| 23 | typedef struct { | ||
| 24 | lobpcg_Tolerance tol; | ||
| 25 | lobpcg_BLASLAPACKFunctions blap_fn; | ||
| 26 | mv_InterfaceInterpreter ii; | ||
| 27 | ST st; | ||
| 28 | Vec w; | ||
| 29 | PetscInt bs; /* block size */ | ||
| 30 | } EPS_BLOPEX; | ||
| 31 | |||
| 32 | 4708 | static void Precond_FnSingleVector(void *data,void *x,void *y) | |
| 33 | { | ||
| 34 | 4708 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)data; | |
| 35 | 4708 | MPI_Comm comm = PetscObjectComm((PetscObject)blopex->st); | |
| 36 | 4708 | KSP ksp; | |
| 37 | |||
| 38 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
4708 | PetscFunctionBegin; |
| 39 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
4708 | PetscCallAbort(comm,STGetKSP(blopex->st,&ksp)); |
| 40 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
4708 | PetscCallAbort(comm,KSPSolve(ksp,(Vec)x,(Vec)y)); |
| 41 |
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.
|
4708 | PetscFunctionReturnVoid(); |
| 42 | } | ||
| 43 | |||
| 44 | 1579 | static void Precond_FnMultiVector(void *data,void *x,void *y) | |
| 45 | { | ||
| 46 | 1579 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)data; | |
| 47 | |||
| 48 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1579 | PetscFunctionBegin; |
| 49 | 1579 | blopex->ii.Eval(Precond_FnSingleVector,data,x,y); | |
| 50 |
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.
|
1579 | PetscFunctionReturnVoid(); |
| 51 | } | ||
| 52 | |||
| 53 | 4918 | static void OperatorASingleVector(void *data,void *x,void *y) | |
| 54 | { | ||
| 55 | 4918 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)data; | |
| 56 | 4918 | MPI_Comm comm = PetscObjectComm((PetscObject)blopex->st); | |
| 57 | 4918 | Mat A,B; | |
| 58 | 4918 | PetscScalar sigma; | |
| 59 | 4918 | PetscInt nmat; | |
| 60 | |||
| 61 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
4918 | PetscFunctionBegin; |
| 62 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
4918 | PetscCallAbort(comm,STGetNumMatrices(blopex->st,&nmat)); |
| 63 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
4918 | PetscCallAbort(comm,STGetMatrix(blopex->st,0,&A)); |
| 64 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
4918 | if (nmat>1) PetscCallAbort(comm,STGetMatrix(blopex->st,1,&B)); |
| 65 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
4918 | PetscCallAbort(comm,MatMult(A,(Vec)x,(Vec)y)); |
| 66 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
4918 | PetscCallAbort(comm,STGetShift(blopex->st,&sigma)); |
| 67 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
4918 | if (sigma != 0.0) { |
| 68 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
1068 | if (nmat>1) PetscCallAbort(comm,MatMult(B,(Vec)x,blopex->w)); |
| 69 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
825 | else PetscCallAbort(comm,VecCopy((Vec)x,blopex->w)); |
| 70 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
1068 | PetscCallAbort(comm,VecAXPY((Vec)y,-sigma,blopex->w)); |
| 71 | } | ||
| 72 |
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.
|
4918 | PetscFunctionReturnVoid(); |
| 73 | } | ||
| 74 | |||
| 75 | 1633 | static void OperatorAMultiVector(void *data,void *x,void *y) | |
| 76 | { | ||
| 77 | 1633 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)data; | |
| 78 | |||
| 79 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1633 | PetscFunctionBegin; |
| 80 | 1633 | blopex->ii.Eval(OperatorASingleVector,data,x,y); | |
| 81 |
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.
|
1633 | PetscFunctionReturnVoid(); |
| 82 | } | ||
| 83 | |||
| 84 | 2796 | static void OperatorBSingleVector(void *data,void *x,void *y) | |
| 85 | { | ||
| 86 | 2796 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)data; | |
| 87 | 2796 | MPI_Comm comm = PetscObjectComm((PetscObject)blopex->st); | |
| 88 | 2796 | Mat B; | |
| 89 | |||
| 90 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
2796 | PetscFunctionBegin; |
| 91 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2796 | PetscCallAbort(comm,STGetMatrix(blopex->st,1,&B)); |
| 92 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
2796 | PetscCallAbort(comm,MatMult(B,(Vec)x,(Vec)y)); |
| 93 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
2796 | PetscFunctionReturnVoid(); |
| 94 | } | ||
| 95 | |||
| 96 | 969 | static void OperatorBMultiVector(void *data,void *x,void *y) | |
| 97 | { | ||
| 98 | 969 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)data; | |
| 99 | |||
| 100 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
969 | PetscFunctionBegin; |
| 101 | 969 | blopex->ii.Eval(OperatorBSingleVector,data,x,y); | |
| 102 |
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.
|
969 | PetscFunctionReturnVoid(); |
| 103 | } | ||
| 104 | |||
| 105 | 42 | static PetscErrorCode EPSSetDimensions_BLOPEX(EPS eps,PetscInt *nev,PetscInt *ncv,PetscInt *mpd) | |
| 106 | { | ||
| 107 | 42 | EPS_BLOPEX *ctx = (EPS_BLOPEX*)eps->data; | |
| 108 | 42 | PetscInt k; | |
| 109 | |||
| 110 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
42 | PetscFunctionBegin; |
| 111 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
|
42 | if (*nev==0) *nev = 1; |
| 112 | 42 | k = ((*nev-1)/ctx->bs+1)*ctx->bs; | |
| 113 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
42 | if (*ncv!=PETSC_DETERMINE) { /* ncv set */ |
| 114 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
6 | PetscCheck(*ncv>=k,PetscObjectComm((PetscObject)eps),PETSC_ERR_USER_INPUT,"The value of ncv is not sufficiently large"); |
| 115 | 36 | } else *ncv = k; | |
| 116 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
42 | if (*mpd==PETSC_DETERMINE) *mpd = *ncv; |
| 117 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6 | else PetscCall(PetscInfo(eps,"Warning: given value of mpd ignored\n")); |
| 118 |
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.
|
14 | PetscFunctionReturn(PETSC_SUCCESS); |
| 119 | } | ||
| 120 | |||
| 121 | 42 | static PetscErrorCode EPSSetUp_BLOPEX(EPS eps) | |
| 122 | { | ||
| 123 | 42 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)eps->data; | |
| 124 | 42 | PetscBool flg; | |
| 125 | 42 | KSP ksp; | |
| 126 | |||
| 127 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
42 | PetscFunctionBegin; |
| 128 |
4/10✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 6 times.
✓ Branch 5 taken 6 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 6 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
|
42 | EPSCheckHermitianDefinite(eps); |
| 129 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
42 | EPSCheckNotStructured(eps); |
| 130 |
4/6✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 6 times.
✗ Branch 5 not taken.
|
66 | if (!blopex->bs) blopex->bs = PetscMin(16,eps->nev?eps->nev:1); |
| 131 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42 | PetscCall(EPSSetDimensions_BLOPEX(eps,&eps->nev,&eps->ncv,&eps->mpd)); |
| 132 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
42 | if (eps->max_it==PETSC_DETERMINE) eps->max_it = PetscMax(100,2*eps->n/eps->ncv); |
| 133 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
42 | if (!eps->which) eps->which = EPS_SMALLEST_REAL; |
| 134 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
42 | PetscCheck(eps->which==EPS_SMALLEST_REAL,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"This solver supports only smallest real eigenvalues"); |
| 135 |
8/18✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 4 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 4 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 6 times.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
|
42 | EPSCheckUnsupported(eps,EPS_FEATURE_ARBITRARY | EPS_FEATURE_REGION | EPS_FEATURE_STOPPING); |
| 136 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 4 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.
|
42 | EPSCheckIgnored(eps,EPS_FEATURE_BALANCE | EPS_FEATURE_EXTRACTION); |
| 137 | |||
| 138 | 42 | blopex->st = eps->st; | |
| 139 | |||
| 140 |
3/6✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 6 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
|
42 | PetscCheck(eps->converged==EPSConvergedRelative || eps->converged==EPSConvergedAbsolute,PetscObjectComm((PetscObject)eps),PETSC_ERR_SUP,"Convergence test not supported in this solver"); |
| 141 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
42 | if (eps->converged == EPSConvergedRelative) { |
| 142 | 24 | blopex->tol.absolute = 0.0; | |
| 143 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
30 | blopex->tol.relative = SlepcDefaultTol(eps->tol); |
| 144 | } else { /* EPSConvergedAbsolute */ | ||
| 145 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
|
18 | blopex->tol.absolute = SlepcDefaultTol(eps->tol); |
| 146 | 18 | blopex->tol.relative = 0.0; | |
| 147 | } | ||
| 148 | |||
| 149 | 42 | SLEPCSetupInterpreter(&blopex->ii); | |
| 150 | |||
| 151 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42 | PetscCall(STGetKSP(eps->st,&ksp)); |
| 152 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42 | PetscCall(PetscObjectTypeCompare((PetscObject)ksp,KSPPREONLY,&flg)); |
| 153 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
42 | if (!flg) PetscCall(PetscInfo(eps,"Warning: ignoring KSP, should use KSPPREONLY\n")); |
| 154 | |||
| 155 | /* allocate memory */ | ||
| 156 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 6 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.
|
42 | if (!eps->V) PetscCall(EPSGetBV(eps,&eps->V)); |
| 157 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42 | PetscCall(PetscObjectTypeCompareAny((PetscObject)eps->V,&flg,BVVECS,BVCONTIGUOUS,"")); |
| 158 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
42 | if (!flg) { /* blopex only works with BVVECS or BVCONTIGUOUS */ |
| 159 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(BVSetType(eps->V,BVCONTIGUOUS)); |
| 160 | } | ||
| 161 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42 | PetscCall(EPSAllocateSolution(eps,0)); |
| 162 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
42 | if (!blopex->w) PetscCall(BVCreateVec(eps->V,&blopex->w)); |
| 163 | |||
| 164 | #if PetscDefined(USE_COMPLEX) | ||
| 165 | 21 | blopex->blap_fn.zpotrf = PETSC_zpotrf_interface; | |
| 166 | 21 | blopex->blap_fn.zhegv = PETSC_zsygv_interface; | |
| 167 | #else | ||
| 168 | 21 | blopex->blap_fn.dpotrf = PETSC_dpotrf_interface; | |
| 169 | 21 | blopex->blap_fn.dsygv = PETSC_dsygv_interface; | |
| 170 | #endif | ||
| 171 |
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.
|
42 | PetscFunctionReturn(PETSC_SUCCESS); |
| 172 | } | ||
| 173 | |||
| 174 | 42 | static PetscErrorCode EPSSolve_BLOPEX(EPS eps) | |
| 175 | { | ||
| 176 | 42 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)eps->data; | |
| 177 | 42 | PetscScalar sigma,*eigr=NULL; | |
| 178 | 42 | PetscReal *errest=NULL; | |
| 179 | 42 | int i,j,info,its,nconv; | |
| 180 | 42 | double *residhist=NULL; | |
| 181 | 42 | mv_MultiVectorPtr eigenvectors,constraints; | |
| 182 | #if PetscDefined(USE_COMPLEX) | ||
| 183 | 21 | komplex *lambda=NULL,*lambdahist=NULL; | |
| 184 | #else | ||
| 185 | 21 | double *lambda=NULL,*lambdahist=NULL; | |
| 186 | #endif | ||
| 187 | |||
| 188 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
42 | PetscFunctionBegin; |
| 189 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42 | PetscCall(STGetShift(eps->st,&sigma)); |
| 190 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
42 | PetscCall(PetscMalloc1(blopex->bs,&lambda)); |
| 191 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 6 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.
|
42 | if (eps->numbermonitors>0) PetscCall(PetscMalloc4(blopex->bs*(eps->max_it+1),&lambdahist,eps->ncv,&eigr,blopex->bs*(eps->max_it+1),&residhist,eps->ncv,&errest)); |
| 192 | |||
| 193 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
|
42 | if (eps->nini<eps->ncv) { |
| 194 | /* Complete the initial basis with random vectors */ | ||
| 195 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
48 | for (i=0;i<eps->nini;i++) { /* in case the initial vectors were also set with VecSetRandom */ |
| 196 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6 | PetscCall(BVSetRandomColumn(eps->V,eps->nini)); |
| 197 | } | ||
| 198 |
7/8✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 6 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
|
246 | for (i=eps->nini;i<eps->ncv;i++) PetscCall(BVSetRandomColumn(eps->V,i)); |
| 199 | } | ||
| 200 | |||
| 201 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
96 | while (eps->reason == EPS_CONVERGED_ITERATING) { |
| 202 | |||
| 203 | /* Create multivector of constraints from leading columns of V */ | ||
| 204 |
7/10✓ Branch 0 taken 6 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 6 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
54 | PetscCall(PetscObjectComposedDataSetInt((PetscObject)eps->V,slepc_blopex_useconstr,1)); |
| 205 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
54 | PetscCall(BVSetActiveColumns(eps->V,0,eps->nconv)); |
| 206 | 54 | constraints = mv_MultiVectorCreateFromSampleVector(&blopex->ii,eps->nds+eps->nconv,eps->V); | |
| 207 | |||
| 208 | /* Create multivector where eigenvectors of this run will be stored */ | ||
| 209 |
5/10✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 6 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
|
54 | PetscCall(PetscObjectComposedDataSetInt((PetscObject)eps->V,slepc_blopex_useconstr,0)); |
| 210 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
54 | PetscCall(BVSetActiveColumns(eps->V,eps->nconv,eps->nconv+blopex->bs)); |
| 211 | 54 | eigenvectors = mv_MultiVectorCreateFromSampleVector(&blopex->ii,blopex->bs,eps->V); | |
| 212 | |||
| 213 | #if PetscDefined(USE_COMPLEX) | ||
| 214 |
2/2✓ Branch 0 taken 3 times.
✓ Branch 1 taken 3 times.
|
45 | info = lobpcg_solve_complex(eigenvectors,blopex,OperatorAMultiVector, |
| 215 | 27 | eps->isgeneralized?blopex:NULL,eps->isgeneralized?OperatorBMultiVector:NULL, | |
| 216 | blopex,Precond_FnMultiVector,constraints, | ||
| 217 | blopex->blap_fn,blopex->tol,eps->max_it,0,&its, | ||
| 218 |
2/2✓ Branch 0 taken 3 times.
✓ Branch 1 taken 3 times.
|
27 | lambda,lambdahist,blopex->bs,eps->errest+eps->nconv,residhist,blopex->bs); |
| 219 | #else | ||
| 220 |
2/2✓ Branch 0 taken 3 times.
✓ Branch 1 taken 3 times.
|
45 | info = lobpcg_solve_double(eigenvectors,blopex,OperatorAMultiVector, |
| 221 | 27 | eps->isgeneralized?blopex:NULL,eps->isgeneralized?OperatorBMultiVector:NULL, | |
| 222 | blopex,Precond_FnMultiVector,constraints, | ||
| 223 | blopex->blap_fn,blopex->tol,eps->max_it,0,&its, | ||
| 224 |
2/2✓ Branch 0 taken 3 times.
✓ Branch 1 taken 3 times.
|
27 | lambda,lambdahist,blopex->bs,eps->errest+eps->nconv,residhist,blopex->bs); |
| 225 | #endif | ||
| 226 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
54 | PetscCheck(info==0,PetscObjectComm((PetscObject)eps),PETSC_ERR_LIB,"BLOPEX failed with exit code=%d",info); |
| 227 | 54 | mv_MultiVectorDestroy(constraints); | |
| 228 | 54 | mv_MultiVectorDestroy(eigenvectors); | |
| 229 | |||
| 230 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
318 | for (j=0;j<blopex->bs;j++) { |
| 231 | #if PetscDefined(USE_COMPLEX) | ||
| 232 | 105 | eps->eigr[eps->nconv+j] = PetscCMPLX(lambda[j].real,lambda[j].imag); | |
| 233 | #else | ||
| 234 | 105 | eps->eigr[eps->nconv+j] = lambda[j]; | |
| 235 | #endif | ||
| 236 | } | ||
| 237 | |||
| 238 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
|
54 | if (eps->numbermonitors>0) { |
| 239 | ✗ | for (i=0;i<its;i++) { | |
| 240 | nconv = 0; | ||
| 241 | ✗ | for (j=0;j<blopex->bs;j++) { | |
| 242 | #if PetscDefined(USE_COMPLEX) | ||
| 243 | ✗ | eigr[eps->nconv+j] = PetscCMPLX(lambdahist[j+i*blopex->bs].real,lambdahist[j+i*blopex->bs].imag); | |
| 244 | #else | ||
| 245 | ✗ | eigr[eps->nconv+j] = lambdahist[j+i*blopex->bs]; | |
| 246 | #endif | ||
| 247 | ✗ | errest[eps->nconv+j] = residhist[j+i*blopex->bs]; | |
| 248 | ✗ | if (residhist[j+i*blopex->bs]<=eps->tol) nconv++; | |
| 249 | } | ||
| 250 | ✗ | PetscCall(EPSMonitor(eps,eps->its+i,eps->nconv+nconv,eigr,eps->eigi,errest,eps->nconv+blopex->bs)); | |
| 251 | } | ||
| 252 | } | ||
| 253 | |||
| 254 | 54 | eps->its += its; | |
| 255 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
|
54 | if (info==-1) { |
| 256 | ✗ | eps->reason = EPS_DIVERGED_ITS; | |
| 257 | ✗ | break; | |
| 258 | } else { | ||
| 259 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
264 | for (i=0;i<blopex->bs;i++) { |
| 260 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
210 | if (sigma != 0.0) eps->eigr[eps->nconv+i] += sigma; |
| 261 | } | ||
| 262 | 54 | eps->nconv += blopex->bs; | |
| 263 |
2/2✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
|
54 | if (eps->nconv>=eps->nev) eps->reason = EPS_CONVERGED_TOL; |
| 264 | } | ||
| 265 | } | ||
| 266 | |||
| 267 |
5/8✓ Branch 0 taken 6 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
42 | PetscCall(PetscFree(lambda)); |
| 268 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 6 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.
|
42 | if (eps->numbermonitors>0) PetscCall(PetscFree4(lambdahist,eigr,residhist,errest)); |
| 269 |
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.
|
14 | PetscFunctionReturn(PETSC_SUCCESS); |
| 270 | } | ||
| 271 | |||
| 272 | 12 | static PetscErrorCode EPSBLOPEXSetBlockSize_BLOPEX(EPS eps,PetscInt bs) | |
| 273 | { | ||
| 274 | 12 | EPS_BLOPEX *ctx = (EPS_BLOPEX*)eps->data; | |
| 275 | |||
| 276 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
12 | PetscFunctionBegin; |
| 277 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
|
12 | if (bs==PETSC_DEFAULT || bs==PETSC_DECIDE) { |
| 278 | ✗ | ctx->bs = 0; | |
| 279 | ✗ | eps->state = EPS_STATE_INITIAL; | |
| 280 | } else { | ||
| 281 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 6 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
12 | PetscCheck(bs>0,PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Block size must be >0"); |
| 282 | 12 | ctx->bs = bs; | |
| 283 | } | ||
| 284 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
4 | PetscFunctionReturn(PETSC_SUCCESS); |
| 285 | } | ||
| 286 | |||
| 287 | /*@ | ||
| 288 | EPSBLOPEXSetBlockSize - Sets the block size of the BLOPEX solver. | ||
| 289 | |||
| 290 | Logically Collective | ||
| 291 | |||
| 292 | Input Parameters: | ||
| 293 | + eps - the linear eigensolver context | ||
| 294 | - bs - the block size | ||
| 295 | |||
| 296 | Options Database Key: | ||
| 297 | . -eps_blopex_blocksize bs - sets the block size | ||
| 298 | |||
| 299 | Level: advanced | ||
| 300 | |||
| 301 | .seealso: [](ch:eps), `EPSBLOPEX`, `EPSBLOPEXGetBlockSize()` | ||
| 302 | @*/ | ||
| 303 | 12 | PetscErrorCode EPSBLOPEXSetBlockSize(EPS eps,PetscInt bs) | |
| 304 | { | ||
| 305 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
12 | PetscFunctionBegin; |
| 306 |
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.
|
12 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 307 |
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.
|
12 | PetscValidLogicalCollectiveInt(eps,bs,2); |
| 308 |
8/14✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 6 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 6 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.
|
12 | PetscTryMethod(eps,"EPSBLOPEXSetBlockSize_C",(EPS,PetscInt),(eps,bs)); |
| 309 |
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.
|
12 | PetscFunctionReturn(PETSC_SUCCESS); |
| 310 | } | ||
| 311 | |||
| 312 | 6 | static PetscErrorCode EPSBLOPEXGetBlockSize_BLOPEX(EPS eps,PetscInt *bs) | |
| 313 | { | ||
| 314 | 6 | EPS_BLOPEX *ctx = (EPS_BLOPEX*)eps->data; | |
| 315 | |||
| 316 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
6 | PetscFunctionBegin; |
| 317 | 6 | *bs = ctx->bs; | |
| 318 |
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.
|
6 | PetscFunctionReturn(PETSC_SUCCESS); |
| 319 | } | ||
| 320 | |||
| 321 | /*@ | ||
| 322 | EPSBLOPEXGetBlockSize - Gets the block size used in the BLOPEX solver. | ||
| 323 | |||
| 324 | Not Collective | ||
| 325 | |||
| 326 | Input Parameter: | ||
| 327 | . eps - the linear eigensolver context | ||
| 328 | |||
| 329 | Output Parameter: | ||
| 330 | . bs - the block size | ||
| 331 | |||
| 332 | Level: advanced | ||
| 333 | |||
| 334 | .seealso: [](ch:eps), `EPSBLOPEX`, `EPSBLOPEXSetBlockSize()` | ||
| 335 | @*/ | ||
| 336 | 6 | PetscErrorCode EPSBLOPEXGetBlockSize(EPS eps,PetscInt *bs) | |
| 337 | { | ||
| 338 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
6 | PetscFunctionBegin; |
| 339 |
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.
|
6 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 340 |
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.
|
6 | PetscAssertPointer(bs,2); |
| 341 |
9/16✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 6 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
|
6 | PetscUseMethod(eps,"EPSBLOPEXGetBlockSize_C",(EPS,PetscInt*),(eps,bs)); |
| 342 |
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.
|
6 | PetscFunctionReturn(PETSC_SUCCESS); |
| 343 | } | ||
| 344 | |||
| 345 | 36 | static PetscErrorCode EPSReset_BLOPEX(EPS eps) | |
| 346 | { | ||
| 347 | 36 | EPS_BLOPEX *blopex = (EPS_BLOPEX*)eps->data; | |
| 348 | |||
| 349 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
36 | PetscFunctionBegin; |
| 350 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(VecDestroy(&blopex->w)); |
| 351 |
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.
|
12 | PetscFunctionReturn(PETSC_SUCCESS); |
| 352 | } | ||
| 353 | |||
| 354 | 36 | static PetscErrorCode EPSDestroy_BLOPEX(EPS eps) | |
| 355 | { | ||
| 356 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
36 | PetscFunctionBegin; |
| 357 |
5/8✓ Branch 0 taken 6 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
36 | PetscCall(PetscFree(eps->data)); |
| 358 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(PetscObjectComposeFunction((PetscObject)eps,"EPSBLOPEXSetBlockSize_C",NULL)); |
| 359 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(PetscObjectComposeFunction((PetscObject)eps,"EPSBLOPEXGetBlockSize_C",NULL)); |
| 360 |
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.
|
12 | PetscFunctionReturn(PETSC_SUCCESS); |
| 361 | } | ||
| 362 | |||
| 363 | 6 | static PetscErrorCode EPSView_BLOPEX(EPS eps,PetscViewer viewer) | |
| 364 | { | ||
| 365 | 6 | EPS_BLOPEX *ctx = (EPS_BLOPEX*)eps->data; | |
| 366 | 6 | PetscBool isascii; | |
| 367 | |||
| 368 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
6 | PetscFunctionBegin; |
| 369 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
6 | PetscCall(PetscObjectTypeCompare((PetscObject)viewer,PETSCVIEWERASCII,&isascii)); |
| 370 |
5/8✓ Branch 0 taken 6 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
6 | if (isascii) PetscCall(PetscViewerASCIIPrintf(viewer," block size %" PetscInt_FMT "\n",ctx->bs)); |
| 371 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
2 | PetscFunctionReturn(PETSC_SUCCESS); |
| 372 | } | ||
| 373 | |||
| 374 | 36 | static PetscErrorCode EPSSetFromOptions_BLOPEX(EPS eps,PetscOptionItems PetscOptionsObject) | |
| 375 | { | ||
| 376 | 36 | PetscBool flg; | |
| 377 | 36 | PetscInt bs; | |
| 378 | |||
| 379 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
36 | PetscFunctionBegin; |
| 380 |
1/12✗ Branch 0 not taken.
✓ Branch 1 taken 6 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.
|
36 | PetscOptionsHeadBegin(PetscOptionsObject,"EPS BLOPEX Options"); |
| 381 | |||
| 382 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(PetscOptionsInt("-eps_blopex_blocksize","Block size","EPSBLOPEXSetBlockSize",20,&bs,&flg)); |
| 383 |
6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
36 | if (flg) PetscCall(EPSBLOPEXSetBlockSize(eps,bs)); |
| 384 | |||
| 385 |
1/14✗ Branch 0 not taken.
✓ 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.
|
36 | PetscOptionsHeadEnd(); |
| 386 |
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.
|
12 | PetscFunctionReturn(PETSC_SUCCESS); |
| 387 | } | ||
| 388 | |||
| 389 | /*MC | ||
| 390 | EPSBLOPEX - EPSBLOPEX = "blopex" - A wrapper to BLOPEX {cite:p}`Kny07`. | ||
| 391 | |||
| 392 | Note: | ||
| 393 | In principle, the method used here is the same as in `EPSLOBPCG`, but in | ||
| 394 | practice convergence may be different. | ||
| 395 | |||
| 396 | Level: beginner | ||
| 397 | |||
| 398 | .seealso: [](ch:eps), `EPS`, `EPSType`, `EPSSetType()` | ||
| 399 | M*/ | ||
| 400 | 36 | SLEPC_EXTERN PetscErrorCode EPSCreate_BLOPEX(EPS eps) | |
| 401 | { | ||
| 402 | 36 | EPS_BLOPEX *ctx; | |
| 403 | |||
| 404 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
36 | PetscFunctionBegin; |
| 405 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(PetscNew(&ctx)); |
| 406 | 36 | eps->data = (void*)ctx; | |
| 407 | |||
| 408 | 36 | eps->categ = EPS_CATEGORY_PRECOND; | |
| 409 | |||
| 410 | 36 | eps->ops->solve = EPSSolve_BLOPEX; | |
| 411 | 36 | eps->ops->setup = EPSSetUp_BLOPEX; | |
| 412 | 36 | eps->ops->setupsort = EPSSetUpSort_Basic; | |
| 413 | 36 | eps->ops->setfromoptions = EPSSetFromOptions_BLOPEX; | |
| 414 | 36 | eps->ops->destroy = EPSDestroy_BLOPEX; | |
| 415 | 36 | eps->ops->reset = EPSReset_BLOPEX; | |
| 416 | 36 | eps->ops->view = EPSView_BLOPEX; | |
| 417 | 36 | eps->ops->backtransform = EPSBackTransform_Default; | |
| 418 | 36 | eps->ops->setdefaultst = EPSSetDefaultST_GMRES; | |
| 419 | |||
| 420 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(PetscObjectComposeFunction((PetscObject)eps,"EPSBLOPEXSetBlockSize_C",EPSBLOPEXSetBlockSize_BLOPEX)); |
| 421 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
36 | PetscCall(PetscObjectComposeFunction((PetscObject)eps,"EPSBLOPEXGetBlockSize_C",EPSBLOPEXGetBlockSize_BLOPEX)); |
| 422 |
5/8✓ Branch 0 taken 6 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
36 | if (slepc_blopex_useconstr < 0) PetscCall(PetscObjectComposedDataRegister(&slepc_blopex_useconstr)); |
| 423 |
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.
|
12 | PetscFunctionReturn(PETSC_SUCCESS); |
| 424 | } | ||
| 425 |