relay_routcos

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++);
          }
        }
    }