[MOAB-dev] r5712 - MOAB/trunk/src/moab/point_locater

[rhl6856 at mcs.anl.gov]

Author: rhl6856
Date: 2012-08-29 13:17:26 -0500 (Wed, 29 Aug 2012)
New Revision: 5712

Added:
   MOAB/trunk/src/moab/point_locater/quadratic_hex_map.hpp
   MOAB/trunk/src/moab/point_locater/spectral_hex_map.hpp
Log:
spectral and quads

Added: MOAB/trunk/src/moab/point_locater/quadratic_hex_map.hpp
===================================================================
--- MOAB/trunk/src/moab/point_locater/quadratic_hex_map.hpp	                        (rev 0)
+++ MOAB/trunk/src/moab/point_locater/quadratic_hex_map.hpp	2012-08-29 18:17:26 UTC (rev 5712)
@@ -0,0 +1,203 @@
+#ifndef MOAB_QUADRATIC_HEX_HPP
+#define MOAB_QUADRATIC_HEX_HPP
+
+#include "moab/Matrix3.hpp"
+#include "moab/CartVect.hpp"
+#include <sstream>
+#include <iomanip>
+#include <iostream>
+
+namespace moab { 
+
+namespace element_utility {
+
+namespace {
+
+double SH(const int i, const double xi)
+{
+  switch (i)
+  {
+  case -1: return (xi*xi-xi)/2;
+  case 0: return 1-xi*xi;
+  case 1: return (xi*xi+xi)/2;
+  default: return 0.;
+  }
+}
+double DSH(const int i, const double xi)
+{
+  switch (i)
+  {
+  case -1: return xi-0.5;
+  case 0: return -2*xi;
+  case 1: return xi+0.5;
+  default: return 0.;
+  }
+}
+
+
+} //non-exported functionality
+
+template< typename _Matrix>
+class Quadratic_hex_map {
+  public:
+	typedef _Matrix Matrix;
+  private:
+	typedef Quadratic_hex_map< Matrix> Self;
+  public: