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

[rhl6856 at mcs.anl.gov]

Author: rhl6856
Date: 2012-08-03 13:37:27 -0500 (Fri, 03 Aug 2012)
New Revision: 5671

Added:
   MOAB/trunk/src/moab/point_locater/element_utility.hpp
Removed:
   MOAB/trunk/src/moab/point_locater/ElemUtil.cpp
   MOAB/trunk/src/moab/point_locater/ElemUtil.hpp
   MOAB/trunk/src/moab/point_locater/ElemUtilTest.cpp
   MOAB/trunk/src/moab/point_locater/nat_param.hpp
Log:
Cleaning up nat param

Deleted: MOAB/trunk/src/moab/point_locater/ElemUtil.cpp
===================================================================
--- MOAB/trunk/src/moab/point_locater/ElemUtil.cpp	2012-08-03 18:37:26 UTC (rev 5670)
+++ MOAB/trunk/src/moab/point_locater/ElemUtil.cpp	2012-08-03 18:37:27 UTC (rev 5671)
@@ -1,789 +0,0 @@
-#include <iostream>
-#include <limits>
-#include <assert.h>
-
-#include "ElemUtil.hpp"
-#include "types.h"
-
-namespace moab { 
-namespace ElemUtil {
-
-  bool debug = false;
-
-/**\brief Class representing a 3-D mapping function (e.g. shape function for volume element) */
-class VolMap {
-  public:
-      /**\brief Return $\vec \xi$ corresponding to logical center of element */
-    virtual CartVect center_xi() const = 0;
-      /**\brief Evaluate mapping function (calculate $\vec x = F($\vec \xi)$ )*/
-    virtual CartVect evaluate( const CartVect& xi ) const = 0;
-      /**\brief Evaluate Jacobian of mapping function */
-    virtual Matrix3 jacobian( const CartVect& xi ) const = 0;
-      /**\brief Evaluate inverse of mapping function (calculate $\vec \xi = F^-1($\vec x)$ )*/
-    bool solve_inverse( const CartVect& x, CartVect& xi, double tol ) const ;
-};
-
-
-bool VolMap::solve_inverse( const CartVect& x, CartVect& xi, double tol ) const
-{
-  const double error_tol_sqr = tol*tol;
-  double det;
-  xi = center_xi();
-  CartVect delta = evaluate(xi) - x;
-  Matrix3 J;
-  while (delta % delta > error_tol_sqr) {
-    J = jacobian(xi);
-    det = J.determinant();
-    if (det < std::numeric_limits<double>::epsilon())
-      return false;
-    xi -= J.inverse(1.0/det) * delta;
-    delta = evaluate( xi ) - x;
-  }
-  return true;
-}
-
-/**\brief Shape function for trilinear hexahedron */
-class LinearHexMap : public VolMap {