Estoy tratando de usar la precisión quad en gfortran, pero parece que elreal*16 No funciona. Después de un poco de pesca alrededor, he encontrado que se puede implementar comoreal*10. Esreal*10 actualmente quad precision?

¿Cómo puedo probar la precisión de mi código? ¿Hay un algoritmo simple estándar para probar la precisión? Por ejemplo, cuando quiero averiguar qué es el cero del equipo, continúo dividiendo por 2.0 hasta que llego a 0.0. Hacer un seguimiento de los valores me permite saber cuándo la computadora "piensa" que mi número distinto de cero es cero, lo que me da cero.

¿Hay una buena manera de averiguar la precisión con un tipo de algoritmo como el que describí?