Discussion :

[Loïc B.]

Bonjour,

j'aurai une question simple en fortran. Est-il possible d'optimiser encore plus le calcul suivant pour éviter de passer par une boucle (en utilisant les capacités de calcul vectoriel du f90...) .

Code :

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

1
2
3
forall (i=1:cut_number-1)
    sum(i)=sum(surfacic_data*fractions(i,:,:))
end forall

Merci !

[FR119492]

Salut!
Tu vas sur www.netlib.org --> Browse ; dans la liste, tu choisis blas, tu descends jusqu'à sscal.f et tu regardes comment ils font pour "increment equal to 1".
Jean-Marc Blanc

[Sylvain Bergeron]

2 commentaires :

• Le forall n'a pas la réputation d'être bien optimisé dans plusieurs compilateurs. La simple boucle do est souvent plus performante.
• Si tu pouvais inverser les indices de ton tableau pour sommer sur (:,:,i) plutôt que sur (i,:,: ), ce serait plus rapide puisque la somme se ferait alors sur des éléments consécutifs en mémoire.

[Loïc B.]

Sylvain, merci !
FR119492 : merci, mais j'ai bêtement oublié de dire que surfacic_data est un tableau à 2 dimensions, pas un scalaire, mais tu pouvais pas le deviner. Merci pour le lien !