use nodal interpolation for inpaint and fine mesh

scripts/run_experiments.jl

@@ -1354,14 +1354,16 @@ function inpaint(ctx)
         if ctx.params.n_refine == 0
             # if we use a grid at image resolution, we wan't to avoid
             # information loss at all cost.
-            # projet_l2_pixel is not suited for inpainting as it is a global
+            # project_l2_pixel is not suited for inpainting as it is a global
             # operator, which may leak data from the inpainting domain
             interpolate!(st.g, x -> evaluate_bilinear(g_arr, x))
+            interpolate!(st.tdata, x -> evaluate_bilinear(mask_arr, x))
         else
             project_image!(st.g, g_arr)
-        end
+            # TODO: maybe a more conservative choice would be good
             project_image!(st.tdata, mask_arr)
         end
+    end
 
     save_step(i) =
         output(st, joinpath(ctx.outdir, "output_$(lpad(i, 5, '0')).vtu"),