Re: Re: Multiplying a vector over multiple vectors
- To: mathgroup at smc.vnet.net
- Subject: [mg91200] Re: [mg91184] Re: Multiplying a vector over multiple vectors
- From: Andrzej Kozlowski <akoz at mimuw.edu.pl>
- Date: Sat, 9 Aug 2008 07:45:46 -0400 (EDT)
- References: <g7ed6u$2pm$1@smc.vnet.net> <200808081117.HAA11833@smc.vnet.net>
On 8 Aug 2008, at 13:17, Jean-Marc Gulliet wrote: > Donald DuBois wrote: > >> How does one get a vector to distribute over >> a list of vectors under multiplication? >> >> For example: >> >> {{1,2}, {3,4}, {5,6}} times {a,b} should >> equal (under this operation) >> >> {{a, 2b}, {3a, 4b}, {5a, 6b}} >> >> If {a,b} is duplicated three times one can >> use MapThread as in: >> >> MapThread[Times, {{{1, 2}, {3, 4}, {5, 6}}, {{a, b}, {a, b}, {a, >> b}}}] >> >> but duplicating {a,b} using Table is very slow >> when the vectors are large (consisting of thousands >> of elemets). > > FWIW, here are some possible ways with their respective timings: > > In[1]:= lst = RandomInteger[{-100, 100}, {10^5, 2}]; > Array[{lst[[#, 1]] a, lst[[#, 2]] b} &, Length@lst]; // > Timing // First > # {a, b} & /@ lst; // Timing // First > MapThread[Times, {lst, Table[{a, b}, {Length[lst]}]}]; // > Timing // First > MapThread[Times, {lst, ConstantArray[{a, b}, Length[lst]]}]; // > Timing // First > Table[{lst[[i, 1]] a, lst[[i, 2]] b}, {i, Length@lst}]; // > Timing // First > Transpose[{a lst[[All, 1]], b lst[[All, 2]]}]; // Timing // First > > Out[2]= 0.555598 > > Out[3]= 0.522737 > > Out[4]= 0.510453 > > Out[5]= 0.474941 > > Out[6]= 0.407999 > > Out[7]= 0.323102 > > Regards, > -- Jean-Marc > Somewhat surprisingly no-one seems to have suggested the approach that seems to me the most obvious: lst /. x_?VectorQ :> x {a, b} This seems to be the second fastest approach, on my machine loosing only slightly to Jean-Marc's last one. Andrzej Kozlowski
- References:
- Re: Multiplying a vector over multiple vectors
- From: Jean-Marc Gulliet <jeanmarc.gulliet@gmail.com>
- Re: Multiplying a vector over multiple vectors