Merge pull request #267 from ckoehn84/master

Added 3D output function
coolfluid · May 27, 2015 · 523e10b · 523e10b
2 parents fdbdaf6 + 402775f
commit 523e10b
Show file tree

Hide file tree

Showing 2 changed files with 41 additions and 5 deletions.
diff --git a/cf3/mesh/LagrangeP1/Tetra.cpp b/cf3/mesh/LagrangeP1/Tetra.cpp
@@ -54,17 +54,51 @@ const RealMatrix& Tetra::local_coordinates()
   static const RealMatrix loc_coord =
       (RealMatrix(nb_nodes, dimensionality) <<
 
-       0.,  0.,  0.,
-       1.,  0.,  0.,
-       0.,  1.,  0.,
-       0.,  0.,  1.
-
+	0., 0., 0.,
+        1., 0., 0.,
+        0., 1., 0.,
+        0., 0., 1.
+       
        ).finished();
   return loc_coord;
 }
 
 ////////////////////////////////////////////////////////////////////////////////
 
+const RealMatrix& Tetra::mononomial_coefficients()
+{
+  static const RealMatrix coeffs=
+      (RealMatrix(nb_nodes, nb_nodes) <<
+
+
+       1.,  -1.,  -1., -1.,
+       0.,   1.,   0.,  0.,
+       0.,   0.,   1.,  0.,
+       0.,   0.,   0.,  1.
+
+       ).finished();
+  return coeffs;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+const RealMatrix& Tetra::mononomial_exponents()
+{
+  static const RealMatrix exponents=
+      (RealMatrix(nb_nodes, dimensionality) <<
+
+       0, 0, 0,
+       1, 0, 0,
+       0, 1, 0,
+       0, 0, 1
+
+       ).finished();
+  return exponents;
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
+
 } // LagrangeP1
 } // mesh
 } // cf3
diff --git a/cf3/mesh/LagrangeP1/Tetra.hpp b/cf3/mesh/LagrangeP1/Tetra.hpp
@@ -45,6 +45,8 @@ struct Mesh_LagrangeP1_API Tetra : public ShapeFunctionBase<Tetra,Tetra_traits>
   static const RealMatrix& local_coordinates();
   static void compute_value(const MappedCoordsT& mapped_coord, ValueT& result);
   static void compute_gradient(const MappedCoordsT& mapped_coord, GradientT& result);
+  static const RealMatrix& mononomial_coefficients();
+  static const RealMatrix& mononomial_exponents();  
 };
 
 ////////////////////////////////////////////////////////////////////////////////