src/HRSC/smoother.jl

@@ -138,7 +138,7 @@ function smooth_C0_2nd_order!(bulk_v, cell_v, mesh::Mesh1d, isperiodic=true)
     v_l, v_r, v_c = cell_v[l], cell_v[r], cell_v[i]
     avg_1, avg_2, avg_3 = avg(v_l, v_c), v_c, avg(v_r, v_c)
     for j = 1:Nx
-      v[j] = lagrange_interpolation_2nd(z[j], avg_1, avg_2, avg_3)
+      v[j] = lagrange_interpolation_1d(z[j], avg_1, avg_2, avg_3)
     end
   end