prod
/* Copyright INRIA */
#include "../machine.h"
/* Element wise product */
void
prod(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[];
integer *nin,outsz[],*nout;
double x[],xd[],z[],tvec[],rpar[];
double *inptr[],*outptr[],*t;
{
int k,i,n;
double *y;
double *u;
n=outsz[0]; /* insz[0]==insz[1] .. ==insz[*nin]== outsz[0] */
y=(double *)outptr[0];
for (i=0;i<n;i++) {
y[i]=1.0;
for (k=0;k<*nin;k++) {
u=(double *)inptr[k];
y[i]=y[i]*u[i];
}
}
}