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