projects
/
physik
/
dfe.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
|
inline
| side by side (from parent 1:
7b6edf3
)
some changes
author
hackbard
<hackbard>
Tue, 16 Nov 2004 17:39:25 +0000
(17:39 +0000)
committer
hackbard
<hackbard>
Tue, 16 Nov 2004 17:39:25 +0000
(17:39 +0000)
dfe.c
patch
|
blob
|
history
diff --git
a/dfe.c
b/dfe.c
index
42e2abf
..
c507fa8
100644
(file)
--- a/
dfe.c
+++ b/
dfe.c
@@
-28,13
+28,13
@@
typedef struct s_gp {
double pp;
} t_gp;
double pp;
} t_gp;
-#define RANGE
10
+#define RANGE
4
-#define X_ 50
0
-#define Y_
5
0
+#define X_ 50
+#define Y_
1
0
-#define STEPS 100
0
-#define PLOT 1
0
+#define STEPS 100
+#define PLOT 1
int main(int argc,char **argv) {
int main(int argc,char **argv) {
@@
-112,7
+112,7
@@
int main(int argc,char **argv) {
if((dx!=0)||(dy!=0)) {
x=(gp_ptr->x+X_+dx)%X_;
y=gp_ptr->y+dy;
if((dx!=0)||(dy!=0)) {
x=(gp_ptr->x+X_+dx)%X_;
y=gp_ptr->y+dy;
- if(y>=0) {
+ if(y>=0
&&y<Y_
) {
gp.x=x; gp.y=y;
if(list_search_data(&list,&gp,2*sizeof(int))==L_SUCCESS)
gp_ptr->pp+=exp(-1.0*sqrt(x*x+y*y)/xi);
gp.x=x; gp.y=y;
if(list_search_data(&list,&gp,2*sizeof(int))==L_SUCCESS)
gp_ptr->pp+=exp(-1.0*sqrt(x*x+y*y)/xi);
@@
-134,7
+134,6
@@
int main(int argc,char **argv) {
sf[0]=0; sf[1]=0;
tmp=list.current;
sf[0]=0; sf[1]=0;
tmp=list.current;
- bf=0;
gp_ptr=(t_gp *)list.current->data;
for(dx=-RANGE;dx<=RANGE;dx++) {
for(dy=-RANGE;dy<=RANGE;dy++) {
gp_ptr=(t_gp *)list.current->data;
for(dx=-RANGE;dx<=RANGE;dx++) {
for(dy=-RANGE;dy<=RANGE;dy++) {
@@
-157,7
+156,7
@@
int main(int argc,char **argv) {
else {
/* und hier wahrscheinlich auch */
if(dy==0) lf[0]+=(-1.0*(dx/abs(dx))*(by+co));
else {
/* und hier wahrscheinlich auch */
if(dy==0) lf[0]+=(-1.0*(dx/abs(dx))*(by+co));
- if(dx==0) lf[1]+=(-1.0*(dy/abs(dy))*(by+co));
+
else
if(dx==0) lf[1]+=(-1.0*(dy/abs(dy))*(by+co));
else{
lf[0]+=(-1.0*dx/abs(dy)*(by+co));
lf[1]+=(-1.0*dy/abs(dx)*(by+co));
else{
lf[0]+=(-1.0*dx/abs(dy)*(by+co));
lf[1]+=(-1.0*dy/abs(dx)*(by+co));
@@
-171,7
+170,7
@@
int main(int argc,char **argv) {
gp_ptr->fx=sf[0]+lf[0]+((1.0*rand()/RAND_MAX)-0.5)/rs;
gp_ptr->fy=sf[1]+lf[1]+((1.0*rand()/RAND_MAX)-0.5)/rs;
if(gp_ptr->status&FEST) {
gp_ptr->fx=sf[0]+lf[0]+((1.0*rand()/RAND_MAX)-0.5)/rs;
gp_ptr->fy=sf[1]+lf[1]+((1.0*rand()/RAND_MAX)-0.5)/rs;
if(gp_ptr->status&FEST) {
- gp_ptr->fy+=(bf*exp(-y/100));
+ gp_ptr->fy+=(bf*exp(-y
*1.0
/100));
if(gp_ptr->fx>0) {
gp_ptr->fx-=gp_ptr->pp;
if(gp_ptr->fx<0) gp_ptr->fx=0;
if(gp_ptr->fx>0) {
gp_ptr->fx-=gp_ptr->pp;
if(gp_ptr->fx<0) gp_ptr->fx=0;
@@
-186,7
+185,7
@@
int main(int argc,char **argv) {
}
else {
gp_ptr->fy+=gp_ptr->pp;
}
else {
gp_ptr->fy+=gp_ptr->pp;
- if(gp_ptr->fy>0) gp_ptr->f
x
=0;
+ if(gp_ptr->fy>0) gp_ptr->f
y
=0;
}
}
if(gp_ptr->y==0) gp_ptr->fy+=by; /* wenn rand by dazu, arschloch !!1 */
}
}
if(gp_ptr->y==0) gp_ptr->fy+=by; /* wenn rand by dazu, arschloch !!1 */