Discussion :

[cpalperou]

Bonjour,
j'ai un vecteur T1 de grandes dimensions (>100 et variable selon l'itération). J'aimerais utiliser ce vecteur si et seulement si il contient plus de 2 valeurs finis.
J'ai écrit ça:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
if Jack==1
      listeSpec1=[];
                for altTest = 21:41
                    for lamSPEC=1:1416
                        T1 = squeeze(Y(altTest,lamSPEC,:));
                             if numel(isfinite(T1)==1)>=2
                                 JackSpec1=Jack(T1);                 
                                 listeSpec1=[listeSpec1,JackSpec1];
                             end
                    end
                end
end

Y a t'il un moyen de simplifier cela?

[Jerome Briot]

La condition

Code :

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

if numel(isfinite(T1)==1)>=2

peut s'écrire plus simplement :

Code :

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

if sum(isfinite(T1))>=2

[cpalperou]

Merci Dut,
c'est vrai que c'est plus joli mais ça ne me fait pas gagner de temps de calcul.
Existe t'il une possibilité d'écrire quelquechose du style "si au moins 2 valeurs sont finis..." avec any par exemple?
Peut être est-il inutile d'essayer de gagner du temps de calcul sur cette condition?

[Jerome Briot]

Tu peux enlever le SQUEEZE comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
if Jack==1
    listeSpec1=[];
    for altTest = 21:41
        for lamSPEC=1:1416
            T1 = Y(altTest,lamSPEC,:);
            if sum(isfinite(T1(:)))>=2
                listeSpec1=
[listeSpec1,Jack(T1(:))];
            end
        end
    end
end

Ensuite, pour optimiser le code, il faut préallouer la mémoire pour listeSpec1 avant la double boucle FOR-END
=> : Qu'est-ce que la préallocation de mémoire ?

[cpalperou]

Merci,
le problème pour la préallocation de mémoire est que je ne connais pas la taille définitive de listeSPEC1. Donc je ne peux pas préallouer!
Tes corrections me conviennent, je t'en remercie

[Jerome Briot]

le problème pour la préallocation de mémoire est que je ne connais pas la taille définitive de listeSPEC1. Donc je ne peux pas préallouer!

Tu ne connais pas sa taille réelle... mais tu peux très certainement déterminer sa taille maximale...

Il te suffit donc de créer listeSPEC1 avec sa taille maximale théorique en le remplissant de 0 et de ne conserver ensuite que les lignes/colonnes contenant des valeurs différentes de 0.