|
From: | Lachlan Andrew |
Subject: | [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps |
Date: | Fri, 20 May 2016 06:44:34 +0000 (UTC) |
User-agent: | Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:46.0) Gecko/20100101 Firefox/46.0 |
Follow-up Comment #19, bug #47800 (project octave): Hmmm... Before this can be applied, it will definitely need to handle vectors, even if only by loops. There are many special cases in this code. I'd imagine that at least some of them can be vectorised efficiently. The iterative components are probably the most expensive, and they seem to be vectorisable, although iteration would need to continue until all had converged, and so the actual FLOP count could be much higher than currently. I'll add it to my to-do list, but will need reminding... _______________________________________________________ Reply to this item at: <http://savannah.gnu.org/bugs/?47800> _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/
[Prev in Thread] | Current Thread | [Next in Thread] |