From 734473b3e78e0cc9d36c92587099f09a7c30d2c9 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?St=C3=A9phane=20Del=20Pino?= <stephane.delpino44@gmail.com>
Date: Tue, 9 Mar 2021 19:57:37 +0100
Subject: [PATCH] Add cell center automatically in VTK outputs

This is useful/necessary to define properly cuts
---
 src/output/VTKWriter.cpp | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/output/VTKWriter.cpp b/src/output/VTKWriter.cpp
index f84d13cb3..f1f87e9cd 100644
--- a/src/output/VTKWriter.cpp
+++ b/src/output/VTKWriter.cpp
@@ -2,6 +2,8 @@
 
 #include <mesh/Connectivity.hpp>
 #include <mesh/Mesh.hpp>
+#include <mesh/MeshData.hpp>
+#include <mesh/MeshDataManager.hpp>
 #include <utils/Messenger.hpp>
 #include <utils/RevisionInfo.hpp>
 
@@ -279,6 +281,7 @@ VTKWriter::_write(const std::shared_ptr<const MeshType>& mesh,
   // Adding basic mesh information
   output_named_item_value_set.add(NamedItemValue{"cell_number", mesh->connectivity().cellNumber()});
   output_named_item_value_set.add(NamedItemValue{"node_number", mesh->connectivity().nodeNumber()});
+  output_named_item_value_set.add(NamedItemValue{"cell_center", MeshDataManager::instance().getMeshData(*mesh).xj()});
 
   if (parallel::rank() == 0) {   // write PVTK file
     std::ofstream fout(_getFilenamePVTU());
-- 
GitLab