Untitled

double t1, t2;
  t1 = omp_get_wtime();
  for (y = 0; y < height; ++y) {
    ncalls_line = 0;
    for (x = 0; x < width; ++x) {
      pixel = &image[y*width+x];
      xx = (2 * ((x + 0.5) * invWidth) - 1) * angle * aspectratio;
      yy = (1 - 2 * ((y + 0.5) * invHeight)) * angle;
      Vec3_new(&raydir, xx, yy, -1);
      Vec3_normalize(&raydir);
      ncalls = trace(pixel, &origin, &raydir, size, spheres, 0);
      ncalls_line += ncalls;
    }
    ncalls_line = ncalls_line/width;
    if (ncalls_line < min_ncalls_line)
      min_ncalls_line = ncalls_line;
    if (ncalls_line > max_ncalls_line)
      max_ncalls_line = ncalls_line;
    update_histogram(histo, ncalls_line);
  }
  t2 = omp_get_wtime();
  printf("Time: %f", t2-t1);