dot_prod.c

二维热能方程的求解

#include <stdio.h>
#include <math.h>
typedef  double ATYPE;
typedef  float DTYPE;

void dot_prod(nxny,level,x,y,sum,ierr)
int nxny,level;
int *ierr;
ATYPE *sum;
ATYPE x[];
ATYPE y[];
{
   int i;

   *sum=0.0;
   for(i=0;i<nxny;i++)
      *sum+=x[i]*y[i];
}