Merge pull request ITHACA-FV#646 from vressegu/mergeForkRedlum

Merge fork from redlum developpers

Author: giovastabile

docs/images/SmagorinskyCylinder.png

@@ -121,4 +121,49 @@ vectorTensorProduct(
     const Eigen::Matrix<float, Eigen::Dynamic, 1>& g,
     const Eigen::Tensor<float, 3 >& c,
     const Eigen::Matrix<float, Eigen::Dynamic, 1>& a);
+
+Eigen::VectorXd repeatElements(Eigen::VectorXd input, int n)
+{
+    if (!(n > 0))
+    {
+        Info << "The integer for the EigenFunctions::repeat method must be positive. Aborting" << endl;
+        abort();
+    }
+    else
+    {
+        Eigen::VectorXd output(input.size() * n);
+        for (int i = 0; i < input.size(); i++)
+        {
+            for (int r = 0; r < n; r++)
+            {
+                output(i * n + r) = input(i);
+            }
+        }
+        return output;
+    }
+}
+
+Eigen::VectorXd reorderVectorFromDim(Eigen::VectorXd input, int dim)
+{
+    if ((input.size() % dim) != 0)
+    {
+        Info << "EigenFunctions::reorderVectorFromDim has incompatible vector size " <<
+                "and reordering dimension. Aborting" << endl;
+        abort();
+    }
+    else
+    {
+        Eigen::VectorXd output(input.size());
+        int n = int(input.size() / dim);
+        for (int i = 0; i < n; i++)
+        {
+            for (int d = 0; d < dim; d++)
+            {
+                output(i * dim + d) = input(i + d * dim);
+            }
+        }
+        return output;
+    }
+}
+
 }

src/ITHACA_CORE/EigenFunctions/EigenFunctions.C

@@ -249,6 +249,29 @@ Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> vectorTensorProduct(
     const Eigen::Tensor<T, 3 >& c,
     const Eigen::Matrix<T, Eigen::Dynamic, 1>& a);
 
+
+//------------------------------------------------------------------------------
+/// @brief      Repeat each element n times after themselves
+///
+/// @details    Example: input = (a, b); n = 3
+///                      output = (a, a, a, b, b, b)
+///
+/// @param      input Eigen::VectorXd
+/// @param      n int 
+Eigen::VectorXd repeatElements(Eigen::VectorXd input, int n);
+
+//------------------------------------------------------------------------------
+/// @brief      Changes the order of the vector's elements given a dimension
+///
+/// @details    Example: input = (a_1, b_1, a_2, b_2, a_3, b_3); dim = 3
+///                      output = (a_1, a_2, a_3, b_1, b_2, b_3)
+///             Typically switch from a component then coordinate order to the opposite
+///             as did the change in the Foam2Eigen mapping 
+///
+/// @param      input Eigen::VectorXd
+/// @param      dim int 
+Eigen::VectorXd reorderVectorFromDim(Eigen::VectorXd input, int dim);
+
 };
 
 template <typename T>