removed debug printouts and updated changelog

Author: Da-Be-Ru

CHANGELOG.md

@@ -1,10 +1,11 @@
 # CHANGES
 
-## v1.5.0 November 03, 2025
+## v1.5.0 November 21, 2025
   - added type parameter to `PointEvaluator` so `lazy_interpolate!` is now accessible for automatic differentiation
   - added new `compute_lazy_interpolation_jacobian` function for computing Jacobians of interpolations between finite element spaces
   - added `H1Pk_to_HDIVRT0_interpolator` function for (faster) computation of interpolation matrix from `H1Pk`-conforming spaces into `HDIVRT0`-conforming spaces on the same grid
   - fixed interpolation of `HCURLN1` finite element
+  - updated `MomentInterpolator` to deliver incident cells to `QPInfos` to speed up `lazy_interpolate!` 
 
 ## v1.4.0 July 17, 2025
   - added new weighted reconstruction operator `WeightedReconstruct` (so far only for H1BR{2})

src/interpolators.jl

@@ -112,7 +112,6 @@ function NodalInterpolator(
                     continue
                 end
                 cell = xNodeCells[1, j]
-                #@show cell
                 QP.item = cell
                 QP.cell = cell
                 QP.region = xCellRegions[cell]
@@ -187,8 +186,6 @@ function MomentInterpolator(
         kwargs...
     ) where {Tv, Ti, FEType <: AbstractFiniteElement, APT}
 
-    #@info "MomentInterpolator $(AT)"
-
     itemvolumes = xgrid[GridComponentVolumes4AssemblyType(AT)]
     itemnodes = xgrid[GridComponentNodes4AssemblyType(AT)]
     itemregions = xgrid[GridComponentRegions4AssemblyType(AT)]
@@ -418,7 +415,6 @@ function MomentInterpolator(
                 eval_trafo!(QP.x, L2G, xref[qp])
 
                 QP.cell = itemcells[item]
-                #@show QP.item QP.cell
 
                 exact_function!(result_f, QP)
                 if (bestapprox)

src/lazy_interpolate.jl

@@ -89,10 +89,8 @@ function lazy_interpolate!(
 
     if same_cells || use_cellparents
         xCellParents::Array{Ti, 1} = same_cells ? (1:num_cells(target.FES.xgrid)) : target.FES.xgrid[CellParents]
-        #@show xCellParents
         function point_evaluation_parentgrid!(result, qpinfo)
             x = xtrafo !== nothing ? xtrafo(x_source, qpinfo.x) : qpinfo.x
-            #@show qpinfo.x qpinfo.cell
             cell = gFindLocal!(xref, CF, x; icellstart = xCellParents[qpinfo.cell], eps = eps, trybrute = !only_localsearch)
             return evaluate_bary!(result, PE, xref, cell)
         end