Add comments

Author: mikeingold

src/specializations/BezierCurve.jl

@@ -50,6 +50,7 @@ end
 #                              Parametric
 ################################################################################
 
+# Wrap (::BezierCurve)(t, ::BezierEvalMethod) into f(t) by embedding second argument
 function _parametric(curve::Meshes.BezierCurve, t, alg::Meshes.BezierEvalMethod)
     return curve(t, alg)
 end

src/specializations/Line.jl

@@ -13,16 +13,18 @@ function integral(
         rule::IntegrationRule;
         kwargs...
 )
+    # Generate a _ParametricGeometry whose parametric function spans the domain [0, 1]
     paramfunction(t) = _parametric(line, t)
-    param_line = _ParametricGeometry(paramfunction, 1)
+
+    # Integrate the _ParametricGeometry using the standard methods
     return _integral(f, param_line, rule; kwargs...)
 end
 
 ################################################################################
 #                              Parametric
 ################################################################################
 
-# Map [0, 1] ↦ (-∞, ∞)
+# Map argument domain from [0, 1] to (-∞, ∞) for (::Line)(t)
 function _parametric(line::Meshes.Line, t)
     # [-1, 1] ↦ (-∞, ∞)
     f1(t) = t / (1 - t^2)

src/specializations/Plane.jl

@@ -14,19 +14,24 @@ function integral(
         rule::IntegrationRule;
         kwargs...
 )
+    # Generate a _ParametricGeometry whose parametric function spans the domain [0, 1]²
     paramfunction(t1, t2) = _parametric(plane, t1, t2)
-    param_plane = _ParametricGeometry(paramfunction, 2)
+
+    # Integrate the _ParametricGeometry using the standard methods
     return _integral(f, param_plane, rule; kwargs...)
 end
 
 ############################################################################################
 #                                       Parametric
 ############################################################################################
 
-# Map [0, 1]² ↦ (-∞, ∞)²
+# Map argument domain from [0, 1]² to (-∞, ∞)² for (::Plane)(t1, t2)
 function _parametric(plane::Meshes.Plane, t1, t2)
+    # [-1, 1] ↦ (-∞, ∞)
     f1(t) = t / (1 - t^2)
+    # [0, 1] ↦ [-1, 1]
     f2(t) = 2t - 1
+    # Compose the two transforms
     f = f1 ∘ f2
     return plane(f(t1), f(t2))
 end

src/specializations/Ray.jl

@@ -14,16 +14,18 @@ function integral(
         rule::IntegrationRule;
         kwargs...
 )
+    # Generate a _ParametricGeometry whose parametric function spans the domain [0, 1]
     paramfunction(t) = _parametric(ray, t)
-    param_ray = _ParametricGeometry(paramfunction, 1)
+
+    # Integrate the _ParametricGeometry using the standard methods
     return _integral(f, param_ray, rule; kwargs...)
 end
 
 ############################################################################################
 #                                     Parametric
 ############################################################################################
 
-# Map [0, 1] ↦ [0, ∞)
+# Map argument domain from [0, 1] to [0, ∞) for (::Ray)(t)
 function _parametric(ray::Meshes.Ray, t)
     f(t) = t / (1 - t^2)
     return ray(f(t))