comp.lang.idl-pvwave archive: archive » vectorization challenge! (help!)

Home » Public Forums » archive » vectorization challenge! (help!)

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

Re: vectorization challenge! (help!) [message #54891 is a reply to message #54844]

Thu, 19 July 2007 10:56

MarioIncandenza
Messages: 231
Registered: February 2005

Senior Member

Let's throw some memory at the problem, shall we?
Note: this is completely untested.

n = 24
npeeps = 50
gn1 = findgen(n,npeeps)
gn2 = findgen(n,npeeps)
cutoff = .01
; here is the array-based version of your task
swap=randomu(seed,n,npeeps)
bswap=swap le cutoff; binary
nswap_per_row=total(bswap,1); number to swap in each row
nswap_mask=rebin(transpose(nswap_per_row),n,npeeps); each cell has N
to swap in that row

swapi=reverse(lindgen(n,npeeps) mod n),1); each row N-1 to 0
bswapi=swapi lt nswap_mask; if NSWAP_PER_ROW=N, swaps entire row

fswap=where(bswapi,nfswap)
if(nswap gt 0) then gn1[fswap] = gn2[fswap]

Report message to a moderator

[Message index]

		vectorization challenge! (help!) By: Conor on Tue, 17 July 2007 09:52
		Re: vectorization challenge! (help!) By: MarioIncandenza on Thu, 19 July 2007 10:56
		Re: vectorization challenge! (help!) By: alvin[1] on Thu, 19 July 2007 09:40
		Re: vectorization challenge! (help!) By: Brian Larsen on Thu, 19 July 2007 06:33
		Re: vectorization challenge! (help!) By: mattf on Thu, 19 July 2007 04:53
		Re: vectorization challenge! (help!) By: Maarten[1] on Thu, 19 July 2007 01:35
		Re: vectorization challenge! (help!) By: alvin[1] on Thu, 19 July 2007 01:01
		Re: vectorization challenge! (help!) By: Brian Larsen on Wed, 18 July 2007 13:37

Previous Topic:	Re: warp_tri()
Next Topic:	Re: how can idl find out information of other softwares

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 14:16:43 PDT 2025

Total time taken to generate the page: 1.36159 seconds