src/MEDMEM/PointLocator.cxx

   1 // Copyright (C) 2007-2013  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19
  20 #include "MEDNormalizedUnstructuredMesh.txx"
  21 #include "MEDMEM_Mesh.hxx"
  22 #include "MEDMEM_Exception.hxx"
  23 #include "PointLocatorAlgos.txx"
  24 #include "PointLocator.hxx"
  25
  26 #include <list>
  27
  28 MEDMEM::PointLocator::PointLocator(const MEDMEM::MESH& mesh)
  29 {
  30   int meshdim=mesh.getMeshDimension();
  31   int spacedim=mesh.getSpaceDimension();
  32   if (meshdim != spacedim) throw MEDMEM::MEDEXCEPTION("Locator is not implemented for meshdim != spacedim");
  33   switch (meshdim)
  34       {
  35       case 2:
  36         _medmesh = new MEDNormalizedUnstructuredMesh<2,2> (&mesh);
  37         _point_locator=new INTERP_KERNEL::PointLocatorAlgos<MEDNormalizedUnstructuredMesh<2,2> >(*(static_cast<MEDNormalizedUnstructuredMesh<2,2>* >(_medmesh)));
  38         break;
  39       case 3:
  40         _medmesh = new MEDNormalizedUnstructuredMesh<3,3> (&mesh);
  41         _point_locator=new INTERP_KERNEL::PointLocatorAlgos<MEDNormalizedUnstructuredMesh<3,3> >(*(static_cast<MEDNormalizedUnstructuredMesh<3,3>* >(_medmesh)));
  42         break;
  43       }
  44 }
  45
  46 MEDMEM::PointLocator::~PointLocator()
  47 {
  48   delete _medmesh;
  49   delete _point_locator;
  50 }
  51
  52 std::list<int> MEDMEM::PointLocator::locate(const double* x)
  53 {
  54   return _point_locator->locates(x,1e-12);
  55 }
  56
  57 MEDMEM::PointLocatorInSimplex::PointLocatorInSimplex(const MEDMEM::MESH& mesh)
  58 {
  59   int meshdim=mesh.getMeshDimension();
  60   int spacedim=mesh.getSpaceDimension();
  61   if (meshdim != spacedim) throw MEDMEM::MEDEXCEPTION("Locator is not implemented for meshdim != spacedim");
  62   switch (meshdim)
  63       {
  64       case 2:
  65         _medmesh = new MEDNormalizedUnstructuredMesh<2,2> (&mesh);
  66         _point_locator=new INTERP_KERNEL::PointLocatorInSimplex<MEDNormalizedUnstructuredMesh<2,2> >(*(static_cast<MEDNormalizedUnstructuredMesh<2,2>* >(_medmesh)));
  67         break;
  68       case 3:
  69         _medmesh = new MEDNormalizedUnstructuredMesh<3,3> (&mesh);
  70         _point_locator=new INTERP_KERNEL::PointLocatorInSimplex<MEDNormalizedUnstructuredMesh<3,3> >(*(static_cast<MEDNormalizedUnstructuredMesh<3,3>* >(_medmesh)));
  71         break;
  72       }
  73 }
  74
  75 void MEDMEM::PointLocatorInSimplex::getNodeWightsInSimplex( const std::vector<const double*>& nodes,
  76                                                             const double* p, double* weghts )
  77 {
  78   INTERP_KERNEL::barycentric_coords( nodes, p, weghts );
  79 }