Computational routine
eng

chuafonc

File content

/* Modnumlib Scicos interfacing function
 * Copyright (C) 2009-2011 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.
 */
 
/* chuafonc Scicos Chua non-linear function block
 * Type 2 simulation function ver 1.0 - scilab-2.6&2.7
 * 8 novembre 2003 - IRCOM GROUP - Author : A.Layec
 */

/* REVISION HISTORY :
 * $Log$
 */

#define Abs(x) ( ( (x) >= 0) ? (x) : -( x) )
#include "machine.h"

/* Cette fonction de simulation réalise la fonction non linéaire de la diode de Chua :
 *
 * y=m1*u+1/2*(m0-m1)*(abs(u+E)-abs(u-E))
 *
 * où m0 et m1 sont des vecteurs de paramètres réels,
 * E=1 et u est le vecteurs des entrées, y le vecteur des sorties
 *
 * Entrée régulière : u[0..nu-1] : vecteur des entrées
 * Sortie régulière : y[0..nu-1] : vecteur des sorties
 * Entrée évènementielle : néant (héritage)
 * sortie évènementielle : néant
 * paramètres : rpar[0..nu-1] vecteur m0
 *              rpar[nu..2*nu-1] vecteur m1
 */

/*prototype*/
void
chuafonc(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 nu,i;

    /*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 u*/
    nu = insz[0];

    /*Calcule du registre de sortie*/
    for(i=0;i<nu;i++)
        y[i]=rpar[nu+i]*u[i]+0.5*(rpar[i]-rpar[nu+i])*(Abs(u[i]+1)-Abs(u[i]-1));
}