comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » looking for sort procedure
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: looking for sort procedure [message #7871 is a reply to message #7778] Tue, 21 January 1997 00:00 Go to previous messageGo to previous message
llobet is currently offline  llobet
Messages: 10
Registered: May 1993
Junior Member
In article <5blhri$ks1@post.gsfc.nasa.gov>, thompson@orpheus.nascom.nasa.gov
(William Thompson) writes:
=kak@sat.ipp-garching.mpg.de (Karl Krieger) writes:
=
=>"R. Bauer" <r.bauer@kfa-juelich.de> writes:
[...]
=>>for this example it needs on may RS6000 AIX more than 2 minutes.
=
=>>a = indgen(10000)
=>>b = [a,a]
=>>print,systime(0)
=>>s = b(sort(b))
=>>print,systime(0)
=
=>>end
=
=>>This is much too long.
=
=>Hi all,
=
=>there seems to be a bug in the implementation of this routine:
=>I tested it on a SUN Ultrasparc and on an IBM RS6000, which
=>has about the same speed (at least for the Monte Carlo simulation
=>coded in F77, which usually runs on these boxes).
=
=>Result for SUN: below 1 second, 6 seconds for a=lindgen(100000L)
=>Result for IBM: about 145 seconds
=
=>Is this a known bug/feature? There seems to be a major problem
=>either in the implementation of sort or in the way it uses
=>the machine's resources (bad optimization?).
=
=I also tried the above example on a DEC AXP 3000/600 where it took about 60
=seconds. I wonder if there's something in the code that is optimized for Sun
=workstations, maybe going back to the days when the first Unix port of IDL was
=called SunIDL?
=
=Bill Thompson

Here are the results on a DEC 3000 Model 400:


IDL. Version 3.6.1 (vms alpha).
[...]

IDL> a = indgen(10000) & b=[a,a]
IDL> print,systime(0) & s = b(sort(b)) & print,systime(0)
Tue Jan 21 09:58:59 1997
Tue Jan 21 09:59:00 1997
IDL>
IDL> a=lindgen(100000L) & b=[a,a]
IDL> print,systime(0) & s = b(sort(b)) & print,systime(0)
Tue Jan 21 10:00:18 1997
Tue Jan 21 10:00:27 1997

Which IDL version are you using?

-xavier
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Speckle
Next Topic: KILL_NOTIFY Programs

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Sat Oct 11 12:20:05 PDT 2025

Total time taken to generate the page: 0.40012 seconds