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

[rhl6856 at mcs.anl.gov]

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

Modified:
   MOAB/trunk/src/moab/point_locater/point_locater.hpp
Log:
point_locater top level 0

Modified: MOAB/trunk/src/moab/point_locater/point_locater.hpp
===================================================================
--- MOAB/trunk/src/moab/point_locater/point_locater.hpp	2012-08-03 18:37:13 UTC (rev 5654)
+++ MOAB/trunk/src/moab/point_locater/point_locater.hpp	2012-08-03 18:37:13 UTC (rev 5655)
@@ -24,403 +24,58 @@
 #ifndef POINT_LOCATER_HPP
 #define POINT_LOCATER_HPP
 
-#include "iBase.h"
-#include "moab/Range.hpp"
-#include "moab/Interface.hpp"
-#include "moab/CartVect.hpp"
-
-#include "moab/ParallelComm.hpp"
-#include "moab/AdaptiveKDTree.hpp"
-#include "moab/GeomUtil.hpp"
-#include "ElemUtil.hpp"
-#include "moab/CN.hpp"
-#include "moab/Types.hpp"
-//#include "iMesh_extensions.h"
-
-#include <stdio.h>
-#include <iostream>
-
-#include "moab/gs.hpp"
-#include "moab/TupleList.hpp"
-
-#include "assert.h"
-
 namespace moab {
 
-namespace _point_search {
-template< typename Tuple_list, typename Vector> 
-void write_data( Tuple_list & target_points, 
-		 const Tuple_list & source_points, 
-		 Vector & located_points){
-	// count non-negatives
-	// TODO: replace with int num_pts = std::count()
-	int num_pts = 0;
-	for (unsigned int i = 0; i < source_points.get_n(); i++){
-	  if (-1 != source_points.vi_rd[3*i+2]) num_pts++;  
-	}
-	target_points.initialize(3, 0, 0, 1, num_pts);
-	target_points.enableWriteAccess();
-	for (unsigned int i = 0; i < source_points.get_n(); i++) {
-	  //Don't send message for discarded points
-	  if (source_points.vi_rd[3*i+2] != -1) { 
-	    const int locIndex = source_points.vi_rd[3*i+1];
-	    //asked 2+ procs if they have point p, they both said 
-	    //yes, we'll keep the one with lowest rank
-	    if(located_points[locIndex]){