diff --git a/src/image.jl b/src/image.jl
index cfb3c9fb2c9afadcc268b31dcbefced541f242ea..a5449d2e2969b82537b4e3e36108941e2ace6a34 100644
--- a/src/image.jl
+++ b/src/image.jl
@@ -160,7 +160,7 @@ function project_l2_lagrange!(u::FeFunction, img::AbstractArray)
     opparams = (; f)
 
     nrdims = prod(space.size)
-    # number of element dofs (i.e. local dofs not counting range dimensions)
+    # number of element dofs (i.e. local dofs not counting codomain dimensions)
     nldofs = ndofs(space.element)
 
     a(xloc, u, du, v, dv; f) = dot(u, v)
@@ -348,7 +348,7 @@ function project_l2_pixel!(u::FeFunction, img)
     opparams = (; f)
 
     nrdims = prod(space.size)
-    # number of element dofs (i.e. local dofs not counting range dimensions)
+    # number of element dofs (i.e. local dofs not counting codomain dimensions)
     nldofs = ndofs(space.element)
 
     # first loop to count number of triangles intersecting pixel evaluation
diff --git a/src/operator.jl b/src/operator.jl
index c8adc04c3e086f49c84c9837469c2defc6a709ae..dfb77946f025585f06fb56059a2a4fcc493dd744 100644
--- a/src/operator.jl
+++ b/src/operator.jl
@@ -52,7 +52,7 @@ function assemble(space::FeSpace, a, l; params...)
 
     d = size(qx, 1) # domain dimension
     nrdims = prod(space.size)
-    nldofs = ndofs(space.element) # number of element dofs (i.e. local dofs not counting range dimensions)
+    nldofs = ndofs(space.element) # number of element dofs (i.e. local dofs not counting codomain dimensions)
     nqpts = length(qw) # number of quadrature points
 
     qphi_ = zeros(nrdims, nrdims, nldofs, nqpts)