| 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 monitors | ||
| 12 | */ | ||
| 13 | |||
| 14 | #include <slepc/private/epsimpl.h> /*I "slepceps.h" I*/ | ||
| 15 | #include <petscdraw.h> | ||
| 16 | |||
| 17 | /* | ||
| 18 | Runs the user provided monitor routines, if any. | ||
| 19 | */ | ||
| 20 | 450740 | PetscErrorCode EPSMonitor(EPS eps,PetscInt it,PetscInt nconv,PetscScalar *eigr,PetscScalar *eigi,PetscReal *errest,PetscInt nest) | |
| 21 | { | ||
| 22 | 450740 | PetscInt i,n = eps->numbermonitors; | |
| 23 | |||
| 24 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
450740 | PetscFunctionBegin; |
| 25 |
7/8✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
|
476876 | for (i=0;i<n;i++) PetscCall((*eps->monitor[i])(eps,it,nconv,eigr,eigi,errest,nest,eps->monitorcontext[i])); |
| 26 |
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.
|
90131 | PetscFunctionReturn(PETSC_SUCCESS); |
| 27 | } | ||
| 28 | |||
| 29 | /*@C | ||
| 30 | EPSMonitorSet - Sets an ADDITIONAL function to be called at every | ||
| 31 | iteration to monitor the error estimates for each requested eigenpair. | ||
| 32 | |||
| 33 | Logically Collective | ||
| 34 | |||
| 35 | Input Parameters: | ||
| 36 | + eps - the linear eigensolver context | ||
| 37 | . monitor - pointer to function (if this is `NULL`, it turns off monitoring), | ||
| 38 | see `EPSMonitorFn` | ||
| 39 | . mctx - [optional] context for private data for the monitor routine | ||
| 40 | (use `NULL` if no context is desired) | ||
| 41 | - monitordestroy - [optional] routine that frees monitor context (may be `NULL`), | ||
| 42 | see `PetscCtxDestroyFn` for the calling sequence | ||
| 43 | |||
| 44 | Options Database Keys: | ||
| 45 | + -eps_monitor - print only the first error estimate | ||
| 46 | . -eps_monitor_all - print error estimates at each iteration | ||
| 47 | . -eps_monitor_conv - print the eigenvalue approximations only when | ||
| 48 | convergence has been reached | ||
| 49 | . -eps_monitor draw::draw_lg - sets line graph monitor for the first unconverged | ||
| 50 | approximate eigenvalue | ||
| 51 | . -eps_monitor_all draw::draw_lg - sets line graph monitor for all unconverged | ||
| 52 | approximate eigenvalues | ||
| 53 | . -eps_monitor_conv draw::draw_lg - sets line graph monitor for convergence history | ||
| 54 | - -eps_monitor_cancel - cancels all monitors that have been hardwired into | ||
| 55 | a code by calls to `EPSMonitorSet()`, but does not cancel | ||
| 56 | those set via the options database. | ||
| 57 | |||
| 58 | Notes: | ||
| 59 | The options database option `-eps_monitor` and related options are the easiest way | ||
| 60 | to turn on `EPS` iteration monitoring. | ||
| 61 | |||
| 62 | `EPSMonitorRegister()` provides a way to associate an options database key with `EPS` | ||
| 63 | monitor function. | ||
| 64 | |||
| 65 | Several different monitoring routines may be set by calling `EPSMonitorSet()` multiple | ||
| 66 | times; all will be called in the order in which they were set. | ||
| 67 | |||
| 68 | Fortran Note: | ||
| 69 | Only a single monitor function can be set for each `EPS` object. | ||
| 70 | |||
| 71 | Level: intermediate | ||
| 72 | |||
| 73 | .seealso: [](ch:eps), `EPSMonitorFirst()`, `EPSMonitorAll()`, `EPSMonitorConverged()`, `EPSMonitorFirstDrawLG()`, `EPSMonitorAllDrawLG()`, `EPSMonitorConvergedDrawLG()`, `EPSMonitorCancel()` | ||
| 74 | @*/ | ||
| 75 | 1317 | PetscErrorCode EPSMonitorSet(EPS eps,EPSMonitorFn *monitor,void *mctx,PetscCtxDestroyFn *monitordestroy) | |
| 76 | { | ||
| 77 | 1317 | PetscInt i; | |
| 78 | 1317 | PetscBool identical; | |
| 79 | |||
| 80 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
1317 | PetscFunctionBegin; |
| 81 |
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.
|
1317 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 82 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
1341 | for (i=0;i<eps->numbermonitors;i++) { |
| 83 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
24 | PetscCall(PetscMonitorCompare((PetscErrorCode(*)(void))(PetscVoidFn*)monitor,mctx,monitordestroy,(PetscErrorCode (*)(void))(PetscVoidFn*)eps->monitor[i],eps->monitorcontext[i],eps->monitordestroy[i],&identical)); |
| 84 |
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.
|
24 | if (identical) PetscFunctionReturn(PETSC_SUCCESS); |
| 85 | } | ||
| 86 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1317 | PetscCheck(eps->numbermonitors<MAXEPSMONITORS,PetscObjectComm((PetscObject)eps),PETSC_ERR_ARG_OUTOFRANGE,"Too many EPS monitors set"); |
| 87 | 1317 | eps->monitor[eps->numbermonitors] = monitor; | |
| 88 | 1317 | eps->monitorcontext[eps->numbermonitors] = mctx; | |
| 89 | 1317 | eps->monitordestroy[eps->numbermonitors++] = monitordestroy; | |
| 90 |
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.
|
1317 | PetscFunctionReturn(PETSC_SUCCESS); |
| 91 | } | ||
| 92 | |||
| 93 | /*@ | ||
| 94 | EPSMonitorCancel - Clears all monitors for an `EPS` object. | ||
| 95 | |||
| 96 | Logically Collective | ||
| 97 | |||
| 98 | Input Parameter: | ||
| 99 | . eps - the linear eigensolver context | ||
| 100 | |||
| 101 | Options Database Key: | ||
| 102 | . -eps_monitor_cancel - cancels all monitors that have been hardwired into a code by calls to | ||
| 103 | `EPSMonitorSet()`, but does not cancel those set via the options database. | ||
| 104 | |||
| 105 | Level: intermediate | ||
| 106 | |||
| 107 | .seealso: [](ch:eps), `EPSMonitorSet()` | ||
| 108 | @*/ | ||
| 109 | 6960 | PetscErrorCode EPSMonitorCancel(EPS eps) | |
| 110 | { | ||
| 111 | 6960 | PetscInt i; | |
| 112 | |||
| 113 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
6960 | PetscFunctionBegin; |
| 114 |
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.
|
6960 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 115 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
8277 | for (i=0; i<eps->numbermonitors; i++) { |
| 116 |
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.
|
1317 | if (eps->monitordestroy[i]) PetscCall((*eps->monitordestroy[i])(&eps->monitorcontext[i])); |
| 117 | } | ||
| 118 | 6960 | eps->numbermonitors = 0; | |
| 119 |
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.
|
6960 | PetscFunctionReturn(PETSC_SUCCESS); |
| 120 | } | ||
| 121 | |||
| 122 | /*@C | ||
| 123 | EPSGetMonitorContext - Gets the monitor context, as set by | ||
| 124 | `EPSMonitorSet()` for the FIRST monitor only. | ||
| 125 | |||
| 126 | Not Collective | ||
| 127 | |||
| 128 | Input Parameter: | ||
| 129 | . eps - the linear eigensolver context | ||
| 130 | |||
| 131 | Output Parameter: | ||
| 132 | . ctx - monitor context | ||
| 133 | |||
| 134 | Level: intermediate | ||
| 135 | |||
| 136 | .seealso: [](ch:eps), `EPSMonitorSet()` | ||
| 137 | @*/ | ||
| 138 | ✗ | PetscErrorCode EPSGetMonitorContext(EPS eps,void *ctx) | |
| 139 | { | ||
| 140 | ✗ | PetscFunctionBegin; | |
| 141 | ✗ | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); | |
| 142 | ✗ | *(void**)ctx = eps->monitorcontext[0]; | |
| 143 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 144 | } | ||
| 145 | |||
| 146 | 240 | static inline PetscErrorCode EPSMonitorPrintEval(EPS eps,PetscViewer viewer,PetscScalar er,PetscScalar ei) | |
| 147 | { | ||
| 148 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
240 | PetscFunctionBegin; |
| 149 |
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.
|
240 | if (eps->problem_type==EPS_HEP || eps->problem_type==EPS_GHEP || eps->problem_type==EPS_BSE) PetscCall(PetscViewerASCIIPrintf(viewer," %g",(double)PetscRealPart(er))); |
| 150 | else { | ||
| 151 | #if defined(PETSC_USE_COMPLEX) | ||
| 152 |
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.
|
130 | PetscCall(PetscViewerASCIIPrintf(viewer," %g%+gi",(double)PetscRealPart(er),(double)PetscImaginaryPart(er))); |
| 153 | #else | ||
| 154 |
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.
|
110 | PetscCall(PetscViewerASCIIPrintf(viewer," %g",(double)er)); |
| 155 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 5 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.
|
110 | if (ei!=0.0) PetscCall(PetscViewerASCIIPrintf(viewer,"%+gi",(double)ei)); |
| 156 | #endif | ||
| 157 | } | ||
| 158 |
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.
|
48 | PetscFunctionReturn(PETSC_SUCCESS); |
| 159 | } | ||
| 160 | |||
| 161 | /*@C | ||
| 162 | EPSMonitorFirst - Print the first unconverged approximate value and | ||
| 163 | error estimate at each iteration of the eigensolver. | ||
| 164 | |||
| 165 | Collective | ||
| 166 | |||
| 167 | Input Parameters: | ||
| 168 | + eps - the linear eigensolver context | ||
| 169 | . its - iteration number | ||
| 170 | . nconv - number of converged eigenpairs so far | ||
| 171 | . eigr - real part of the eigenvalues | ||
| 172 | . eigi - imaginary part of the eigenvalues | ||
| 173 | . errest - error estimates | ||
| 174 | . nest - number of error estimates to display | ||
| 175 | - vf - viewer and format for monitoring | ||
| 176 | |||
| 177 | Options Database Key: | ||
| 178 | . -eps_monitor - activates `EPSMonitorFirst()` | ||
| 179 | |||
| 180 | Note: | ||
| 181 | This is not called directly by users, rather one calls `EPSMonitorSet()`, with this | ||
| 182 | function as an argument, to cause the monitor to be used during the `EPS` solve. | ||
| 183 | |||
| 184 | Level: intermediate | ||
| 185 | |||
| 186 | .seealso: [](ch:eps), `EPSMonitorSet()`, `EPSMonitorAll()`, `EPSMonitorConverged()` | ||
| 187 | @*/ | ||
| 188 | 20 | PetscErrorCode EPSMonitorFirst(EPS eps,PetscInt its,PetscInt nconv,PetscScalar eigr[],PetscScalar eigi[],PetscReal errest[],PetscInt nest,PetscViewerAndFormat *vf) | |
| 189 | { | ||
| 190 | 20 | PetscScalar er,ei; | |
| 191 | 20 | PetscViewer viewer = vf->viewer; | |
| 192 | |||
| 193 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
20 | PetscFunctionBegin; |
| 194 |
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.
|
20 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 195 |
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.
|
20 | PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,8); |
| 196 |
3/10✓ 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.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
|
20 | if (its==1 && ((PetscObject)eps)->prefix) PetscCall(PetscViewerASCIIPrintf(viewer," Eigenvalue approximations and residual norms for %s solve.\n",((PetscObject)eps)->prefix)); |
| 197 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
20 | if (nconv<nest) { |
| 198 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerPushFormat(viewer,vf->format)); |
| 199 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerASCIIAddTab(viewer,((PetscObject)eps)->tablevel)); |
| 200 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerASCIIPrintf(viewer,"%3" PetscInt_FMT " EPS nconv=%" PetscInt_FMT " first unconverged value (error)",its,nconv)); |
| 201 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerASCIIUseTabs(viewer,PETSC_FALSE)); |
| 202 | 20 | er = eigr[nconv]; ei = eigi[nconv]; | |
| 203 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(STBackTransform(eps->st,1,&er,&ei)); |
| 204 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(EPSMonitorPrintEval(eps,viewer,er,ei)); |
| 205 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerASCIIPrintf(viewer," (%10.8e)\n",(double)errest[nconv])); |
| 206 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerASCIIUseTabs(viewer,PETSC_TRUE)); |
| 207 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerASCIISubtractTab(viewer,((PetscObject)eps)->tablevel)); |
| 208 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
20 | PetscCall(PetscViewerPopFormat(viewer)); |
| 209 | } | ||
| 210 |
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); |
| 211 | } | ||
| 212 | |||
| 213 | /*@C | ||
| 214 | EPSMonitorAll - Print the current approximate values and | ||
| 215 | error estimates at each iteration of the eigensolver. | ||
| 216 | |||
| 217 | Collective | ||
| 218 | |||
| 219 | Input Parameters: | ||
| 220 | + eps - the linear eigensolver context | ||
| 221 | . its - iteration number | ||
| 222 | . nconv - number of converged eigenpairs so far | ||
| 223 | . eigr - real part of the eigenvalues | ||
| 224 | . eigi - imaginary part of the eigenvalues | ||
| 225 | . errest - error estimates | ||
| 226 | . nest - number of error estimates to display | ||
| 227 | - vf - viewer and format for monitoring | ||
| 228 | |||
| 229 | Options Database Key: | ||
| 230 | . -eps_monitor_all - activates `EPSMonitorAll()` | ||
| 231 | |||
| 232 | Note: | ||
| 233 | This is not called directly by users, rather one calls `EPSMonitorSet()`, with this | ||
| 234 | function as an argument, to cause the monitor to be used during the `EPS` solve. | ||
| 235 | |||
| 236 | Level: intermediate | ||
| 237 | |||
| 238 | .seealso: [](ch:eps), `EPSMonitorSet()`, `EPSMonitorFirst()`, `EPSMonitorConverged()` | ||
| 239 | @*/ | ||
| 240 | 110 | PetscErrorCode EPSMonitorAll(EPS eps,PetscInt its,PetscInt nconv,PetscScalar eigr[],PetscScalar eigi[],PetscReal errest[],PetscInt nest,PetscViewerAndFormat *vf) | |
| 241 | { | ||
| 242 | 110 | PetscInt i; | |
| 243 | 110 | PetscScalar er,ei; | |
| 244 | 110 | PetscViewer viewer = vf->viewer; | |
| 245 | |||
| 246 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
110 | PetscFunctionBegin; |
| 247 |
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.
|
110 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 248 |
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.
|
110 | PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,8); |
| 249 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerPushFormat(viewer,vf->format)); |
| 250 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerASCIIAddTab(viewer,((PetscObject)eps)->tablevel)); |
| 251 |
3/10✓ 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.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
|
110 | if (its==1 && ((PetscObject)eps)->prefix) PetscCall(PetscViewerASCIIPrintf(viewer," Eigenvalue approximations and residual norms for %s solve.\n",((PetscObject)eps)->prefix)); |
| 252 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerASCIIPrintf(viewer,"%3" PetscInt_FMT " EPS nconv=%" PetscInt_FMT " Values (Errors)",its,nconv)); |
| 253 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerASCIIUseTabs(viewer,PETSC_FALSE)); |
| 254 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
270 | for (i=0;i<nest;i++) { |
| 255 | 160 | er = eigr[i]; ei = eigi[i]; | |
| 256 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
160 | PetscCall(STBackTransform(eps->st,1,&er,&ei)); |
| 257 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
160 | PetscCall(EPSMonitorPrintEval(eps,viewer,er,ei)); |
| 258 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
160 | PetscCall(PetscViewerASCIIPrintf(viewer," (%10.8e)",(double)errest[i])); |
| 259 | } | ||
| 260 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerASCIIPrintf(viewer,"\n")); |
| 261 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerASCIIUseTabs(viewer,PETSC_TRUE)); |
| 262 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerASCIISubtractTab(viewer,((PetscObject)eps)->tablevel)); |
| 263 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
110 | PetscCall(PetscViewerPopFormat(viewer)); |
| 264 |
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.
|
22 | PetscFunctionReturn(PETSC_SUCCESS); |
| 265 | } | ||
| 266 | |||
| 267 | /*@C | ||
| 268 | EPSMonitorConverged - Print the approximate values and | ||
| 269 | error estimates as they converge. | ||
| 270 | |||
| 271 | Collective | ||
| 272 | |||
| 273 | Input Parameters: | ||
| 274 | + eps - the linear eigensolver context | ||
| 275 | . its - iteration number | ||
| 276 | . nconv - number of converged eigenpairs so far | ||
| 277 | . eigr - real part of the eigenvalues | ||
| 278 | . eigi - imaginary part of the eigenvalues | ||
| 279 | . errest - error estimates | ||
| 280 | . nest - number of error estimates to display | ||
| 281 | - vf - viewer and format for monitoring | ||
| 282 | |||
| 283 | Options Database Key: | ||
| 284 | . -eps_monitor_conv - activates `EPSMonitorConverged()` | ||
| 285 | |||
| 286 | Note: | ||
| 287 | This is not called directly by users, rather one calls `EPSMonitorSet()`, with this | ||
| 288 | function as an argument, to cause the monitor to be used during the `EPS` solve. | ||
| 289 | |||
| 290 | Level: intermediate | ||
| 291 | |||
| 292 | .seealso: [](ch:eps), `EPSMonitorSet()`, `EPSMonitorFirst()`, `EPSMonitorAll()` | ||
| 293 | @*/ | ||
| 294 | 10 | PetscErrorCode EPSMonitorConverged(EPS eps,PetscInt its,PetscInt nconv,PetscScalar eigr[],PetscScalar eigi[],PetscReal errest[],PetscInt nest,PetscViewerAndFormat *vf) | |
| 295 | { | ||
| 296 | 10 | PetscInt i; | |
| 297 | 10 | PetscScalar er,ei; | |
| 298 | 10 | PetscViewer viewer = vf->viewer; | |
| 299 | 10 | SlepcConvMon ctx; | |
| 300 | |||
| 301 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
10 | PetscFunctionBegin; |
| 302 |
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.
|
10 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 303 |
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.
|
10 | PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,8); |
| 304 | 10 | ctx = (SlepcConvMon)vf->data; | |
| 305 |
2/10✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 10 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
|
10 | if (its==1 && ((PetscObject)eps)->prefix) PetscCall(PetscViewerASCIIPrintf(viewer," Convergence history for %s solve.\n",((PetscObject)eps)->prefix)); |
| 306 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
10 | if (its==1) ctx->oldnconv = 0; |
| 307 |
1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
|
10 | if (ctx->oldnconv!=nconv) { |
| 308 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10 | PetscCall(PetscViewerPushFormat(viewer,vf->format)); |
| 309 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10 | PetscCall(PetscViewerASCIIAddTab(viewer,((PetscObject)eps)->tablevel)); |
| 310 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
70 | for (i=ctx->oldnconv;i<nconv;i++) { |
| 311 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
60 | PetscCall(PetscViewerASCIIPrintf(viewer,"%3" PetscInt_FMT " EPS converged value (error) #%" PetscInt_FMT,its,i)); |
| 312 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
60 | PetscCall(PetscViewerASCIIUseTabs(viewer,PETSC_FALSE)); |
| 313 | 60 | er = eigr[i]; ei = eigi[i]; | |
| 314 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
60 | PetscCall(STBackTransform(eps->st,1,&er,&ei)); |
| 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.
|
60 | PetscCall(EPSMonitorPrintEval(eps,viewer,er,ei)); |
| 316 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
60 | PetscCall(PetscViewerASCIIPrintf(viewer," (%10.8e)\n",(double)errest[i])); |
| 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.
|
60 | PetscCall(PetscViewerASCIIUseTabs(viewer,PETSC_TRUE)); |
| 318 | } | ||
| 319 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10 | PetscCall(PetscViewerASCIISubtractTab(viewer,((PetscObject)eps)->tablevel)); |
| 320 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
10 | PetscCall(PetscViewerPopFormat(viewer)); |
| 321 | 10 | ctx->oldnconv = nconv; | |
| 322 | } | ||
| 323 |
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); |
| 324 | } | ||
| 325 | |||
| 326 | 16 | PetscErrorCode EPSMonitorConvergedCreate(PetscViewer viewer,PetscViewerFormat format,void *ctx,PetscViewerAndFormat **vf) | |
| 327 | { | ||
| 328 | 16 | SlepcConvMon mctx; | |
| 329 | |||
| 330 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
16 | PetscFunctionBegin; |
| 331 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(PetscViewerAndFormatCreate(viewer,format,vf)); |
| 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.
|
16 | PetscCall(PetscNew(&mctx)); |
| 333 | 16 | mctx->ctx = ctx; | |
| 334 | 16 | (*vf)->data = (void*)mctx; | |
| 335 |
6/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
|
16 | PetscFunctionReturn(PETSC_SUCCESS); |
| 336 | } | ||
| 337 | |||
| 338 | 16 | PetscErrorCode EPSMonitorConvergedDestroy(PetscViewerAndFormat **vf) | |
| 339 | { | ||
| 340 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
16 | PetscFunctionBegin; |
| 341 |
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.
|
16 | if (!*vf) PetscFunctionReturn(PETSC_SUCCESS); |
| 342 |
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.
|
16 | PetscCall(PetscFree((*vf)->data)); |
| 343 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
16 | PetscCall(PetscViewerDestroy(&(*vf)->viewer)); |
| 344 |
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.
|
16 | PetscCall(PetscFree(*vf)); |
| 345 |
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); |
| 346 | } | ||
| 347 | |||
| 348 | /*@C | ||
| 349 | EPSMonitorFirstDrawLG - Plots the error estimate of the first unconverged | ||
| 350 | approximation at each iteration of the eigensolver. | ||
| 351 | |||
| 352 | Collective | ||
| 353 | |||
| 354 | Input Parameters: | ||
| 355 | + eps - the linear eigensolver context | ||
| 356 | . its - iteration number | ||
| 357 | . nconv - number of converged eigenpairs so far | ||
| 358 | . eigr - real part of the eigenvalues | ||
| 359 | . eigi - imaginary part of the eigenvalues | ||
| 360 | . errest - error estimates | ||
| 361 | . nest - number of error estimates to display | ||
| 362 | - vf - viewer and format for monitoring | ||
| 363 | |||
| 364 | Options Database Key: | ||
| 365 | . -eps_monitor draw::draw_lg - activates `EPSMonitorFirstDrawLG()` | ||
| 366 | |||
| 367 | Notes: | ||
| 368 | This is not called directly by users, rather one calls `EPSMonitorSet()`, with this | ||
| 369 | function as an argument, to cause the monitor to be used during the `EPS` solve. | ||
| 370 | |||
| 371 | Call `EPSMonitorFirstDrawLGCreate()` to create the context used with this monitor. | ||
| 372 | |||
| 373 | Level: intermediate | ||
| 374 | |||
| 375 | .seealso: [](ch:eps), `EPSMonitorSet()`, `EPSMonitorFirstDrawLGCreate()` | ||
| 376 | @*/ | ||
| 377 | 344 | PetscErrorCode EPSMonitorFirstDrawLG(EPS eps,PetscInt its,PetscInt nconv,PetscScalar eigr[],PetscScalar eigi[],PetscReal errest[],PetscInt nest,PetscViewerAndFormat *vf) | |
| 378 | { | ||
| 379 | 344 | PetscViewer viewer = vf->viewer; | |
| 380 | 344 | PetscDrawLG lg; | |
| 381 | 344 | PetscReal x,y; | |
| 382 | |||
| 383 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
344 | PetscFunctionBegin; |
| 384 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 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.
|
344 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 385 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 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.
|
344 | PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,8); |
| 386 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
344 | PetscCall(PetscViewerPushFormat(viewer,vf->format)); |
| 387 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
344 | PetscCall(PetscViewerDrawGetDrawLG(viewer,0,&lg)); |
| 388 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
344 | if (its==1) { |
| 389 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscDrawLGReset(lg)); |
| 390 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscDrawLGSetDimension(lg,1)); |
| 391 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscDrawLGSetLimits(lg,1,1,PetscLog10Real(eps->tol)-2,0.0)); |
| 392 | } | ||
| 393 |
1/2✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
|
344 | if (nconv<nest) { |
| 394 | 344 | x = (PetscReal)its; | |
| 395 |
1/2✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
|
344 | if (errest[nconv] > 0.0) y = PetscLog10Real(errest[nconv]); |
| 396 | ✗ | else y = 0.0; | |
| 397 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
344 | PetscCall(PetscDrawLGAddPoint(lg,&x,&y)); |
| 398 |
6/6✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 8 times.
✓ Branch 5 taken 8 times.
|
344 | if (its <= 20 || !(its % 5) || eps->reason) { |
| 399 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
216 | PetscCall(PetscDrawLGDraw(lg)); |
| 400 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
216 | PetscCall(PetscDrawLGSave(lg)); |
| 401 | } | ||
| 402 | } | ||
| 403 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
344 | PetscCall(PetscViewerPopFormat(viewer)); |
| 404 |
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.
|
43 | PetscFunctionReturn(PETSC_SUCCESS); |
| 405 | } | ||
| 406 | |||
| 407 | /*@C | ||
| 408 | EPSMonitorFirstDrawLGCreate - Creates the plotter for the first error estimate. | ||
| 409 | |||
| 410 | Collective | ||
| 411 | |||
| 412 | Input Parameters: | ||
| 413 | + viewer - the viewer | ||
| 414 | . format - the viewer format | ||
| 415 | - ctx - an optional user context | ||
| 416 | |||
| 417 | Output Parameter: | ||
| 418 | . vf - the viewer and format context | ||
| 419 | |||
| 420 | Level: intermediate | ||
| 421 | |||
| 422 | .seealso: [](ch:eps), `EPSMonitorSet()` | ||
| 423 | @*/ | ||
| 424 | 16 | PetscErrorCode EPSMonitorFirstDrawLGCreate(PetscViewer viewer,PetscViewerFormat format,void *ctx,PetscViewerAndFormat **vf) | |
| 425 | { | ||
| 426 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
16 | PetscFunctionBegin; |
| 427 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscViewerAndFormatCreate(viewer,format,vf)); |
| 428 | 16 | (*vf)->data = ctx; | |
| 429 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscViewerMonitorLGSetUp(viewer,NULL,"First Error Estimate","Log Error Estimate",1,NULL,PETSC_DECIDE,PETSC_DECIDE,400,300)); |
| 430 |
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.
|
2 | PetscFunctionReturn(PETSC_SUCCESS); |
| 431 | } | ||
| 432 | |||
| 433 | /*@C | ||
| 434 | EPSMonitorAllDrawLG - Plots the error estimate of all unconverged | ||
| 435 | approximations at each iteration of the eigensolver. | ||
| 436 | |||
| 437 | Collective | ||
| 438 | |||
| 439 | Input Parameters: | ||
| 440 | + eps - the linear eigensolver context | ||
| 441 | . its - iteration number | ||
| 442 | . nconv - number of converged eigenpairs so far | ||
| 443 | . eigr - real part of the eigenvalues | ||
| 444 | . eigi - imaginary part of the eigenvalues | ||
| 445 | . errest - error estimates | ||
| 446 | . nest - number of error estimates to display | ||
| 447 | - vf - viewer and format for monitoring | ||
| 448 | |||
| 449 | Options Database Key: | ||
| 450 | . -eps_monitor_all draw::draw_lg - activates `EPSMonitorAllDrawLG()` | ||
| 451 | |||
| 452 | Notes: | ||
| 453 | This is not called directly by users, rather one calls `EPSMonitorSet()`, with this | ||
| 454 | function as an argument, to cause the monitor to be used during the `EPS` solve. | ||
| 455 | |||
| 456 | Call `EPSMonitorAllDrawLGCreate()` to create the context used with this monitor. | ||
| 457 | |||
| 458 | Level: intermediate | ||
| 459 | |||
| 460 | .seealso: [](ch:eps), `EPSMonitorSet()`, `EPSMonitorAllDrawLGCreate()` | ||
| 461 | @*/ | ||
| 462 | 16 | PetscErrorCode EPSMonitorAllDrawLG(EPS eps,PetscInt its,PetscInt nconv,PetscScalar eigr[],PetscScalar eigi[],PetscReal errest[],PetscInt nest,PetscViewerAndFormat *vf) | |
| 463 | { | ||
| 464 | 16 | PetscViewer viewer = vf->viewer; | |
| 465 | 16 | PetscDrawLG lg; | |
| 466 | 16 | PetscInt i,n = PetscMin(eps->nev,255); | |
| 467 | 16 | PetscReal *x,*y; | |
| 468 | |||
| 469 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
16 | PetscFunctionBegin; |
| 470 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 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.
|
16 | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); |
| 471 |
3/16✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 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.
|
16 | PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,8); |
| 472 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscViewerPushFormat(viewer,vf->format)); |
| 473 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscViewerDrawGetDrawLG(viewer,0,&lg)); |
| 474 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
16 | if (its==1) { |
| 475 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
8 | PetscCall(PetscDrawLGReset(lg)); |
| 476 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
8 | PetscCall(PetscDrawLGSetDimension(lg,n)); |
| 477 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
8 | PetscCall(PetscDrawLGSetLimits(lg,1,1,PetscLog10Real(eps->tol)-2,0.0)); |
| 478 | } | ||
| 479 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscMalloc2(n,&x,n,&y)); |
| 480 |
2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
|
32 | for (i=0;i<n;i++) { |
| 481 | 16 | x[i] = (PetscReal)its; | |
| 482 |
2/4✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 8 times.
✗ Branch 3 not taken.
|
16 | if (i < nest && errest[i] > 0.0) y[i] = PetscLog10Real(errest[i]); |
| 483 | ✗ | else y[i] = 0.0; | |
| 484 | } | ||
| 485 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscDrawLGAddPoint(lg,x,y)); |
| 486 |
1/6✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
|
16 | if (its <= 20 || !(its % 5) || eps->reason) { |
| 487 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscDrawLGDraw(lg)); |
| 488 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscDrawLGSave(lg)); |
| 489 | } | ||
| 490 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscFree2(x,y)); |
| 491 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
16 | PetscCall(PetscViewerPopFormat(viewer)); |
| 492 |
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.
|
2 | PetscFunctionReturn(PETSC_SUCCESS); |
| 493 | } | ||
| 494 | |||
| 495 | /*@C | ||
| 496 | EPSMonitorAllDrawLGCreate - Creates the plotter for all the error estimates. | ||
| 497 | |||
| 498 | Collective | ||
| 499 | |||
| 500 | Input Parameters: | ||
| 501 | + viewer - the viewer | ||
| 502 | . format - the viewer format | ||
| 503 | - ctx - an optional user context | ||
| 504 | |||
| 505 | Output Parameter: | ||
| 506 | . vf - the viewer and format context | ||
| 507 | |||
| 508 | Level: intermediate | ||
| 509 | |||
| 510 | .seealso: [](ch:eps), `EPSMonitorSet()` | ||
| 511 | @*/ | ||
| 512 | 8 | PetscErrorCode EPSMonitorAllDrawLGCreate(PetscViewer viewer,PetscViewerFormat format,void *ctx,PetscViewerAndFormat **vf) | |
| 513 | { | ||
| 514 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
8 | PetscFunctionBegin; |
| 515 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
8 | PetscCall(PetscViewerAndFormatCreate(viewer,format,vf)); |
| 516 | 8 | (*vf)->data = ctx; | |
| 517 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 7 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
8 | PetscCall(PetscViewerMonitorLGSetUp(viewer,NULL,"All Error Estimates","Log Error Estimate",1,NULL,PETSC_DECIDE,PETSC_DECIDE,400,300)); |
| 518 |
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.
|
1 | PetscFunctionReturn(PETSC_SUCCESS); |
| 519 | } | ||
| 520 | |||
| 521 | /*@C | ||
| 522 | EPSMonitorConvergedDrawLG - Plots the number of converged eigenvalues | ||
| 523 | at each iteration of the eigensolver. | ||
| 524 | |||
| 525 | Collective | ||
| 526 | |||
| 527 | Input Parameters: | ||
| 528 | + eps - the linear eigensolver context | ||
| 529 | . its - iteration number | ||
| 530 | . nconv - number of converged eigenpairs so far | ||
| 531 | . eigr - real part of the eigenvalues | ||
| 532 | . eigi - imaginary part of the eigenvalues | ||
| 533 | . errest - error estimates | ||
| 534 | . nest - number of error estimates to display | ||
| 535 | - vf - viewer and format for monitoring | ||
| 536 | |||
| 537 | Options Database Key: | ||
| 538 | . -eps_monitor_conv draw::draw_lg - activates `EPSMonitorConvergedDrawLG()` | ||
| 539 | |||
| 540 | Notes: | ||
| 541 | This is not called directly by users, rather one calls `EPSMonitorSet()`, with this | ||
| 542 | function as an argument, to cause the monitor to be used during the `EPS` solve. | ||
| 543 | |||
| 544 | Call `EPSMonitorConvergedDrawLGCreate()` to create the context used with this monitor. | ||
| 545 | |||
| 546 | Level: intermediate | ||
| 547 | |||
| 548 | .seealso: [](ch:eps), `EPSMonitorSet()`, `EPSMonitorConvergedDrawLGCreate()` | ||
| 549 | @*/ | ||
| 550 | ✗ | PetscErrorCode EPSMonitorConvergedDrawLG(EPS eps,PetscInt its,PetscInt nconv,PetscScalar eigr[],PetscScalar eigi[],PetscReal errest[],PetscInt nest,PetscViewerAndFormat *vf) | |
| 551 | { | ||
| 552 | ✗ | PetscViewer viewer = vf->viewer; | |
| 553 | ✗ | PetscDrawLG lg; | |
| 554 | ✗ | PetscReal x,y; | |
| 555 | |||
| 556 | ✗ | PetscFunctionBegin; | |
| 557 | ✗ | PetscValidHeaderSpecific(eps,EPS_CLASSID,1); | |
| 558 | ✗ | PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,8); | |
| 559 | ✗ | PetscCall(PetscViewerPushFormat(viewer,vf->format)); | |
| 560 | ✗ | PetscCall(PetscViewerDrawGetDrawLG(viewer,0,&lg)); | |
| 561 | ✗ | if (its==1) { | |
| 562 | ✗ | PetscCall(PetscDrawLGReset(lg)); | |
| 563 | ✗ | PetscCall(PetscDrawLGSetDimension(lg,1)); | |
| 564 | ✗ | PetscCall(PetscDrawLGSetLimits(lg,1,1,0,eps->nev)); | |
| 565 | } | ||
| 566 | ✗ | x = (PetscReal)its; | |
| 567 | ✗ | y = (PetscReal)eps->nconv; | |
| 568 | ✗ | PetscCall(PetscDrawLGAddPoint(lg,&x,&y)); | |
| 569 | ✗ | if (its <= 20 || !(its % 5) || eps->reason) { | |
| 570 | ✗ | PetscCall(PetscDrawLGDraw(lg)); | |
| 571 | ✗ | PetscCall(PetscDrawLGSave(lg)); | |
| 572 | } | ||
| 573 | ✗ | PetscCall(PetscViewerPopFormat(viewer)); | |
| 574 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 575 | } | ||
| 576 | |||
| 577 | /*@C | ||
| 578 | EPSMonitorConvergedDrawLGCreate - Creates the plotter for the convergence history. | ||
| 579 | |||
| 580 | Collective | ||
| 581 | |||
| 582 | Input Parameters: | ||
| 583 | + viewer - the viewer | ||
| 584 | . format - the viewer format | ||
| 585 | - ctx - an optional user context | ||
| 586 | |||
| 587 | Output Parameter: | ||
| 588 | . vf - the viewer and format context | ||
| 589 | |||
| 590 | Level: intermediate | ||
| 591 | |||
| 592 | .seealso: [](ch:eps), `EPSMonitorSet()` | ||
| 593 | @*/ | ||
| 594 | ✗ | PetscErrorCode EPSMonitorConvergedDrawLGCreate(PetscViewer viewer,PetscViewerFormat format,void *ctx,PetscViewerAndFormat **vf) | |
| 595 | { | ||
| 596 | ✗ | SlepcConvMon mctx; | |
| 597 | |||
| 598 | ✗ | PetscFunctionBegin; | |
| 599 | ✗ | PetscCall(PetscViewerAndFormatCreate(viewer,format,vf)); | |
| 600 | ✗ | PetscCall(PetscNew(&mctx)); | |
| 601 | ✗ | mctx->ctx = ctx; | |
| 602 | ✗ | (*vf)->data = (void*)mctx; | |
| 603 | ✗ | PetscCall(PetscViewerMonitorLGSetUp(viewer,NULL,"Convergence History","Number Converged",1,NULL,PETSC_DECIDE,PETSC_DECIDE,400,300)); | |
| 604 | ✗ | PetscFunctionReturn(PETSC_SUCCESS); | |
| 605 | } | ||
| 606 |