Computational routine
eng

nintegr

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.
 */
 
/* nintegr Scicos vector intgerator
 * Type 2 simulation function ver 1.0 - scilab-2.6&2.7
 * 16 novembre 2003 - IRCOM GROUP - Author : A.Layec
 */

/* REVISION HISTORY :
 * $Log$
 */

#include "machine.h"

/* Cette fonction de simulation réalise l'intégration simple
 * du vecteur u[] par le solveur ODE intégré et place le résultat dans y[]
 *
 * 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 : néant (temps dépendant)
 * sortie évènementielle : néant.
 * état continu : u[0..in_size-1] vecteur des états continus
 */

/*Prototype*/
void nintegr(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*/
 double *y;
 double *u;
 int i,nu;

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

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

 /*Le flag 1 place le vecteur des états continus x[]*/
 /*dans le registre de sortie y[]*/
 if(*flag == 1)
    for(i=0;i<nu;i++) y[i]=x[i];

 /*Le flag 2 place les valeurs d'entrée u[]*/
 /*dans le vecteur des dérivées xd[]*/
 else if(*flag ==0)
    for(i=0;i<nu;i++) xd[i]=u[i];
}