projects
/
physik
/
nlsop.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
new run
[physik/nlsop.git]
/
nlsop.c
diff --git
a/nlsop.c
b/nlsop.c
index
ae62925
..
683803e
100644
(file)
--- a/
nlsop.c
+++ b/
nlsop.c
@@
-10,7
+10,7
@@
* selforganizing lamella structure in the observed behaviour.
*
* refs:
* selforganizing lamella structure in the observed behaviour.
*
* refs:
- * - J. K. N. Lindner. Habil
ationss
chrift, Universitaet Augsburg.
+ * - J. K. N. Lindner. Habil
.S
chrift, Universitaet Augsburg.
* - Maik Haeberlen. Diplomarbeit, Universitaet Augsburg.
*
* Copyright (C) 2004 Frank Zirkelbach
* - Maik Haeberlen. Diplomarbeit, Universitaet Augsburg.
*
* Copyright (C) 2004 Frank Zirkelbach
@@
-119,6
+119,7
@@
int process_cell(d3_lattice *d3_l,u32 x,u32 y,u32 z,info *my_info,u32 nel_z)
j+=(*(d3_l->status+x+((y-1+d3_l->max_y)%d3_l->max_y)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
j+=(*(d3_l->status+x+y*d3_l->max_x+((z+1+d3_l->max_z)%d3_l->max_z)*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
j+=(*(d3_l->status+x+y*d3_l->max_x+((z-1+d3_l->max_z)%d3_l->max_z)*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
j+=(*(d3_l->status+x+((y-1+d3_l->max_y)%d3_l->max_y)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
j+=(*(d3_l->status+x+y*d3_l->max_x+((z+1+d3_l->max_z)%d3_l->max_z)*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
j+=(*(d3_l->status+x+y*d3_l->max_x+((z-1+d3_l->max_z)%d3_l->max_z)*d3_l->max_x*d3_l->max_y)&AMORPH)?1:0;
+
p+=((q/6)*j);
if(get_rand(URAND_MAX)>p) MAKE_CRYST(thiz);
}
p+=((q/6)*j);
if(get_rand(URAND_MAX)>p) MAKE_CRYST(thiz);
}
@@
-1168,7
+1169,7
@@
int main(int argc,char **argv)
y_c=get_rand(d3_l.max_y);
// z_c=get_rand_reject(d3_l.max_z,ne_max,n_e_loss);
z_c=get_rand(d3_l.max_z);
y_c=get_rand(d3_l.max_y);
// z_c=get_rand_reject(d3_l.max_z,ne_max,n_e_loss);
z_c=get_rand(d3_l.max_z);
- nel_z=
(URAND_MAX*n_e_loss[z_c])/ne_max
;
+ nel_z=
URAND_MAX*(1.0*n_e_loss[z_c]/ne_max)
;
process_cell(&d3_l,x_c,y_c,z_c,&my_info,nel_z);
}
distrib_c(&d3_l,&my_info,i,c_ratio,ip_max,c_profile);
process_cell(&d3_l,x_c,y_c,z_c,&my_info,nel_z);
}
distrib_c(&d3_l,&my_info,i,c_ratio,ip_max,c_profile);