You could try changing the variable order in the CRE form using
rat( expression, var1, var2, ...)
var's can be not just variables, but also 'kernels' like f(x,y) or exp(z).
-s
On Thu, Jun 14, 2012 at 11:33 AM, Dmitry Shkirmanov <piminusmeson at bk.ru>wrote:
> Thanks, with this improvements maple still wins, but now it is only 3
> times faster than maxima. ( 66 seconds against 196 seconds for N= 6 in my
> loop). Is there possibility to improve this result more?
>
> Another speed up: load(lrats) an replace subst with lratsubst. Since
>> ratsubst knows a bit more about the meaning of expressions than
>> does subst, you'll need to decide it that's OK.
>>
>> --Barton
>>
>> ______________________________**__________
>>
>> Using CRE expressions reduces the run time for the big loop in your code
>> by a factor of about 3. To do this, insert something
>> like
>>
>> f : rat(f)
>> C[0] : rat(1/(2*pi^(3/2)*sqrt(rho[0, 0]*d)*m*sigma));
>>
>> into your code
>> --Barton
>>
>> ______________________________**__________
>> F
>>
>>
>>
>>
>
> ______________________________**_________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/**mailman/listinfo/maxima<http://www.math.utexas.edu/mailman/listinfo/maxima>;
>