Cette thèse présente les limites du mode s.i.m.d. Dans le cadre de la programmation parallèle d'algorithmes d'algèbre linéaire. Plus précisément, celles de la règle d'or du parallélisme massif: un élément de la matrice par processeur, sont développées. Des expérimentations sont effectuées sur une connection machine 2. Néanmoins, la première partie montre comment la création de procédures de communications écrites a partir d'un nouvel algorithme de construction d'arbres équilibres, et un placement de données judicieux permettent d'atteindre des performances proches de la puissance crête. Mais ce type de travail ne peut pas être effectue sur n'importe quel algorithme, et tout ne s'adapte pas aussi bien. Dans la deuxième partie, nous présentons les avantages de la décomposition en blas pour la construction d'algorithmes massivement parallèles. Elle met, dans le chapitre 4, en évidence la barrière de synchronisation pour la methode du gradient conjugue. Nous proposons dans ce cas particulier comme solution, une ancienne methode qui bien qu'elle soit, en séquentiel, de convergence plus lente, est plus rapide en parallèle. De plus, la structure des matrices est un facteur important. Elle permet d'accélérer les calculs et d'augmenter la dimension des problèmes a résoudre. L'architecture des machines actuelles en limite encore trop l'utilisation. La dernière partie est entièrement consacrée aux permutations, et aux communications qu'elles entrainent. Dans le cadre de l'algorithme de Burg, nous proposons une solution qui calcule a la fois les coefficients de réflexion et ceux d'autoregression sans cout supplémentaire

Pour le calcul d'ombres en informatique graphique, il est courant de s'intéresser à la vue qu'un observateur a d'une scène géométrique. En particulier, il est important de caractériser les changements structurels, appelés événements visuels, qui se produisent dans cette vue lorsque l'observateur se déplace. En se basant sur la définition combinatoire de la vue proposée par Gigus et Malik et la classification des événements visuels qui en découle, de nombreux travaux se heurtent à des problèmes de complexité en temps et en espace. C'est notamment le cas de la méthode du maillage de discontinuités. Nous suggérons donc une approche nouvelle qui repose sur la remise en cause de cette notion de vue.<br /><br />Pour un ensemble d'objets convexes disjoints, nous proposons une définition topologique de la vue qui fait la part belle aux silhouettes visibles des objets de la scène et nous caractérisons géométriquement les lieux où se produisent les événements visuels. Nous utilisons cette caractérisation pour proposer une méthode qui permet d'extraire les limites entre lumière et pénombre et entre ombre et pénombre dans une scène éclairée par des sources surfaciques. Nous arrivons ainsi à réduire considérablement la taille des objets intermédiaires utilisés pour la construction des limites entre les régions.<br /><br />De plus, nous démontrons les premières bornes théoriques non triviales sur la complexité des limites entre lumière et pénombre ainsi qu'entre ombre et pénombre.

Dans la fabrication des produits pétroliers en raffinerie, les lois linéaires de mélange permettent de représenter les ensembles de mélanges faisables par des zonotopes. La faisabilité d'un mélange est un probleme important qui est résolu par des méthodes d'optimisation convexe. Le but du travail présente est de montrer que, dans le cas de la dimension trois, la géométrie algorithmique apporte d'autres solutions a ce probleme. La spécificité des zonoedres et l'utilisation d'une structure de données de type arête-ailée permettent la mise en œuvre d'algorithmes de géométrie optimaux pour les représenter, puis des algorithmes de manipulation et visualisation rapides et robustes destines a être utilises de manière concrète. Le logiciel développe a partir de ces outils apporte une aide efficace dans la décision de la fabrication des gazoles. Dans le cadre plus vaste de la séparation, l'état de séparation d'un système physico-chimique est représente par un zonoide. Les Zonodres fournissent une approche géométrique pour l'étude de tels objets

Cette thèse a été menée en collaboration avec la der-edf. Il s'agit d'étudier ici la parallélisation d'un probleme de répartition de charges dans les réseaux électriques. Ce probleme correspond mathématiquement a la resolution successive de systèmes linéaires dont les matrices sont proches les unes des autres. Une methode originale de resolution est tout d'abord présentée dans un cadre séquentiel, puis une parallélisation sur plusieurs types d'architectures (mimd vectoriel, simd massivement parallèle et mimd a topologie reconfigurable) est proposée. Les machines cibles ont ete étudiées en profondeur et modélisées théoriquement. Plusieurs idées de parallélisation ont été envisagées. Il est intéressant de constater que les meilleures méthodes de resolution de ce probleme concret modèle (c'est-a-dire les plus rapides) sont différentes suivant le type de machine parallèle considéré

Cette dernière décennie a vu apparaître et se développer toute une théorie sur l'échantillonnage des surfaces lisses. L'objectif était de trouver des conditions d'échantillonnage qui assurent une bonne reconstruction d'une surface lisse S à partir d'un sous-ensemble fini E de points de S. Parmi ces conditions, l'une des plus importantes est sans conteste la condition d'e-échantillonnage, introduite par Amenta et Bern, qui stipule que tout point p de S doit être à distance de E au plus e fois lfs(p), où lfs(p) désigne la distance de p à l'axe médian de S. Amenta et Bern ont montré qu'il est possible d'extraire de la triangulation de Delaunay d'un e-échantillon E une surface affine par morceaux qui approxime S du point de vue topologique (isotopie) et géométrique (distance de Hausdorff). Néanmoins restaient ouvertes les questions cruciales de pouvoir vérifier si un ensemble de points donné est un e-échantillon d'une part, et de construire des e-échantillons d'une surface lisse donnée d'autre part. De plus, les conditions d'échantillonnage proposées jusque là n'offraient des garanties que dans le cas lisse, puisque lfs s'annule aux points où la surface n'est pas différentiable. Dans cette thèse, nous introduisons le concept d'e-échantillon lâche, qui peut être vu comme une version faible de la notion d'e-échantillon. L'avantage majeur des e-échantillons lâches sur les e-échantillons classiques est qu'ils sont plus faciles à vérifier et à construire. Plus précisément, vérifier si un ensemble fini de points est un e-échantillon lâche revient à regarder si les rayons d'un nombre fini de boules sont suffisamment petits. Quand la surface S est lisse, nous montrons que les e-échantillons sont des e-échantillons lâches et réciproquement, à condition que e soit suffisamment petit. Il s'ensuit que les e-échantillons lâches offrent les mêmes garanties topologiques et géométriques que les e-échantillons. Nous étendons ensuite nos résultats au cas où la surface échantillonnée est non lisse en introduisant une nouvelle grandeur, appelée rayon Lipschitzien, qui joue un rôle similaire à lfs dans le cas lisse, mais qui s'avère être bien défini et positif sur une plus large classe d'objets. Plus précisément, il caractérise la classe des surfaces Lipschitziennes, qui inclut entre autres toutes les surfaces lisses par morceaux pour lesquelles la variation des normales aux abords des points singuliers n'est pas trop forte. Notre résultat principal est que, si S est une surface Lipschitzienne et E un ensemble fini de points de S tel que tout point de S est à distance de E au plus une fraction du rayon Lipschitzien de S, alors nous obtenons le même type de garanties que dans le cas lisse, à savoir : la triangulation de Delaunay de E restreinte à S est une variété isotope à S et à distance de Hausdorff O(e) de S, à condition que ses facettes ne soient pas trop aplaties. Nous étendons également ce résultat aux échantillons lâches. Enfin, nous donnons des bornes optimales sur la taille de ces échantillons. Afin de montrer l'intérêt pratique des échantillons lâches, nous présentons ensuite un algorithme très simple capable de construire des maillages certifiés de surfaces. Etant donné une surface S compacte, Lipschitzienne et sans bord, et un paramètre positif e, l'algorithme génère un e-échantillon lâche E de S de taille optimale, ainsi qu'un maillage triangulaire extrait de la triangulation de Delaunay de E. Grâce à nos résultats théoriques, nous pouvons garantir que ce maillage triangulaire est une bonne approximation de S, tant sur le plan topologique que géométrique, et ce sous des hypothèses raisonnables sur le paramètre d'entrée e. Un aspect remarquable de l'algorithme est que S n'a besoin d'être connue qu'à travers un oracle capable de détecter les points d'intersection de n'importe quel segment avec la surface. Ceci rend l'algorithme assez générique pour être utilisé dans de nombreux contextes pratiques et sur une large gamme de surfaces. Nous illustrons cette généricité à travers une série d'applications : maillage de surfaces implicites, remaillage de polyèdres, sondage de surfaces inconnues, maillage de volumes.

L'objet de cette thèse est l'étude de la parallélisation d'algorithmes du calcul scientifique et<br />leur implémentation sur des ordinateurs parallèles à mémoire partagée et sur des réseaux systoliques.<br /> Un accent particulier est mis sur l'obtention de résultats de complexité. La thèse est organisée autour<br /> d'articles et textes de conférences qui sont analysés et discutés dans une première partie de façon à <br />permettre de replacer les problèmes traités dans leur contexte.<br />Dans le premier chapitre, nous présentons les principaux résultats théoriques concernant <br />l'étude de complexité des algorithmes parallèles, ainsi qu'une description critique de l'architecture <br />de référence, qui est une machine de type MIMD à mémoire partagée. Le chapitre suivant est dédie" à <br />l'ensemble des résultats de complexité concernant les algorithmes de diagonalisation et <br />l'élimination de Gauss, il a pour but d'illustrer la méthodologie. Il existe en tout dix écritures possibles de la méthode de Gauss, qui conduisent principalement à deux grandes classes de graphes de précédente, conceptuellement différents : les graphes de type "glouton" et ceux du type "2 pas".<br />Ces types de graphes se rencontrent d'une manière plus générale dans d'autres problèmes <br />d'algèbre linéaire et même dans certaines méthodes non numériques de la théorie des graphes. <br />Nous développons les résultats de complexité concernant ces deux types de graphes sur les exemples<br /> les plus courant (versions kji et kij de Gauss en parallèle), puis nous montrons comment adapter<br /> l'étude en prenant en compte t'es temps de communication entre tes processeurs, ce qui rend le modèle<br /> théorique plus réaliste.<br />Le chapitre 6 est consacré aux architectures systoliques. Le problème du chemin algébrique permet <br />d'unifier plusieurs problèmes informatiques. Nous présentons un réseau résolvant ce problème en Sn-2 <br />pas sur un réseau de taille n(n+l ). De plus, quelques modifications permettent de calculer des projections<br /> en filtrage adaptatif en vu d'obtenir une solution en temps réel pour le traitement numérique des signaux.<br />Avant de conclure, nous présentons des résultats complémentaires de parallélisation effective sur d'autres<br /> types d'architectures : l'étude de l'algorithme du gradient conjugué sur des super calculateurs <br />(CRAY-XMP et IBM 3090-VF).

Les quatres chapîtres de cette thèse aborde quatre thèmes de la théorie des itérations: 1) nous élaborons un algorithme de vérification de l'attraction d'un point fixe d'une itération discrète dans son voisinage second. Cet algorithme est comparé aux conditions nécessaires et suffisantes énoncées par F. Robert avant d'être généralisé à d'autres attracteurs et d'autres bassins d'attraction. 2) Après un tour d'horizon des méthodes de calcul de racines pième de matrices réelles nous proposons un algorithme de calcul de racines carrées de matrices booléennes quelconques. 3) Nous utilisons un opérateur monotone pour étudier les itérations bloc-séquentielles de réseaux à seuil: on caractérise ainsi leurs dynamiques. Nous étendons ces méthodes aux fonctions majorité et verres de spin généralisés. 4) Après avoir comparé les différents outils d'observation des dynamiques des réseaux booléens aléatoires d'interconnectivité 2, nous proposons une approche basée sur le calcul d'une approximation de chacune des 3 composantes: le coeur stable du réseau, le coeur oscillant, les paliers (notion introduite ici). En application nous nous intéressons au problème de la reconnaissance de séquences booléennes par ce type de réseaux

.

Les travaux présentés dans cette thèse ont pour but d'améliorer la couche transport de l'architecture réseau de l'OSI. La couche transport est de nos jour dominée par l'utilisation de TCP et son contrôle de congestion. Récemment de nouveaux mécanismes de contrôle de congestion ont été proposés. Parmi eux TCP Friendly Rate Control (TFRC) semble être le plus abouti. Cependant, tout comme TCP, ce mécanisme ne prend pas en compte ni les évolutions du réseau ni les nouveaux besoins des applications. La première contribution de cette thèse consiste en une spécialisation de TFRC afin d'obtenir un protocole de transport avisé de la Qualité de Service (QdS) spécialement défini pour des réseaux à QdS offrant une garantie de bande passante. Ce protocole combine un mécanisme de contrôle de congestion orienté QdS qui prend en compte la réservation de bande passante au niveau réseau, avec un service de fiabilité totale afin de proposer un service similaire à TCP. Le résultat de cette composition constitue le premier protocole de transport adapté à des réseau à garantie de bande passante. En même temps que cette expansion de service au niveau réseau, de nouvelles technologies ont été proposées et déployées au niveau physique. Ces nouvelles technologies sont caractérisées par leur affranchissement de support filaire et la mobilité des systèmes terminaux. De plus, elles sont généralement déployées sur des entités où la puissance de calcul et la disponibilité mémoire sont inférieures à celles des ordinateurs personnels. La deuxième contribution de cette thèse est la proposition d'une adaptation de TFRC à ces entités via la proposition d'une version allégée du récepteur. Cette version a été implémentée, évaluée quantitativement et ses nombreux avantages et contributions ont été démontrés par rapport à TFRC. Enfin, nous proposons une optimisation des implémentations actuelles de TFRC. Cette optimisation propose tout d'abord un nouvel algorithme pour l'initialisation du récepteur basé s ur l'utilisation de l'algorithme de Newton. Nous proposons aussi l'introduction d'un outil nous permettant d'étudier plus en détails la manière dont est calculé le taux de perte du côté récepteur.

Développement d'une notation algorithmique adaptée à une démarche de programmation systématique et validation de l'emploi de méthodes précises d'analyse et de techniques de programmation au travers d'une expérience significative. La première partie de la thèse propose une réflexion sur les concepts sous-jacents à la maitrise des diverses constructions algorithmiques de base. Tandis que la deuxième partie aborde la spécification des divers composants d'un environnement de programmation