Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
Drexel University or The Math Forum.



Re: deterministic integration
Posted:
Feb 27, 2013 11:43 PM


On Wednesday, February 27, 2013 2:03:06 AM UTC6, Alex Krasnov wrote: > This issue has already been discussed on this list, but the previous > > solutions are unsatisfactory. Integrate returns nondeterministic results > > depending on machine speed and kernel cache state. Example (Mathematica > > 8.0.4): > > > > In: Assuming[Element[z, Reals], Integrate[1/Sqrt[x^2+y^2+z^2], {x, 1, 1}, {y, 1, 1}]] > > Out: 4*(z*ArcCot[z*Sqrt[2 + z^2]] + Log[1  I*z] + Log[(1 + I*z)/(3 + z^2 + 2*Sqrt[2 + z^2])]) > > > > In: Assuming[Element[z, Reals], Integrate[1/Sqrt[x^2+y^2+z^2], {x, 1, 1}, {y, 1, 1}]] > > Out: 4*(z*ArcCot[z*Sqrt[2 + z^2]] + Log[(1 + z^2)/(3 + z^2 + 2*Sqrt[2 + z^2])]) > > > > The previous solutions involve adjusting the machine speed and clearing > > the kernel cache before each evaluation. This strategy attempts to achieve > > the least transformed result, though this is presumably not guaranteed. > > Instead, how can one achieve the most transformed result? The behavior of > > Integrate appears to be similar to that of Refine, Simplify, FullSimplify. > > However, unlike the latter, Integrate does not expose a TimeConstraint > > option. How can one achieve the effect of TimeConstraint > Infinity? Is > > the absence of this option related to the undecidability of Risch's > > algorithm or is the nondeterminism entirely in a subsequent > > simplification phase? > > > > Alex
Integrate uses a very large number of timeconstrained computations, and some of them also use time constraints under the hood.
There is some amount of discussion of this here:
http://library.wolfram.com/infocenter/Conferences/5832/
A revised/updated version may be found here:
http://www.sigsam.org/cca/issues/issue175.html
See in particular section 5.7.
As for effectively setting all time constraints to infinity, this might by and large work (except there are also several places where Integrate code bypasses the evaluator for timeconstraining).
Unprotect[TimeConstrained]; TimeConstrained[expr_,args___] := expr
This alters toe speed of your example above from 10 seconds to a minute; plan accordingly.
Daniel Lichtblau Wolfram Research



