comp.lang.idl-pvwave archive: archive » arithmetic operation on array

Home » Public Forums » archive » arithmetic operation on array

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

Re: arithmetic operation on array [message #85481 is a reply to message #85479]

Mon, 12 August 2013 16:08

Phillip Miller
Messages: 4
Registered: August 2013

Junior Member

On 2013-08-12 22:22:14 +0000, Phillip Bitzer said:

> OK, I'll bite. There are three ways I can think to do this off the top
> of my head:
> 1) Do a loop, like Phillip said (what a fantastic name :-) )
> 2) Rebin, like David said
> 3) Use the mysterious "add an extra dimension" method
> ( https://groups.google.com/d/msg/comp.lang.idl-pvwave/Vu9rzqc kBNQ/HvkK_QnJrsgJ
> and more recently
> https://groups.google.com/d/msg/comp.lang.idl-pvwave/dM8XXas Eio0/d3__pvX7svMJ)
>
>
> Here are the sample code I used:
>
> data = RANDOMU(1L, 360, 180, 456)
> avg = MEAN(data, DIM=3)
> mData1 = FLTARR(360, 180, 456)
>
> tic & for i=0, 455 do mData1[*,*,i] = data[*,*,i] - avg & toc
>
> mData2 = FLTARR(360, 180, 456)
>
> tic & mData2 = data - Rebin(avg, 360, 180, 456) & toc
> tic & data = TEMPORARY(data) - Rebin(avg, 360, 180, 456) & toc ;just
> for completeness
>
> data = RANDOMU(1L, 360, 180, 456) ;redefine data - we changed it above
> mData3 = FLTARR(360, 180, 456)
>
> tic & mData3 = data - avg[*, *, 0] & toc
>
> I used mData (modified data) arrays so I can check I get the same
> answer, regardless of the method.
>
> The four times I get are, in order relative to the above:
> % Time elapsed: 0.12749481 seconds.
> % Time elapsed: 0.33231401 seconds.
> % Time elapsed: 0.33304191 seconds.
> % Time elapsed: 0.019619942 seconds.
>
> So, it seems the mysterious extra dimension method is the fastest, by
> an order of magnitude. Whoa.
>
> For prosperity,
> IDL> print, !VERSION
> { x86_64 darwin unix Mac OS X 8.2.2 Jan 23 2013 64 64}

Excellent, thanks for the tip.
I tried it on my 360x180x1000 double precision array and got the
similar results as you, though not quite an order of magnitude
difference between it and the "improved" loop method, the mysterious
extra dimension method was indeed the fastest

Report message to a moderator

[Message index]

		arithmetic operation on array By: Phillip Miller on Mon, 12 August 2013 13:59
		Re: arithmetic operation on array By: David Fanning on Mon, 12 August 2013 14:24
		Re: arithmetic operation on array By: Phillip Bitzer on Mon, 12 August 2013 15:22
		Re: arithmetic operation on array By: Phillip Miller on Mon, 12 August 2013 16:08
		Re: arithmetic operation on array By: Phillip Miller on Mon, 12 August 2013 16:17
		Re: arithmetic operation on array By: David Fanning on Mon, 12 August 2013 16:18
		Re: arithmetic operation on array By: Phillip Bitzer on Mon, 12 August 2013 16:37
		Re: arithmetic operation on array By: Phillip Miller on Mon, 12 August 2013 15:51
		Re: arithmetic operation on array By: John Correira on Wed, 14 August 2013 12:24
		Re: arithmetic operation on array By: wlandsman on Wed, 14 August 2013 12:54

Previous Topic:	Coyote Graphics PS/PDF output size/orientation
Next Topic:	Not IDL Related Question

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 20:22:07 PDT 2025

Total time taken to generate the page: 0.31950 seconds