Improvement parallel computations, start of including data reduction

Author: FvanSteen

Common/include/grid_movement/CRadialBasisFunctionInterpolation.hpp

@@ -40,23 +40,40 @@ class CRadialBasisFunctionInterpolation : public CVolumetricMovement {
 protected:
 
   vector<CRadialBasisFunctionNode*> boundaryNodes;  /*!< \brief Vector with boundary nodes.*/
-  vector<CRadialBasisFunctionNode*> boundaryNodes_local;  /*!< \brief Vector with boundary nodes.*/
-  vector<unsigned long> boundNodes, boundNodes_local;
   vector<unsigned long> internalNodes;              /*!< \brief Vector with internal nodes.*/
   unsigned long nBoundaryNodes,                     /*!< \brief Number of boundary nodes*/
   nInternalNodes;                                   /*!< \brief Number of internal nodes*/
 
   vector<CRadialBasisFunctionNode*>* controlNodes;  /*!< \brief Vector with control nodes*/
 
   vector<passivedouble> deformationVector;  /*!< \brief Deformation vector.*/
-  vector<passivedouble> deformationVector_local;  /*!< \brief Deformation vector.*/
 
   vector<passivedouble> coefficients;       /*!< \brief Control node interpolation coefficients.*/
   CSymmetricMatrix interpMat;               /*!< \brief Interpolation matrix.*/  
   // su2activematrix interpMat;
 
   RADIAL_BASIS kindRBF; /*!< \brief Type of Radial Basis Function.*/
   su2double radius;     /*!< \brief Support radius of compact Radial Basis Function.*/
+
+
+  /*--- data reduction parameters ---*/
+  vector<CRadialBasisFunctionNode*> greedyNodes;  /*!< \brief Vector with selected control nodes in greedy algorithm. */
+  bool dataReduction;                             /*!< \brief Determines whether data reduction is used. */
+  unsigned long MaxErrorNode;
+  su2double MaxError;
+
+
+  unsigned long Global_nControlNodes{0};
+  /*--- mpi related*/
+  
+  unsigned long Local_nControlNodes;
+  vector<unsigned long> Local_nControlNodesVec;
+
+  vector<su2double> LocalCoords;
+  vector<su2double> GlobalCoords;
+
+
+
 
 public:
 
@@ -149,4 +166,10 @@ class CRadialBasisFunctionInterpolation : public CVolumetricMovement {
   inline static bool Equal2(unsigned long a, unsigned long b){
       return a == b;    
   };
+
+
+  /*--- Data reduction functions ---*/
+  void GreedyIteration(CGeometry* geometry, CConfig* config);
+
+  void GetInitMaxErrorNode(CGeometry* geometry, CConfig* config);
 };