point.C

Go to the documentation of this file.

// The libMesh Finite Element Library.
// Copyright (C) 2002-2025 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner

// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.

// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.

// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA




// Local includes
#include "libmesh/point.h"

// C++ includes
#include <type_traits> // std::is_trivially_copyable


// Boost doesn't satisfy `is_trivially_copyable` for its float128
// wrapper, but surely it's memcpyable anyway?
#if !LIBMESH_DEFAULT_QUADRUPLE_PRECISION
static_assert(std::is_trivially_copyable<libMesh::TypeVector<libMesh::Point>>::value,
              "Someone made Point non-TriviallyCopyable");
#endif // !LIBMESH_DEFAULT_QUADRUPLE_PRECISION

// ------------------------------------------------------------
// Point class member functions
// unsigned int Point::key() const
// {
//   unsigned int tempx,tempy,tempz;

//   int i,j=2,cnt=0;
//   unsigned int index[3];
//   const Real deg = 1.e12;

//   tempx = static_cast<unsigned int>(((*this)(0)*deg));
//   tempy = static_cast<unsigned int>(((*this)(1)*deg));
//   tempz = static_cast<unsigned int>(((*t!his)(2)*deg));

//   index[0]=0;
//   index[1]=0;
//   index[2]=0;

//   for (i=sizeof(unsigned int)*8-1;i>=0;i--)
//     {
//       index[j] += (tempx >> i) & 01;
//       index[j]  = index[j] << 01;

//       if (( cnt % (sizeof(unsigned int)*8) == 0) && (cnt !=0 ) )
// {
//   cnt = 0;
//   j--;
// }
//       else
// cnt++;

//       index[j] += (tempy >> i) & 01;
//       index[j]  = index[j] << 01;

//       if (( cnt % (sizeof(unsigned int)*8) == 0) && (cnt !=0 ) )
// {
//   cnt = 0;
//   j--;
// }
//       else
// cnt++;

//       index[j] += (tempz >> i) & 01;
//       index[j]  = index[j] << 01;

//       if (( cnt % (sizeof(unsigned int)*8) == 0) && (cnt !=0 ) )
// {
//   cnt = 0;
//   j--;
// }
//       else
// cnt++;
//     }

//   return index[2];
// }