diff --git a/python/dune/perftool/sumfact/vectorization.py b/python/dune/perftool/sumfact/vectorization.py
index da1926783b994a77d65b579dab05eca7a389ad26..4dab904154d80072585453383e482725bfa5da97 100644
--- a/python/dune/perftool/sumfact/vectorization.py
+++ b/python/dune/perftool/sumfact/vectorization.py
@@ -189,6 +189,11 @@ def decide_vectorization_strategy():
             for qp in quad_points:
                 for strat in fixed_quad_vectorization_opportunity_generator(frozenset(stage1_sumfacts), width, qp):
                     if strategy == int(get_option("vectorization_list_index")):
+                        set_quadrature_points(qp)
+                        # Output the strategy and its cost into a separate file
+                        if get_global_context_value("form_type") == "jacobian_apply":
+                            with open("strategycosts.csv", "a") as f:
+                                f.write("{} {}\n".format(strategy, strategy_cost(strat)))
                         return qp, strat
                     strategy = strategy + 1