[Débutant] Mex files

**brel380** · 08/06/2007, 11h12

bonjour,

j'ai un programme matlab qui marche bien (tant mieux), mais vu que je fais des calculs ca me prend beaucoup de temps surtout quand je fais un assemblage de matrice (pour info : je modelise par elements finis)...
vous pouvez m'aide pour ecrire cette partie en Mex files, j'y arrive pas...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
for i=1:5
        if i~=5
            M1 (1, LMC(i)) = M1( 1,LMC(i) ) + k1(1,i) ;
        end
        for j=1:5
            M2( LMC(i) , LMC(j) ) = M2( LMC(i) , LMC(j) ) + k2(i,j);
 
            M3( LMC(i) , LMC(j) ) = M3( LMC(i) , LMC(j) ) + k(i,j) +k3(i,j) +k4(i,j) +  0.5 * k5(i,j) ;
        end
 
    end

merci

**Caro-Line** · 08/06/2007, 11h35

1 => mets les balises code (#) et indente ton code (CTRL+A,CTRL+I sous l'éditeur de Matlab), c'est plus agréable à lire

2 => avant de vouloir faire des mex-files (et là j'avoue ne pas pouvoir aider), il faut voir si une simple optimisation de code sous Matlab n'arrangerait pas les choses.

Par exemple : pourquoi faire 5 fois dans ta boucle le test "if i~=5" alors que tu pourrais simplement faire une boucle de i=1:4 puis ton traitement pour i=5 en dehors :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
for i=1:4%5
%     if i~=5
%         M1 (1, LMC(i)) = M1( 1,LMC(i) ) + k1(1,i) ;
%     end
    for j=1:5
        M2( LMC(i) , LMC(j) ) = M2( LMC(i) , LMC(j) ) + k2(i,j);
        M3( LMC(i) , LMC(j) ) = M3( LMC(i) , LMC(j) ) + k(i,j) +k3(i,j) +k4(i,j) + 0.5 * k5(i,j) ;
    end
 
end
%sortir i=5 de la boucle evite de faire 5 fois le test pour rien
%i=5
M1 (1, LMC(5)) = M1( 1,LMC(5) ) + k1(1,5) ;

Après il y a sûrement moyen de vectoriser tout ça mais c'est difficile sans avoir les données.

Autre remarque : il faut mieux éviter les variables M1, M2, M3, k1, k2 ...etc.
Tu peux utiliser des cell array :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part