tcsltj4.c

Go to the documentation of this file.

#include <math.h>
#include <memory.h>
#include "machine.h"
#include "scicos_block.h"
#include "sci_mem_alloc.h"  /* malloc */

extern int C2F(dmmul)();
extern int C2F(dmmul1)();

void tcsltj4(scicos_block *block,int flag)
{
  /* Copyright INRIA

     Scicos block simulator
     continuous state space linear system simulator
     rpar(1:nx*nx)=A
     rpar(nx*nx+1:nx*nx+nx*nu)=B
     rpar(nx*nx+nx*nu+1:nx*nx+nx*nu+nx*ny)=C
     rpar(nx*nx+nx*nu+nx*ny+1:nx*nx+nx*nu+nx*ny+ny*nu)=D */

  int un=1,lb,lc;
  int nx=block->nx;
  double* x=block->x;
  double* xd=block->xd;
  double* rpar=block->rpar;
  double* y=block->outptr[0];
  double* u1=block->inptr[0];
  double* u2=block->inptr[1];
  int* outsz=block->outsz;
  int* insz=block->insz;
    
  lb=nx*nx;
  lc=lb;
  
  if (flag ==1 || flag ==6){
    /* y=c*x*/  
    C2F(dmmul)(&rpar[lc],outsz,x,&nx,y,outsz,outsz,&nx,&un);    
  }else if (flag == 2 && block->nevprt == 1){
    
    /* x+=u2 */
    memcpy(x,u2 ,nx*sizeof(double));
  }else if (flag ==0 && block->nevprt == 0){
    /* xd=a*x*/
    C2F(dmmul)(&rpar[0],&nx,x,&nx,xd,&nx,&nx,&nx,&un);
  }
}

---