@@ -81,0 +95,4 @@
+ * \{ */
+
+template<typename T>
+static Span<T> gaussian_blur_1D_ex(const bool is_cyclic,

@@ -81,0 +143,4 @@
+   * Note these weights only work because the averaging is done in relative coordinates.
+   */
+
+  /* Avoid computation if the mask is empty */

@@ -81,0 +159,4 @@
+
+  const int64_t first_pt = curve_points.first();
+  const int64_t last_pt = curve_points.last();
+  const int64_t nb_pts = curve_points.size();

@@ -81,0 +172,4 @@
+    }
+  });
+
+  for (int step = iterations; step > 0; step--) {

@@ -81,0 +180,4 @@
+        return;
+      }
+
+      float w_before = float(w - w2);

@@ -81,0 +263,4 @@
+static int grease_pencil_stroke_smooth_exec(bContext *C, wmOperator *op)
+{
+  using namespace blender;
+  Scene *scene = CTX_data_scene(C);

@@ -81,0 +291,4 @@
+        /* Position. */
+        Array<float3> curves_positions_copy;
+        if (smooth_position && !curves.positions().is_empty()) {
+          curves_positions_copy = curves.positions();

@@ -81,0 +293,4 @@
+        if (smooth_position && !curves.positions().is_empty()) {
+          curves_positions_copy = curves.positions();
+        }
+        /* Opacity. */

@@ -81,0 +294,4 @@
+          curves_positions_copy = curves.positions();
+        }
+        /* Opacity. */
+        Array<float> curves_opacities_copy;

@@ -81,0 +304,4 @@
+          curves_radii_copy = drawing.radii().get_internal_span();
+        }
+        /* Cyclic. */
+        const offset_indices::OffsetIndices<int> points_by_curve = curves.points_by_curve();

@@ -81,0 +307,4 @@
+        const offset_indices::OffsetIndices<int> points_by_curve = curves.points_by_curve();
+        const VArray<bool> cyclic = curves.cyclic();
+        /* Selection. */
+        bke::AttributeAccessor curves_attributes = curves.attributes();

@@ -81,0 +337,4 @@
+                                      smooth_ends,
+                                      keep_shape,
+                                      is_cyclic,
+                                      curves.positions_for_write(),

@@ -62,5 +63,25 @@ void create_blank(Main &bmain, Object &object, int frame_number);
 void create_stroke(Main &bmain, Object &object, float4x4 matrix, int frame_number);
 void create_suzanne(Main &bmain, Object &object, float4x4 matrix, const int frame_number);
 
+Span<float3> gaussian_blur_1D_float3(const IndexRange curve_points,