#include "stdio.h"
Go to the source code of this file.
|
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) |
|
t_real | sum (int n1, int n2, int n3, t_real u[n1][n2][n3]) |
|
int | main (int argc, char *argv[]) |
|
◆ t_real
◆ init()
Definition at line 44 of file stencil.c.
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++) {
Referenced by main().
◆ main()
int main |
( |
int |
argc, |
|
|
char * |
argv[] |
|
) |
| |
Definition at line 70 of file stencil.c.
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)
=======================================================================
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)
References init(), is2(), printf(), stencil8(), and sum().
◆ stencil8()
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 |
|
) |
| |
Definition at line 10 of file stencil.c.
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]);
References is2().
Referenced by main().
◆ sum()
Definition at line 57 of file stencil.c.
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++) {
Referenced by bounds_equal_p(), lu_back_substitution(), lu_decomposition(), main(), sc_elim_double_constraints(), size_of_actual_array(), vect_sum(), xsum(), and zeroing_multiple_edges().