11 int is1,
int ie1,
int is2,
int ie2,
16 t_real c_4,c_3,c_1, c_2, c0, c1, c2,c3,c4;
18 c_4 = c[
L-4]; c_3 = c[
L-3]; c_2 = c[
L-2]; c_1 = c[
L-1];
20 c4 = c[
L+4]; c3 = c[
L+3]; c2 = c[
L+2]; c1 = c[
L+1];
23 for (i1=is1+
L; i1<ie1-
L ; i1++) {
25 for (i2=
is2+
L; i2<ie2-
L ; i2++) {
27 for (i3=is3+
L; i3<ie3-
L ; i3++) {
29 c_4 * (v[i1-4][i2][i3] + v[i1][i2-4][i3] + v[i1][i2][i3-4])
30 + c_3 * (v[i1-3][i2][i3] + v[i1][i2-3][i3] + v[i1][i2][i3-3])
31 + c_2 * (v[i1-2][i2][i3] + v[i1][i2-2][i3] + v[i1][i2][i3-2])
32 + c_1 * (v[i1-1][i2][i3] + v[i1][i2-1][i3] + v[i1][i2][i3-1])
33 + c0 * v[i1][ i2][i3] * 3.f
34 + c1 * (v[i1+1][i2][i3] + v[i1][i2+1][i3] + v[i1][i2][i3+1])
35 + c2 * (v[i1+2][i2][i3] + v[i1][i2+2][i3] + v[i1][i2][i3+2])
36 + c3 * (v[i1+3][i2][i3] + v[i1][i2+3][i3] + v[i1][i2][i3+3])
37 + c4 * (v[i1+4][i2][i3] + v[i1][i2+4][i3] + v[i1][i2][i3+4]);
45 int i = 0, j = 0, k = 0;
46 for (i=0; i<n1 ; i++) {
47 for (j=0; j<n2 ; j++) {
48 for (k=0; k<n3 ; k++) {
59 int i = 0, j = 0, k = 0;
60 for (i=0; i<n1 ; i++) {
61 for (j=0; j<n2 ; j++) {
62 for (k=0; k<n3 ; k++) {
70 int main (
int argc,
char * argv[]) {
72 int is1,ie1,
is2,ie2,is3,ie3,i;
77 if(argc >100000) n1=n2=n3=
L=78;
86 for (i=0; i<2*
L+1; i++) {
91 init (n1,n2,n3,u , 1.0f);
92 init (n1,n2,n3,v , 1.0f);
93 stencil8(
L,n1,n2,n3,u,v,c,is1,ie1,
is2,ie2,is3,ie3);
95 printf (
"the sum is : %f\n",
sum (n1,n2,n3,u));
static int is2(Pproblem XX, Pproblem VV, struct rproblem *RR)
=======================================================================
int main(int argc, char *argv[])
t_real sum(int n1, int n2, int n3, t_real u[n1][n2][n3])
void stencil8(int L, int n1, int n2, int n3, t_real u[n1][n2][n3], t_real v[n1][n2][n3], t_real c[2L+1], int is1, int ie1, int is2, int ie2, int is3, int ie3)
void init(int n1, int n2, int n3, t_real u[n1][n2][n3], t_real val)