Optimisation de code - somme

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:

---

1 2 3

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

---

Merci !

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

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.

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 !