diff --git a/ProcessLib/LocalAssemblerInterface.h b/ProcessLib/LocalAssemblerInterface.h
index 355bb5d8f1cdb0aecf2a2f15b05e377a7b06b3a7..f1f887c4f497e90c59480ea01fb2deec13fcd5e6 100644
--- a/ProcessLib/LocalAssemblerInterface.h
+++ b/ProcessLib/LocalAssemblerInterface.h
@@ -11,6 +11,7 @@
 #define PROCESSLIB_LOCALASSEMBLERINTERFACE_H
 
 #include "NumLib/NumericsConfig.h"
+#include "MathLib/Point3d.h"
 
 namespace NumLib
 {
@@ -53,6 +54,15 @@ public:
                               NumLib::LocalToGlobalIndexMap const& dof_table,
                               GlobalVector const& x);
 
+    /// Computes the flux in the point \c p_local_coords that is given in local
+    /// coordinates using the values from \c local_x.
+    virtual std::vector<double> getFlux(
+        MathLib::Point3d const& p_local_coords,
+        std::vector<double> const& local_x) const
+    {
+        return std::vector<double>();
+    }
+
 private:
     virtual void preTimestepConcrete(std::vector<double> const& /*local_x*/,
                                      double const /*t*/, double const /*dt*/)