Re: element selection
- To: mathgroup at smc.vnet.net
- Subject: [mg53621] Re: [mg53605] element selection
- From: yehuda ben-shimol <benshimo at bgumail.bgu.ac.il>
- Date: Fri, 21 Jan 2005 06:35:25 -0500 (EST)
- References: <200501200847.DAA04095@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Bruyndonckx P. wrote:
>I have two long sorted lists of Integers. If the difference between a number in list 1 and a
>number in list 2 is less than a certain value, both are selected to create a couple. I am
>looking for a fast way to go through both lists and extract these elements to create a new list
>of couples. I tried a number of things, but all of them were rather slow so I need some expert
>advice here.
>
>Thx,
>
>Peter
>
>
>
Hi,
Call the first list a and the second b and the threshold th. In case
that you want to check elements of the same positions in both lists say
a[[1]] -b[[1]]<th, a[[2]]-b[[2]]<th etc use the following code(note that
both lists need to be of identical length)
sel[l1_List, l2_List, th_] := Select[Thread[{l1, l2}], #[[
1]] - #[[2]] < th &];
examining its performance on a modest laptop gives with two lists of
100000 elements each
a=Table[Random[Integer,{1,10}],{100000}];
b=Table[Ransom[Integer,{1,10}],{100000}];
pairs=sel[a,b,3];//Timing
results with
{0.751 Second, Null}
In case you want to compare EACH element of l1 agianst each element of
l2 (here l1 and l2 need not be of identical length)
you can use the following code
selAll[l1_List,l2_List,th_]:=Select[Tuples[{l1,l2}],#[[1]]-#[[2]]<th&];
here the number of pairs to be generated and checked is Length[l1] *
Length[l2] so the performance of two list with 1000 elements each
(1000000 pairs are generated and checked) is
a=Table[Random[Integer,{1,10}],{1000}];
b=Table[Ransom[Integer,{1,10}],{1000}];
pairsAll=selAll[a,b,3];//Timing
{9.153 Second, Null}
yehuda
- References:
- element selection
- From: "Bruyndonckx P." <pbruynd@vub.ac.be>
- element selection