float3 intersectSphere(float3 rayOrigin, float3 rayDir, float radius)
{
    float b = dot(rayDir, rayOrigin);
    float c = dot(rayOrigin, rayOrigin) - radius*radius;
    float t = -b - sqrt(b*b-c);
    return rayOrigin + t * rayDir;
}