[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" f
From: |
Brian Gladman |
Subject: |
[Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess |
Date: |
Wed, 02 Oct 2013 13:36:20 +0000 |
User-agent: |
Mozilla/5.0 (Windows NT 6.2; WOW64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/29.0.1547.76 Safari/537.36 |
Follow-up Comment #4, bug #39713 (project gsl):
When the secant method is used to find a root with a linear function, its root
prediction is essentially perfect and it will give a function value that is
zero (or _very_ close). In this situation, when gsl_root_test_delta() is used
to test for termination without also using gsl_root_test_residual(), they will
get a return which indicates the need to continue the iteration even though
this is not necessary. And this WILL lead to the "derivative value is not
finite" error.
So, in my view, the real question raised by this issue is what should the
solver do if it is entered at a position where the previous iteration has
already found the root and the user has not terminated the iteration? To do
nothing risks an infinite loop. But to do what is done now indicates failure
for what is really success! Maybe a warning return value such as "root
already found'?
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?39713>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Brian Gladman, 2013/10/01
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Brian Gladman, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess,
Brian Gladman <=
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Rhys Ulerich, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Brian Gladman, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Rhys Ulerich, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Rhys Ulerich, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Brian Gladman, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Rhys Ulerich, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Brian Gladman, 2013/10/02
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Brian Gladman, 2013/10/03
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Brian Gladman, 2013/10/03
- [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess, Max, 2013/10/05