/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        E T U D E   D E   L A   C O M P L E X I T E   S T R U C T U R E L L E                                                      */
/*        D ' U N E   T R A N S F O R M A T I O N   1 / Z   S U R   U N   T O R E  :                                                 */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrC/ObjetComplexe.91$vv$c' :                                                                                   */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 20130506111247).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#ifdef    COMMENTS
#         undef     COMMENTS
#else
#endif

#define   COMMENTS  "Transformation 1/Z des coordonnees sur un tore (format octet) -noir ou blanc-"

#include  "images_1octet.01.vv.I"

#ifndef   FACTEUR
#         define    FACTEUR                                                                                                             \
                              gMUL2(64000,gDIVI(dimXY,512))
#else
#endif

int       facteur=FACTEUR;

MAIN(
          {
          for       (y=Ymin ; y <= Ymax ; y=ADD2(y,1))
                    {
                    for       (x=Xmin ; x <= Xmax ; x=ADD2(x,1))
                              {
                              int       xc=SOUS(x,dimXs2);
                              int       yc=SOUS(y,dimYs2);
                              if        ((xc!=0) || (yc!=0))
                                        {
                                        STORE_IMAGE_TORE(ADD2(DIVI(MUL2(facteur,(+xc)),ADD2(MUL2(xc,xc),MUL2(yc,yc))),dimXs2)
                                                  ,ADD2(DIVI(MUL2(facteur,(-yc)),ADD2(MUL2(xc,xc),MUL2(yc,yc))),dimYs2)
                                                  ,BLANC
                                                   );
                                        }
                              else
                                        {
                                        }
                              }
                    }

          SORTIE_DU_FICHIER_IMAGE;
          EDITER_LES_COMPTEURS_DES_FONCTIONS;
          }
     )