Arrondir les grands flottants à l'entier

Question

Dans Fortran, je souhaite pouvoir arrondir un grand flottant double précision à l'entier le plus proche. Je tenterai cela comme suit,

PROGRAM rounding 
IMPLICIT NONE 
INTEGER, PARAMETER :: DP = 8 
REAL(KIND=DP) :: value 

value = 12345678987.123456 

print *, CEILING(value) 

END PROGRAM ROUNDING 

Cependant, cela ne sort pas la bonne réponse, mais donne lieu -2147483648. J'apprécie que cela soit lié à la nature et à la précision des points flottants, mais arrondir un nombre comme celui-ci semble un objectif raisonnable.

Quelqu'un peut-il me diriger dans la bonne direction?

Answer 1

La variable entière renvoyée par ceiling() est de type trop petit pour contenir la valeur. Vous devez indiquer au plafond pour produire un entier d'un type plus grand. Sinon, il utilise l'entier par défaut.

Par exemple par

print *, CEILING(value, kind=dp) 

(en supposant qu'il existe ce genre, mais il devrait exister, si DP est sélectionné judicieusement)

Notez que la kind=8 est laid et non portable. Vous ne devez pas utiliser de constantes littérales telles que 8. Il n'y a aucune garantie qu'une telle sorte existe.

J'utiliserais:

INTEGER, PARAMETER :: DP = kind(1.d0) 

Pour l'entier genre vous pouvez utiliser SELECTED_INT_KIND pour sélectionner le type avec des chiffres assez comme

INTEGER, PARAMETER :: IP = SELECTED_INT_KIND(15) 

print *, CEILING(value, kind=IP) 

ou utiliser toute autre méthode de Fortran: integer*4 vs integer(4) vs integer(kind=4)