int main(int argc, char ** argv) { offset_momentum(); printf ("%.9f\n", energy()); advance(atoi(argv[1])); printf ("%.9f\n", energy()); return 0; }
int main(int argc, char ** argv) { long double N; scanf("%ld",&N); offset_momentum(); printf ("%.9f\n", energy()); advance(N); printf ("%.9f\n", energy()); return 0; }
int main() { int n = 5000000; int i; offset_momentum(NBODIES, bodies); printf ("%.9f\n", energy(NBODIES, bodies)); for (i = 1; i <= n; i++) advance(NBODIES, bodies, 0.01); printf ("%.9f\n", energy(NBODIES, bodies)); return 0; }
int main(int argc, char ** argv) { int n = atoi(argv[1]); int i; offset_momentum(NBODIES, bodies); printf ("%.9f\n", energy(NBODIES, bodies)); for (i = 1; i <= n; i++) advance(NBODIES, bodies, 0.01); printf ("%.9f\n", energy(NBODIES, bodies)); return 0; }
int main(int argc, char ** argv) { int n,i; // = atoi(argv[1]); //int i; scanf("%d",&n); offset_momentum(NBODIES, bodies); printf ("%.9f\n", energy(NBODIES, bodies)); v2df dtdt = {0.01, 0.01}; for (i = 1; i <= n; i++) advance(NBODIES, bodies, dtdt); printf ("%.9f\n", energy(NBODIES, bodies)); return 0; }
int main(int argc, char ** argv) { int n = atoi(argv[1]); int i; offset_momentum(NBODIES, bodies); printf ("%.9f\n", energy(NBODIES, bodies)); scale_bodies(NBODIES, bodies, DT); for (i = 1; i <= n; i++) { advance(NBODIES, bodies); } scale_bodies(NBODIES, bodies, RECIP_DT); printf ("%.9f\n", energy(NBODIES, bodies)); return 0; }
void inner_iter(int n) { int i; offset_momentum(NBODIES, bodies); checksum += energy(NBODIES, bodies); for (i = 1; i <= n; i++) n_advance(NBODIES, bodies, 0.01); checksum += energy(NBODIES, bodies); if (abs(checksum - EXPECT_CHECKSUM) >= EPSILON) { errx(EXIT_FAILURE, "bad checksum: %.52f vs %.52f", checksum, EXPECT_CHECKSUM); } }
int main(int argc, char ** argv) { int n = 50000000; int i; clock_t t0, t1; double r0, r1; t0 = clock(); offset_momentum(NBODIES, bodies); r0 = energy(NBODIES, bodies); for (i = 1; i <= n; i++) advance(NBODIES, bodies, 0.01); r1 = energy(NBODIES, bodies); t1 = clock(); printf("%.9f\n", r0); printf("%.9f\n", r1); printf("%f\n", ((double) (t1 - t0)) / CLOCKS_PER_SEC); return 0; }