Computational routine
eng

relay

File content

/* Copyright INRIA */
#include "../machine.h"

void 
relay(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;

{

    int k;
    double *y;
    double *u;
    int nev,ic;
    ic=(int)z[0];
    if ((*flag)<3) {
      if ((*nevprt)>0) {
	ic=-1;
	nev=*nevprt;
	while (nev>=1) {
	  ic=ic+1;
	  nev=nev/2;
	}
      }
      if ((*flag)==2) {z[0]=ic;return ;}
      if (*nin>1) {
	y=(double *)outptr[0];
	u=(double *)inptr[ic];
	for (k=0;k<outsz[0];k++)
	  *(y++)=*(u++);  
      }
      else {
	y=(double *)outptr[ic];
	u=(double *)inptr[0];
	for (k=0;k<outsz[0];k++)
	  *(y++)=*(u++);  
      }
    }
}