Branch data     Line data    Source code

       1                 :            : /**
       2                 :            :  * MOAB, a Mesh-Oriented datABase, is a software component for creating,
       3                 :            :  * storing and accessing finite element mesh data.
       4                 :            :  *
       5                 :            :  * Copyright 2004 Sandia Corporation.  Under the terms of Contract
       6                 :            :  * DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government
       7                 :            :  * retains certain rights in this software.
       8                 :            :  *
       9                 :            :  * This library is free software; you can redistribute it and/or
      10                 :            :  * modify it under the terms of the GNU Lesser General Public
      11                 :            :  * License as published by the Free Software Foundation; either
      12                 :            :  * version 2.1 of the License, or (at your option) any later version.
      13                 :            :  *
      14                 :            :  */
      15                 :            : 
      16                 :            : #ifndef MB_UTIL_HPP
      17                 :            : #define MB_UTIL_HPP
      18                 :            : 
      19                 :            : #include "moab/MOABConfig.h"
      20                 :            : #include "moab/Forward.hpp"
      21                 :            : #include "moab/CartVect.hpp"
      22                 :            : 
      23                 :            : #include <math.h>
      24                 :            : #if defined MOAB_HAVE_ISFINITE
      25                 :            : #define moab_isfinite( f ) isfinite( f )
      26                 :            : #elif defined MOAB_HAVE_STDISFINITE
      27                 :            : #include <cmath>
      28                 :            : #define moab_isfinite( f ) std::isfinite( f )
      29                 :            : #elif defined MOAB_HAVE_FINITE
      30                 :            : #define moab_isfinite( f ) finite( f )
      31                 :            : #else
      32                 :            : #define moab_isfinite( f ) ( !isinf( f ) && !isnan( f ) )
      33                 :            : #endif
      34                 :            : 
      35                 :            : namespace moab
      36                 :            : {
      37                 :            : 
      38                 :            : /** \class Util
      39                 :            :  *
      40                 :            :  * \brief Utility functions for computational geometry and mathematical calculations
      41                 :            :  */
      42                 :            : class Util
      43                 :            : {
      44                 :            :   public:
      45                 :            :     template < typename T >
      46                 :            :     static bool is_finite( T value );
      47                 :            : 
      48                 :            :     static void normal( Interface* MB, EntityHandle handle, double& x, double& y, double& z );
      49                 :            : 
      50                 :            :     static void centroid( Interface* MB, EntityHandle handle, CartVect& coord );
      51                 :            : 
      52                 :            :     // static void edge_centers(Interface *MB, EntityHandle handle, std::vector<CartVect>
      53                 :            :     // &coords_list);
      54                 :            : 
      55                 :            :     // static void face_centers(Interface *MB, EntityHandle handle, std::vector<CartVect>
      56                 :            :     // &coords_list);
      57                 :            : 
      58                 :            :   private:
      59                 :            :     Util() {}
      60                 :            : };
      61                 :            : 
      62                 :            : template < typename T >
      63                 :    1954558 : inline bool Util::is_finite( T value )
      64                 :            : {
      65                 :    1954558 :     return moab_isfinite( (double)value );
      66                 :            : }
      67                 :            : 
      68                 :            : }  // namespace moab
      69                 :            : 
      70                 :            : #endif