comp.lang.idl-pvwave archive: archive » Optimizing code for faster calculation

Home » Public Forums » archive » Optimizing code for faster calculation

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

Optimizing code for faster calculation [message #88032]

Wed, 12 March 2014 23:33

Kenneth D
Messages: 3
Registered: March 2014

Junior Member

I've been looking at this block of code now for... ever.

I've been editing a program created by my Adviser to reduce run time wherever possible. So far I've reduced the run time by nearly half, and I'm trying to juice any performance I can get from absolutely anywhere. My final project will use an array roughly 17,000 by 17,000. And I have to iterate through the program at least 17,000*10 times. If I'm lucky it won't take a month to process my data-sets now. This code is about all I have left to work with:

exceed_subs = where(min_rmse GT rmse_threshold, counter)
if counter GT 0 then modeled_class(exceed_subs) = "unmodeled"

min_rmse is an array Float[200], such as [0.347272, 0.312437, 0.360164,...]
rmse_threshold = 0.025
modeled_class is an array String[200], such as ["soil","quag","soil","grass",...]

The code find the locations where min_rmse is greater than a threshold value, and replaces those index locations in the string array (modeled_class) with "unmodeled".

This may well be the most efficient way to do this (this code will run a minimum of 17,000 times) but a look at the histograms page at Exelis:
http://www.exelisvis.com/docs/HISTOGRAM.html

shows:
For example, make the histogram of array A:
H = HISTOGRAM(A, REVERSE_INDICES = R)
;Set all elements of A that are in the ith bin of H to 0.
IF R[i] NE R[i+1] THEN A[R[R[I] : R[i+1]-1]] = 0

;The above is usually more efficient than the following:
bini = WHERE(A EQ i, count)
IF count NE 0 THEN A[bini] = 0

Which looks so similar to what I'm trying to do. I tried to implement this with no luck (maybe because strings?). Is there anything else I can do? That is, besides taking out iterations, they simply must be there to do what I need.

Report message to a moderator

[Message index]

		Optimizing code for faster calculation By: Kenneth D on Wed, 12 March 2014 23:33
		Re: Optimizing code for faster calculation By: Helder Marchetto on Thu, 13 March 2014 00:55
		Re: Optimizing code for faster calculation By: David Fanning on Thu, 13 March 2014 05:28
		Re: Optimizing code for faster calculation By: Heinz Stege on Thu, 13 March 2014 06:34
		Re: Optimizing code for faster calculation By: chris_torrence@NOSPAM on Thu, 13 March 2014 06:46
		Re: Optimizing code for faster calculation By: Kenneth D on Thu, 13 March 2014 17:46

Previous Topic:	Difficulty using cgimage2kml for swath data plotted in a cylindrical (Map_set) projection
Next Topic:	Force to print exponential to e-6

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 20:01:50 PDT 2025

Total time taken to generate the page: 0.04150 seconds