[ET Trac] #2617: AEILocalInterp can produce GBs of level 1 warnings if interpolation points are outside of the domain

Roland Haas trac-noreply at einsteintoolkit.org
Wed Jul 6 17:06:58 CDT 2022


#2617: AEILocalInterp can produce GBs of level 1 warnings if interpolation points are outside of the domain

 Reporter: Roland Haas
   Status: new
Milestone: 
  Version: 
     Type: bug
 Priority: minor
Component: EinsteinToolkit thorn

AEILocalInterp will produce warning about each point that it is asked to interpolate to that is outside of the domain. This can easily lead to many many essentially identical warnings if eg the center of a spherical surface is set to NaN but a failure of shift tracking.

Instead of producing warnings all the time I would propose the following change in behaviour:

* by default output no more than \(say\) 10 points, then be silent and report “and NN more warnnings” at the end
* add an option to either make the warning level 0 right away \(an error\) or once 10 pints have been reached 

the limits would apply to each interpolator call \(LocalInterpUniform call really so technically per grid component\) and thus reset after each call.

Note: AEILocalInterp uses a level 1 warning since it actually does follow documented \(but mostly ignored\) Cactus design that states that a level 1 warning is almost certainly going to lead to incorrect result but that the code can continue \([http://einsteintoolkit.org/referencemanual/ReferenceManual.html#x1-235000A2](http://einsteintoolkit.org/referencemanual/ReferenceManual.html#x1-235000A2)\). So it _is_ indeed using the warning level as documented, just not like everyone else does.

> #define CCTK\_WARN\_ABORT    0    /\* abort the Cactus run \*/    
> \#define CCTK\_WARN\_ALERT    1    /\* the results of this run will probably \*/    
>                                 /\* be wrong, but this isn’t quite certain, \*/    
>                                 /\* so we’re not going to abort the run \*/    
> \#define CCTK\_WARN\_COMPLAIN 2    /\* the user should know about this, but \*/    
>                                 /\* the results of this run are probably ok \*/    
> \#define CCTK\_WARN\_PICKY    3    /\* this is for small problems that can \*/    
>                                 /\* probably be ignored, but that careful \*/    
>                                 /\* people may want to know about \*/    
> \#define CCTK\_WARN\_DEBUG    4    /\* these messages are probably useful \*/    
>                                 /\* only for debugging purposes \*/

‌

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2617/aeilocalinterp-can-produce-gbs-of-level-1
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/trac/attachments/20220706/ea886854/attachment.html 


More information about the Trac mailing list