/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        R E C H E R C H E   D U   P L U S   G R A N D   N O M B R E   N O N   S I G N I F I C A T I F                              */
/*        P A R   R A P P O R T   A   1  :                                                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/epsilon.01$c' :                                                                                            */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#define   ITER                                                                                                                          \
                    20
#define   UN                                                                                                                            \
                    1.0

main()
     {
     int       n;
     float     f=UN;
     double    d=UN;

     for       (n=1 ; n<=ITER ; n++)
               {
               float     xf=UN-f;
               double    xd=UN-d;
               if        (xf == UN)
                         {
                         printf("\n float=%g",f);
                                        /* On n'edite que les valeurs de 'f' qui ne sont pas significatives devant UN.               */
                         }
               else
                         {
                         }

               if        (xd == UN)
                         {
                         printf("\n double=%g",d);
                                        /* On n'edite que les valeurs de 'd' qui ne sont pas significatives devant UN.               */
                         }
               else
                         {
                         }
               f = f/10.0;
               d = d/10.0;
               }
     }



Copyright © Jean-François Colonna, 2021-2023.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2021-2023.