projects
/
physik
/
morpheus.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
z relaxation, fixed bug in random.c
[physik/morpheus.git]
/
main.c
diff --git
a/main.c
b/main.c
index
2ac8d6a
..
d171490
100644
(file)
--- a/
main.c
+++ b/
main.c
@@
-134,14
+134,15
@@
int process_cell(cell *cell_p,u32 x,u32 y,u32 z,u32 x_max,u32 y_max,u32 z_max,in
{
for(j=-range;j<=range;j++)
{
{
for(j=-range;j<=range;j++)
{
- for(k=-range;k<=range;k++)
- {
- if(!(i==0 && j==0 && k==0))
+ // z relaxation, no pressure in z direction
+ // for(k=-range;k<=range;k++)
+ // {
+ if(!(i==0 && j==0)) // && k==0))
{
{
- // if((cell_p+((x+x_max+i)%x_max)+((y+j+y_max)%y_max)*x_max+((z+k+z_max)%z_max)*x_max*y_max)->status&AMORPH) pressure+=
(cell_p+((x+x_max+i)%x_max)+((y+j+y_max)%y_max)*x_max+((z+k+z_max)%z_max)*x_max*y_max)->conc*
faktor*URAND_2BYTE_MAX/(i*i+j*j+k*k);
- if((cell_p+((x+x_max+i)%x_max)+((y+j+y_max)%y_max)*x_max+
((z+k+z_max)%z_max)*x_max*y_max)->status&AMORPH) pressure+=faktor*URAND_2BYTE_MAX/(i*i+j*j+k*k
);
+ // if((cell_p+((x+x_max+i)%x_max)+((y+j+y_max)%y_max)*x_max+((z+k+z_max)%z_max)*x_max*y_max)->status&AMORPH) pressure+=faktor*URAND_2BYTE_MAX/(i*i+j*j+k*k);
+ if((cell_p+((x+x_max+i)%x_max)+((y+j+y_max)%y_max)*x_max+
z*x_max*y_max)->status&AMORPH) pressure+=faktor*URAND_2BYTE_MAX/(i*i+j*j
);
}
}
- }
+
//
}
}
}
pressure*=this_cell->conc;
}
}
pressure*=this_cell->conc;