comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » Vectorizing Code

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Vectorizing Code [message #29642 is a reply to message #29539]

Fri, 01 March 2002 05:40

robert.dimeo
Messages: 42
Registered: November 2001

Member

Hi Steve,

I vectorized it as follows:

u = 1+bytarr(nstate)
v = 1+bytarr(nvec)
d = u#(indgen(nvec)) - (indgen(nstate))#v
sa = 10.^2*exp(-abs(d)*dz/h)

Your way took 0.125 seconds on my pc for nstate = 500 and nvec = 120.
The vectorized way above took 0.03 seconds.

By the way, the rebin method proposed by Kenneth Bowman also took 0.03 seconds.

Hope this helps,

Rob

Steve Jones <no@email.com> wrote in message news:<270220022244527274%no@email.com>...
> Is it possible to vectorize a simple double for-loop?
>
> for i=0,nstate-1 do begin
> for j=0,nvec-1 do begin
> sa(i,j)=10.^2*exp(-abs(i-j)*dz/h)
> endfor
> endfor
>
> I tend to write a large number of such loops and my indexes have been
> steadily increasing of late... Is there a faster alternative? Thanks
> in advance

Report message to a moderator

[Message index]

		Vectorizing Code By: Steve Jones on Wed, 27 February 2002 14:44
		Re: Vectorizing Code By: robert.dimeo on Fri, 01 March 2002 05:40
		Re: Vectorizing Code By: Pavel A. Romashkin on Thu, 28 February 2002 09:26
		Re: Vectorizing Code By: David Fanning on Thu, 28 February 2002 08:23

Previous Topic:	Re: IDL to MATLAB binary files
Next Topic:	Best way to move a window

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 08:57:05 PDT 2025

Total time taken to generate the page: 0.00639 seconds