Computational routine
eng


fmodulo

File content


/* fmodulo Scicos float modulo function
 * Type 2 simulation function ver 1.0 - scilab-2.6&2.7
 * 19 octobre 2003 - IRCOM GROUP - Author : A.Layec
 */

/* REVISION HISTORY :
 * $Log$
 */

#include<math.h>

#include "machine.h"

/* Cette fonction de simulation réalise la fonction non-lineaire
 * modulo sur des nombres définis en virgule flottante.
 * y=u-0.5*m-m*floor( (u-0.5m) / m)-0.5*m
 * entrées régulières : vecteur des entrées u[0..nu-1]
 *
 * sorties régulières : vecteur des sorties y[nu..2nu-1]
 * paramètres : rpar[0..nu-1] : m : vecteur des amplitudes
 * entrées d'évenement : néant
 * sorties d'évenement : néant
 *
 * rmq : la sortie est soit positive soit négative
 */

/*prototype*/
void fmodulo(flag,nevprt,t,xd,x,nx,z,nz,tvec,ntvec,rpar,nrpar,
             ipar,nipar,inptr,insz,nin,outptr,outsz,nout)
integer *flag,*nevprt,*nx,*nz,*ntvec,*nrpar,ipar[],*nipar,insz[],*nin,outsz[],*nout;
double x[],xd[],z[],tvec[],rpar[];
double *inptr[],*outptr[],*t;
{
 /*déclaration des variables*/
 int i,nu;
 double *y;
 double *u;

 /*Récupération des adresses des ports réguliers*/
 y=(double *)outptr[0];
 u=(double *)inptr[0];

 /*Récupération de la taille du port d'entrée*/
 nu=insz[0];

 /*Calcul de la valeur modulo et place dans le registre de sortie*/
 for(i=0;i<nu;i++)
      y[i] = (u[i]-0.5*rpar[i])-rpar[i]*floor((u[i]-0.5*rpar[i])/rpar[i])-0.5*rpar[i];
}