diff --git a/AssemblerLib/SerialExecutor.h b/AssemblerLib/SerialExecutor.h
index 2bb5c1ac3c36a41175f25af5aa1b75cbbf2ac0cf..039fea7d2a7f5bcafb149ed8633f37d5cc37b64a 100644
--- a/AssemblerLib/SerialExecutor.h
+++ b/AssemblerLib/SerialExecutor.h
@@ -40,6 +40,32 @@ struct SerialExecutor
             f(i, c[i], std::forward<Args_>(args)...);
     }
 
+    /// Same as execute(f, c), but with two containers, where the second one is
+    /// modified.
+    ///
+    /// \tparam F    \c f type.
+    /// \tparam C    input container type.
+    /// \tparam Data input/output container type.
+    ///
+    /// \param f     a function that accepts a pointer to container's elements,
+    ///              an index, and a second container element as arguments, which
+    ///              is modified.
+    /// \param c     a container supporting const access over operator[] and size().
+    /// \param data  a container supporting non-const access over operator[] and size().
+    template <typename F, typename C, typename Data, typename ...Args_>
+    static
+    void
+#if defined(_MSC_VER) && (_MSC_VER >= 1700)
+    execute(F& f, C const& c, Data& data, Args_&&... args)
+#else
+    execute(F const& f, C const& c, Data& data, Args_&&... args)
+#endif
+    {
+        assert(c.size() == data.size());
+
+        for (std::size_t i = 0; i < c.size(); i++)
+            f(i, c[i], data[i], std::forward<Args_>(args)...);
+    }
 };
 
 }   // namespace AssemblerLib