/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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 E   P O I N T S   P S E U D O - A L E A T O I R E S   A V E C   P E R I O D I C I T E  :                                 */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xrC/ObjetComplexe.B1$vv$c' :                                                                                   */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 20130507135241).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#ifdef    COMMENTS
#         undef     COMMENTS
#else
#endif

#define   COMMENTS  "Champ pseudo-aleatoire type 'Fsinus(...)' (format octet) -noir ou blanc-"

#include  "images_1octet.01.vv.I"

#include  "Fmodulo.01.vv.I"
#include  "Fsinus.01.vv.I"

#ifndef   INCREMENT
#         define    INCREMENT                                                                                                           \
                              (1.033351371227)
#else
#endif

double    increment=INCREMENT;

MAIN(
          {
          double    Xi=1;
          double    i=1.277231953;

          for       (y=Ymin ; y <= Ymax ; y=ADD2(y,1))
                    {
                    for       (x=Xmin ; x <= Xmax ; x=ADD2(x,1))
                              {
                              Xi=ADD2f(Xi,i);
                              i=ADD2f(i,increment);
                              STORE_IMAGE_TORE(x,y,COND(IFLT(FonctionPrivee_Fsinus(Xi),0),NOIR,BLANC));
                              }
                    }

          SORTIE_DU_FICHIER_IMAGE;
          EDITER_LES_COMPTEURS_DES_FONCTIONS;
          }
     )