Re: Slow Show/Graphics in v6.0
- To: mathgroup at smc.vnet.net
- Subject: [mg81053] Re: Slow Show/Graphics in v6.0
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Sun, 9 Sep 2007 06:17:30 -0400 (EDT)
- Organization: The Open University, Milton Keynes, UK
- References: <fbljvd$qf5$1@smc.vnet.net> <fbqpr3$6hv$1@smc.vnet.net> <fbtks5$8ih$1@smc.vnet.net>
Alex Shvydky wrote: <snip> > Here is another purely numerical example: > > profr[r_, s_] := E^(Log[0.05]*r^s); > fonebeam[th_, s_] := profr[Sin[th], s]*Max[0, Cos[th]]; > lmax = 10; > Table[ > Interpolation[({#1, fonebeam[#1, s]} &) /@ Range[0, Pi, > Pi/(30.*lmax)]]; > NIntegrate[finterp[th]*LegendreP[#1, Cos[th]]*Sin[th], {th, 0, Pi}, > MaxRecursion -> 16] & /@ Range[0, lmax], > {s, 1.9, 3., 0.02}]; > > took 3.78sec in Mathematica V5.2 and 39.97sec in V6.0. Alex, Could you, please, post the definition of finterp (and of any other required definitions if any) so one can experiment on one's system? Regards, -- Jean-Marc