Add basisType to RCSSurfaceGeometry and update related file handling

Author: AlejandroMunozManterola

src/components/RCSSurfaceExporter.cpp

@@ -54,11 +54,13 @@ RCSSurfaceExporter::RCSSurfaceExporter(
         totalQuadPts += ir->GetNPoints();
     }
 
-    int32_t header[4] = { 
+    int basisType = fec->GetBasisType();
+    int32_t header[5] = { 
         static_cast<int32_t>(spaceDim_), 
         static_cast<int32_t>(numDofs_), 
         static_cast<int32_t>(numBdr),
-        static_cast<int32_t>(totalQuadPts) 
+        static_cast<int32_t>(totalQuadPts),
+        static_cast<int32_t>(basisType)
     };
     dataFile_.write(reinterpret_cast<const char*>(header), sizeof(header));
 

src/components/RCSSurfacePostProcessor.cpp

@@ -79,12 +79,13 @@ RCSSurfacePostProcessor::readRankData(const std::string& rankPath) const
     std::ifstream f(rankPath + "/surface_data.bin", std::ios::binary);
     if (!f) throw std::runtime_error("Cannot open " + rankPath + "/surface_data.bin");
 
-    int32_t hdr[4];
+    int32_t hdr[5];
     f.read(reinterpret_cast<char*>(hdr), sizeof(hdr));
     rd.geometry.spaceDimension = hdr[0];
     rd.geometry.numDofs        = hdr[1];
     rd.geometry.numBdrElements = hdr[2];
     rd.geometry.numQuadPoints  = hdr[3];
+    rd.geometry.basisType      = hdr[4];
 
     const int nqp = rd.geometry.numQuadPoints;
     const int sdim = rd.geometry.spaceDimension;
@@ -309,7 +310,7 @@ void RCSSurfacePostProcessor::computeAndWriteResults(
         auto mesh = Mesh::LoadFromFile(rp + "/mesh", 1, 0);
         auto pmesh = ParMesh(MPI_COMM_WORLD, mesh);
         int order = determineFECOrder(pmesh, nDofs);
-        DG_FECollection fec(order, pmesh.Dimension());
+        DG_FECollection fec(order, pmesh.Dimension(), rd.geometry.basisType);
         ParFiniteElementSpace fes(&pmesh, &fec);
         std::cout << " done. (FEC order: " << order << ")\n";
 

src/components/RCSSurfacePostProcessor.h

@@ -13,6 +13,7 @@ struct RCSSurfaceGeometry {
     int numDofs;
     int numBdrElements;
     int numQuadPoints;
+    int basisType;
 
     std::vector<double> positions;  // [numQuadPoints * spaceDimension]
     std::vector<double> normals;    // [numQuadPoints * 3]