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