annotate main_GL.cbc @ 0:249965d0a68f

GL is no finish.
author kent
date Thu, 29 May 2008 19:35:24 +0900
parents
children 09e774f4433f
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
249965d0a68f GL is no finish.
kent
parents:
diff changeset
1 #include<stdio.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
2 #include<stdlib.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
3 #include<unistd.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
4 #include<time.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
5 #include<float.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
6 #include<SDL.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
7 #include<OpenGL/gl.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
8 #include<OpenGL/glu.h>
249965d0a68f GL is no finish.
kent
parents:
diff changeset
9
249965d0a68f GL is no finish.
kent
parents:
diff changeset
10 #define DEBUGlog(f, args...) \
249965d0a68f GL is no finish.
kent
parents:
diff changeset
11 fprintf(stderr, "in %s: "f, __FUNCTION__, ## args)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
12
249965d0a68f GL is no finish.
kent
parents:
diff changeset
13 #define W_HEIGHT 480
249965d0a68f GL is no finish.
kent
parents:
diff changeset
14 #define W_WIDTH 640
249965d0a68f GL is no finish.
kent
parents:
diff changeset
15
249965d0a68f GL is no finish.
kent
parents:
diff changeset
16 /*
249965d0a68f GL is no finish.
kent
parents:
diff changeset
17 N body problem
249965d0a68f GL is no finish.
kent
parents:
diff changeset
18 */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
19
249965d0a68f GL is no finish.
kent
parents:
diff changeset
20 /* parameters */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
21 static int NUM_BODY = 3;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
22 //static float Gravitation = 6.67e-11 ;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
23 //static float delta = 100;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
24 //static float FIELD = 2e11;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
25 static float Gravitation = 0.2f; // ?
249965d0a68f GL is no finish.
kent
parents:
diff changeset
26 static float delta = 100.0f; // 0.01 ~ 100 ?
249965d0a68f GL is no finish.
kent
parents:
diff changeset
27 static float FIELD = 400.0f; // -100 ~ 100
249965d0a68f GL is no finish.
kent
parents:
diff changeset
28 static const float eps = 0.0f;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
29
249965d0a68f GL is no finish.
kent
parents:
diff changeset
30 /* for OpenGL Utility. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
31 GLUquadricObj **sphere;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
32
249965d0a68f GL is no finish.
kent
parents:
diff changeset
33 typedef struct
249965d0a68f GL is no finish.
kent
parents:
diff changeset
34 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
35 /* star's parameter. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
36 float weight;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
37 float a[3]; /* acceleration */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
38 float v[3]; /* velocity */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
39 float r[3]; /* position */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
40 /* for SDL. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
41 //SDL_Rect rect;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
42 } body;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
43 body *stars_old, *stars_new;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
44
249965d0a68f GL is no finish.
kent
parents:
diff changeset
45
249965d0a68f GL is no finish.
kent
parents:
diff changeset
46 void start(void);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
47 __code initialize(int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
48 __code randomInit(SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
49 __code starsInit(SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
50 __code loop(int count, SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
51 __code compute(int count, SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
52 __code nextTurn(int count, SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
53 __code moveCenter(int count, SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
54 __code CenteringVelocity(int count, SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
55 __code drawStars(int count, SDL_Surface *screen, int num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
56
249965d0a68f GL is no finish.
kent
parents:
diff changeset
57 int main(int argc, char **argv)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
58 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
59 int ch;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
60 while ((ch = getopt(argc, argv, "s:g:")) != -1) {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
61 switch (ch) {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
62 case 's':
249965d0a68f GL is no finish.
kent
parents:
diff changeset
63 delta = atof(optarg);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
64 break;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
65 case 'g':
249965d0a68f GL is no finish.
kent
parents:
diff changeset
66 Gravitation = atof(optarg);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
67 break;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
68 case '?':
249965d0a68f GL is no finish.
kent
parents:
diff changeset
69 default:
249965d0a68f GL is no finish.
kent
parents:
diff changeset
70 break;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
71 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
72 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
73 start();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
74 return 0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
75 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
76 void start()
249965d0a68f GL is no finish.
kent
parents:
diff changeset
77 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
78 goto initialize(NUM_BODY);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
79 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
80
249965d0a68f GL is no finish.
kent
parents:
diff changeset
81 __code finish(int ret, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
82 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
83 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
84 DEBUGlog("Gravity = %e\n", Gravitation);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
85 DEBUGlog("FLT_MAX = %e\n", FLT_MAX);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
86 DEBUGlog("FLT_MAX_EXP = %d\n", FLT_MAX_EXP);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
87 DEBUGlog("FLT_MIN = %e\n", FLT_MIN);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
88 DEBUGlog("FLT_MIN_EXP = %d\n", FLT_MIN_EXP);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
89 DEBUGlog("FLT_EPSILON = %e\n", FLT_EPSILON);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
90 free(stars_old);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
91 free(stars_new);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
92 for( i=0; i<NUM_BODY; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
93 gluDeleteQuadric(sphere[i]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
94 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
95 free(sphere);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
96 exit(ret);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
97 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
98
249965d0a68f GL is no finish.
kent
parents:
diff changeset
99
249965d0a68f GL is no finish.
kent
parents:
diff changeset
100 __code initialize(int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
101 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
102 SDL_Surface *screen;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
103 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
104
249965d0a68f GL is no finish.
kent
parents:
diff changeset
105 /* malloc. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
106 stars_old = malloc( sizeof(body)*num );
249965d0a68f GL is no finish.
kent
parents:
diff changeset
107 stars_new = malloc( sizeof(body)*num );
249965d0a68f GL is no finish.
kent
parents:
diff changeset
108 sphere = malloc( sizeof(GLUquadricObj*)*num );
249965d0a68f GL is no finish.
kent
parents:
diff changeset
109 if (stars_old==NULL||stars_new==NULL||sphere==NULL){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
110 perror("malloc");
249965d0a68f GL is no finish.
kent
parents:
diff changeset
111 goto finish(1, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
112 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
113
249965d0a68f GL is no finish.
kent
parents:
diff changeset
114 /* SDL initialization. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
115 if(SDL_Init(SDL_INIT_VIDEO) < 0){ //Could we start SDL_VIDEO?
249965d0a68f GL is no finish.
kent
parents:
diff changeset
116 fprintf(stderr,"Couldn't init SDL"); //Nope, output to stderr and quit
249965d0a68f GL is no finish.
kent
parents:
diff changeset
117 exit(1);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
118 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
119 SDL_GL_SetAttribute( SDL_GL_RED_SIZE, 5);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
120 SDL_GL_SetAttribute( SDL_GL_GREEN_SIZE, 5);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
121 SDL_GL_SetAttribute( SDL_GL_BLUE_SIZE, 5);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
122 SDL_GL_SetAttribute( SDL_GL_DEPTH_SIZE, 16);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
123 SDL_GL_SetAttribute( SDL_GL_DOUBLEBUFFER, 1);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
124
249965d0a68f GL is no finish.
kent
parents:
diff changeset
125 screen = SDL_SetVideoMode(W_WIDTH, W_HEIGHT, 32, SDL_HWSURFACE | SDL_RESIZABLE | SDL_OPENGL ); //Create a 640x480x32 resizable window
249965d0a68f GL is no finish.
kent
parents:
diff changeset
126 atexit(SDL_Quit); //Now that we're enabled, make sure we cleanup
249965d0a68f GL is no finish.
kent
parents:
diff changeset
127
249965d0a68f GL is no finish.
kent
parents:
diff changeset
128 /* initialize OpenGL. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
129 glViewport(0,0,screen->w,screen->h);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
130 glMatrixMode(GL_PROJECTION);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
131 glLoadIdentity();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
132 gluPerspective( 60.0, (float)screen->w/(float)screen->h, 1.0, 1000.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
133 gluLookAt( 500.0,500.0,500.0, 0.0,0.0,0.0, 1.0,0.0,0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
134 glClearColor(0.0, 0.0, 0.0, 0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
135 glMatrixMode(GL_MODELVIEW);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
136
249965d0a68f GL is no finish.
kent
parents:
diff changeset
137 //DEBUGlog("scr->w=%d\n", screen->w);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
138 //DEBUGlog("scr->h=%d\n", screen->h);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
139 for( i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
140 sphere[i] = gluNewQuadric();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
141 gluQuadricDrawStyle(sphere[i], GLU_FILL);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
142 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
143
249965d0a68f GL is no finish.
kent
parents:
diff changeset
144 goto starsInit(screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
145 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
146
249965d0a68f GL is no finish.
kent
parents:
diff changeset
147 __code starsInit0(SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
148 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
149 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
150 srandom(time(NULL));
249965d0a68f GL is no finish.
kent
parents:
diff changeset
151 for (i=0; i<num; i++){ // this loop should be split into few code segment..
249965d0a68f GL is no finish.
kent
parents:
diff changeset
152 stars_old[i].weight = random()/(RAND_MAX+1.0)*5+5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
153 stars_old[i].v[0] = random()/(RAND_MAX+1.0)*5+5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
154 stars_old[i].v[1] = random()/(RAND_MAX+1.0)*5+5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
155 stars_old[i].v[2] = random()/(RAND_MAX+1.0)*5+5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
156 stars_old[i].r[0] = random()/(RAND_MAX+1.0)*5+5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
157 stars_old[i].r[1] = random()/(RAND_MAX+1.0)*5+5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
158 stars_old[i].r[2] = random()/(RAND_MAX+1.0)*5+5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
159 stars_new[i].weight = stars_old[i].weight;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
160 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
161
249965d0a68f GL is no finish.
kent
parents:
diff changeset
162 goto loop(0, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
163 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
164 __code starsInit(SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
165 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
166 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
167 /* */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
168 stars_old[0].weight = 110;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
169 stars_old[0].v[0] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
170 stars_old[0].v[1] = -1.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
171 stars_old[0].v[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
172 stars_old[0].r[0] = 100.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
173 stars_old[0].r[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
174 stars_old[0].r[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
175 /* */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
176 stars_old[1].weight = 110;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
177 stars_old[1].v[0] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
178 stars_old[1].v[1] = -1.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
179 stars_old[1].v[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
180 stars_old[1].r[0] = -100.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
181 stars_old[1].r[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
182 stars_old[1].r[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
183 /* */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
184 stars_old[2].weight = 110;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
185 stars_old[2].v[0] = -1.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
186 stars_old[2].v[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
187 stars_old[2].v[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
188 stars_old[2].r[0] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
189 stars_old[2].r[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
190 stars_old[2].r[2] = -70.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
191
249965d0a68f GL is no finish.
kent
parents:
diff changeset
192 for( i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
193 stars_new[i].weight = stars_old[i].weight;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
194 //stars_new[i].rect.h = 5, stars_new[i].rect.w = 5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
195 //stars_old[i].rect.h = 5, stars_old[i].rect.w = 5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
196 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
197
249965d0a68f GL is no finish.
kent
parents:
diff changeset
198 goto loop(0, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
199 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
200
249965d0a68f GL is no finish.
kent
parents:
diff changeset
201 __code starsInit1(SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
202 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
203 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
204 /* Sun */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
205 stars_old[0].weight = 1.9891e30; // 1.9891*10^30
249965d0a68f GL is no finish.
kent
parents:
diff changeset
206 stars_old[0].v[0] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
207 stars_old[0].v[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
208 stars_old[0].v[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
209 stars_old[0].r[0] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
210 stars_old[0].r[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
211 stars_old[0].r[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
212 /* Venus */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
213 stars_old[1].weight = 4.869e24; // 4.869*10^24
249965d0a68f GL is no finish.
kent
parents:
diff changeset
214 stars_old[1].v[0] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
215 stars_old[1].v[1] = 3.50214e4; // 35.0214 km/s
249965d0a68f GL is no finish.
kent
parents:
diff changeset
216 stars_old[1].v[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
217 stars_old[1].r[0] = 1.08e11; // 108,208,930 km
249965d0a68f GL is no finish.
kent
parents:
diff changeset
218 stars_old[1].r[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
219 stars_old[1].r[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
220 /* Earth */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
221 stars_old[2].weight = 5.9742e24; // 5.9742*10^24
249965d0a68f GL is no finish.
kent
parents:
diff changeset
222 stars_old[2].v[0] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
223 stars_old[2].v[1] = 2.97859e4; // 29.7859 km/s
249965d0a68f GL is no finish.
kent
parents:
diff changeset
224 stars_old[2].v[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
225 stars_old[2].r[0] = 1.49e11; // 149,597,870km
249965d0a68f GL is no finish.
kent
parents:
diff changeset
226 stars_old[2].r[1] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
227 stars_old[2].r[2] = 0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
228
249965d0a68f GL is no finish.
kent
parents:
diff changeset
229 for( i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
230 stars_new[i].weight = stars_old[i].weight;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
231 //stars_new[i].rect.h = 5, stars_new[i].rect.w = 5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
232 //stars_old[i].rect.h = 5, stars_old[i].rect.w = 5;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
233 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
234
249965d0a68f GL is no finish.
kent
parents:
diff changeset
235 goto loop(0, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
236 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
237
249965d0a68f GL is no finish.
kent
parents:
diff changeset
238 __code loop(int count, SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
239 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
240 SDL_Event event;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
241
249965d0a68f GL is no finish.
kent
parents:
diff changeset
242 /* check SDL event. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
243 while(SDL_PollEvent(&event)){ //Poll events
249965d0a68f GL is no finish.
kent
parents:
diff changeset
244 switch(event.type){ //Check event type
249965d0a68f GL is no finish.
kent
parents:
diff changeset
245 case SDL_QUIT: //User hit the X (or equivelent)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
246 goto finish(1, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
247 //case SDL_VIDEORESIZE: //User resized window
249965d0a68f GL is no finish.
kent
parents:
diff changeset
248 //screen = SDL_SetVideoMode(event.resize.w, event.resize.h, 32,
249965d0a68f GL is no finish.
kent
parents:
diff changeset
249 //SDL_HWSURFACE | SDL_RESIZABLE); // Create new window
249965d0a68f GL is no finish.
kent
parents:
diff changeset
250 //break; //Event handled, fetch next :)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
251 case SDL_KEYDOWN:
249965d0a68f GL is no finish.
kent
parents:
diff changeset
252 if (event.key.keysym.sym==SDLK_UP && event.key.state==SDL_PRESSED)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
253 FIELD *= 2.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
254 else if (event.key.keysym.sym==SDLK_DOWN && event.key.state==SDL_PRESSED)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
255 FIELD /= 2.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
256 else if (event.key.keysym.sym==SDLK_r && event.key.state==SDL_PRESSED)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
257 goto moveCenter(count, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
258 else if (event.key.keysym.sym==SDLK_v && event.key.state==SDL_PRESSED)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
259 goto CenteringVelocity(count, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
260 else if (event.key.keysym.sym==SDLK_l && event.key.state==SDL_PRESSED)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
261 goto AdjustLook(count, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
262 break;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
263 default:
249965d0a68f GL is no finish.
kent
parents:
diff changeset
264 break;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
265 } //Finished with current event
249965d0a68f GL is no finish.
kent
parents:
diff changeset
266 } //Done with all events for now
249965d0a68f GL is no finish.
kent
parents:
diff changeset
267
249965d0a68f GL is no finish.
kent
parents:
diff changeset
268 if ( count<num ){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
269 DEBUGlog("count %d, goto commpute().\n", count);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
270 goto compute(count, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
271 }else{
249965d0a68f GL is no finish.
kent
parents:
diff changeset
272 DEBUGlog("count %d, goto nextTurn()\n", count);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
273 goto nextTurn(count, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
274 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
275 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
276
249965d0a68f GL is no finish.
kent
parents:
diff changeset
277 __code AdjustLook(int count, SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
278 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
279 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
280 float c0,c1,c2;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
281 float v0,v1,v2,v;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
282 c0=c1=c2=v0=v1=v2=0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
283
249965d0a68f GL is no finish.
kent
parents:
diff changeset
284 for (i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
285 c0 += stars_new[i].r[0];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
286 c1 += stars_new[i].r[1];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
287 c2 += stars_new[i].r[2];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
288 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
289 c0/=3.0; c1/=3.0; c2/=3.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
290 for (i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
291 v0 += (stars_new[i].r[0]-c0)*(stars_new[i].r[0]-c0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
292 v1 += (stars_new[i].r[1]-c1)*(stars_new[i].r[1]-c1);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
293 v2 += (stars_new[i].r[2]-c2)*(stars_new[i].r[2]-c2);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
294 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
295 v = v0+v1+v2;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
296 v0 = sqrt(v0); v1 = sqrt(v1); v2 = sqrt(v2);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
297 v = sqrt(v);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
298
249965d0a68f GL is no finish.
kent
parents:
diff changeset
299 glMatrixMode(GL_PROJECTION);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
300 glLoadIdentity();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
301 gluPerspective( 60.0, (float)screen->w/(float)screen->h, 1.0, v);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
302 gluLookAt( v,v,v, 0.0,0.0,0.0, 1.0,0.0,0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
303 glClearColor(0.0, 0.0, 0.0, 0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
304 glMatrixMode(GL_MODELVIEW);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
305
249965d0a68f GL is no finish.
kent
parents:
diff changeset
306 goto loop(0, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
307 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
308
249965d0a68f GL is no finish.
kent
parents:
diff changeset
309 __code CenteringVelocity(int count, SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
310 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
311 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
312 float v0,v1,v2;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
313 v0=v1=v2=0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
314
249965d0a68f GL is no finish.
kent
parents:
diff changeset
315 for (i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
316 v0 += stars_old[i].v[0];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
317 v1 += stars_old[i].v[1];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
318 v2 += stars_old[i].v[2];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
319 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
320 v0/=(float)num; v1/=(float)num; v2/=(float)num;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
321 for (i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
322 stars_old[i].v[0] -= v0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
323 stars_old[i].v[1] -= v1;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
324 stars_old[i].v[2] -= v2;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
325 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
326
249965d0a68f GL is no finish.
kent
parents:
diff changeset
327 goto loop(0, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
328 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
329 __code moveCenter(int count, SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
330 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
331 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
332 float m0,m1,m2;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
333 m0=m1=m2=0.0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
334
249965d0a68f GL is no finish.
kent
parents:
diff changeset
335 for (i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
336 m0 += stars_old[i].r[0];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
337 m1 += stars_old[i].r[1];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
338 m2 += stars_old[i].r[2];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
339 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
340 m0/=(float)num; m1/=(float)num; m2/=(float)num;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
341 for (i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
342 stars_old[i].r[0] -= m0;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
343 stars_old[i].r[1] -= m1;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
344 stars_old[i].r[2] -= m2;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
345 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
346
249965d0a68f GL is no finish.
kent
parents:
diff changeset
347 goto loop(0, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
348 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
349
249965d0a68f GL is no finish.
kent
parents:
diff changeset
350 /*
249965d0a68f GL is no finish.
kent
parents:
diff changeset
351 * next x = x+dx
249965d0a68f GL is no finish.
kent
parents:
diff changeset
352 * dx = v*dt => dx/dt = v
249965d0a68f GL is no finish.
kent
parents:
diff changeset
353 * dv = a*dt => dv/dt = a
249965d0a68f GL is no finish.
kent
parents:
diff changeset
354 * a = F/m;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
355 * F = F1 + F2 + ... + Fn
249965d0a68f GL is no finish.
kent
parents:
diff changeset
356 * Fi = G m1m2/r^2
249965d0a68f GL is no finish.
kent
parents:
diff changeset
357 *
249965d0a68f GL is no finish.
kent
parents:
diff changeset
358 */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
359 __code compute(int count, SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
360 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
361 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
362 /* a is accel this planet receive now. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
363 stars_old[count].a[0]=stars_old[count].a[1]=stars_old[count].a[2]=0.0f;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
364 DEBUGlog("count=%d\n", count);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
365 for (i=0; i<num; i++){ // this loop should be split into few code segment..
249965d0a68f GL is no finish.
kent
parents:
diff changeset
366 //float F;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
367 float d0, d1, d2, d;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
368 float a;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
369 //body *o, *m;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
370 //o = &stars_old[i]; m = &stars_old[count];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
371
249965d0a68f GL is no finish.
kent
parents:
diff changeset
372 /* skip compute with itself. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
373 if ( i==count ) continue;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
374 /* compute distance between two i-th planet and itself. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
375 d0 = stars_old[i].r[0] - stars_old[count].r[0];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
376 d1 = stars_old[i].r[1] - stars_old[count].r[1];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
377 d2 = stars_old[i].r[2] - stars_old[count].r[2];
249965d0a68f GL is no finish.
kent
parents:
diff changeset
378 d = ( d0*d0+d1*d1+d2*d2 );
249965d0a68f GL is no finish.
kent
parents:
diff changeset
379 /* compute force it receive from i-th planet. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
380 //F = Gravitation * stars_old[i].weight * stars_old[count].weight / d;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
381 /* and accel. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
382 //a = F/stars_old[count].weight;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
383 a = Gravitation/d * stars_old[i].weight ;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
384 stars_old[count].a[0] += a*d0/sqrt(d);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
385 stars_old[count].a[1] += a*d1/sqrt(d);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
386 stars_old[count].a[2] += a*d2/sqrt(d);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
387 DEBUGlog("a=%e, d=%e, d0=%e\n", a, d, d0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
388 DEBUGlog("g*w=%e\n", Gravitation*stars_old[i].weight);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
389 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
390
249965d0a68f GL is no finish.
kent
parents:
diff changeset
391 stars_new[count].v[0] = stars_old[count].v[0]+stars_old[count].a[0]*delta;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
392 stars_new[count].v[1] = stars_old[count].v[1]+stars_old[count].a[1]*delta;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
393 stars_new[count].v[2] = stars_old[count].v[2]+stars_old[count].a[2]*delta;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
394 stars_new[count].r[0] = stars_old[count].r[0]+stars_new[count].v[0]*delta/*+stars_old[count].a[0]*delta*delta/2.0*/;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
395 stars_new[count].r[1] = stars_old[count].r[1]+stars_new[count].v[1]*delta/*+stars_old[count].a[1]*delta*delta/2.0*/;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
396 stars_new[count].r[2] = stars_old[count].r[2]+stars_new[count].v[2]*delta/*+stars_old[count].a[2]*delta*delta/2.0*/;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
397
249965d0a68f GL is no finish.
kent
parents:
diff changeset
398 DEBUGlog("x = %e\n", stars_new[count].r[0]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
399 DEBUGlog("y = %e\n", stars_new[count].r[1]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
400 DEBUGlog("z = %e\n", stars_new[count].r[2]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
401 DEBUGlog("vx = %e\n", stars_new[count].v[0]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
402 DEBUGlog("vy = %e\n", stars_new[count].v[1]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
403 DEBUGlog("vz = %e\n", stars_new[count].v[2]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
404 DEBUGlog("a0=%e,a1=%e,a2=%e\n", stars_new[count].a[0], stars_new[count].a[1], stars_new[count].a[2]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
405
249965d0a68f GL is no finish.
kent
parents:
diff changeset
406 goto loop(++count, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
407 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
408
249965d0a68f GL is no finish.
kent
parents:
diff changeset
409 /* it is executed after all bodies parameter is computed. */
249965d0a68f GL is no finish.
kent
parents:
diff changeset
410 __code nextTurn(int count, SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
411 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
412 body *tmp;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
413 tmp = stars_old;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
414 stars_old = stars_new;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
415 stars_new = tmp;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
416 SDL_UpdateRect(screen, 0, 0, 0, 0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
417
249965d0a68f GL is no finish.
kent
parents:
diff changeset
418 goto drawStars(count, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
419 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
420
249965d0a68f GL is no finish.
kent
parents:
diff changeset
421 __code drawStars(int count, SDL_Surface *screen, int num)
249965d0a68f GL is no finish.
kent
parents:
diff changeset
422 {
249965d0a68f GL is no finish.
kent
parents:
diff changeset
423 int i;
249965d0a68f GL is no finish.
kent
parents:
diff changeset
424 DEBUGlog("draw\n");
249965d0a68f GL is no finish.
kent
parents:
diff changeset
425 glClear( GL_COLOR_BUFFER_BIT );
249965d0a68f GL is no finish.
kent
parents:
diff changeset
426 for (i=0; i<num; i++){
249965d0a68f GL is no finish.
kent
parents:
diff changeset
427 glPushMatrix();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
428 glTranslatef( stars_new[i].r[0], stars_new[i].r[1], stars_new[i].r[2]);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
429 gluSphere( sphere[i], 200.0, 8.0, 8.0 );
249965d0a68f GL is no finish.
kent
parents:
diff changeset
430 glPopMatrix();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
431 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
432
249965d0a68f GL is no finish.
kent
parents:
diff changeset
433 glColor3d( 1.0, 1.0, 1.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
434 glBegin(GL_LINES);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
435 glVertex3d( 0.0, 0.0, 0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
436 glVertex3d( 10.0, 0.0, 0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
437 glEnd();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
438 glBegin(GL_LINES);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
439 glVertex3d( 0.0, 0.0, 0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
440 glVertex3d( 0.0, 10.0, 0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
441 glEnd();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
442 glBegin(GL_LINES);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
443 glVertex3d( 0.0, 0.0, 0.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
444 glVertex3d( 0.0, 0.0, 10.0);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
445 glEnd();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
446
249965d0a68f GL is no finish.
kent
parents:
diff changeset
447 SDL_GL_SwapBuffers();
249965d0a68f GL is no finish.
kent
parents:
diff changeset
448 goto loop(0, screen, num);
249965d0a68f GL is no finish.
kent
parents:
diff changeset
449 }
249965d0a68f GL is no finish.
kent
parents:
diff changeset
450