diff --git a/MathLib/FormattingUtils.cpp b/MathLib/FormattingUtils.cpp
new file mode 100644
index 0000000000000000000000000000000000000000..7123ce449da96ed9298164d50f81f98f765fc16a
--- /dev/null
+++ b/MathLib/FormattingUtils.cpp
@@ -0,0 +1,17 @@
+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2024, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#include "FormattingUtils.h"
+
+namespace MathLib
+{
+Eigen::IOFormat const EigenIOFormat::full_precision{
+    Eigen::FullPrecision, 0, " ", "\n", "", "", "", "", ' '};
+}  // namespace MathLib
diff --git a/MathLib/FormattingUtils.h b/MathLib/FormattingUtils.h
index ba75fc4e6afb3d19c61e681ab75ce3c062e8f42d..8c346d64121d44b484af99024bbcb56e0123b27e 100644
--- a/MathLib/FormattingUtils.h
+++ b/MathLib/FormattingUtils.h
@@ -15,8 +15,21 @@
 #include <Eigen/Core>
 #include <concepts>
 
+namespace MathLib
+{
+struct EigenIOFormat
+{
+    static const Eigen::IOFormat full_precision;
+};
+}  // namespace MathLib
+
 template <typename T>
 requires std::derived_from<T, Eigen::DenseBase<T>>
 struct fmt::formatter<T> : fmt::ostream_formatter
 {
+    auto format(T const& value, fmt::format_context& ctx) const
+    {
+        return fmt::ostream_formatter::format(
+            value.format(MathLib::EigenIOFormat::full_precision), ctx);
+    }
 };