#include <math.h>
#include "timer.h"

#define NUMBALLS 1000

struct Ball {float x, y, z, Size;} Balls [NUMBALLS];

void Setup (void)
{
    int i;

    for (i=0; i<NUMBALLS; i++) {
        Balls [i].Size = i * 10.0f/NUMBALLS;
        Balls [i].x = (i % 31) * 4.0f;
        Balls [i].y = (i / 31) * 2.0f;
        Balls [i].z = (i % 5) * 3.0f;
    }
}

int DoTest (void)
{
    int i, j, NumCollisions;
    float dx, dy, dz, d;

    NumCollisions = 0;
    for (i=0; i<NUMBALLS-1; i++) {
        for (j=i+1; j<NUMBALLS; j++) {
            dx = Balls [i].x - Balls [j].x;
            dy = Balls [i].y - Balls [j].y;
            dz = Balls [i].z - Balls [j].z;
            d = Balls [i].Size + Balls [j].Size;
            if (dx * dx + dy * dy + dz * dz < d * d) {
                NumCollisions++;
            }
        }
    }
    return NumCollisions;
}

int main () {
	unsigned a;

        Setup ( );

	START_TIMER();

        for ( a = 0; a < 200; a++ )
                DoTest ( );

	STOP_TIMER ( );	
	return 0;
}