Low level routine
fr - eng

genint_c - random integer number generator computational routine

Module

• Modnum

Parameters

• n : size of vectors
• m : address of length vector of the binary words in number of bits (scalar)
• typ : address of type of binary generators (scalar)
  • 0 : Generate random Non Return Zero (NRZ) single bits.
  • 1 : Generate random Return Zero (RZ) single bits or unsigned integers.
  • 2 : Generate random signed integers.
• y : address of output vector

File content

/* genint_c subroutine
 * Random Integer Number Generator
 *
 * Copyright (C) 2007-2009 Alan Layec
 *
 * This file is part of modnumlib.
 *
 * modnumlib is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * modnumlib is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with modnumlib; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 *
 */

/* REVISION HISTORY :
 * $Log$
 */

#include <stdlib.h>

#include "modnum_lib.h"

/* genint_c routine de calcul de mots entiers aléatoires
 *
 * Entrées :
 *  n     : dimension 1 de la matrice de sortie (scalaire)
 *  m     : dimension 2 de la matrice de sortie (scalaire)
 *  nbit  : longueur des mots binaires des éléments du vecteur de sorties (scalaire)
 *  typ   : type des générateurs binaires (scalaire)
 *          0 : génére 1 seul bit codé NRZ
 *          1 : génère 1 seul bit codé RZ, ou un mot non signé
 *          2 : génère un mot codé code complément à 2
 *
 * Sorties :
 *  y : pointeur sur double de sortie
 *
 * dépendances
 * math.h
 */

void genint_c(int *n,int *m,int *nbit,int *typ,double *y)
{
 /*déclaration des variables compteurs*/
 int i,j,k,l;

 switch (*typ)
 {
   /*Type 0 -> dédié signal NRZ*/
   case 0 :
   {
    for(l=0;l<(*m);l++) {
      for(i=0;i<(*n);i++) {
        y[l*(*n)+i] = (double) ((rand()&1)*2-1);
      }
    }
    break;
   }

   /*Type 1 -> mot non signé*/
   case 1 :
   {
    k = 1;
    k = k<<(*nbit);
    for(l=0;l<(*m);l++) {
      for(i=0;i<(*n);i++) {
        j = rand();
        y[l*(*n)+i] = (double) (j&(k-1));
      }
    }
    break;
   }

   /*Type 2 -> mot codé code complémént à 2*/
   case 2 :
   {
    for(l=0;l<(*m);l++) {
      for(i=0;i<(*n);i++) {
        j  = rand();
        j -= 2<<((*nbit)-2);
        j &= (2<<((*nbit)-1)) - 1;
        j -= 2<<((*nbit)-2);
        y[l*(*n)+i] = (double) j;
      }
    }
   }
   break;
 }
 return;
}

void geninti_c(int *n,int *m,int *nbit,int *typ,int *y)
{
 /*déclaration des variables compteurs*/
 int i,j,k,l;

 switch (*typ)
 {
   /*Type 0 -> dédié signal NRZ*/
   case 0 :
   {
    for(l=0;l<(*m);l++) {
      for(i=0;i<(*n);i++) {
        y[l*(*n)+i] = ((rand()&1)*2-1);
      }
    }
    break;
   }

   /*Type 1 -> mot non signé*/
   case 1 :
   {
    k = 1;
    k = k<<(*nbit);
    for(l=0;l<(*m);l++) {
      for(i=0;i<(*n);i++) {
        j = rand();
        y[l*(*n)+i] = (j&(k-1));
      }
    }
    break;
   }

   /*Type 2 -> mot codé code complémént à 2*/
   case 2 :
   {
    for(l=0;l<(*m);l++) {
      for(i=0;i<(*n);i++) {
        j  = rand();
        j -= 2<<((*nbit)-2);
        j &= (2<<((*nbit)-1)) - 1;
        j -= 2<<((*nbit)-2);
        y[l*(*n)+i] = j;
      }
    }
   }
   break;
 }
 return;
}

/* genintv_c routine de calcul de mots entiers aléatoires
 *
 * Entrées :
 * n   : taille du  vecteur de sortie (scalaire)
 * m   : longueur des mots binaires des éléments du vecteur de sorties (vecteur)
 * typ : type des générateurs binaires (vecteur)
 * (typ=0: génére 1 seul bit codé NRZ
 *    =1: génère 1 seul bit codé RZ, ou un mot non signé
 *    =2: génère un mot codé code complément à 2)
 * Sorties :
 * y : vecteur de sortie
 *
 * dépendances
 * math.h
 */
void genintv_c(int *n,int *m,int *typ,double *y)
{
 /*déclaration des variables compteurs*/
 int k,i,j;

 for(i=0;i<(*n);i++)
 {
  switch (typ[i])
  {
   /*Type 0 -> dédié signal NRZ*/
   case 0:
   {
    y[i]=(int)((rand()&1)*2-1);
    break;
   }
   /*Type 1 -> entier non signé*/
   case 1 :
   {
    k=1;
    k=k<<m[i];
    j=rand();
    y[i]=(j&(k-1));
    break;
   }
   /*Type 2 -> entier code complémént à 2*/
   case 2 :
   {
    j=rand();
    j -= 2<<(m[i]-2);
    j &= (2<<(m[i]-1)) - 1;
    j -= 2<<(m[i]-2);
    y[i]=j;
    break;
   }
   break;
  }
 }
 return;
}

Authors

IRCOM Group Alan Layec