Function sample_unit_disk_concentric 

pub fn sample_unit_disk_concentric(u1: Real, u2: Real) -> (Real, Real)

Uniformly sample a point on the surface of a unit disk.

Uses Shirley’s concentric-disk mapping for low distortion. Returns (x, y) with x² + y² ≤ 1.