Add here a paragraph of the function description.
/* int_euler Scicos discrete integral by euler method block * Type 2 simulation function ver 1.0 - scilab-2.6&2.7 * 17 novembre 2003 - IRCOM GROUP - Author : A.Layec */ /* REVISION HISTORY : * $Log$ */ #include "machine.h" /* Cette fonction de simulation réalise une intégration discrète * du vecteur d'entrée par la méthode d'Euler arrière (dérivée à gauche) : * y[k]=h*u[k]+y[k-1] * où h est le pas d'intégration, y[k] le vecteur de sortie, * et u[k] le vecteur d'entrée * * entrée régulières : u[0..nu-1] vecteur d'entrée * sorties régulières : y[0..nu-1] registre de sortie * entrée évènementielle : Instants de déclenchement * sortie évènementielle : néant. * état discret : z[0..nu-1] vecteur des états de sortie discrets précédents * * paramètre réel : rpar[0] pas d'intégration */ /*prototype*/ void int_euler(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*/ double *y; double *u; int i,nu; /*récupération de l'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]; /*Le flag 1 calcule l'integrale et place le résultat dans le registre y[]*/ if(*flag==1) { for(i=0;i<nu;i++) y[i]=rpar[0]*u[i]+z[i]; } /*le flag 2 place l'état de sortie y dans z*/ else if(*flag==2) { for(i=0;i<nu;i++) z[i]=y[i]; } }