How can I prealloate a comlex-valued matrix?

3 Ansichten (letzte 30 Tage)
Hanna Liu
Hanna Liu am 12 Sep. 2019
Bearbeitet: Matthew Bayer am 22 Sep. 2021
The following code seems really slow, which I do allocating a matrix and then write some complex numbers in it:
tic,for kk = 1:1000
aa = (zeros(10000,500)); aa(:,300) = 1i * 3;
end;toc
The code is really slow:
Elapsed time is 42.463659 seconds.
Is there any way to accerlate this?
  4 Kommentare
Hanna Liu
Hanna Liu am 12 Sep. 2019
I think I didn't really explain my problem clear enough in the original description. I am currently working on a function, in which I need to compute a large complex-valued matrix column by column. When I run the matlab profiler, I find out that the code I first write values in my matrix is super slow.
The code I provide in the orignal description is just a simplified example of my current problem. The second line of the following is much slower:
aa = (zeros(10000,500)); aa(:,300) = 3;
bb = (zeros(10000,500)); bb(:,300) = 1i * 3;
I wonder if there is some faster way to work with large complex-valued matrix.
Many thanks for your comments!
Bruno Luong
Bruno Luong am 12 Sep. 2019
Bearbeitet: Bruno Luong am 12 Sep. 2019
"When I run the matlab profiler, I find out that the code I first write values in my matrix is super slow. "
MATLAB has complicated memory management than you would think. It can allocate and copy a big array when you assign something even tiny (copy-on-write). It can do nothing when you think you allocate a big array, but then the first time you read something it will allocated (allocate-on-use).
It is crucial for you to describe how you allocate, if your matrix is "shared" by something else, etc... or the data flow of your algorithm so we can advide.
You cannot make a simplify your code and ask for the advise based on the simplification, especially "how to allocate faster than zeros(...)", no there is NO simple answer to such question, because the time depends what your program did before.

Melden Sie sich an, um zu kommentieren.

Akzeptierte Antwort

Steven Lord
Steven Lord am 12 Sep. 2019
I would probably use the 'like' syntax for the zeros function.
A = zeros(10000, 500, 'like', 1i);
  1 Kommentar
Matthew Bayer
Matthew Bayer am 22 Sep. 2021
Bearbeitet: Matthew Bayer am 22 Sep. 2021
This is better than the other answer, A = complex(zeros(10000,500)). I recently tried both and found the 'like' syntax was considerably faster. The complex() statement is presumably first creating the real array and then converting to a complex one, while the 'like' syntax just creates the complex array right away.

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

Mehr zu Historical Contests finden Sie in Help Center und File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by