MPI Point-to-Point Refactoring + New Periodic BC Implementation

Common/include/config_structure.hpp

@@ -1076,8 +1076,8 @@ class CConfig {
   bool Compute_Entropy;                      /*!< \brief Whether or not to compute the entropy in the fluid model. */
   bool Use_Lumped_MassMatrix_DGFEM;          /*!< \brief Whether or not to use the lumped mass matrix for DGFEM. */
   bool Jacobian_Spatial_Discretization_Only; /*!< \brief Flag to know if only the exact Jacobian of the spatial discretization must be computed. */
-  bool Compute_Average; /*!< \brief Whether or not to compute averages for unsteady simulations in FV or DG solver. */
-
+  bool Compute_Average;                      /*!< \brief Whether or not to compute averages for unsteady simulations in FV or DG solver. */
+  unsigned short Comm_Level;                 /*!< \brief Level of MPI communications to be performed. */
 
   ofstream *ConvHistFile;       /*!< \brief Store the pointer to each history file */
   bool Time_Domain;             /*!< \brief Determines if the multizone problem is solved in time-domain */
@@ -1399,7 +1399,7 @@ class CConfig {
   /*!
    * \brief Constructor of the class which reads the input file.
    */
-  CConfig(char case_filename[MAX_STRING_SIZE], unsigned short val_software, unsigned short val_iZone, unsigned short val_nZone, unsigned short val_nDim, unsigned short verb_level);
+  CConfig(char case_filename[MAX_STRING_SIZE], unsigned short val_software, unsigned short val_iZone, unsigned short val_nZone, unsigned short val_nDim, bool verb_high);
 
   /*!
    * \brief Constructor of the class which reads the input file.
@@ -1444,15 +1444,6 @@ class CConfig {
    * \return Total number of domains in the grid file.
    */
   static unsigned short GetnDim(string val_mesh_filename, unsigned short val_format);
-
-  /*!
-   * \brief Determine whether there are periodic BCs in the grid.
-   * \param[in] val_mesh_filename - Name of the file with the grid information.
-   * \param[in] val_format - Format of the file with the grid information.
-   * \param[in] config - Definition of the particular problem.
-   * \return Boolean for whether or not there are periodic BCs in the grid.
-   */
-  static bool GetPeriodic(string val_mesh_filename, unsigned short val_format, CConfig *config);
 
   /*!
    * \brief Initializes pointers to null
@@ -5985,7 +5976,7 @@ class CConfig {
    * \return <code>TRUE</code> if there is a rotational frame; otherwise <code>FALSE</code>.
    */
   bool GetAxisymmetric(void);
-
+    
   /*!
    * \brief Get information about the axisymmetric frame.
    * \return <code>TRUE</code> if there is a rotational frame; otherwise <code>FALSE</code>.
@@ -6560,7 +6551,7 @@ class CConfig {
    * \param[in] translate - Pointer to a vector containing the coordinate of the center.
    */
   void SetPeriodicTranslate(unsigned short val_index, su2double* translate);
-  
+
   /*!
    * \brief Get the translation vector for a periodic transformation.
    * \param[in] val_index - Index corresponding to the periodic transformation.
@@ -9140,6 +9131,12 @@ class CConfig {
   su2double GetMax_Time(void);
 
   /*!
+   * \brief Get the level of MPI communications to be performed.
+   * \return Level of MPI communications.
+   */
+  unsigned short GetComm_Level(void);
+
+  /*
    * \brief Check if the mesh read supports multiple zones.
    * \return YES if multiple zones can be contained in the mesh file.
    */
@@ -9168,6 +9165,7 @@ class CConfig {
    * \return YES if the forces breakdown file is written.
    */
   bool GetWrt_ForcesBreakdown(void);
+
 };
 
 #include "config_structure.inl"

Common/include/config_structure.inl

@@ -2014,6 +2014,8 @@ inline ofstream* CConfig::GetHistFile(void) { return ConvHistFile; }
 
 inline void CConfig::SetHistFile(ofstream *HistFile) { ConvHistFile = HistFile; }
 
+inline unsigned short CConfig::GetComm_Level(void) { return Comm_Level; }
+
 inline bool CConfig::GetTopology_Optimization(void) const { return topology_optimization; }
 
 inline string CConfig::GetTopology_Optim_FileName(void) const { return top_optim_output_file; }

Common/include/dual_grid_structure.hpp

@@ -147,10 +147,12 @@ class CPoint : public CDualGrid {
   vector<unsigned long> Point;        /*!< \brief Points surrounding the central node of the control volume. */
   vector<long> Edge;                  /*!< \brief Edges that set up a control volume. */
   su2double *Volume;                  /*!< \brief Volume or Area of the control volume in 3D and 2D. */
+  su2double Periodic_Volume;          /*!< \brief Missing component of volume or area of a control volume on a periodic marker in 3D and 2D. */
   bool Domain,                        /*!< \brief Indicates if a point must be computed or belong to another boundary */
   Boundary,                           /*!< \brief To see if a point belong to the boundary (including MPI). */
   PhysicalBoundary,                   /*!< \brief To see if a point belong to the physical boundary (without includin MPI). */
-  SolidBoundary;                      /*!< \brief To see if a point belong to the physical boundary (without includin MPI). */
+  SolidBoundary,                      /*!< \brief To see if a point belong to the physical boundary (without includin MPI). */
+  PeriodicBoundary;                      /*!< \brief To see if a point belongs to a periodic boundary (without including MPI). */
   long *Vertex;                       /*!< \brief Index of the vertex that correspond which the control volume (we need one for each marker in the same node). */
   su2double *Coord,                   /*!< \brief vector with the coordinates of the node. */
   *Coord_Old,                         /*!< \brief Old coordinates vector for geometry smoothing. */
@@ -403,6 +405,18 @@ class CPoint : public CDualGrid {
 	 */
 	su2double GetVolume(void);
 
+  /*!
+   * \brief Get the missing component of area or volume for a control volume on a periodic marker.
+   * \return Periodic component of area or volume for a control volume on a periodic marker.
+   */
+  su2double GetPeriodicVolume(void);
+
+  /*!
+   * \brief Set the missing component of area or volume for a control volume on a periodic marker.
+   * \param[in] val_volume - Value of the volume from the missing components of the CV on the periodic marker.
+   */
+  void SetPeriodicVolume(su2double val_volume);
+
 	/*!
 	 * \brief Get the maximum cell-center to cell-center length.
 	 * \return The maximum cell-center to cell-center length.
@@ -451,6 +465,12 @@ class CPoint : public CDualGrid {
 	 */
 	void SetSolidBoundary(bool val_boundary);
 
+  /*!
+   * \brief Set if a point belongs to a periodic boundary.
+   * \param[in] val_boundary - <code>TRUE</code> if the point belongs to a periodic boundary; otherwise <code>FALSE</code>.
+   */
+  void SetPeriodicBoundary(bool val_boundary);
+
   /*!
 	 * \brief Provides information about if a point belong to the physical boundaries (without MPI).
 	 * \return <code>TRUE</code> if the point belong to the boundary; otherwise <code>FALSE</code>.
@@ -463,6 +483,12 @@ class CPoint : public CDualGrid {
 	 */
 	bool GetSolidBoundary(void);
 
+  /*!
+   * \brief Provides information about if a point belongs to a periodic boundary (without MPI).
+   * \return <code>TRUE</code> if the point belongs to a periodic boundary; otherwise <code>FALSE</code>.
+   */
+  bool GetPeriodicBoundary(void);
+
 	/*! 
 	 * \brief Set a color to the point that comes from the grid partitioning.
 	 * \note Each domain has a different color.

Common/include/dual_grid_structure.inl

@@ -80,6 +80,10 @@ inline unsigned long CPoint::GetPoint(unsigned short val_point) { return Point[v
 
 inline su2double CPoint::GetVolume (void) { return Volume[0]; }
 
+inline su2double CPoint::GetPeriodicVolume (void) { return Periodic_Volume; }
+
+inline void CPoint::SetPeriodicVolume (su2double val_volume) { Periodic_Volume = val_volume; }
+
 inline su2double CPoint::GetMaxLength(void) {return MaxLength;}
 
 inline bool CPoint::GetMove (void) { return Move; }
@@ -96,6 +100,10 @@ inline void CPoint::SetSolidBoundary(bool val_boundary) { SolidBoundary = val_bo
 
 inline bool CPoint::GetSolidBoundary(void) { return SolidBoundary; }
 
+inline void CPoint::SetPeriodicBoundary(bool val_boundary) { PeriodicBoundary = val_boundary; }
+
+inline bool CPoint::GetPeriodicBoundary(void) { return PeriodicBoundary; }
+
 inline void CPoint::AddVolume (su2double val_Volume) { Volume[0] += val_Volume; }
 
 inline void CPoint::SetVolume (su2double val_Volume) { Volume[0] = val_Volume; }