[ET Trac] [Einstein Toolkit] #734: McLachlan triggers arithmetic exception (but results are fine)

Einstein Toolkit trac-noreply at einsteintoolkit.org
Sun Feb 5 20:56:01 CST 2012


#734: McLachlan triggers arithmetic exception (but results are fine)
------------------------------------------+---------------------------------
  Reporter:  wolfgang.kastaun@…           |       Owner:           
      Type:  defect                       |      Status:  closed   
  Priority:  minor                        |   Milestone:           
 Component:  Other                        |     Version:           
Resolution:  fixed                        |    Keywords:  McLachlan
------------------------------------------+---------------------------------
Changes (by eschnett):

  * status:  new => closed
  * resolution:  => fixed


Comment:

 Thank you. This problem was caused by the expression

 etaExpr = Min [SpatialBetaDriverRadius / r, 1];

 which divides by r, and then (luckily?) ignores the result. I have changed
 this to

 etaExpr = SpatialBetaDriverRadius / Max [r, SpatialBetaDriverRadius];

 which is equivalent, but avoids the division by r if r is too small.

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/734#comment:1>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit


More information about the Trac mailing list