[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to
From: |
Marco Caliari |
Subject: |
[Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps |
Date: |
Thu, 2 Jun 2016 07:51:55 +0000 (UTC) |
User-agent: |
Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:46.0) Gecko/20100101 Firefox/46.0 |
Follow-up Comment #37, bug #47800 (project octave):
Excellent, Nir. I see that you do not remove converged values in
__gammainc_lentz as you did for gammainc_s (jj subindex). Isn't it dangerous?
For instance, with your previous file #37329, if I do
gammainc(rand(4),rand(4))
quite often I have no convergence. The problem seems to me that at a certain
iteration, some values have converged, but at the next iteration their error
is a bit above the tolerance and other converged, and so on. However, I cannot
reproduce this behavior with your new file #37342.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?47800>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, anonymous, 2016/06/01
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, Marco Caliari, 2016/06/01
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, Lachlan Andrew, 2016/06/01
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, Nir Krakauer, 2016/06/01
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, Marco Caliari, 2016/06/01
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, Nir Krakauer, 2016/06/01
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, Nir Krakauer, 2016/06/01
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps,
Marco Caliari <=
- [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps, Nir Krakauer, 2016/06/02