Voir la source de Unum

{{Redirect|Posit|homonymie=Posit Software, PBC|autres=la société informatique}}
{{ébauche|informatique}}

Les « nombres universels » (''{{lang|en|universal numbers}}''), dits '''unum''', sont une famille de [[Quantification (signal)|quantifications]] de la [[droite projective]] proposées depuis 2015 par John L. Gustafson. Elles sont conçues comme une alternative au standard [[IEEE 754]], relatif aux [[nombres flottants]].

Gustafson a révisé plusieurs fois son modèle, de telle sorte qu'il existe plusieurs versions : l'unum de type I, l'unum de type II et l'unum de type III, dit '''posit'''. Le nom ''posit'' provient du verbe anglais ''{{lang|en|to posit}}'' qui lorsque utilisé comme un nom peut signifier « assertion » ou bien « [[conjecture]] ». Gustafson définit aussi deux formats dits '''valid''' et '''quire''' qui complémentent le posit en apportant respectivement des fonctionnalités différentes et une aide pour les calculs intermédiaires.

Toutes ces variantes ont en commun les propriétés suivantes :

* le codage en binaire est à longueur fixe, typiquement choisie comme un nombre entier d'[[octet]]s.
* il existe deux valeurs dites exceptionnelles: 0 et infini (<math>\infty</math>). Elles sont encodées différemment des autres valeurs.
* le décodage est [[injection (mathématiques)|injectif]] : deux éléments distincts de l'ensemble discret correspondent toujours à deux points distincts de la droite projective.

Depuis 2017, les efforts de standardisation du système des nombres universels portent essentiellement sur le posit. Cet effort a abouti en 2022 à la publication d'un document normatif<ref name="standard">{{en}} [https://posithub.org/docs/posit_standard-2.pdf Posit Standard Document].</ref>.

== Unum ==
=== Type I ===
L'unum de type I est assez proche de la norme IEEE 754. La principale différence consiste dans l'utilisation du [[bit de poids faible]] pour implémenter une [[arithmétique d'intervalles]].

=== Type II ===
L'unum de type II est radicalement différent de la norme IEEE 754. En plus du [[bit de signe]] et du bit d'intervalle mentionné précédemment, l'unum de type II utilise un bit pour signifier l'[[inverse|inversion]]. La donnée de ces trois opérations permet, à partir d'un [[ensemble fini]] de points entre l'unité et l'infini, de quantifier toute la droite projective hormis quatre points : les deux exceptions 0 et <math>\infty</math>, puis 1 et -1. Cet ensemble de points est choisi arbitrairement et les calculs arithmétiques les concernant sont effectués non pas de façon logique, mais en faisant appel à une [[table de correspondance]]. La taille d'une telle table devient rédhibitoire pour un format d'encodage de plusieurs octets. Cette difficulté a rendu nécessaire le développement de l'unum de type III, dit ''posit'', discuté ci-après.

== Posit ==
Le posit est plus proche de IEEE 754 que ne l'est l'unum de type II. En substance, la principale différence avec IEEE 754 est que l'exposant est [[Code à longueur variable|encodé en taille variable]]. Il en résulte que la [[précision arithmétique|précision]] est plus élevée auprès de l'unité qu'elle ne l'est auprès de zéro ou de l'infini.

Le posit n'implémente pas d'arithmétique d'intervalles. Gustafson a cependant conçu un format companion, le ''valid'', qui lui permet ce genre de calculs. Gustafson semble considérer ''posit'' et ''valid'' comme deux modes de fonctionnement du format ''unum''.

=== Principe général ===
Le format posit de précision <math>n</math> est un format d'encodage de [[fraction dyadique|fractions dyadiques]] sur <math>n</math> bits. Toute fraction dyadique qui peut être encodée de façon exacte sous un tel format sera ici dite ''représentable''.
Le format nécessite un processus d'[[arrondi (mathématiques)|arrondissage]] pour encoder les réels qui ne sont pas des fractions dyadiques, ou les fractions dyadiques qui ne sont pas représentables.

==== Fractions dyadiques représentables ====
Soit <math>r</math> une fraction dyadique représentable que Gustafson écrit sous la forme :

:<math>r = K 2^M</math>

où <math>K</math> et <math>M</math> sont deux [[entiers relatifs]] qu'il faut encoder ensemble dans <math>n</math> bits, la précision du format. On se réduit ici au cas <math>r > 0</math> et donc <math>K > 0</math>.  L'encodage des [[Nombre négatif|nombres négatifs]] est discuté plus loin.

Comme mentionné plus haut, <math>M</math> doit être encodé en longueur variable : les bits restants, s'il y en a, doivent être utilisés pour encoder <math>K</math>. Comme <math>K</math> est nécessairement non nul puisque zéro est un point exceptionnel, son [[bit de poids fort]] est nécessairement 1 et on peut le considérer implicite. Considérons donc <math>m</math> tel que K s'écrit sur <math>m+1</math> bits. <math>K</math> peut alors s'écrire <math>K = 2^m + F</math>, où <math>0 \leq F < 2^m </math>.  On doit donc encoder non pas <math>K</math> mais <math>F = K-2^m</math>.  <math>F</math> contrairement à <math>K</math>, peut être nul. <math>m</math> doit coïncider avec le nombre de bits restant après l'encodage de <math>M</math>. C'est aussi le nombre de bits nécessaires pour encoder <math>F</math>.

On a:

:<math>r = K 2^M = (2^m + F) 2^M = (1 + \frac{F}{2^m}) 2^{M+m}</math>

Notons <math>f=1 + \frac{F}{2^m}</math> et <math>N = M+m</math>. On a alors: <math>r = f 2^N</math> et, par construction :

:<math>1 \le f < 2 </math>

Il suit:

:<math>\begin{align}
0 & \le & log_2(f) & < & 1\\
N & \le & log_2(f) + N & < & N + 1\\
N & \le & log_2(f 2^N) & < & N + 1\\
N & \le & log_2(r) & < & N + 1
\end{align}</math>

Et donc, par définition de la [[partie entière]]:

:<math>N = \lfloor log_2(r) \rfloor</math>

Il suffit alors d'encoder <math>N</math>. Le nombre de bits restants est <math>m</math>, ce qui nous donne <math>M = N - m</math>, et donc <math>K = r/2^M</math>.

==== Arrondissage des nombres réels et des fractions dyadiques non-représentables ====
{{...}}

=== Codage de l'exposant ===
L'exposant N est encodé avec un [[codage de Rice]] adapté pour encoder des [[entiers relatifs]]. Le paramètre utilisé est <math>\text{es}</math>.  Il peut être nul et est typiquement choisi petit (inférieur à trois). Si les versions initiales du format posit faisaient de <math>\text{es}</math> un paramètre au même titre que <math>n</math>, le standard de 2022 semble avoir finalisé la valeur 2.

Le signe du quotient <math>q</math> est encodé en faisant usage des deux conventions possibles pour un [[codage unaire]], selon que le bit répété est 0 ou 1. Plus précisément, <math>q</math> est encodé avec une séquence commençant par un bit <math>b</math> égal à 0 si <math>q</math> est strictement négatif, et 1 sinon :

:<math>b = \sgn(1+\sgn(q))</math>

Ce bit est répété <math>m</math> fois{{Efn|Il ne faut pas ici confondre <math>m</math> avec celui de la section précédente.}}, avec : 

:<math>m = |q| + b</math>

La séquence de bits encodant le quotient q peut alors s'écrire:

:<math>\underbrace{b,\ldots,b}_{m},\overline{b}</math>

Le reste de la division euclidienne est quant à lui encodé comme dans le codage Rice proprement dit.

Seuls les <math>m</math> bits sont nécessaires. Le codage de Rice est tronqué lorsque <math>m = n - 1</math>.  Les bits absents sont alors implicitement égaux à zéro.

=== Exceptions ===
Zéro est codé par un champ de bits tous égaux à 0. Il s'agit d'une exception car s'il devait être décodé par l'algorithme décrit plus haut, {{passage à vérifier|un tel champ de bits correspondrait à la valeur <math>2^{(n-2)2^\text{es}}</math>}}.

L'infini est codé par un champ de bits dont le premier est 1 et tous les autres sont 0. S'il devait être décodé par l'algorithme décrit plus haut, {{passage à vérifier|un tel champ de bits correspondrait à la valeur <math>2^{-(n-1)2^\text{es}}</math>}}.

=== Complément à deux ===
Afin d'assurer un décodage injectif, les valeurs négatives (marquées comme telles par le [[bit de signe]]) doivent être [[complément à deux|complémentée à deux]] en première étape de décodage et en dernière étape de codage.

== Notes et références ==
=== Notes ===
{{Notes}}
=== Références ===
{{References}}

{{Portail|informatique}}

[[Catégorie:Virgule flottante]]