Parallel Cgns reader/writer

.gitignore

@@ -1 +1,2 @@
 build
+settings.json

CMakeLists.txt

@@ -22,14 +22,26 @@ if(USE_XSDK_DEFAULTS)
   xsdk_compiler_flags()
 endif()
 
+# require c++14
+option(ENABLE_CGNS "Enable the CGNS reader: requires c++14 extensions" OFF)
+message(STATUS "ENABLE_CGNS: ${ENABLE_CGNS}")
+
 # Set some default compiler flags that should always be used
 if(NOT USE_XSDK_DEFAULTS)
   bob_set_shared_libs()
   bob_begin_cxx_flags()
   bob_end_cxx_flags()
   set(CMAKE_C_FLAGS "${CMAKE_CXX_FLAGS}")
   if(SCOREC_ENABLE_CXX11)
-    bob_cxx11_flags()
+    if(ENABLE_CGNS)
+      bob_cxx14_flags()
+    else()
+      bob_cxx14_flags()
+    endif()
+  else()
+    if(ENABLE_CGNS)
+      bob_cxx14_flags()
+    endif()
   endif()
 endif()
 message(STATUS "CMAKE_CXX_FLAGS = ${CMAKE_CXX_FLAGS}")
@@ -91,6 +103,7 @@ option(ENABLE_SIMMETRIX "Build with Simmetrix support" OFF)
 message(STATUS "ENABLE_SIMMETRIX: ${ENABLE_SIMMETRIX}")
 option(ENABLE_OMEGA_H "Enable the Omega_h interface" OFF)
 message(STATUS "ENABLE_OMEGA_H: ${ENABLE_OMEGA_H}")
+
 if(ENABLE_SIMMETRIX)
   add_definitions(-DHAVE_SIMMETRIX)
 endif()
@@ -117,6 +130,15 @@ if(ENABLE_OMEGA_H)
   bob_public_dep(Omega_h)
 endif()
 
+if(ENABLE_CGNS)
+  find_package(CGNS)
+  include_directories(SYSTEM ${CGNS_INCLUDE_DIR})
+  #
+  find_package(HDF5) 
+  #
+  add_definitions(-DHAVE_CGNS)
+endif()
+
 # Include the SCOREC project packages
 add_subdirectory(lion)
 add_subdirectory(pcu)
@@ -142,6 +164,9 @@ add_subdirectory(omega_h)
 # this INTERFACE target bundles all the enabled libraries together
 add_library(core INTERFACE)
 target_link_libraries(core INTERFACE ${SCOREC_EXPORTED_TARGETS})
+if(ENABLE_CGNS)
+  target_link_libraries(core INTERFACE ${CGNS_LIBRARIES} ${HDF5_LIBRARIES} ${HDF5_HL_LIBRARIES})
+endif()
 scorec_export_library(core)
 
 if(BUILD_EXES)

apf/CMakeLists.txt

@@ -49,6 +49,10 @@ set(SOURCES
   apfMIS.cc
 )
 
+if(ENABLE_CGNS)
+  set(SOURCES ${SOURCES} apfCGNS.cc)
+endif(ENABLE_CGNS)
+
 # Package headers
 set(HEADERS
   apf.h
@@ -94,6 +98,12 @@ target_link_libraries(apf
      mth
    )
 
+
+if(ENABLE_CGNS)
+  message(STATUS ${CGNS_LIBRARIES})
+  target_link_libraries(apf PRIVATE ${CGNS_LIBRARIES} ${HDF5_LIBRARIES} ${HDF5_HL_LIBRARIES})
+endif(ENABLE_CGNS)
+
 scorec_export_library(apf)
 
 bob_end_subdir()

apf/apf.cc

@@ -69,7 +69,7 @@ Field* makeField(
     FieldShape* shape,
     FieldData* data)
 {
-  PCU_ALWAYS_ASSERT( ! m->findField(name));
+  PCU_ALWAYS_ASSERT_VERBOSE( ! m->findField(name), name);
   Field* f = 0;
   if (valueType == SCALAR)
     f = new ScalarField();

apf/apf.h

@@ -13,6 +13,7 @@
 #include "apfDynamicArray.h"
 
 #include <vector>
+#include <map>
 #include <limits>
 
 /** \file apf.h
@@ -34,6 +35,7 @@ class Field;
 class Element;
 class Mesh;
 class MeshEntity;
+class MeshTag;
 class VectorElement;
 /** \brief Mesh Elements represent the mesh coordinate vector field. */
 typedef VectorElement MeshElement;
@@ -629,6 +631,28 @@ for (t::iterator i = (w).begin(); \
 for (t::const_iterator i = (w).begin(); \
      (i) != (w).end(); ++(i))
 
+struct CGNSInfo
+{
+  // cgns_bc_name
+  std::string cgnsBCSName;
+  /* tag value
+
+    Tag value holds [0, 1] as a
+    marker to indicate mesh_entities 
+    within bc group. 1="in group", 0="not in group"
+    Tags set on vertices, edges, faces, and cells
+  */
+  apf::MeshTag* bcsMarkerTag = nullptr;
+  // model dimension
+  int mdlDim = -1;
+  // model id
+  int mdlId = -1;
+};
+
+//using CGNSBCMap = std::map<std::string, std::vector<std::tuple<std::string, apf::MeshTag *, int>>>;
+using CGNSBCMap = std::map<std::string, std::vector<CGNSInfo>>;
+void writeCGNS(const char *prefix, Mesh *m, const CGNSBCMap &cgnsBCMap);
+
 /** \brief Write a set of parallel VTK Unstructured Mesh files from an apf::Mesh
   * with binary (base64) encoding and zlib compression (if LION_COMPRESS=ON)
   * \details Nodal fields whose shape differs from the mesh shape will