Raytracing неправильная мягкая выборка теней

Здравствуйте, я изучаю алгоритм трассировки лучей и застрял в алгоритме Монте-Карло. При рендеринге без подсветки области мой вывод был правильным, но когда я добавил реализацию подсветки области в исходный код для генерации мягкой тени, я столкнулся с проблемой.

Вот изображение до и после.

Когда я перемещаю синюю сферу вниз, проблема продолжается (обратите внимание, что артефакт продолжается, когда сфера вдоль белой пунктирной линии).
Обратите внимание, что эта сфера и arealight имеют одинаковое смещение по оси z. Когда я подношу синюю сферу к экрану, артефакт исчезнет. Я думаю, что проблема вызвана равномерным конусом выборки или функцией сферы выборки, но не уверен.

Вот функция:

template <typename T>
CVector3<T> UConeSample(T u1, T u2, T costhetamax,
const CVector3<T>& x, const CVector3<T>& y, const CVector3<T>& z) {
T costheta = Math::Lerp(u1, costhetamax, T(1));
T sintheta = sqrtf(T(1) - costheta*costheta);
T phi = u2 * T(2) * T(M_PI);

return cosf(phi) * sintheta * x +
sinf(phi) * sintheta * y +
costheta * z;
}

Я генерирую случайное значение с плавающей точкой u1, u2 из последовательности Ван Дер Корпута.
Это метод сферной выборки

CPoint3<float> CSphere::Sample(const CLightSample& ls, const CPoint3<float>& p, CVector3<float> *n) const {
// translate object to world space
CPoint3<float> pCentre = o2w(CPoint3<float>(0.0f));
CVector3<float> wc = Vector::Normalize(pCentre - p);
CVector3<float> wcx, wcy;
//create local coordinate system from wc for uniform sample cone
Vector::CoordinateSystem(wc, &wcx, &wcy);

//check if inside, epsilon val. this is true?
if (Point::DistSquare(p, pCentre) - radius*radius < 1e-4f)
return Sample(ls, n);

// Else outside evaluate cosinus theta value
float sinthetamax2 = radius * radius / Point::DistSquare(p, pCentre);
float costhetamax = sqrtf(Math::Max(0.0f, 1.0f - sinthetamax2));

// Surface properties
CSurfaceProps dg_sphere;
float thit, ray_epsilon;
CPoint3<float> ps;

//create ray direction from sampled point then send ray to sphere
CRay ray(p, Vector::UConeSample(ls.u1, ls.u2, costhetamax, wcx, wcy, wc), 1e-3f);
// Check intersection against sphere, fill surface properties and calculate hit point
if (!Intersect(ray, &thit, &ray_epsilon, &dg_sphere))
thit = Vector::Dot(pCentre - p, Vector::Normalize(ray.d));

// Evaluate surface normal
ps = ray(thit);
*n = CVector3<float>(Vector::Normalize(ps - pCentre));

//return sample point
return ps;
}

У кого-нибудь есть предложения? Благодарю.