range -> codomain

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

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)