sdfCylindricalRepeat

fn sdfCylindricalRepeat(position: vec3<f32>, angleRepeat: f32, heightRepeat: f32, axis: vec3<f32>) -> vec3<f32> {
  let n = normalize(axis);
  
  // Project onto axis for height
  let h = dot(position, n);
  let radial = position - h * n;
  
  // Repeat in height
  let newH = h - heightRepeat * round(h / heightRepeat);
  
  // Repeat in angle
  let radius = length(radial);
  if (radius < 0.001) {
    return newH * n;
  }
  
  let angle = atan2(radial.y, radial.x);
  let sectorAngle = 6.28318530718 / angleRepeat;
  let newAngle = angle - sectorAngle * round(angle / sectorAngle);
  
  let newRadial = radius * vec2<f32>(cos(newAngle), sin(newAngle));
  
  // This assumes axis is along Z - for general axis, need proper basis vectors
  return newH * n + vec3<f32>(newRadial.x, newRadial.y, 0.0);
}