/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        I N C R E M E N T A T I O N / D E C R E M E N T A T I O N  :                                                               */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/increment.01$c' :                                                                                          */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#define   N         11                                                                                                                  \
                                        /* Definition du nombre d'incrementations/decrementations.                                   */
#define   Xd        -1.0
#define   Xa        1.0
                                        /* Definition des extrema de 'x'.                                                            */
#define   PAS       0.2
                                        /* Pas des incrementations/decrementations.                                                  */
main()
{
int       n;
                                        /* Numero de l'iteration.                                                                    */
double    xm=Xd,xd=Xa;
                                        /* 'x' "montant" et "descendant"...                                                          */
double    dx=PAS;
                                        /* Pas des incrementations/decrementations.                                                  */
for       (n=1 ; n<=N ; n++)
          {
          printf("xm=%.16g   xd=%.16g\n",xm,xd);
                                        /* Les suites "montante" et "descendante" devraient etre symetriques, or elles ne le         */
                                        /* sont pas pour certaines valeurs (et suivant les SYSTEMEs...) :                            */
                                        /*                                                                                           */
                                        /* 1-sur 'SYSTEME_SG...', 'SYSTEME_HP...' :                                                  */
                                        /*                                                                                           */
                                        /*                  xm=-1                                   xd=1                             */
                                        /*                  xm=-0.8                                 xd=0.8                           */
                                        /*                  xm=-0.6000000000000001                  xd=0.6000000000000001            */
                                        /*                  xm=-0.4000000000000001   <----          xd=0.4000000000000001            */
                                        /*                  xm=-0.2000000000000001        |         xd=0.2000000000000001            */
                                        /*                  xm=-5.551115123125783e-17     |         xd=5.551115123125783e-17         */
                                        /*                  xm=0.2                        |         xd=-0.2                          */
                                        /*                  xm=0.4                         ---->    xd=-0.4                          */
                                        /*                  xm=0.6                                  xd=-0.6                          */
                                        /*                  xm=0.8                                  xd=-0.8                          */
                                        /*                  xm=1                                    xd=-1                            */
                                        /*                                                                                           */
                                        /*                                                                                           */
                                        /* 2-sur 'SYSTEME_VAX9000_ULTRIX_GC' :                                                       */
                                        /*                                                                                           */
                                        /*                  xm=-1                                   xd=1                             */
                                        /*                  xm=-0.8                                 xd=0.8                           */
                                        /*                  xm=-0.6                                 xd=0.6                           */
                                        /*                  xm=-0.4                                 xd=0.4                           */
                                        /*                  xm=-0.2                                 xd=0.2                           */
                                        /*                  xm=-6.938893903907228e-18               xd=6.938893903907228e-18         */
                                        /*                  xm=0.2                                  xd=-0.2                          */
                                        /*                  xm=0.4                                  xd=-0.4                          */
                                        /*                  xm=0.6                                  xd=-0.6                          */
                                        /*                  xm=0.8                                  xd=-0.8                          */
                                        /*                  xm=1                                    xd=-1                            */
                                        /*                                                                                           */
          xm = xm + dx;
          xd = xd - dx;
          }
}



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.