next up previous contents index
suivant: III.D.3.2 Vectorisation du calcul de I2. monter: III.D.3 Recettes d'implémentation sur machine vectorielle. précédent: III.D.3 Recettes d'implémentation sur machine vectorielle.   Table des matières   Index


III.D.3.1 Vectorisation du calcul de I1.

Rappelons l'algorithme optimal de programmation de I1: Cet algorithme, à partir de $\forall \alpha \in {\mathbb{C}}$, peut se transformer en sans rien changer à la précision du calcul, de l'ordre de pm. Sur calculateur Cray, nous proposons de calculer la fonction indicatrice à l'aide de la formule
(III.D.22) \begin{displaymath}
{{\mathcal{I}}}_{\left[\vert\alpha\vert\le \epsilon_1 \right]}=\max(0,nint(sign(1.,\epsilon_1-\vert\alpha\vert)))
\end{displaymath}

qui n'utilise que des fonctions intrinsèques au calculateur, fonctions dont la vectorisation est ``totale''.

Remarque 53   On peut aussi calculer ${{\mathcal{I}}}_{\left[\vert\alpha\vert\le \epsilon_1 \right]}$ par
(III.D.23) \begin{displaymath}%
{{\mathcal{I}}}_{\left[\vert\alpha\vert\le \epsilon_1 \right]}=nint(0.5+sign(0.5,\vert\alpha\vert-\epsilon_1))
\end{displaymath}

et aussi la fonction $I_1(\alpha)$ par
(III.D.24) \begin{displaymath}%
I_1(\alpha)= L e^{i{{\mathbf k}}\vec{G}} e^{i\alpha} \left(...
...}}}_{\left[\vert\alpha\vert\le \epsilon_1 \right]} \right) \ .
\end{displaymath}


next up previous contents index
suivant: III.D.3.2 Vectorisation du calcul de I2. monter: III.D.3 Recettes d'implémentation sur machine vectorielle. précédent: III.D.3 Recettes d'implémentation sur machine vectorielle.   Table des matières   Index
Cessenat Olivier 2007-04-21