Brevet FR2647914A1 — publié le 19901207

L'invention a pour objet un calculateur optoélectronique systolique constitué par un réseau pipeline de multiplieurs ootiques mis en cascade à l'intérieur desquels un faisceau lumineux se propage pour ezfec- tuer une pLur-alite d'additions et de multiplications. Elle trouve une application dans les calculs ultra-rapides parallèles de type algèbre linéaire (produits matrice x vecteur, inversion de matrice, etc...) tels que ceux habituellement utilisés en traitement numérique d'images, en traitement du signal, dans les opérations de détection ou d'identification L'un des objectifs recherchés en calcul numérique lorsque le nombre de données à traiter est important est l'augmentation des cadences de calcul, celLes-ci pouvant atteindre plusieurs gigaopérations par seconde. Les calculateurs classiques de type von Neuman sont trop lents du fait de leur fonctionnement séquentiel dans le temps qui empeche tout paraltètisme de calcul. Une solution est constituée par tes architectures systoliques dans lesquelles l'unité centrale de calcul est remplacée par un réseau d'unités de calcut élémentaires pouvant effectuer en paral- Lèle des opérations simples et répétitives sur un flot de données qui s'écoulent dans le réseau à chaque cycle de calcul. Le gain en temps de calcul est alors proportionnel au nombre d'unités de calcul élémentaires formant le réseau. Les réseaux systoliques sont des processeurs spécial i- ses dans les calculs d'algebre Linéaire, tels par exemple Les calculs de produits matrice x vecteur A.x = y. Ceux-ci font intervenir une succession de multiplications étémentaires répétitives entre les éléments aij de la matrice A et les composantes x. du vecteur x, suivies d'additions élémentaires donnant les composantes yi du vecteur Les produits matrice x vecteur sont à la base de calculs pLus élaborés, tels que produits matrice x matrice, calcul des va Leurs et des vecteurs propres d'une matrice, inversion de matrices et constituent une part importante des calculs en traitement du signal ou d'images. Ainsi l'accès au contour d'un objet dans les systèmes de vision pour La robotique s'effectue par le produit de convolution de l'objet qui s'exprime, en calcul numérique, sous la forme d'un produit matrice x vecteur. En détection radar, L'extraction d'un signal noyé dans du bruit et l'améLioration du rapport signal-sur-bruit du signal détecté peut être effectuée par une opération de corrélation. Ce(Le-ci s'effectue ici encore sous forme d'un produit entre une matrice de Toeplitz et un vecteur colonne. De la même manière, l'évaluation de structures mécaniques en avionique par exemple fait souvent intervenir des systèmes d'équations linéaires de grande dimension dont la résolu- tion nécessite des calculs matriciels qui impliquent des durées de calcul très Longues. Les calculateurs systoliques permettent de réduire les durées de calcul de plusieurs ordres de grandeur. De nouveaux algorithmes de calcul adaptés au calcul parallèle font actuellement l'objet de recherches dans les laboratoires afin d'accroitre encore davantage la rapidité de calcul et le champ d'application de ces calculateurs. A coté des calculateurs systoliques tout électroniques, des cal cul teurs systoliques opto-électroniques sont actuellement développés. Les calculs s'effectuent alors optiquement au moyen de faisceaux lumineux. Ceci représente l'avantage potentiel de réduire Les consommations d'énergie, de permettre des cadences de calcul au moins comparables à celles des calculateurs systoliques tout électroniques (106 à 108 additions/multiplications par seconde, en fonction du nombre de bits de précision) et d'obtenir un prix de revient beaucoup plus faible. Ce dernier point est d'importance par exemple dans les systèmes embarqués. Les processeurs systoliques optoélectroniques connus de l'homme de L'art font généralement intervenir des effets acousto-optiques pour effectuer des multiplications éLémentaires. Le flot de données est alors constitué par des ondes ultrasonores se propageant dans un cristal acoustooptique éclairé par une pluralité de sources lumineuses (diodes laser ou diodes électroluminescentes). Par exemple, Caulfield et al (Optics Communicstions, 40, 86, 1981), Casasent (App2ied Optics,21, 1859-1865, 1982), Tsai et al (App2ied Physics Letters,47,6, 549-551, 1985) décrivent par exemple des processeurs analogiques dans Lesquels le nombre de sources lumineuses est égal au nombre de colonnes de la matrice A, ou au nombre de diagonales non nulles lorsque la matrice A est une matrice à structure de bande. Rhodes et al (Proc,IEEE,72,7, 820-829, 1984), et Guitfoyle(Optical Engineer*ng,23, 23, 1, 20-25, 1984) décrivent d'autres architectures, de type numérique, dans lesquelles une pluralité de pixels lumineux générés par exemple par un premier cristal acousto-optique éclairent un second cristal acousto-optique pour effectuer des produits matrice x vecteur seLon L'algorithme de convolution imaginé par Whitehouse et Speiser (Aspects cf Signez Processirg, vol.2, G.Tacconi Ed, Reidel, Hingham, Mass. 1977). Un désavantage important de ces procédés de calcul est de nécessiter une pluralité de sources dont le nombre est d'autant plus élevé que La dimension des matrices est grandie . Ceci conduit à La mise en oeuvre de circuits de contrôle et de régulation de La puissance d'émission de chaque source et donc a une complexité du système. De plus ceci conduit à un manque de fiabilité, les sources pouvant présenter des dérives importantes en fonction de la température ou du vieillissement, en particulier pour les diodes laser. Un autre désavantage réside en ce que ces procédés nécessitent une synchronisation précise de chaque source avec les données se oropageant dans le cristal acousto-optique. De plus, celui-ci introduit une atténuation acoustique sur le flot de données qui s'y propage et doit donc être pris en compte pour calibrer La puissance d'émission de chacune des sources. Ceci peut expliquer que, bien que plusieurs études de principe aient été rapportées dans La littérature,seuls quelques rares prototypes sont actuellement expérimentés dans le monde (par exemple, Casasent et a!, AppZied Optics, 25 14, 2258-2262, 1986)) La présente invention porte sur une nouvel Le architecture de calculateur optoélectronique systolique. L'un de ses objets est de permettre l'utilisation d'une source lumineuse unique tout en gardant la possibilité d'un parallèlisme de calcul, et ainsi de simplifier de manière considérable les conditions de synchronisation, les circuits d'entrée des données et, de façon générale, l'ensemble du hardware, avec, pour avantage, une augmentation de la fiabilité de fonctionnement. Un autre objet de l'invention est de constituer un moyen d'obtenir des vitesses de calcul très élevées en codant l'ensemble des données à traiter sur un faisceau lumineux unique le long duquel les opérations s'effectuent à La vitesse de La Lumière à mesure que la lumière se propage dans Le réseau de calcul. Un autre objet de l'invention est de constituer un moyen d'augmenter Les possibilités de calcul parallèle actuelles, l'invention étant compatible avec l'ensemble des autres systèmes systoliques optiques faisant intervenir par exemple des cellules acousto-optiques ou magnéto-optiques, ou de façon générale des valves optiques fonctionnant en lumière monochromatique ou polychromatique. Le domaine d'application de l'invention recouvre de façon générale le calcul rapide faisant intervenir des algorithmes matriciels, avec des cadences d'horloge pouvant atteindre plusieurs centaines de MHz. Durant un cycle d'horloge, la cadence de calcul peut atteindre 108 à 10 additions/multiplications par seconde. La précision de calcul peut atteindre une dizaine de bits ou plus selon le type de sources lumineuses Le principe de base del'invention repose sur L'utilisation d'un réseau de calcul formé 4'une pluralité de multiplieurs optiques mis en cascade et à L'intérieur desquels se propage un faisceau lumineux préalablement modulé en intensité Chaque multiplieur effectue en temps réel (égal au temps de propagation de la lumière dans le multiplieur, de l'ordre de 10-10 seconde) une multiplication entre son signal de commande et L'énergie du faisceau lumineux le traversant. Plus précisément, chaque aultiplieur code son signal de commande sous forme d'un retard optique D supérieur à la Longueur de cohérence L = AZ/dS du faisceau lumineux ( A: longueur d'onde, dA : largeur spectrale du rayonnement). Il en découle que le résultat de la multiplication est codé sous forme d'une modulation de la coherence temporelle du faisceau Lumineux, non détectable de façon directe. A mesure que le faisceau lumineux se propage dans le réseau de calcul, les différentes multiplications sont ainsi effectuées de proche en proche et codées sur Le faisceau lumineux. En sortie du réseau de multirLieurs, L'accès aux résultats s'effectue au moyen d'une pluralité de décodeurs mis en parallèle et introduisant des retards optiques D sensiblement équiva Lents à ceux des multiplieurs. Le résultat de chaque multiplication élé- mèntaire est alors obtenu à La sortie de chaque décodeur sous forme d'une modulation d'intensité Lumineuse. Ceci a pour effet de rendre possible une détection directe de chaque multiplication éLémentaire au moyen d'une photodiode. Les multiplieurs évoqués plus haut, et dont il sera question dans la suite, recouvrent tout moyen optique dont le fonctionnement peut se ramener à celui d'un modulateur (ou d'un filtre spectral accordable ou d'un interféromètre) introduisant un retard optique supérieur à la longueur de cohérence de la Lumière et tel qu'aucune modulation de L'intensité Lumineuse de sortie n'apparaisse sous l'effet de l'accord du filtre ou de l'interféromètre. Les multiplieurs peuvent par exemple etre des interféromètres à deux ondes ou à ondes multiples, des boucles optiques, des coupleurs directifs, des modulateurs électro-, acousto-, magnéto-, piézo- ou élasto-optiques ou tout autre moyen faisant intervenir des effets électriques, acoustiques, magnétiques ou mécaniques pour induire et moduler un retard optique ou induire un filtrage spectral des longueurs d'onde. Les différents objets et caractéristiques de l'invention apparaitront mieux à la lumière de la description qui va suivre, faite à titre d'exemple non limitatif. Cette description se préfère à des dessins annexés sur lesquels : - la figure 1 montre le principe de base d'un multiplieur fonctionnant par modulation de cohérence. - La figure 2 représente un produit matrice x vecteur tel qu'il peut etre calculé par le réseau systolique de la figure 3, et dans lequel Les composantes sont analogiques. - la figure 3 illustre le principe de fonctionnement du cal cul teur systolique et la méthode utilisée pour calculer le produit matrice x vecteur de la figure 2. - La figure 4 montre les séquences de calcul et la position des différentes multiplications élémentaires à chaque étape du calcul du produit matrice x vecteur de la figure 2. - la figure 5 représente un exemple de réalisation du réseau systolique de la figure 3 à partir de composants discrets disponibles sur le marché. - La figure 6 montre une variante réalisée à partir de composants intégrés et de fibres optiques. - la figure 7 représente un produit de matrice x vecteur dans lequel les composantes sont binaires. - la figure 8 illustre un système systolique numérique selon la presente invention permettant le calcul du produit matrice x vecteur binaire de la figure 7 avec une précision digitale. - La figure 9 représente Les séquences de calcul conduisant au calcul du produit matrice x vecteur de la figure 7 au moyen du processeur numérique de la figure 8. - la figure 10 décrit un mode de réalisation du processeur de la figure 8. - la figure Il montre un processeur numérique permettant de calculer le produit matrice x vecteur de La figure 7. - La figure 12 représente un processeur d'architecture parallèle. - la figure 13 donne un exemple des séqueoces de calcul du processeur de la figure 12 pour le calcul du produit matrice x vecteur de La figure 2. La figure 1 donne Le principe d'un multiplieur permettant d'effectuer des multiplications par moduLation de cohérence d'un faisceau Lumineux La figure la représente La densité spectrale d'énergie B(o) de La source Lumineuse 1 8(a) = 1 pour a 1 = 0 ailleurs La variable o est un nombre d'ondes. La longueur d'onde moyenne émise par La source est A . Cet exemple n'est pris qu'à titre o illustratif pour simplifier les notations qui vont suivre, des sources présentant d'autres densités spectrales d'énergie (gaussiennes, lorentziennes, à structure de modes, etc...) pouvant être utilisées sans inconvénient. La figure lb montre te cas simple du calcul d'un produit unique a1 b1. Le système comprend une source lumineuse 1 modulée en intensité par a1, un moyen de modulation i A sur lequel est adressée La donnée b1, un moyen de démodulation 2A, un photodétecteur 4 et un circuit de mise en forme 5. L'énergie Lumineuse émise par la source 1 est proportionnelle à a1 I1 - a1 (1) Les moyens de modulation 1A et de démodulation 2A sont constitués par tout moyen introduisant une loi de filtrage spectral à maximums multiples sur Le faisceau 2 Les traversant. Par exemple, les éléments 1A et 2A peuvent etre des dispositifs dont Le comportement se ramène à celui d'un interféromètre de Fabry-Pérot, ou à celui d'un interfero- mètre de Michelson sans que cela puisse exclure tout autre système dont la loi de filtrage spectral en fonction de la longueur d'onde ou des fréquences temporel Les de la lumière est sinusoidale ou, de façon générale, présente des maximums multiples. A titre d'illustration l'ensemble des descriptions dans la demande d'invention est fait à partir de dispositifs dont Le fonctionnement peut etre ramené à celui d'interféromètres à deux ondes introduisant des retards optiques, c'està-dire pour Lesquels La Loi de filtrage spectral est de type sinusoidal. Suivant cette description, le modulateur 1A introduit sur Le faisceau lumineux 2 le traversant, un retard optique D(A) tel que D(A) = DA + b1 ( L) (2) ou DA est une composante continue modulée par une seconde composante proportionneLle à b1. Le retard optique D(A) ainsi introduit sur Le faisceau 3 est toujours supérieur à la longueur de cohérence L : #2/ | #2 - #1| = 1/| c2 - c1| de la source 1, L'énergie lumineuse du faisceau 3 résulte alors d'un filtrage spectral sinusoidal, et a pour expression IL en résulte qu'à la sortie du modulateur 1A, L'énergie Lumineuse I2 n'est pas modulée par b Le retard optique D(A) étant suoérieur à la Longueur de cohérence et ne pouvant donc donner lieu à un phénomène d'interférence. A La sortie du modulateur 1A, la donnée al est alors codée sous forme d'une modulation d'intensité lumineuse du faisceau 3 tandis que La donnée b2 est codée sous forme d'une modulation de retard optique à l'intérieur du même faisceau. Toujours suivant La même description, Le démodulateur 2 A est constitué par un moyen introduisant un retard optique D'(A) sensiblement équivalent à D(A). De façon plus précise et dans une configuration préférentielle adaptée à des données ana Logiques ou multiniveaux, on a sensiblement : D'(A) = D - A /4 (4) o La condition (4) assure Les meilleures conditions de linéarité du démodulateur 2A lorsque les données sont analogiques. Dans le cas de données binaires, la Linéarité est moins critique. La dynamique est alors optimale quand on a sensiblement D'(A) = DA (5) L'énergie lumineuse 13 détectée par Le photodétecteur 4 est alors : En introduisant la relation (4) dans (6), l'expression (6) se réduit à : I3 = 1/4 {a1 + K a1b1} (7) où K est une constante (K = 2##c). Le deuxième terme de l'expression (7) donne le produit a1b1. Celui-ci est obtenu en envoyant Le signal électrique délivré par Le photodétecteur 4, proportionnel à I3, dans le circuit de mise en forme 5 sur lequel est appliquée la donnée a1. Ce circuit est constitué par tout moyen permettant de soustraire le permier terme de l'expression (7) de manière à donner, à sa sortie, le produit a1b1, La figure 1 c étend Le principe de calcul précédent au cas du calcul simultané de deux produits a1b1 et a1b2. La donnée a1 est adressée sur La source lumineuse 1 et sur les deux circuits de mise en forme 5 et 8. Les données b1 et b2 sont adressées sur les modulateurs 1 A et 1 a respectivement. Ceux-ci sont constitués par des moyens introduisant des retards optiques différents D(A) = DA + b1 ( L) (8) D(9) = DB + b2 ( L) (9) Les retards optiques D(A) et D(B) sont supérieurs à la longueur de cohérence de la source 1. L'intensité lumineuse 12 du faisceau lumineux 3 est donnée par l'expression (3). L'intensité lumineuse 14 du faisceau lumineux 6 s'écrit = 1/4 a1 (11) le faisceau 6 est divisé en deux faisceaux qui traversent Les démodulateurs 2A et 28. Ceux-ci introduisent des retards optiques D'(A) et D'(8) sensiblement équivalents à ceux introduits par Les modulateurs 1A et 1B. En reprenant un raisonnement sembLable à celui donne à la figure lb, on peut montrer que L'énergie lumineuse détectée par Les photodétecteurs 4 et 7 est proportionnelle à -sur Le photodétecteur 4 : 16 { a1 + K a1b1 J (12) - sur le photodétecteur7 : 16 { a1 + K a1b2} (13) Après mise en forme des signaux par les éléments 5 et 8, le processeur délivre les produits a1b1 et a1b2 sous forme électrique sur chacune de ses deux sorties. Les produits a1b1 et a1b2 ainsi obtenus peuvent ensuite être stockés dans des mémoires par les moyens habituels connus de l'homme de l'art. il est à noter que la durée de calcul est déterminée par le temps de propagation de La Lumière dans les éléments 1A, 1B, 2A et 2B, Celui-ci est généralement négligeable devant les constantes de temps des circuits d'adressage des données et devant la durée des cycles d'horloge. L'adressage des données a b1,b2 s'effectue donc en mode parallèle sur Les éléments 1,1A, 18, 5 et 8 sans avoir à introduire d'éventuels retards de temporisation tant que le temps de propagation de la lumière peut être considéré comme négligeable. Le résultat des produits éLémentaires est obtenu lui aussi en mode parallèle. Un tel système peut effectuer des calculs de produits matrice x vecteur quand une pluralité de multiplieurs sont connectés entre eux. La figure 2 représente le produit d'une matrice 9 par un vecteur 10. Le résultat du produit matrice x vecteur est Le vecteur 11. La matrice 9 comprend une pluralité de lignes (p) et de colonnes (Q) formées de composantes a... Le vecteur d'entrée 10 est forme des don il nées (x1,x2,...xQ). Le vecteur de sortie 11 comprend les composantes (y1,y2,...yp). Celles-ci sont obtenues nar addition de multiplications élémentaires aijxi entre les composantes de la matrice 9 et celles du vecteur d'entrée 10 Le produit matrice x vecteur décrit par La figure 2 peut etre calculé en adressant séquentiellement dans le temps les éléments a.. et x. à une pluralité de multiplieurs mis en série selon L'architecture J de La figure 3. La figure 3 montre L'architecture et le principe de fonctionnement d'un calculateur systolique adapté au calcul du produit matrice x vecteur de La figure 2. Le système est formé d'une pluralité de modulateurs 1A, 18,...1P mis en série et dont Le nombre est égal au nombre P de lignes de la matrice. Une pluralité (P) de démodulateurs 2A, 2B... 2P sont mis en parallèle et fonctionnent suivant Le principe exposé à la figure 1. Le réseau de calcul est éclairé par une source Lumineuse 12 de longueur de cohérence L. A chaque séquence de calcuL, une composante xj du vecteur d'entrée 10 et les P composantes aij de la colonne j de la matrice 9 sont appliquéesen mode parallèle sur la source 12 et sur chacun des modulateurs. Les différentes séquences de calcul sont indiquées à la figure 3, ainsi que L'entrée des données. Une séquence de calcul n j permet de calculer simultanément P produits élèmentaires aijxj, a2jxj, a3jxj,... apjxj. Ceux-ci sont obtenus sur chacune des sorties des circuits de mise en forme 3A, 3B,... 3P et sont stockés et additionnéS dans Les registres 4A, 4B,...4P, A La séquence de calcul suivante, on obtient P nouveaux produits élémentaires qui s'additionnent aux résultats précédents. Aprés Q séquences de calcul, on obtient finalement dans Les registres 4A, 4B,..4P Les composantes y1,y2,...yp du vecteur de sortie 1'. Une meilleure compréhension des séquences de calcul de la figure 3 peut être obtenue en se référant en même temps à la figure 4. La figure 4 montre les séquences de calcul et la position des données d'entrée et des multiplications élémentaires à chaque étape du calcul du produit matrice x vecteur de la figure 2. Dans La figure 4a, qui représente La séquence de calcul n 1, La composante x1 est envoyée dans le réseau de calcul (sous forme d'une modulation d'intensité lumineuse) en même temps que les composantes a11, a21,..., a11,...ap1 de la colonne n 1 de la matrice 9. A la sortie du processeur, on obtient simultanément les multiplications élémentaires a11x1,a21x1,...aj1x1,..., aP1x1, chacune étant alors stockée dans des registres additionneurs 4A, 48, 4C,.... 4P. Dans La figure 4b correspondant à la séquence de calcul n 2, Les composantes a12,a22, a32 , ai2,.. a P2 de la colonne n 2 de la matrice 9 et la composante x2 du vecteur d'entrée 10 sont adressées en mode parallèle sur Le réseau de calcul. Le système calcule alors en parallèle les produits élémentaires a12x2, a22x2,....ai2x2,....aP2x2. Ceux-ci s'additionnent aux produits élémentaires du cycle précédent, de sorte que Le contenu des circuits d'addition et de stockage devient : y1 = a11x1 + a12x2 y2 = a21x1 + a22x2 (15) yi = ai1x1 + ai2x2 yP = aP1x1 + aP2x2 La troisième séquence de calcul, représentée à la figure 4c, fournit un autre ensemble de produits élémentaires qui s'ajoutent, dans chaque circuit de stockage et d'addition, aux résultats obtenus dans Les cycles de calcul précédents. Les P circuits de stockage contiennent alors, respectivement : y1 = a11x1 + a12x2 + a13x3 y2 = a21x1 + a22x2 + a23x3 (16) yi = ai1x1 + ai2x2 + ai3x3 yP = aP1x1 + aP2x2 + aP3x3 Les cycles de calcul continuent de cette manière jusqu'à ce que le produit matrice x vecteur soit terminé. Les données stockées dans Les P additionneurs sont alors y1 = a11x1 + a12x2 + .... + a1QxQ y2 = a21x1 + a22x2 + .....+ a2QxQ (17) y@ = a@@x@ + a@@x@ + .....+ a@ @ Q Les données (y1,y2,...yP) exprimées par Les relations (17) sont les conposantes du produit matrice x vecteur décrites par ta relation (14). On constate que le nombre de cycles de calcul est égal au nombre Q de coLonnes de la matrice et que te nombre de modulateurs est égal au nombre P de Lignes de la matrice. Au cours d'une séquence de calcul donnée, les produits élémentaires sont obtenus par croisement des données appliquées sur les modulateurs 1A, 1B, 1C.... 1P avec cette appliquée sur La source 12 et se propageant sous forme d'un faisceau lumineux dans les modulateurs. Théoriquement l'adressage de ces données doit être synchronisé pour prendre en compte le temps de propagation du faisceau Lumineux d'un modulateur au suivant. En pratique ce temps de propagation étant négligeable (10 12 à seconde) devant les constantes de temps des circuits d'adressage, les données peuvent être appliquées simultanément à ta source 12 et aux modulateurs 1A, 1B, iC. ..en mode parallèle. Les conditions de synchronisation assurant le fonctionnement du calculateur sont donc en général très faciles à remplir. De la même manière, l'utilisation d'une source lumineuse 12 unique pour effectuer en parallèle plusieurs multiplications élémentaires procure une fiabilité élevée et évite les problèmes de calibrage et d'étalonnage de sources. Le fait d'utiliser une architecture série mettant en oeuvre une pluralité de modulateurs mis en cascade évite également les problèmes de calibrage des modulateurs dûs à une éventuelle dispersion de leurs caractéristiques de fonctionnement et de leurs pertes d'insertion. Ainsi l'atte- nuation de la lumière à mesure qu'elle se propage dans la pluralité de modulateurs n'est pcs à être prise en compte au niveau de chaque démodu- lateur et ne nécessite donc pas de correction, l'intensité lumineuse du faisceau à La sortie des modulateurs étant déterminée par tes pertes de propagation dans l'ensemble du réseau de modulateurs. Ces particularités du système faisant l'objet de l'invention conduisent à une simplification des circuits de commande et d'horloge en réduisant au naximum les circuits d'asservissement et de régulation généralement nécessaires dans les processeurs optoélectroniques pour le contra le des sources et des systèmes de détection. Les conditions clefs assurant le fonctionnement de ce type de calculateur optoélectronique sont définies par la valeur des retards optiques D(A), D(B), D(C),.... caractérisant les modulateurs 1A, lB, 1C, La première condition est que les retards optiques doivent etre supérieurs à la longueur de cohérence L de La source lumineuse 12. Ainsi pour une source lumineuse de Longueur d'onde À et de largeur de raie d#, cette condition est remplie si Les retards optiques sont supérieurs à L = #2/d# . La deuxième condition est que L'écart entre les différents retards optiques soit choisi pour éviter toute diaphonie entre les étages de calcul fonctionnant en parallèle. Ainsi une séquence possible de retards optiques est donnée par la progression pL, 3pL, 8 pL,...connue dans ce domaine (Journal of Light Wawe Technology Vol.LT-3 n 5 Oct.85 p. 1065) A titre d'exemple, pour une diode superluminescente tel Le que A = 800 nm et ## = 40 nm et un taux de diaphonie de -20dB optique obtenu avec p = 2, les retards optiques sont de L'ordre D(A) = 32 microns, D(B) = 96 microns, D(C) = 256 microns, etc... La source 12 peut être une diode Laser monomode ou multimode, une diode électroluminescente ou une diode superluminescente modulée en intensité. En principe une source à Large spectre d'émission est recommandée, de manière à obtenir une faible longueur de cohérence L et par conséquent desretards optiques D(A), D(8), D(C),.... facilement réalisables. Les modulateurs 1A, 18, 1C, ... sont constitués par tout moyen conçu pour introduire des retards optiques D(A), DtB), D(C).... supérieurs à la longueur de cohérence L. La modulation du retard optique peut Les séquences de calcut données dans les figures 3 et 4 s'appliquent au calcul de produits matrice x vecteur dans lesquels la matrice est de dimension (P x Q) et dont Les éléments sont des données analogiques. Dans ce cas, Le réseau de calcul comporte P modulateurs asscciés à P démodulateurs et le nombre de cycles de calcul est Q. Le processeur calcule donc Les P composantes du produit y en 9 unités de temps déterminées par La période d'horloge à laquelle les données sont adressées dans le système. Par conséquent la cadence de calcul est de P additions-multiplications par période d'horloge. La figure 5 montre un exemple de réalisation du calculateur de la figure 3 dont te fonctionnement a été vérifié expérimentalement. Chaque modulateur 1A, 18 1C comprend un modulateur électro-optique 15, 18, 21 à champ électrique transverse, une Lame biréfringente 14, 17, 20 et un polariseur 13, 16, 19, 22. Les directions de polarisation des poLariseurs sont parallèles entre elles et orientées à 450 des lignes neutres des lames biréfringentes et des modulateurs électro-optiques. Les lames biréfringentes 14, 17, 20 sont des lames de quartz ou de spath de coupe 0 Z et d'épaisseurs différentes. Elles introduisent des retards optiques DA,DB et DC supérieurs à la longueur de cohérence L de la source. Celle ci est une diode Laser multimode émettant à A = 800 nm, et de longueur o de cohérence L = 80 microns. Les retards optiques ont été choisis de manière à vérifier la condition de fonctionnement évoquée auparavant : DA = 160 microns, D8 = 480 microns et DC = 1280 microns. Après traversée des trois modulateurs, Le faisceau Lumineux est envoye sur trois démodulateurs au moyen des lames séparatrices 23, 24 et du miroir 25. Chaque démodulateur comprend une Lame biréfringente 26, 28, 30 associée à un polariseur 27, 29, 31. Les polariseurs 27, 29, 31 sont orientés parallèlement aux polariseurs 13, 16, 19, 22 alors que les lignes neutres des Lames biréfringentes 26, 28, 30 sont orientees à 450 des directions de polarisation des polariseurs 27, 29, 31. Les lames biréfringentes 26, 28, 30 introduisent des retards optiques DA1, DB et DC1 identiques à O, Da et O . L'introduction d'un écart de retard optique de - # A0/4 peut etre réalisé électriquement en appliquant une ten- sion quart d'onde sur chacun des moduLateurs 15, 18, 21. Les signaux lumineux obtenus à La sortie des décodeurs sont détectés par les photodétecteurs 32, 33, 34 reliés aux circuits 35, 36, 37 de mise en forme des signaux. Les éléments 38, 39, 40 sont des circuits de maintien et d'addition dans Lesquels les produits élémentaires sont stockés et additionnés à chaque séquence de calcul. Les expériences de contrâle ont été menées avec des cadences d'horloge de 2 MHz, sur des matrices à 3 lignes et 2800 colonnes. La précision de calcul fournie par le processeur expérimenté était de 3%. La figure 6 donne une autre version du processeur adapté au calcul du produit matrice x vecteur de la figure 2 réalisée par la technologie de L'optique intégrée. La source 12 est couplée à des modulateurs électro-optiques intégrés 42, 44, 46 reliés par Les fibres 41, 43, 45, 47. La Lumière transmise par Les moduLateurs est injectée par l'intermédiaire d'un coupLeur 48 et des fibres 49, 50, 51 dans Les démodulateurs. Ceux-ci peuvent etre constitués par des modulateurs 52, 53, 54 analogues aux composants 42, 44, 46 c'est-à-dire introduisant des retards optiques sensiblement équivalents. L'architecture donnée dans Les figures 3, 4, 5 et 6 est adaptée au calcul de produits matrice x vecteur dans lesquels Les données sont analogiques. La cadence de calcul peut être très élevée mais La précision est relativement faible (quelques pour cent) par rapport a celle obtenue avec un calculateur numérique. On va maintenant décrire deux modes de réalisation du procédé de calcuL adapté à des données numériques, et en particulier binaires bien que le procédé de calcul s'applique tout aussi bien à des données codées dans des bases supérieures à 2. La figure 7 représente un produit matrice (55) x vecteur (56) dans lequel Les composantes sont codées sur n = 5 bits. Le résultat est mis, dans Le vecteur 57, sous forme binaire mixte (les composantes sont alors écrites dans une base supérieure à 2 ; par exemple a11x1 = 111 x 11 = 1221). Le vecteur 58 donne Le résultat en format binaire t par exemple 1221 = 1û1O1 en binaire). La matrice 55 est une matrice à P Lignes et Q coLonnes. La figure 8 montre l'architecture d'un processeur optonumérique adapté au caLcul du produit matrice x vecteur de la figure 7. Le processeur comprend P moduLateurs 1A, 18,... mis en série et P démodulateurs 2A, 28,... dont le fonctionnement a déjà été décrit. Les modulateurs sont éclairés par la source Lumineuse 59 qui fonctionne en régime continu. L'élément 60 est constitué.par tout moyen pouvant se comporter comme un imageur, c'est-à.dire qui permet d'afficher sous forme optique plusieurs données électriques qui lui sont appliquées et de Les transférer d'un point à un autre le long de la direction 64. Par exemple l'élément 60 peut être un modulateur acoustooptique, ou une batterie de cellules électro- ou magnéto-optiques mises cote à côte et commandées par des registres à décalage assurant le transfert des données d'une cellule à la voisine. Les éléments 62 et 63 sont des barrettes de photodétecteurs reliés à des circuits de mise en forme et de stockage des signaux 65,66. L'élément 61 est constitué par tout système optique permettant de former L'image du composant 60 sur Les eléments 62 et 63. Par exemple, il peut être constitué par des Lentilles sphérocylindriques associées à des lames sénaratrices, ou ave une nappe de coupleurs à fibres optiques. L'algorithme utilisé pour calculer le produit matrice x vecteur de la figure 7 est l'algorithme de multiplication par convolution décrit par Whitehouse et Speiser ("Aspects of Signal Processing, Editor : D.Reidel, Boston MA, 1976, pp. 669-702). La compréhension des séquences de calcul de la figure 8 peut être obtenue en se référant en même temps à la figure 9. La figure 9 montre la position des bits lors des n = 5 séquences de calcul permes tant d'obtenir les produits élémentaires a11x1 et a21x1 du produit matrice x vecteur de la figure 8. Dans la figure 9a, qui représente la séquence de caLul n 1, Les bits maximum (MSB) de a11, a21, a31 sont appliqués sur les modulateurs 1A,1B, 1C,... tandis que Le bit minimum CLSB) de x1 est adressé sur l'élément 64. Le photodétecteur MSB de L'élément 62 détecte Le produit entre le bit MSB de a11 et le bit LSB appliqué sur L'élément 60. De La même manière, Le photodétecteur MSB de l'élément 63 détecte le produit entre Le bit MSB de a21 et le bit LSS de x1 Ces produits sont stockés dans les mémoires 65 et 66. A la séquence de calcul n02 (figure 9b) l'élément 60 contient Les deux premiers bits 1 et 1 de x1, le bit LS8 ayant été décalé vers le haut, Le long de l'axe 64. Les deux photodétecteurs correspon- dants de l'élément 62 détectent Le produit entre le bit O appliqué sur Le modulateur lA et chacun des deux bits présents dans l'élément 60. La situation est semblable pour les deux photodétecteurs de L'élément 63. Ceux-ci détectent les produits entre le bit appliqué sur L'élément 1B et chacun des deux bits présents dans L'élément 60. Le même processus se répète à chaque séquence de calcuL, avec stockage des résultats. La figure 9e montre Le résultat final obtenu à l'issue de la première procédure de calcul. Les n = 5 mémoires éLémentaires de chacun des registres 65 et 66 contiennent Les digits des produits a11x1 et a21 x1 codés en binaire mixte. La conversion en base 2 est faite par des convertisseurs A/D, ou les résultats restent stockés dans Les mémoires en binaire mixte pour La procédure de calcul suivante. Celle-ci donne les produits a12x2 et a12x2 codés en binaire mixte. Le contenu des mémoires devient : mémoires 65 : a11 x1 + a12x2 Cen binaire mixte) mémoires 66 : a21 x1 + a22x2 (en binaire mixte) Lorsque Les procédures de calcul sont terminées, les mémoires élémentaires constituant chacun des registres 65 et 66 contiennent respectivement Les digits de La composante y1 = 2231 et ceux de y2 = 1311 du produit matrice x vecteur de la figure 7. La conversion analogique/ digitale de ces digits donne le résultat final en base 2. La différence entre Le procédé de calcul de l'invention, illustré par La figure 8 et Les procédés de calcul optonumérique connus de L'homme de L'art (par exemple décrits par D.Psaltis et al, "Accurcte numerical computation by opticaZ convoZutzon " dans International Optical Computing Conference, W.T.Rhodes Ed., Proc. SPIE, vol.232, pp.151-156) réside dans le fait que Le parallélisme de calcul est réalisé nar codage des données sous forme d'une modulation de retards optiques introduits sur un faisceau lumineux unique se propageant dans une succession de moyens de modulation mis en série, et non par codage des données sous forme d'une modulation d'intensité de plusieurs faisceaux lumineux. Ceci a des conséquences importantes sur Les moyens utilisables pour mettre en oeuvre le procédé de caleut optonumérique selon l'inven- tion. On a déjà eu L'occasion de souligner qu'un faisceau lumineux unique augmente La fiabilité et simplifie énormément Les problèmes de calibrage et de dérive de la puissance Lumineuse. Par ailleurs, Le procédé de calcul est compatible avec les autres procédés optonumériques systoliques cités ci-dessus, et conduit donc à un accroissement du parallèlisme de calcul. La figure 10a décrit un mode de réalisation du processeur de La figure 8 qui fonctionne avec les séquences de calcul de La figure 9. La source 67 est par exemple une diode super luminescente de faible lon- gueur de cohérence L = 80 microns. Le mode de réalisation des éléments 1A, 1B, 2A, 2B a déjà été donné à La figure 5. La lentille 68 élargit le faisceau lumineux transmis par Les éléments 1A et 1B de manière à éclairer L'élément 69 de préférence en faisceau parallèle. Celui-ci est un déflecteur acousto-optique fonctionnant en régime de Bragg. L'élément 70 est un masque où chaque pixel délimite une région éclairée du déflecteur acousto-optique. L'élément 61 de la figure 8, est, dans la figure 10, constitué par la lentille 71 et La arille nultiplicatrice 72, de sorte que L'image du masque 70 est formée sur les photodétecteurs 73 et 74, à travers les décodeurs 2A et 2B. Le nombre de cellutes composant les barrettes de photodétecteurs 73 et 74 est égal au nombre d'ouvertures dans le masque ; les cellules de détection sont reliées à des circuits de mise en forme 75,76 et à des circuits de maintien et d'addition 77,78 connus de l'homme de l'art. La grille multiplicatrice 72 est constituee par un réseau d'amplitude ou de phase unidimensionnel ou bidimensionnel dont un mode de réalisation est donné par Viénot et al dans Applied Optics, vol. 12, 5, pp.950-960, 1973 ou par N.Aebischer et al dans La Nouvelle Rev. d'Optique, vol.6, 1, pp.37-47, 1975. Une grille multiplicatrice réalisée par ce procédé permet d'obtenir jusqu'à 30 x 30 images identiques. La figure 10 b décrit une variante de L'élément 61 de la figure 8. Suivant cette variante, un réseau de coupleurs à fibres optiques 90 est utilisé pour former une pluralité d'images de La valve optique 69. La fibre d'entrée 91 de chaque coupleur est couplée à un pixel particulier de La valve optique 69 tandis que les fibres de Sortie 92 sont assignées à des pixels de détection particuliers des barrettes 73 et 74, selon La configuration illustrée à La figure 10 b. Dans ce mode de réalisation, Les coupleurs 90 peuvent etre des coupleurs 100/140 microns avec une fibre d'entrée unique et P fibres de sortie, P étant Le nombre de Lignes de La matrice 55. La figure 11 est une variante du calcul numérique du produit matrice x vecteur de la figure 7 mettant en oeuvre Le processeur déjà décrit à la figure 3 mais avec des séquences de calcul différentes. Le nombre n de modulateurs 1A, 1B, 1C... et de démodulateur 2A, 2B, 2C, est ici déterminé par Le nombre n de bits. Les séquences de calcul conduisant au calcul du produit éLémentaire a11x1 sont indiquées. Les bits formant x1 sont appliqués en mode série sur la source 79 qui est modulée en intensité. Les bits formant a11 sont appliqués en mode paralLèLe sur chaque modulateur. A chaque séquence de calcul, un registre à décalage 80 permet de décaler un bit de a11 vers Le précédent. Les séquences de calcul se répètent ainsi n fois (n = 5 dans l'exemple donné à la figure 11) pour obtenir Le résultat a11x1. La procédure de calcul se répète pour obtenir a12x2, etc... La figure 12 donne une variante d'un processeur mettant en oeuvre une architecture parallèle. Le processeur est formé d'une pluralité de sources lumineuses 81, 82, 83, d'une pluralité de modulateurs 1A, 1B, lC et d'un démodulateur 84. Les retards optiques D(A), D(B), D(C) des modulateurs sont supérieurs à la longueur de cohérence des sources. Une séquence possible de retards optiques est donnée par La progression pL, 2pL, 3pL, 4pL,... Le facteur p(p > l) détermine Le taux de diaphonie. Par exemple une isolation optique de L'ordre de -2ûdB (équivalente à - 40 dB électrique) entre Les multiplieurs est obtenue en prenant p = 2. Ce type d'architecture caractérisée en ce que Les multiplieurs 1A, 1B, 1C sont mis en parallèle peut être utilisée dans l'ensemble des exemples de réalisation présentés auparavant. L'inconvénient par rapport aux systèmes précédents est de nécessiter plusieurs sources Lumineuses 81, 82, 83. Toutefois Le fait que Leur Longueur d'onde puisse être quelconque,reduit en partie cet inconvénient. L'avantage est d'autoriser un nombre plus grand de multiplieurs et donc des dimensions de matrice plus élevées. La figure 12 utilise ce procédé dans le cas du calcul du produit matrice x vecteur de la figure 7. Les éléments 60, 61, 2A, 2B, 2C, 62, 63, 65, 66 sont ceux déjà décrits à propos de La figure 8 ; L'élément 84 est un coupleur à fibre. Les séquences de calcul indiquées à La figure 12 sont identiques à celles déjà discutées à La figure 9. D'autres versions plus élaborées de processeurs peuvent etre facilement imaginées pour traiter des matrices de grande dimension ou des tenseurs d'ordre n en associant à La dite architecture parallèle des sources multiplexaes en longueur d'onde. L'ensemble des architectures précédentes peut être associé au multiplexage en longueur d'onde pour augmenter les cadences de calcul, la précision de calcul et La dimension des matrices à traiter. Les processe@rs peuvent être également bouclés sur eux-mêmes pour Le calcul d'algorithmes particuliers connus de l'homme del'art (inversion de matrices, calcul de valeurs propres,etc...) Il est bien évident que Les descriptions qui précédent n'ont été données qu'à titre d'exemple non limitatifs et que de nombreuses variantes peuvent être envisagées sans pour autant sortir du cadre de l'invention. REVENDICATIONS 1 - Multiplieur optoélectronique rour multiplier deux signaux et caractérisé Par le fait qu'il comprend: - une source lumineuse unique (1) emettant un faisceau lumineux (2) dont l'énergie est modulée par un signal (a,), - un moyen de modulation (1A) qui introduit une ioi de filtrage prédétermlnee maximums multiples sur les composantes spectrales du faisceau lumineux sortant (3) et dans lequel la lumière subi= une modulation spectrale accordable sous l'effet dun autre signal (b@) sans qu 'intervienne une modulation de l'énergie du faisceau lumineux qui s'y propage, - un moyen de modulation (2A) qui introduit une loi de filtrage spectral prédéterminée sensiblement identique à la loi de filtrage spectral du moyen de modulation (1A) et dans lequel l'énergie lumineuse du faisceau (3) qui s'y propage est rendue proportionnelle au produit des deux signaux, - un moyen (3A) assurant la détection et la mise e forme du résultat de la multiplication entre les deux signaux et composé d'un photodétecteur (4) et d'un circuit de mise en forme (5). 2 - Multiplieur selon la revendication 1 et caractérisé par le fait que: - le moyen de modulation est un modulateur introduisant un retard optique prédéterminé supérieur à la longueur de cohérence de la source et modulable autour de cette valeur prédéterminée sous l'effet d'un des deux signaux, - le moyen de démodulation est un démodulateur introduisant un retard optique sensiblement équivalent à la valeur du retard optique prédéterminé caractérisant le modulateur (1A), - le fonctionnement du multiplieur est indépendant du fait que le modulateur (1A) est situé avant ou après le démodulateur (2A). - le fonctionnement du multiplieur est indépendant du fait que la modulation d'énergie se fait au niveau de la source ou en tout autre endroit du circuit lumineux par un modulateur externe. 3 - Multiplieur seion les revendications 1 et 2 et caractérisé par le fait que le dit moyen ce modulation est un modulateur électrooptique introduisent un retard optique prédéterminé 4 - Multiplieur selon les reventdications 2 et 3 dns lequel les retards optiques sont obtenus par des moyens faisant intervenir la biréfringence naturelle des cris- taux électrooptiques ou non. 5 - Multiplieur selon la revendication 1 et caractérisé par le fait que les moyens de modulation et de démodula- tion sont constitues Par des coup leurs directifs qui induisent une loi de filtrage spectral 8 maximums multiples équivalente d celle qu'introduirait un retard optique supérieur à la longueur de cohérence de la source. 6 - Processeur optoelectronique pour le calcul de produits matrice vecteur entre un vecteur d'entrée et une matrice, comprenant une pluralité de multiplieurs sels les revendications 1 é 6 et caractérisé par le fait qu'il comprend - une source lumineuse unique (12) servant de source commune émettant un faisceau lumineux (3) dont l'énergie est modulée par les éléments du vecteur d'entrée, - une pluralité de moyens de modulation (1A. 1B,... 1P) mis en série et dont chacun assigne au faisceau lumineux (3) modulé en intensité une loi de filtrage spectral a maximums multiples modulables sous l'effet de signaux représentant les éléments de la matrice sans que ceci n'introduise une! modification de la modulation d'intensité lumineuse du faisceau (3), - une pluralité de moyens de démodulation (2A, 2B,.. 2P) mis en parallèle et dont chacun assigne au faisceau lumineux qui le traverse une loi de filtrage prédéterminée sensiblement équivalente å la loi de filtrage des moyens de modulation et ainsi modifie l'énergie de sortie du faisceau lumineux pour la rendre proportionnelle au produit élémentaire entre un élément du vecteur d'entrée et les éléments de la matrice appliqués sur les moyens de modulation, - une pluralité de moyens de detection et de mise en forme (3A, 3B.. .3P) des signaux délivrés par les moyens de modulation qui fournissent les produits élémentaires, - une pluralité de moyens de maintien et de stockage (4A. 4B. . .4P) des produits élémentaires pour en effectuer la somme et obtenir ainsi les composantes du vecteur de sortie qui constituent le resultat du produit matrice vecteur. 7 - Processeur selon la revendication 6 et caractérisé par Le fait que * les dits moyens de modulation sont constitués par des modulateurs introduisant une série de retards optiques prédéterminés, supérieurs à la longueur de cohérence de La source et modulables autour de ces valeur; prédéterminées sous L'effet de signaux représentant les éléments de La matrice ; * les dits moyens de démodulation sont constitués par des démodu Lateurs introduisant une série de retards optiques sensiblement équiva Lents aux retards optiques caractérisant Les dits moyens de modulation ; * le fonctionnement du processeur est indépendant du fait que Les dits moyens de modulation sont situés avant ou après les dits moyens de démodulation ; 8 - Processeur selon la revendication 7 et caractérisé par le fait que les retards optiques assignés aux dits moyens de modulation suivent sensiblement la progression conne pt, 3*, 8c3, p ét2nt un facteur enter ou non, supérieur à 1, déterminant l'isolation optique entre chaque étage de calcul fonctionnant en paraLLèle 9 - Processeur selon La revendication 6 et caractérisé par Le fait que les dits moyens de modulation et de démodulation sont des coupleurs directifs introduisant une Loi de filtrage spectral à maximums multiples équivalente à La loi de filtrage qu'introduirait un retard optique supérieur à La longueur de cohérence ; 10 - Processeur optoélectronique selon L'une quelconque des revendications 6,7,8,9 et caractérisé en ce qu'il comprend un ou une pluralité de registres à décalage -(80) permettant de transférer les digits des nombres à multiplier d'un modulateur à un autre, et une pluralité de convertisseurs analogiques/digitaux permettant d'obtenir le résultat du produit matrice par vecteur sous forme numérique. 11 - Processeur pour Le calcul numérique de produits matrice x vecteur entre un vecteur d'entrée et une matrice dont les éléments sont sous forme digitale, et caractérisé par Le fait qu'il comprend : * une source lumineuse (59) unique fonctionnant en régime continu * une pluralité de moyens de modulation (1A, 1B,...) mis en série dont chacun assigne à la lumière une Loi de filtrage spectral prédéterminée à maximums multiples modulable sous l'effet de signaux représentant Les bits des éléments de La matrice sans qu'intervienne une modulation de L'intensité de La Lumière ; * une valve optique (60) permettant d'afficher sous forme optique simultanément plusieurs bits d'un élément vecteur d'entrée ; * une pluralité de moyens de démodulation (2A, 28,...) dont chacun assigne à La lumière une loi de filtrage prédéterminée sensiblement équivalente à La loi de filtrage de chacun des dits moyens de modulatin et ainsi modifie L'énergie de sortie de la lumière pour La rendre proportionnelle au produit entre Les bits du vecteur d'entrée et Les bits des éléments de la matrice ; * une pluralité de moyens (62,63) de détection, de mise en forme et de stockage des signaux optiques dé Livrés par chacun des dits moyens de démodulation ; * un moyen (61) permettant de former l'image de la dite valve optique (60) sur Les dits moyens de détection (62,63) 12 - Processeur selon La revendication 11 caractérisé par le fait que * les dits moyens de modulation sont des modulateurs introduisant une série de retards optiques prédéterminés, supérieurs à La Longueur de cohérence de La source, et modulables autour de ces valeurs prédéterminées sous L'effet de signaux représentant les bits des éléments de La matrice ; * les dits moyens de démodulation sont constitués par des démodu Lateurs introduisant une série de retards optiques sensiblement équiva Lents aux retards optiques entrant dans Les dits modulateurs ; *.le fonctionnement du processeur est indépendant du fait que Les dits moyens de modulation sont si tués avant ou après les dits moyens de démodulation ; 13 - Processeur selon la revendication 12 et caractérisé par le fait que le ditsretards optiques suivent sensiblement la progression pL, 3pL, 8pL,..., p etant un facteur entier ou non, supérieur à 1, déterminant l'isolation optique entre Les étages de calcul fonctionnant en parallèle ; 14 - Processeur selon La revendication 11 et caractérisé par Le fait que les dits moyens de modulation et de démodulation sont des coupleurs directifs introduisant une loi de filtrage spectraL à maximums multiples équivalente à La loi de filtrage qu'introduirait un retard optique suoé- rieur à la longueur de cohérence de la source ; 15 - Processeur selon l'une quelconque des revendications 11, 12, 13, 14, caractérisé par le fait cue La dite valve optiaue (60) est un modulateur acousto-optique ; 16 - Processeur selon l'une quelconque des revendications 11, 12, 13, 14, caractérisé par le fait que la dite valve optique (60) est un modulateur magnéto-optique 17 - Processeur selon l'une quelconque des revendications 11, 12, 13 14 caractérisé par le fait que La dite valve optique (60) est un afficheur électro-optique ; 18 - Processeur selon l'une quelconque des revendications 11, 12, 13, 14, caractérisé par le fait que le moyen (61) est un système d'imagerie faisant intervenir des Lentilles sphérocylindriques ou une matrice de coupleurs à fibres optiques ou encore une grille multiplicatrice. 19 - Processeur optoélectronique pour Le calcul de produit matrice x vecteur entre un vecteur d'entrée et une matrice, et caractérisé par le fait qu'il comprend * une pluralité de sources lumineuses (81,82,83) fonctionnant en parallèle t une pluralité de moyens de modulation (1A, 1B, 1C) dont chacun assigne au faisceau lumineux émis par chaque source une Loi de filtrage s@ectral ceterminee a maximums multiples qui est modulable sous l'effet ce sianau@ 20 - Processeur selon la revendication 19 caractérisé par le fait que * Les dits moyens de modulation sont des modulateurs introduit sant des retards optioues prédéterminés supérieurs å la longueur de coherence de La source, qui sont modulables autour de ces va Leurs prédéterminées sous L'effet de signaux représentant Les bits des éléments de la matrice, * les dits moyens de démodulation sont des démodulateurs intro@ulsan un retard optique sensiblement équivalent é L'un tes retards optiques caractérisant Les dits modulateJrs. 21 - Processeur selon La revendication 20 et caractérisé par le fait que les retards optiques des dits moduLateurs suivent sensiblement la progression pL, 2pL, 3pL,...D étant un nombre supérieur à 1 déterminant l'isolation optique entre chaque étage de calcul. 22 - Processeur associant des étages de calcul selon L'une des revendications 1 à 18 a des étages de calcul selon l'une des revendications 19, 20, 21. 23 - Processeur pour le calcul de produits matrice x vecteur ou pour des calculs tensoriels, caractérisé par Le fait qu'il comprend une pluralité d'étages de calcul selon L'une quelconque des revendications 1 à 22 multiplexés en Longueur d'onde.