Optimisation de code MATLAB

Bonjour,

Je travaille sur MATLAB en ce moment et j'ai un problème de temps d'exécution. La cause : 5 boucles imbriquées 11x128x128x128x128 = 3 milliards d'itérations. Je sais que c'est mal mais cela fait plusieurs jours / semaines que je réfléchis à comment optimiser mon code, en passant par l'écriture matricielle, en vain. Ci-dessous un copier-coller de mon code pour voir si vous avez des idées d'optimisation.

J'ai également pensé à convertir ce bout de code en C, C++ ou Fortran, via les MEX-files de MATLAB. J'ai demandé à Mathworks une version d'évaluation de leur toolbox "MATLAB Coder" (http://www.mathworks.fr/products/matlab-coder/index.html), qui m'a l'air intéressant. Je me trompe ?

Merci de me tenir au courant de vos idées. Et merci de votre aide !

-------------

for klat_s = 1 : 11

disp( [ ' ANALYSE DE LA LATITUDE : ' , num2str( lat_s(klat_s) ) ] ) ;

fabcd = zeros( nt , nx ) ;
fst = zeros( nt , nx ) ;
fwr = zeros( nt , nx ) ; fwl = zeros( nt ,nx ) ;
fer = zeros( nt , nx ) ; fel = zeros( nt ,nx ) ;
fwerl = zeros( nt , nx ) ;

% Balayage en temps
for jt = 1 : 128
tic
disp(['jt = ', num2str(jt)])

% Balayage en espace
for ix = 1 : 128
disp([' ix = ', num2str(ix)])

dfabcd = 0 ;
dfss = 0 ;
dfww = 0 ; dfwr = 0 ; dfwl = 0 ;
dfee = 0 ; dfer = 0 ; dfel = 0 ;

% Balayage des pulsations
for jl = 1 : 128

if( pe(jl) <= 10 )
iw_r = 1 ; iw_l = 0 ;
ie_r = 1 ; ie_l = 0 ;
elseif( pe(jl) <= 20 )
iw_r = 0 ; iw_l = 1 ;
ie_r = 1 ; ie_l = 0 ;
else
iw_r = 0 ; iw_l = 1 ;
ie_r = 0 ; ie_l = 1 ;
end

% Balayage des vecteurs d'onde
for ik = 1 : 128

dfabcd = dfabcd + cos(k(ik)*x(ix)) * (a_u( jl , ik )*cos(l(jl)*t(jt)) + b_u( jl , ik )*sin(l(jl)*t(jt))) ...
+ sin(k(ik)*x(ix)) * (c_v( jl , ik )*cos(l(jl)*t(jt)) + d_v( jl , ik)*sin(l(jl)*t(jt))) ;

dfss = dfss + ( sc_f( jl , ik ) * cos(k(ik)*x(ix)) ...
+ ss_f( jl , ik ) * sin(k(ik)*x(ix))) ...
* cos( l(jl)*t(jt) + sph_f( jl , ik ) ) ;

dfw = wc_f( jl , ik ) * cos( k(ik)*x(ix) + l(jl)*t(jt) ) ...
+ ws_f( jl , ik ) * sin( k(ik)*x(ix) + l(jl)*t(jt) ) ;
dfww = dfww + dfw ;
dfwr = dfwr + iw_r * dfw ;
dfwl = dfwl + iw_l * dfw ;

dfe = ec_f( jl , ik ) * cos( k(ik)*x(ix) - l(jl)*t(jt) ) ...
+ es_f( jl , ik ) * sin( k(ik)*x(ix) - l(jl)*t(jt) ) ;
dfee = dfee + dfe ;
dfer = dfer + ie_r * dfe ;
dfel = dfel + ie_l * dfe ;

end
end

fabcd( jt , ix ) = dfabcd ;
fst( jt , ix ) = dfss ;
fwr( jt , ix ) = dfwr ;
fwl( jt , ix ) = dfwl ;
fer( jt , ix ) = dfer ;
fel( jt , ix ) = dfel ;
fwerl( jt , ix ) = dfwr + dfwl + dfer + dfel ;

end
t_toc=toc;
disp(['Temps d''execution : ', num2str(toc)]); % donne 223 secondes en moyenne
end
end