integral::evalf(): don't attempt to ignore problems.

Don't ignore exceptions thrown by numerical integration routine.
In general, the code like this

try {
	// blah-blah
} catch (std::exception& err) { }

is just plain evil. Case in the point:

fsolve((1/(sqrt(2*Pi)))*integral(t,0,x,exp(-1/2*t^2))==0.5,x,0,100)
(cherry picked from commit 515171f0bcd42099c266713c3d605cd92cedd2e2)

diff --git a/ginac/integral.cpp b/ginac/integral.cpp
index 07670f3c9eaefad92831409915bf249506fd1248..f45d1ab04b74a30c1e42d31c4f1d7777097f5cab 100644 (file)
--- a/ginac/integral.cpp
+++ b/ginac/integral.cpp
@@ -188,9 +188,7 @@ ex integral::evalf(int level) const
        // results after subsituting a number for the integration variable.
        if (is_exactly_a<numeric>(ea) && is_exactly_a<numeric>(eb) 
                        && is_exactly_a<numeric>(ef.subs(x==12.34).evalf())) {
        // results after subsituting a number for the integration variable.
        if (is_exactly_a<numeric>(ea) && is_exactly_a<numeric>(eb) 
                        && is_exactly_a<numeric>(ef.subs(x==12.34).evalf())) {

-               try {

                        return adaptivesimpson(x, ea, eb, ef);
                        return adaptivesimpson(x, ea, eb, ef);

-               } catch (runtime_error &rte) {}

        }
 
        if (are_ex_trivially_equal(a, ea) && are_ex_trivially_equal(b, eb)
        }
 
        if (are_ex_trivially_equal(a, ea) && are_ex_trivially_equal(b, eb)