Re: Re: Re: Timing runs for the last part of my previous post

• To: mathgroup at smc.vnet.net
• Subject: [mg62206] Re: [mg62196] Re: Re: Timing runs for the last part of my previous post
• From: "Carl K. Woll" <carl at woll2woll.com>
• Date: Tue, 15 Nov 2005 04:16:05 -0500 (EST)
• References: <dkshq9\$jei\$1@smc.vnet.net> <200511100750.CAA07305@smc.vnet.net> <dl1jjb\$t8d\$1@smc.vnet.net> <200511140538.AAA27648@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```Maxim wrote:
> On Fri, 11 Nov 2005 08:09:15 +0000 (UTC), Oyvind Tafjord
> <tafjord at wolfram.com> wrote:
>
>
>>Here's yet another way which is about 30% faster then testFunc4:
>>
>>Pick[values, Plus @@ Sign[#] & /@ values, -2]
>>
>>At some point you'll need to process all pairs in the list, so it's
>>faster
>>to do this up front, in an effective way, and then use Pick to get the
>>elements you want.
>>
>>Oyvind Tafjord
>>Wolfram Research
>>
>>
>
>
> An even more efficient (and contrived) way is to make all the loops over
> the list elements implicit:
>
> In[1]:=
> L = Array[Random[Integer, {-10, 10}]&, {10^6, 2}];
> Cases[L, {_?Negative, _?Negative}]; // Timing
> Pick[L, Plus@@ Sign@ #& /@ L, -2]; // Timing
> Pick[L, Total@ Transpose@ Sign@ L, -2]; // Timing
>
> Out[2]= {1.563*Second, Null}
>
> Out[3]= {1.062*Second, Null}
>
> Out[4]= {0.468*Second, Null}
>
> The second version uses Map, which is fairly efficient but still not as
> good as applying listable functions to vectors.
>
> Maxim Rytin
> m.r at inbox.ru

It is possible to use SparseArrays to speed this up even further. First
convert the original list into a vector of 1s and 0s, where the 1s
indicate where a pair of negatives were. Then convert to SparseArrays
and extract the position of the nonzero elements:

In[15]:=
r1=L[[SparseArray[
1-Sign[Total[Transpose[1+Sign[L]]]]
]/.SparseArray[_,_,_,{_,{_,x_},_}]:>Flatten[x]]];//Timing
r2=Pick[L,Total@Transpose@Sign@L,-2];//Timing
r1===r2

Out[15]=
{0.188 Second,Null}

Out[16]=
{1.516 Second,Null}

Out[17]=
True

Carl Woll
Wolfram Research

```

• Prev by Date: Re: Defining variables in \$UserBaseDirectory/init.m
• Next by Date: NDSolve in 5.2
• Previous by thread: Re: Re: Timing runs for the last part of my previous post
• Next by thread: Re: Timing runs for the last part of my previous post