From 7b6edf3294e792bd7cef93f6ecfc71726b430296 Mon Sep 17 00:00:00 2001 From: hackbard Date: Tue, 16 Nov 2004 16:40:28 +0000 Subject: [PATCH 1/1] initial checkin --- .cvsignore | 4 + dfe.c | 242 +++++++++++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 246 insertions(+) create mode 100644 .cvsignore create mode 100644 dfe.c diff --git a/.cvsignore b/.cvsignore new file mode 100644 index 0000000..a8b9ca8 --- /dev/null +++ b/.cvsignore @@ -0,0 +1,4 @@ +list.* +bmp.* +dfe +*.bmp diff --git a/dfe.c b/dfe.c new file mode 100644 index 0000000..42e2abf --- /dev/null +++ b/dfe.c @@ -0,0 +1,242 @@ +#define _GNU_SOURCE +#include +#include +#include +#include + +/* + * compile: gcc -Wall -o dfe dfe.c bmp.c list.c -O3 -lm + * + * usage: ./dfe 100 50 2 1 1 10 + * + * author: christian leirer, frank zirkelbach + * + * need api: cvs -d:pserver:anonymous@hackdaworld.dyndns.org:/my-code co api + * + */ + +#include "list.h" +#include "bmp.h" + +typedef struct s_gp { + int x,y; + unsigned char status; +#define FLUESSIG 0 +#define FEST 1 + double fx,fy; + double op; + double pp; +} t_gp; + +#define RANGE 10 + +#define X_ 500 +#define Y_ 50 + +#define STEPS 1000 +#define PLOT 10 + +int main(int argc,char **argv) { + + int x,y; + int dx,dy; + int mx,my; + double lf[2]; + double sf[2]; + double mf; + double lpp; + int i; + int size; + t_gp gp; + t_gp *gp_ptr; + t_list list; + t_list_element *tmp; + t_bmp bmp; + char string[128]; + int fd; + + double by; + double pp; + double xi; + double rs; + double co; + double bf; + + /* parse argv */ + + if(argc!=7) { + printf("usage: %s \n", + argv[0]); + return -1; + } + + by=atof(argv[1]); + pp=atof(argv[2]); + xi=atof(argv[3]); + rs=atof(argv[4]); + co=atof(argv[5]); + bf=atof(argv[6]); + + /* init */ + + fd=open("/dev/null",O_WRONLY); + + list_init(&list,fd); + + size=X_*Y_; + bmp_init(&bmp,1); + bmp.width=X_; + bmp.height=Y_; + bmp.mode=WRITE; + bmp_alloc_map(&bmp); + + srandom(time(NULL)); + + for(x=0;xdata; + tmp=list.current; + for(dx=-RANGE;dx<=RANGE;dx++) { + for(dy=-RANGE;dy<=RANGE;dy++) { + if((dx!=0)||(dy!=0)) { + x=(gp_ptr->x+X_+dx)%X_; + y=gp_ptr->y+dy; + if(y>=0) { + 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); + + } + } + } + } + if(lpp>(0.8*pp)) gp_ptr->status=FLUESSIG; + else gp_ptr->status=FEST; + list.current=tmp; + } + while(list_next(&list)!=L_NO_NEXT_ELEMENT); + + /* force on vortex */ + list_reset(&list); + do { + lf[0]=0; lf[1]=0; + 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++) { + if((dx!=0)||(dy!=0)) { + x=(gp_ptr->x+X_+dx)%X_; + y=gp_ptr->y+dy; + if(y>=0&&ydata,sizeof(t_gp)); + if(gp.status&FEST) { + if(dy==0) sf[0]+=(-1.0/(dx*dx+dy*dy)*co*(dx/abs(dx))); + else if (dx==0) sf[1]+=(-1.0/(dx*dx+dy*dy)*co*(dy/abs(dy))); + /* hier ist die scheisse drin!! */ + else { + sf[0]+=(-1.0/(dx*dx+dy*dy)*co*dx/abs(dy)); + sf[1]+=(-1.0/(dx*dx+dy*dy)*co*dy/abs(dx)); + } + } + 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{ + lf[0]+=(-1.0*dx/abs(dy)*(by+co)); + lf[1]+=(-1.0*dy/abs(dx)*(by+co)); + } + } + } + } + } + } + } + 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)); + if(gp_ptr->fx>0) { + gp_ptr->fx-=gp_ptr->pp; + if(gp_ptr->fx<0) gp_ptr->fx=0; + } + else { + gp_ptr->fx+=gp_ptr->pp; + if(gp_ptr->fx>0) gp_ptr->fx=0; + } + if(gp_ptr->fy>0) { + gp_ptr->fy-=gp_ptr->pp; + if(gp_ptr->fy<0) gp_ptr->fy=0; + } + else { + gp_ptr->fy+=gp_ptr->pp; + if(gp_ptr->fy>0) gp_ptr->fx=0; + } + } + if(gp_ptr->y==0) gp_ptr->fy+=by; /* wenn rand by dazu, arschloch !!1 */ + if((mffx)||(mffy)) { + mf=gp_ptr->fx>gp_ptr->fy?gp_ptr->fx:gp_ptr->fy; + mx=gp_ptr->x; + my=gp_ptr->y; + } + list.current=tmp; + } + while(list_next(&list)!=L_NO_NEXT_ELEMENT); + + /* move vortex with highest force */ + gp.x=mx; gp.y=my; + printf("step %d: move vortex %d %d,",i,mx,my); + list_search_data(&list,&gp,2*sizeof(int)); + gp_ptr=(t_gp *)list.current->data; + dx=0; dy=0; + if(fabs(gp_ptr->fx)>fabs(gp_ptr->fy)) { + if(gp_ptr->fx>0) dx=1; + else dx=-1; + } + else { + if(gp_ptr->fy>0) dy=1; + else dy=-1; + } + printf(" with direction dx=%d dy=%d | force: %f %f\n",dx,dy,gp_ptr->fx,gp_ptr->fy); + gp_ptr->x+=dx; + gp_ptr->y+=dy; + if(gp.y==0) list_add_element(&list,&gp,sizeof(t_gp)); + + /* plot every PLOT steps */ + if(i%PLOT==0) { + memset(bmp.map,0,3*size*sizeof(unsigned char)); + list_reset(&list); + do { + gp_ptr=(t_gp *)list.current->data; + if(gp_ptr->status&FEST) + memset(bmp.map+gp_ptr->y*X_+gp_ptr->x,0xff,1); + else + memset(bmp.map+gp_ptr->y*X_+gp_ptr->x,0xff,3); + } while(list_next(&list)!=L_NO_NEXT_ELEMENT); + sprintf(string,"dfe_%d_of_%d.bmp",i,STEPS); + strcpy(bmp.file,string); + bmp_write_file(&bmp); + } + + } + + close(fd); + + return 1; +} -- 2.39.2