[WIP] working version of parallel RBF computation.

Requires some clean up work

Author: FvanSteen

Common/include/grid_movement/CRadialBasisFunctionInterpolation.hpp

@@ -40,15 +40,20 @@ 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.*/
@@ -108,7 +113,7 @@ class CRadialBasisFunctionInterpolation : public CVolumetricMovement {
   * \brief Selecting internal nodes for the volumetric deformation.
   * \param[in] geometry - Geometrical definition of the problem.
   */
-  void SetInternalNodes(CGeometry* geometry);
+  void SetInternalNodes(CConfig* config, CGeometry* geometry);
 
   /*!
   * \brief Solving of the Radial Basis Function interpolation system, yielding the interpolation coefficients
@@ -140,4 +145,8 @@ class CRadialBasisFunctionInterpolation : public CVolumetricMovement {
   inline static bool Equal(CRadialBasisFunctionNode* a, CRadialBasisFunctionNode* b){
     return a->GetIndex() == b->GetIndex();
   }
+
+  inline static bool Equal2(unsigned long a, unsigned long b){
+      return a == b;    
+  };
 };

Common/include/grid_movement/CRadialBasisFunctionNode.hpp

@@ -47,7 +47,7 @@ class CRadialBasisFunctionNode{
   /*!
   * \brief Constructor of the class.
   */
-  CRadialBasisFunctionNode(CGeometry* geometry, unsigned short marker_val, unsigned long vertex_val);
+  CRadialBasisFunctionNode(unsigned long idx_val, unsigned short marker_val, unsigned long vertex_val);
 
   /*!
   * \brief Returns global index.