Avoiding doble loops to gain efficience
1 回表示 (過去 30 日間)
古いコメントを表示
Hi.
I have solved a problem by using a doble loop but it is taking too long and I am pretty sure there are more efficient ways to solve it, but I am quite a begginer with the program. I would be glad if you could help me with this.
Basically, I have two arrays of different size and I want to multiply each element of the first array for all the elements of the second one. The second vector also changes size in every iteration and that's why I have ended with a doble loop. I am doing something like this.
for j=1:length(s)
t=T-s(j):0.00025:T-0.00025;
for i=1:length(x)
aux1=exp(-x(i)^2/2./(T-t))./((t.*(T-t)).^(3/2));
aux2=sum(aux1)*0.00025;
omega(i)=exp(-mu^2/2*T+mu.*x(i)).*x(i)/2/pi.*aux2;
end
aux3=sum(omega)*0.00025;
Psnj(j)=aux3*exp(-lambda*s(j));
end
How could I do this more efficiently?
Thank you in advance Javier
0 件のコメント
回答 (1 件)
Arjun
2024 年 10 月 28 日
I see that you have solved a problem using double loop, but you are looking for more efficient ways to do this.
You can achieve this by utilizing MATLAB's vectorization capabilities. Instead of iterating over each element of "x" within a loop, you can vectorize the inner loop. Additionally, perform pre-computation to calculate any constant values outside the loop, ensuring that these operations aren't repeated with each loop execution.
Kindly refer to the code below for better understanding of the implementation:
% Assume some sample data for testing the code
T = 1.0;
mu = 0.5;
lambda = 0.1;
s = linspace(0.1, 0.5, 5);
x = linspace(0.1, 1.0, 10);
% Initialize result array at beginning to avoid resizing
Psnj_original = zeros(1, length(s));
% Precompute constants to avoid duplicate computation
dt = 0.00025;
constantFactor = exp(-mu^2/2*T) / (2*pi);
for j = 1:length(s)
t = T - s(j):dt:T - dt;
commonTerm1 = 1 ./ ((t .* (T - t)).^(3/2));
% Vectorized calculation of aux1 for all x
aux1 = exp(-x'.^2 / 2 ./ (T - t)) .* commonTerm1;
aux2 = sum(aux1, 2) * dt;
% Calculate omega using vectorization
omega = constantFactor * exp(mu .* x) .* x .* aux2';
% Sum omega and calculate Psnj
aux3 = sum(omega) * dt;
Psnj_optimized(j) = aux3 * exp(-lambda * s(j));
end
disp('Optimized Approach Results:');
disp(Psnj_optimized);
I hope this will help!
0 件のコメント
参考
カテゴリ
Help Center および File Exchange で Loops and Conditional Statements についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!