std::array<bool, 4> IntersectBoxWithSphere(std::array<const Math::AABB*, 4> boxes, const Physics::Sphere& sphere) const {
                float radiusf = sphere.radius * sphere.radius;
                __m128 zero = _mm_setzero_ps();
                __m128 radius = _mm_load1_ps(&radiusf);
                __m128 centerx = _mm_load1_ps(&sphere.origin.x);
                __m128 centery = _mm_load1_ps(&sphere.origin.y);
                __m128 centerz = _mm_load1_ps(&sphere.origin.z);
                __m128 boxminx = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[0]->BottomLeftClosest));
                __m128 boxmaxx = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[0]->TopRightFurthest));
                __m128 boxminy = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[1]->BottomLeftClosest));
                __m128 boxmaxy = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[1]->TopRightFurthest));
                __m128 boxminz = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[2]->BottomLeftClosest));
                __m128 boxmaxz = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[2]->TopRightFurthest));
                __m128 e       = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[3]->BottomLeftClosest));
                _MM_TRANSPOSE4_PS(boxminx, boxminy, boxminz, e);
                e              = _mm_loadu_ps(reinterpret_cast<const float*>(&boxes[3]->TopRightFurthest));
                _MM_TRANSPOSE4_PS(boxmaxx, boxmaxy, boxmaxz, e);

                // _mm_max_ps(_mm_sub_ps(boxmin, center), zero)
                boxminx = _mm_max_ps(_mm_sub_ps(boxminx, centerx), zero);
                boxminy = _mm_max_ps(_mm_sub_ps(boxminy, centery), zero);
                boxminz = _mm_max_ps(_mm_sub_ps(boxminz, centerz), zero);

                // _mm_max_ps(_mm_sub_ps(center, boxmin), zero)
                boxmaxx = _mm_max_ps(_mm_sub_ps(centerx, boxmaxx), zero);
                boxmaxy = _mm_max_ps(_mm_sub_ps(centery, boxmaxy), zero);
                boxmaxz = _mm_max_ps(_mm_sub_ps(centerz, boxmaxz), zero);

                e = _mm_add_ps(_mm_add_ps(boxminx, boxmaxx), _mm_add_ps(_mm_add_ps(boxminy, boxmaxy), _mm_add_ps(boxminz, boxmaxz)));
                e = _mm_cmple_ps(_mm_mul_ps(e, e), radius);
                float output[4];
                _mm_store_ps(output, e);
                std::array<bool, 4> ret;
                for(int i = 0; i < 4; i++) {
                    ret[i] = output[i];
                }
                return ret;
            }