Skip to main content

list manipulation - How to use Reap and Sow instead of Append to


I have matrix and and i want to do increment in a do loop and and i want to store in list using append to. It works fine for less values, But when i want to do it for values around 1 million the program is slow.


 x = {{0.1},{0.2}};xm = {{0.4},{0.5}};
Do[datax[i] = {}; dataxm[i] = {};, {i, 1, 2}]


Do[
 x = x + t;
xm = xm + t;
 Do[AppendTo[datax[i], Flatten[{t, x[[i]]}]];, {i, 1, 2}];
Do[AppendTo[dataxm[i], Flatten[{t, xm[[i]]}]];, {i, 1, 2}], {t, 0, 1, 0.1}]


 datax[1]


 {{0., 0.1}, {0.1, 0.2}, {0.2, 0.4}, {0.3, 0.7}, {0.4, 1.1}, {0.5, 1.6}, {0.6, 2.2}, {0.7, 2.9}, {0.8, 3.7}, {0.9, 4.6}, {1., 5.6}}

datax[2]

{{0., 0.2}, {0.1, 0.3}, {0.2, 0.5}, {0.3, 0.8}, {0.4, 1.2}, {0.5, 
  1.7}, {0.6, 2.3}, {0.7, 3.}, {0.8, 3.8}, {0.9, 4.7}, {1., 5.7}}

dataxm[1]

 {{0., 0.4}, {0.1, 0.5}, {0.2, 0.7}, {0.3, 1.}, {0.4, 

  1.4}, {0.5, 1.9}, {0.6, 2.5}, {0.7, 3.2}, {0.8, 4.}, {0.9, 
  4.9}, {1., 5.9}}

 dataxm[2]

 {{0., 0.5}, {0.1, 0.6}, {0.2, 0.8}, {0.3, 1.1}, {0.4, 
  1.5}, {0.5, 2.}, {0.6, 2.6}, {0.7, 3.3}, {0.8, 4.1}, {0.9, 5.}, {1.,
   6.}}

Similar to the above I want to use reap and sow functions to speed it up but i can store only the last value. Why?



 x = {{0.1},{0.2}};xm = {{0.4},{0.5}};


Do[
x = x + t;xm = xm + t;

Do[datax[i] = Reap[Sow[{t, x[[i]]}]][[2, 1]], {i, 1, 2}];

Do[dataxm[i] = Reap[Sow[{t, xm[[i]]}]][[2, 1]], {i, 1, 2}];



, {t, 0, 1, 0.1}]

datax[1]
{{1., {5.6}}}

 datax[2]
{{1., {5.7}}}



datax[1]
{{1., {5.9}}}

datax[2]
{{1., {6.0}}

Can anyone help in fixing this issue? Thanks in advance



Answer



You need olny one Reap for all Sows -- that's primarily the advantage of using them.


Moreover, you may use tags to obtain datax[1] and datax[2]:



x = {{0.1}, {0.2}};
ClearAll[datax];
datax = Association@Reap[
    Do[
     x = x + t;
     Do[
      Sow[
       Flatten[{t, x[[tag]]}], 
       tag
       ],

      {tag, 1, 2}],
     {t, 0, 1, 0.1}],
    _, Rule][[2]]


<|1 -> {{0., 0.1}, {0.1, 0.2}, {0.2, 0.4}, {0.3, 0.7}, {0.4, 1.1}, {0.5, 1.6}, {0.6, 2.2}, {0.7, 2.9}, {0.8, 3.7}, {0.9, 4.6}, {1., 5.6}}, 2 -> {{0., 0.2}, {0.1, 0.3}, {0.2, 0.5}, {0.3, 0.8}, {0.4, 1.2}, {0.5, 1.7}, {0.6, 2.3}, {0.7, 3.}, {0.8, 3.8}, {0.9, 4.7}, {1., 5.7}}|>



Now you can access datax[1] and datax[2] as before.



Because the question changed: This should allow you to assemble more than one list in one go:



data = Merge[
   Reap[
     Do[
      x = x + t;
      xm = xm + t;
      Do[
       Sow[Flatten[{t, x[[i]]}], "datax" -> i];
       Sow[Flatten[{t, xm[[i]]}], "dataxm" -> i]
       ,
       {i, 1, 2}], {t, 0, 1, 0.1}],

     _, #1[[1]] -> Association[#1[[2]] -> #2] &][[2]]
   , Association
   ];

Now data["datax"] and data["dataxm"] should behave like datax and dataxm.



You can obtain essentially the same result much more efficiently with arrays produced by Table:


x = {0.1, 0.2};
xm = {0.4, 0.5};
{datax, dataxm} = Transpose[

   Table[{{{t, t}, x += t}, {{t, t}, xm += t}}, {t, 0, 1, 0.1}], 
   {3, 1, 4, 2}];

Now you have (notice the double brackets):


datax[[1]]
datax[[2]]
dataxm[[1]]
dataxm[[2]]



Even faster for this particular problem is to avoid cursive addition and to use vectorized operations.


The Reap/Sow method from above:


n = 1000000;

x = {0.1, 0.2};
xm = {0.4, 0.5};
First@AbsoluteTiming[
  data = Merge[
     Reap[
       Do[

        x = x + t;
        xm = xm + t;
        Do[
         Sow[Flatten[{t, x[[i]]}], "datax" -> i];
         Sow[Flatten[{t, xm[[i]]}], "dataxm" -> i]
         ,
         {i, 1, 2}], {t, 0., 1., 1./n}],
       _, #1[[1]] -> Association[#1[[2]] -> #2] &][[2]]
     , Association
     ];

  ]


11.4845



The Table method from above:


x = {0.1, 0.2};
xm = {0.4, 0.5};
First@AbsoluteTiming[
  {datax, dataxm} = 

    Transpose[
     Table[{{{t, t}, x += t}, {{t, t}, xm += t}}, {t, 0., 1., 
       1./n}], {3, 1, 4, 2}];
  ]


2.74062



Vectorized version:


x = {0.1, 0.2};

xm = {0.4, 0.5};
First@AbsoluteTiming[
  tlist = Subdivide[0., 1., n];
  slist = Accumulate[tlist];
  datax2 = {
    Transpose[{tlist, x[[1]] + slist}],
    Transpose[{tlist, x[[2]] + slist}]
    };
  dataxm2 = {
    Transpose[{tlist, xm[[1]] + slist}],

    Transpose[{tlist, xm[[2]] + slist}]
    };
  ]


0.061355



So, the latter is 44/187 time faster than the other methods. So you should take these considerations into account when using Reap and Sow.


Errors:


Max[Abs[data["datax"][1] - datax[[1]]]]

Max[Abs[data["datax"][2] - datax[[2]]]]
Max[Abs[data["dataxm"][1] - dataxm[[1]]]]
Max[Abs[data["dataxm"][2] - dataxm[[2]]]]
Max[Abs[datax - datax2]]
Max[Abs[dataxm - dataxm2]]


0.


0.


0.



0.


5.82077*10^-11


5.82077*10^-11



Comments

Post a Comment

Popular posts from this blog

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more

How to thread a list

I have data in format data = {{a1, a2}, {b1, b2}, {c1, c2}, {d1, d2}} Tableform: I want to thread it to : tdata = {{{a1, b1}, {a2, b2}}, {{a1, c1}, {a2, c2}}, {{a1, d1}, {a2, d2}}} Tableform: And I would like to do better then pseudofunction[n_] := Transpose[{data2[[1]], data2[[n]]}]; SetAttributes[pseudofunction, Listable]; Range[2, 4] // pseudofunction Here is my benchmark data, where data3 is normal sample of real data. data3 = Drop[ExcelWorkBook[[Column1 ;; Column4]], None, 1]; data2 = {a #, b #, c #, d #} & /@ Range[1, 10^5]; data = RandomReal[{0, 1}, {10^6, 4}]; Here is my benchmark code kptnw[list_] := Transpose[{Table[First@#, {Length@# - 1}], Rest@#}, {3, 1, 2}] &@list kptnw2[list_] := Transpose[{ConstantArray[First@#, Length@# - 1], Rest@#}, {3, 1, 2}] &@list OleksandrR[list_] := Flatten[Outer[List, List@First[list], Rest[list], 1], {{2}, {1, 4}}] paradox2[list_] := Partition[Riffle[list[[1]], #], 2] & /@ Drop[list, 1] RM[list_] := FoldList[Transpose[{First@li...
Post a Comment
Read more

front end - keyboard shortcut to invoke Insert new matrix

I frequently need to type in some matrices, and the menu command Insert > Table/Matrix > New... allows matrices with lines drawn between columns and rows, which is very helpful. I would like to make a keyboard shortcut for it, but cannot find the relevant frontend token command (4209405) for it. Since the FullForm[] and InputForm[] of matrices with lines drawn between rows and columns is the same as those without lines, it's hard to do this via 3rd party system-wide text expanders (e.g. autohotkey or atext on mac). How does one assign a keyboard shortcut for the menu item Insert > Table/Matrix > New... , preferably using only mathematica? Thanks! Answer In the MenuSetup.tr (for linux located in the $InstallationDirectory/SystemFiles/FrontEnd/TextResources/X/ directory), I changed the line MenuItem["&New...", "CreateGridBoxDialog"] to read MenuItem["&New...", "CreateGridBoxDialog", MenuKey["m", Modifiers-...
Post a Comment
Read more