Comment utiliser et vérifier la précision du quad en gfortran?

Je suis en train d'utiliser quad précision dans gfortran, mais il semble que le real*16 ne fonctionne pas. Après un peu de pêche autour, j'ai trouvé qu'il peut être mis en œuvre comme real*10. Est real*10 fait quad précision?

Comment puis-je tester la précision de mon code? Est-il un standard de l'algorithme simple pour les tests de précision? Par exemple, quand je veux comprendre ce que l'ordinateur zéro est, je continue à diviser par 2.0 jusqu'à ce que je atteindre 0.0. Garder la trace des valeurs qui me permet de savoir lorsque l'ordinateur pense que mon nombre différent de zéro est zéro en me donnant des nuls en informatique.

Est-il un bon moyen de déterminer la précision avec un type d'algorithme comme je l'ai décrit?

source d'informationauteur drjrm3