Discussion :

[utawalpa]

Bonjour,

je travaille avec Matlab et Fortran. En gros, petit test sous Matlab pour voir si çà tourne bien puis Fortran pour passer la seconde dans les calculs, puis re-matlab pour visualiser les résultats.

Sous Matlab je vectorise assez facilement mes boucles FOR END et ne fait que du code Fortran pour des combinaisons de FOR IF END END.

ma question est surement très naïve mais, étant donné que graphiquement c'est quand même un peu plus pratique sous matlab, existe-t-il un moyen de vectoriser des boucles contenant de structures conditionnelles ?

On m'a laissé entendre que si je bouclais sous matlab c'était la honte mais je ne vois pas comment passer outre un IF END avec des vecteurs et des matrices ?

Merci !

[Jerome Briot]

Aurais-tu un exemple de code à vectoriser ?

[utawalpa]

Par exemple,

Code Fortran :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
 
CD = 0
DO N=1;NMAX 
     DO WHILE CD==0     
 
     CALL MYFUNCTION(VAL) 
 
     IF  (VAL>0.5) THEN
 
     VAL = 2.*VAL
     CD = 0
 
     ELSE
 
     VAL = VAL + 4.
     VECT(N) = VAL
     CD = 1
 
     END IF
 
     END DO
END DO

Pour l'exemple on peut considérer que "myfunction" est vectorisable

[Jerome Briot]

Prenons le cas de la condition IF sous MATLAB :

Code :

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

1
2
3
4
5
6
7
if VAL>.5
    VAL = 2*VAL
    CD = 0
else
    VAL = VAL+4
    CD = 1
end

La version sans la condition serait :

Code :

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

1
2
CD = VAL<=.5
VAL = 2*VAL*(~CD)+(VAL+4)*CD

C'est l'idée... maintenant pour le gain de temps... tout dépend de ton cas réel.

[utawalpa]

Merci Dut,

J'ai quelques cas un peu plus complexes où l'info doit évoluer en temps réel (fortran quand même donc...) mais ta proposition a le mérite de bien simplifier 90% de mon programme.