Azzera filtri
Azzera filtri

Optimization of code with double for-loops

6 visualizzazioni (ultimi 30 giorni)
I am looking to optimize the following code specifically, and more generally, I am trying to learn optimization techniques in MATLAB for a larger project. The code below is just a snippet example which I am looking to learn techniques and then apply to a larger problem. I am potentially looking at values of L being as high as with four nested loops as opposed to just the two shown below.
Any pointers (or links to decent YouTube video on MATLAB optimization) would be welcomed.
clc; clearvars;
T = 1; Ts = 1e-2; t = 0:Ts:L*T;
L = 2^5;
start_time = clock();
for l = 0:L-1
for q = 0:L-1
flag = (t>=(l-q-1)*T) & (t<=(l-q+1)*T);
partA = (T-abs(t-(l-q)*T));
end
end
elapsed_time = clock();
time = elapsed_time(4)*3600+elapsed_time(5)*60+elapsed_time(6)
  4 Commenti
Howard Wilton
Howard Wilton il 20 Nov 2022
The questions didn't really hit the answer so i thought i would cast the questions differently.
However, your advice is duly noted and i shall be careful in the future. Thanks.
Howard Wilton
Howard Wilton il 20 Nov 2022
Modificato: Howard Wilton il 20 Nov 2022
I have taken the code a bit further and managed to dispense with the nested loops . However, I still need a 1D loop. Performance is substantially improved.
clc; clearvars;
L = 2^8; T = 1; t = 0:0.2:L*T;
start_time = clock();
l = 1:L; q = 1:L;
X = (l(:)-q(:).')*T;
for loop = 1:length(t)
t2(:,:,loop) = t(loop) - X;
end
stop_time = elapsedTime(start_time)
function time = elapsedTime(start_time)
elapsed_time = clock-start_time;
time = elapsed_time(4)*3600+elapsed_time(5)*60+elapsed_time(6);
end
Keen to figure out a way to eliminate the 1D loop and generate a matrix .

Accedi per commentare.

Risposta accettata

Walter Roberson
Walter Roberson il 20 Nov 2022
X = (l(:)-q(:).')*T;
for loop = 1:length(t)
t2(:,:,loop) = t(loop) - X;
end
replace with
t2 = reshape(t, 1, 1, []) - X;

Più risposte (0)

Categorie

Scopri di più su Programming in Help Center e File Exchange

Prodotti


Release

R2020b

Community Treasure Hunt

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

Start Hunting!

Translated by