Example #1
0
int main(int argc, char ** argv)
{
  offset_momentum();
  printf ("%.9f\n", energy());
  advance(atoi(argv[1]));
  printf ("%.9f\n", energy());
  return 0;
}
Example #2
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;
}
Example #3
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;
}
Example #4
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;
}
Example #5
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;
}
Example #6
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;
}
Example #7
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);
  }
}
Example #8
0
File: nbody.c Project: arowM/felix
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;
}