Computational routine
eng

convsp

File content

/* Modnumlib Scicos interfacing function
 * Copyright (C) 2009 Alan Layec
 *
 * This library 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.
 *
 * This library 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 this library; if not, write to the
 * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
 * Boston, MA 02111-1307, USA.
 */
 
/* convsp Scicos Convertisseur série parralèle
 * Type 2 simulation function ver 1.0 - scilab-2.6&2.7&3.0
 * 21 décembre 2003 - IRCOM GROUP - Author : A.Layec
 */

/* REVISION HISTORY :
 * $Log$
 */

#include "modnum_lib.h"

/* Cette fonction réalise la conversion "série-parallèle". Le valeur du vecteur d'entrée
 * U[] est stockée dans Z[] et est délivrée lorsque le nombre d'échantillons stocké est N,
 * avec N la taille du port de sortie.
 *
 * entrées régulières : u[0] : valeur d'entrée
 * sorties régulières : y[0..N-1] : vecteur de sortie
 * états discrets : z[0..N-1] : valeur u[k-N,..,k-2,k-1,k]
 *                 z[N] : compteur échantillons
 */

/*prototype*/
void convsp(flag,nevprt,t,xd,x,nx,z,nz,tvec,ntvec,rpar,nrpar,
            ipar,nipar,inptr,insz,nin,outptr,outsz,nout)

int *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 ny,k;

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

  /*récupération de la taille du port de sortie*/
  ny=outsz[0];

  /*Uniquement sur le flag 1*/
  if(*flag==1)
  {
   /*Test valeur compteur échantillons*/
   if((int)z[ny]==(outsz[0]-1))
   {
    /*Recopie z dans y*/
    copyd_c(&ny,&z[0],&y[0]);
    y[ny-1]=u[0];
   }
  }
  else if(*flag==2)
  {
   /*Test valeur compteur échantillons*/
   if((int)z[ny]==outsz[0])
   {
    /*Ajout du nouveau point*/
    z[0]=u[0];
    /*Incremente compteur échantillon*/
    z[ny]=1.;
   }
   else
   {
    /*ajout du nouveau point*/
    z[(int)z[ny]]=u[0];
    /*Incremente compteur échantillon*/
    z[ny]++;
   }
  }
}