Algorithmes de la morphologie mathématique pour - Pastel - HAL

Jaromír BRAMBOR7.5. NOTES SUR L’IMPLÉMENTATION, RÉSULTATS EXPÉRIMENTAUXau résultat de la propagation verticale précédente et nous l’utilisons en tant que valeur initiale pour lapropagation (au bord de l’image, nous utiliserons le schéma de propagation décrit par la fonction mfoldl1qui ne travaille pas avec un tel élément).Après avoir terminé la propagation, nous gardons la valeur du dernier élément B N du stream résultantdans une mémoire temporaire V 2 pour pouvoir la réutiliser lors de la prochaine propagation dans la mêmedirection.Sur le stream B i , nous appliquons la transposition par diagonale, exprimée par le kernel du calcultr2DDiagNxPVecNbf (cf. l’algorithme 6.4) qui utilise, rappelons-le, les fonctions shuffle pour exécuterla transposition en Nlog 2 N opérations. Le stream résultat C i de cette transposition est mis en relationavec la valeur résultante H de la propagation horizontale 1 précédente pour ensuite appliquer unedeuxième phase de propagation verticale des valeurs à l’aide de la fonction mfoldl (les mêmes suppositionsse mettent en place pour le travail sur les bords de l’image).Les résultats de cette propagation C i sont les résultat semi-finaux. Nous prenons la valeur du dernierélément D N et nous la sauvegardons comme H 2 pour l’évaluation suivante du macro bloc voisin. Lesvaleurs de ce stream résultant sont écrites dans la mémoire. Ainsi, nous avons obtenu le résultat de deuxphases de propagation au niveau du macro bloc.Macro bloc d’entréeE 1E 2E 3 E 1E 2E 3E 1E 2E 3 E 1E 2E 3InterconnexionsÉlément exécutifValeur intermédiairestockée localementValeur intermédiaire transmisepar le réseaux d’interconnexionsMacro bloc résultantItération 0 Itération 1Itération 2 Itération 3FIG. 7.11 : Propagation à l’échelle des macro blocs lors du calcul avec plusieurs unités exécutives. Une desvaleurs intermédiaires est stockée localement dans l’unité, la deuxième est transmise par le réseau d’interconnexionà l’unité suivante.Il faut noter que l’implémentation physique de cette approche demande les connaissance de l’architecturecible car les valeurs intermédiaires de la propagation doivent être gérées différemment sur lesarchitectures à un seul fil d’exécution où nous avons besoin d’une mémoire temporaire pour les sauvegardermais qu’elles peuvent être passées à l’unité d’exécution suivante (e.g. sur les architectures àplusieurs fils d’exécution), comme l’illustre la fig. 7.11. Si nous n’avons pas d’architecture capable detravailler dans une telle configuration, l’utilisation de la mémoire temporaire pour stocker les résultatsintermédiaires de la propagation doit être envisagée.7.5 Notes sur l’implémentation, résultats expérimentauxLes implémentations de la fonction distance et des nivellements que nous avons conçues ont été initialementdestinées au processeur SuperH SH-5 qui bénéficie de 64 registres de 64 bits. Ces implémentationsont été incluses dans l’outil MorphoMedia que nous avons développé afin d’obtenir l’indépendance1Après avoir effectué la transposition par diagonale, la propagation est effectuée, à l’échelle du macro bloc, verticalementmais elle correspond à une propagation horizontale dans l’image159