First let me explain what I am trying to do in my code.

Assume that I send a very short chirp (2,000 ms long) into a domain filled with 20,000 particles randomly scattered in the domain.

enter image description here

The backscattered sounds are received by 8 receivers. Because the position of each particle varies in space relative to each receiver, the arrival time of each backscattered chirp varies in time.

Adding up the backscattered chirps from each of the 20,000 particles gives the backscattered signal. The first backscattered chirp arrives at t= 0 ms and the last one at t= 100,000 ms.

enter image description here

In the loop below I want to add up 20,000 signals with length of 100,000. Each of the backscattered signal include only 2,000 non-zero elements and the zeros before and after are added to mark the position of chirp through the time (t= 0 to 100,000).

This process should be repeated for each of the receivers (8 times). I tried to break the 20,000 signals into 20 x 1000 clusters to speed up the process, but still it is very slow.

   bsAll =

{bsi = ConstantArray[0, maxN]},

 Block[{bsCluster, bsClusTot},

  bsCluster =


        sstart ;; sfinish[[itr, jj]]]],
      ifinish[[itr, jj]]],
    {jj, ii - (cluster - 1), ii}];

  bsClusTot = MapThread[Plus, bsCluster];
  bsi = MapThread[Plus, {bsi, bsClusTot}];
 {ii, cluster, 20000 , cluster}];
{itr, ntr}];

ntr = number of receivers 
sstart = 1
sfinish = is a 8x20000 matrix. most of the components are =2000
ifinish = is a 8x20000 matrix. components are integers between 1 and 98000
maxN = 100000



bsAll for cluster=10, and {ii,cluster, 2cluster, cluster}








  • $\begingroup$ To me, at least, your question is unanswerable without a detailed description of the data you are attempting to process. At a minimum I would like to see a 100 data points that (statistically) model your experimental data. I would also want to see what you would consider a correct summation of the model data. $\endgroup$ – m_goldberg Sep 14 '15 at 23:24
  • $\begingroup$ the data is the chirp (the first pic, but can be any sinusoidal wave). The final summation is displayed in the second pic. $\endgroup$ – Mahdi Razaz Sep 14 '15 at 23:46
  • $\begingroup$ Seems to me populating a sparse array directly and just totaling would be vastly more efficient - iterating and doing two paddings each time spells doom for performance. $\endgroup$ – ciao Sep 14 '15 at 23:56
  • $\begingroup$ I don't have enough RAM to populate the sparse array and then do a total. 20,000 particles are not nearly enough for my work. $\endgroup$ – Mahdi Razaz Sep 15 '15 at 0:18
  • $\begingroup$ Are the 20,000 signals sorted somehow? $\endgroup$ – Marius Ladegård Meyer Sep 15 '15 at 8:26

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.