Forums before death by AOL, social media and spammers... "We can't have nice things"
|    sci.math.symbolic    |    Symbolic algebra discussion    |    10,432 messages    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
|    Message 9,051 of 10,432    |
|    Albert Rich to clicl...@freenet.de    |
|    Re: Who is fastest AND optimal?    |
|    30 Apr 16 16:18:08    |
      From: Albert_Rich@msn.com              On Friday, April 29, 2016 at 11:43:40 AM UTC-10, clicl...@freenet.de wrote:       > And the winner is ... Rubi! FriCAS and Mathematica come second in my       > view, as the former times out on the last integral but supposedly       > returns reasonably compact and real elementary expressions otherwise,       > whereas the latter resorts to complex antiderivatives but doesn't slow       > down seriously. Maple's return of root objects for the first group and       > MuPAD's failure to evaluate the second group should place those systems       > last.       >        > I suppose that allowing not just binomials but also trinomials as       > numerator factors accounts for most of the increased Rubi rule count. On       > the other hand, the new two term-relations work without exceptions, and       > just three rules (of which one merely serves to enhance efficency)       > suffice to reduce any pair of exponents in any nondegenerate integrand       >        > (A + B*x + C*x^2)*(a + b*x + c*x^2)^m*(d + e*x + f*x^2)^n       >        > to the standard interval [-1, 0]; compare the (blazingly fast) Derive       > code appended below. A number of old Rubi rules (integrand expansion for       > integer exponents, special-casing the exponent -2) should therefore have       > become obsolete.       >        > Rubi will presumably suffer a similar size explosion for the P1^m*Q1^n*       > S2^p integrands treated in section 11 of the paper you linked to.       >        > I was actually assuming you had made your original decision in awareness       > of the consequences, and were therefore going to defend it!       >        > Martin.       > [...]              Thanks for announcing Rubi the winner of your product of powers of quadratics       integration challenge. However it should be noted that Rubi won only after       improvements were made in response to the challenge.              You comment that only 3 rules are required to reduce the exponents of       nondegenerate integrands to [-1,0]. Unfortunately, their implementation in       the real-world is not so simple. Understandably, Mathematica's pattern       matcher does not recognize        expressions of the form A, A+B x and A+C x^2 as matching expressions of the       form A+B x+C x^2 (that is B and C cannot be 0 in order to match). So the       new Rubi has to have 3 rules for integrating each of the following forms:                (a+b x+c x^2)^p (d+e x+f x^2)^q        (a+b x+c x^2)^p (d+e x+f x^2)^q (A+B x)        (a+b x+c x^2)^p (d+e x+f x^2)^q (A+C x^2)        (a+b x+c x^2)^p (d+e x+f x^2)^q (A+B x+C x^2)                (a+b x+c x^2)^p (d+f x^2)^q        (a+b x+c x^2)^p (d+f x^2)^q (A+B x)        (a+b x+c x^2)^p (d+f x^2)^q (A+C x^2)        (a+b x+c x^2)^p (d+f x^2)^q (A+B x+C x^2)                (a+c x^2)^p (d+e x+f x^2)^q        (a+c x^2)^p (d+e x+f x^2)^q (A+B x)        (a+c x^2)^p (d+e x+f x^2)^q (A+C x^2)        (a+c x^2)^p (d+e x+f x^2)^q (A+B x+C x^2)                (a+c x^2)^p (d+f x^2)^q        (a+c x^2)^p (d+f x^2)^q (A+B x)        (a+c x^2)^p (d+f x^2)^q (A+C x^2)        (a+c x^2)^p (d+f x^2)^q (A+B x+C x^2)               Clearly as the number of parameters to rules increases, the impracticality of       using pattern matching to implement a rule-based system like Rubi 4 becomes       obvious. Happily however, using a decision-tree to implement a rule-based       system like Rubi 5 avoids        this fan-out problem. This is because like your int22 function, a       decision-tree's if-then-else function can accept arguments having a value of 0.               So I need to resist the urge to add more knowledge to Rubi 4, and instead       focus my efforts on getting Rubi 5 out-the-door...               Albert              --- SoupGate-Win32 v1.05        * Origin: you cannot sedate... all the things you hate (1:229/2)    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
(c) 1994, bbs@darkrealms.ca