Discussion :

[nikopol82]

Bonjour,

La question que je pose porte sur l'éventuelle spécification dans la norme du c++ (et surtout la prochaine norme) de l'implémentation des nombres

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

std::complex<T>

.

En effet, beaucoup de codes numériques issus du fortran pré-supposent que un nombre complexe est implémenté en mémoire comme étant la suite de deux nombres de type T, le premier représentant la partie réelle et l'autre la partie imaginaire.

Une interprétation "objet" de std::complex, qui serait défini par son interface laisserait la possibilité de le représenter d'une manière quelconque. (par exemple module et angle ou encore autre chose).

En revanche le type complex du C99 spécifie lui aussi une représentation en mémoire comme la succession de la partie réelle et imaginaire.

Un intérêt d'une telle représentation est qu'un tableau de double de taille 2*n peut être utilisé comme un container n éléments de type complex à l'aide d'un

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

reinterpret_cast

. Une référence sur un type complexe peut donc facilement être obtenue en changeant les parties imaginaires et réelles "en place". Par exemple, c'est l'hypothèse qui est faite par dans l'ouvrage numerical recipes, dans le chapitre sur la FFT.

Ce problème fait l'objet d'un post de Grosse-Kunstleveet David Abrahams ici:

http://www.open-std.org/jtc1/sc22/wg...002/n1356.html


A la connaissance des auteurs du posts, même si celà n'est pas spécifié dans la norme de 98, tous le compilateurs font une représentation des std::complex autorisant cette confusion.

D'autres liens à ce propos:

http://www.open-std.org/jtc1/sc22/wg...004/n1589.html

http://anubis.dkuug.dk/jtc1/sc22/wg2...ctive.html#387

Si vous savez ce qu'il en est pour la prochaine norme du c++ finalement, ça m'intéresse beaucoup.

Par ailleurs, que pensez-vous d'un code C++ qui ferait déjà cette supposition sur la représentation du std::complex, au sens ou ça fonctionne sous C99 et sur les compilateurs connus ?

S.

[Lavock]

Rien de tout ça, en tout cas dans le n3000. Apparemment, on ne pourra pas construire un conteneur de complex avec un tableau de dimension. Pas plus qu'il n'existe un opérateur de conversion pour faire l'inverse.

[nikopol82]

En fait, ce qui est intéressant surtout c'est la validité du reinterpret_cast entre un std::complex<double> et en tableau double[2], et inversement.
(Pour le problème de la fft en particulier...)

[Alp]

Le std::complex ne "transporte", niveau données, que les deux floats/doubles, rien d'autre. Après c'est que des fonctions. Ce n'est donc pas *si* étonnant que ça.

[JolyLoic]

On peut lire dans le n3000 (26.4/4) :

If z is an lvalue expression of type cv std::complex<T> then:
— the expression reinterpret_cast<cv T(&)[2]>(z) shall be well-formed,
— reinterpret_cast<cv T(&)[2]>(z)[0] shall designate the real part of z, and
— reinterpret_cast<cv T(&)[2]>(z)[1] shall designate the imaginary part of z.
Moreover, if a is an expression of type cv std::complex<T>* and the expression a[i] is well-defined for an
integer expression i, then:
— reinterpret_cast<cv T*>(a)[2*i] shall designate the real part of a[i], and
— reinterpret_cast<cv T*>(a)[2*i + 1] shall designate the imaginary part of a[i].

[nikopol82]

Yes, merci JolyLoic ! C'est exactement ce que je cherchais.

[Lavock]

... et inversement.

Apparemment non pour le "inversement" par contre.