[Octave-bug-tracker] [bug #56884] mean can overflow

[Marco Caliari]

URL:
  <https://savannah.gnu.org/bugs/?56884>

                 Summary: mean can overflow
                 Project: GNU Octave
            Submitted by: caliari
            Submitted on: Thu 12 Sep 2019 09:07:22 AM UTC
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Inaccurate Result
                  Status: None
             Assigned to: None
         Originator Name: 
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: dev
        Operating System: GNU/Linux

    _______________________________________________________

Details:

The function mean can overflow


> mean([1.7e308,1.5e308])
ans =  Inf


even if the result is clearly 1.6e308. It appear to me that the simplest
solution is to replace


y = sum (x, dim) / n;


with


y = sum (x / n, dim);


around line 153 of mean.m. By the way, matlab overflows too.





    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?56884>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/