Hey Gabriel,
I could give you some hints to speed up your code:
- I noticed that some computations are being repeated, for example:
h(i,j)*th(i,j) % This is being computed multiple times for the variable PE
The repeated computations could be stored and reused.
- Another instance of similar repeated computation is
h(i-1,j)*th(i-1,j) % Notice the indices i-1
% if h(i,j)*th(i,j) is already computed then h(i-1,j)*th(i-1,j) would have
% been computed in a prior iteration
- MATLAB is optimized for matrix and vector operations. Try vectorizing the operations instead of computing within nested for loops.
- Simplify and combine mathematical operations where possible to reduce computational overhead.
- Ensure that all matrices (p, th, etc.) are preallocated before the loop to avoid dynamic resizing during each iteration.
- If you have the Parallel Computing Toolbox, you can use parfor to parallelize the loop.
- Use MATLAB’s profile function to identify bottlenecks in your code.
For further information on improving performance, refer to the following documentation: https://www.mathworks.com/help/matlab/matlab_prog/techniques-for-improving-performance.html#buwj1l5
I hope the above hints were useful!
