void axmb(int m, int n, int ns1, int ns2, double *uu, double *ww, double *xg, double *yg,

  int jsta2, int jend1){

  int i,j;

  double pi, hh, hh2;

  hh = xg[1] - xg[0]; hh2 = hh*hh;

  for(i=0;i<m*(ns2-ns1+1);i++) ww[i] = 0.L;

  

  for(j=1;j<ns2-ns1;j++){

    for(i=1;i<m-1;i++){

      ww[i+m*j] += (-2.L*uu[i+m*j]+uu[i-1+m*j] + uu[i+1+m*j])/hh2;

    }

  }

  hh = yg[1] - yg[0]; hh2 = hh*hh;

  for(j=1;j<ns2-ns1;j++){

    for(i=1;i<m-1;i++){

      ww[i+m*j] += (-2.L*uu[i+m*j]+uu[i+m*(j-1)] + uu[i+m*(j+1)])/hh2;

    }

  }

  pi = M_PI;

  for(j=0;j<ns2-ns1+1;j++){

    for(i=0;i<m;i++){

      ww[i+m*j] -= (-2.l*pi*pi*cos(pi*xg[i])*cos(pi*yg[j+ns1-1]));

    }

  }

}



머리가 빠가라서


슬라이싱이 안되네


자살..