sdfCappedTorus

fn sdfCappedTorus(position: vec3<f32>, majorRadius: f32, minorRadius: f32, angle: f32) -> f32 {
  let sc = vec2<f32>(sin(angle), cos(angle));
  let q = vec3<f32>(abs(position.x), position.y, position.z);
  let k = select(
    length(q.xy), 
    dot(q.xy, sc),
    sc.y * q.x > sc.x * q.y
  );
  return sqrt(dot(q, q) + 
              majorRadius * majorRadius - 
              2.0 * majorRadius * k) - 
              minorRadius;
}