Mercurial > hg > Members > kent > N-BodyProblem

comparison main.cbc @ 1:09e774f4433f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
main_GL is runnable.
author kent
date Fri, 30 May 2008 00:00:49 +0900
parents 249965d0a68f
children
comparison
equal deleted inserted replaced
0:249965d0a68f 1:09e774f4433f
4 #include<time.h> 4 #include<time.h>
5 #include<SDL.h> 5 #include<SDL.h>
6 #include<float.h> 6 #include<float.h>
7 7
8 #define DEBUGlog(f, args...) \ 8 #define DEBUGlog(f, args...) \
9 ; 9 fprintf(stderr, "in %s: "f, __FUNCTION__, ## args)
10 //fprintf(stderr, "in %s: "f, __FUNCTION__, ## args)
11 10
12 #define W_HEIGHT 480 11 #define W_HEIGHT 480
13 #define W_WIDTH 640 12 #define W_WIDTH 640
14 13
15 /* 14 /*
21 //static float delta = 100; 20 //static float delta = 100;
22 //static float FIELD = 2e11; 21 //static float FIELD = 2e11;
23 static float Gravitation = 1.0f; // ? 22 static float Gravitation = 1.0f; // ?
24 static float delta = 0.05f; // 23 static float delta = 0.05f; //
25 static float FIELD = 400.0f; // -100 ~ 100 24 static float FIELD = 400.0f; // -100 ~ 100
26 static const float eps = 1.5f; 25 static const float eps = 15.0f;
27 26
28 typedef struct 27 typedef struct
29 { 28 {
30 /* star's parameter. */ 29 /* star's parameter. */
31 float weight; 30 float weight;
151 stars_old[2].v[1] = 0.0; 150 stars_old[2].v[1] = 0.0;
152 stars_old[2].v[2] = 0.0; 151 stars_old[2].v[2] = 0.0;
153 stars_old[2].r[0] = 0.0; 152 stars_old[2].r[0] = 0.0;
154 stars_old[2].r[1] = 0.0; 153 stars_old[2].r[1] = 0.0;
155 stars_old[2].r[2] = -70.0; 154 stars_old[2].r[2] = -70.0;
156 #else 155 #elif 0
157 /* */ 156 /* */
158 stars_old[0].weight = 100; 157 stars_old[0].weight = 1000;
159 stars_old[0].v[0] = -1.0; 158 stars_old[0].v[0] = 0.0;
160 stars_old[0].v[1] = 0.0; 159 stars_old[0].v[1] = 0.0;
161 stars_old[0].v[2] = 0.5; 160 stars_old[0].v[2] = 0.0;
162 stars_old[0].r[0] = 50.0; 161 stars_old[0].r[0] = 0.0;
163 stars_old[0].r[1] = 0.0; 162 stars_old[0].r[1] = 0.0;
164 stars_old[0].r[2] = 0.0; 163 stars_old[0].r[2] = 0.0;
165 /* */ 164 /* */
166 stars_old[1].weight = 100; 165 stars_old[1].weight = 5;
167 stars_old[1].v[0] = -0.5; 166 stars_old[1].v[0] = 0.1;
168 stars_old[1].v[1] = 0.1; 167 stars_old[1].v[1] = 5.0;
169 stars_old[1].v[2] = 1.732; 168 stars_old[1].v[2] = 0.0;
170 stars_old[1].r[0] = -50.0; 169 stars_old[1].r[0] = 10.0;
171 stars_old[1].r[1] = 0.0; 170 stars_old[1].r[1] = 0.0;
172 stars_old[1].r[2] = 0.0; 171 stars_old[1].r[2] = 0.0;
173 /* */ 172 /* */
174 stars_old[2].weight = 100; 173 stars_old[2].weight = 5;
175 stars_old[2].v[0] = 0.5; 174 stars_old[2].v[0] = 0.0;
176 stars_old[2].v[1] = -0.1; 175 stars_old[2].v[1] = -5.0;
177 stars_old[2].v[2] = -1.732; 176 stars_old[2].v[2] = 0.1;
178 stars_old[2].r[0] = 0.0; 177 stars_old[2].r[0] = -10.0;
179 stars_old[2].r[1] = 0.0; 178 stars_old[2].r[1] = 0.0;
180 stars_old[2].r[2] = 86.60; 179 stars_old[2].r[2] = 0.0;
180 #elif 1
181 /* */
182 stars_old[0].weight = 1000;
183 stars_old[0].v[0] = 0.0;
184 stars_old[0].v[1] = 0.0;
185 stars_old[0].v[2] = 0.0;
186 stars_old[0].r[0] = 0.0;
187 stars_old[0].r[1] = 0.0;
188 stars_old[0].r[2] = 0.0;
189 /* */
190 stars_old[1].weight = 5;
191 stars_old[1].v[0] = 0.1;
192 stars_old[1].v[1] = 5.0;
193 stars_old[1].v[2] = 0.0;
194 stars_old[1].r[0] = 100.0;
195 stars_old[1].r[1] = 0.0;
196 stars_old[1].r[2] = 0.0;
197 /* */
198 stars_old[2].weight = 5;
199 stars_old[2].v[0] = 0.0;
200 stars_old[2].v[1] = -5.0;
201 stars_old[2].v[2] = 0.1;
202 stars_old[2].r[0] = -100.0;
203 stars_old[2].r[1] = 0.0;
204 stars_old[2].r[2] = 0.0;
181 #endif 205 #endif
182 206
183 for( i=0; i<num; i++){ 207 for( i=0; i<num; i++){
184 stars_new[i].weight = stars_old[i].weight; 208 stars_new[i].weight = stars_old[i].weight;
185 stars_new[i].rect.h = 5, stars_new[i].rect.w = 5; 209 stars_new[i].rect.h = 5, stars_new[i].rect.w = 5;