float a = Vector3f::dot(d,d);
float b = 2 * (Vector3f::dot(o,d));
float c = Vector3f::dot(o,o) - pow(radius, 2);

float discr = pow(b,2) - 4 * a * c;