Util.hpp

Go to the documentation of this file.

/**
 * MOAB, a Mesh-Oriented datABase, is a software component for creating,
 * storing and accessing finite element mesh data.
 *
 * Copyright 2004 Sandia Corporation.  Under the terms of Contract
 * DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government
 * retains certain rights in this software.
 *
 * 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.
 *
 */

#ifndef MB_UTIL_HPP
#define MB_UTIL_HPP

#include "moab/MOABConfig.h"
#include "moab/Forward.hpp"
#include "moab/CartVect.hpp"

#include <cmath>
#if defined MOAB_HAVE_ISFINITE
#define moab_isfinite( f ) isfinite( f )
#elif defined MOAB_HAVE_STDISFINITE
#include <cmath>
#define moab_isfinite( f ) std::isfinite( f )
#elif defined MOAB_HAVE_FINITE
#define moab_isfinite( f ) finite( f )
#else
#define moab_isfinite( f ) ( !std::isinf( double( f ) ) && !std::isnan( double( f ) ) )
#endif

namespace moab
{

/** \class Util
 *
 * \brief Utility functions for computational geometry and mathematical calculations
 */
class Util
{
  public:
    template < typename T >
    static bool is_finite( T value );

    static void normal( Interface* MB, EntityHandle handle, double& x, double& y, double& z );

    static void centroid( Interface* MB, EntityHandle handle, CartVect& coord );

    // static void edge_centers(Interface *MB, EntityHandle handle, std::vector<CartVect>
    // &coords_list);

    // static void face_centers(Interface *MB, EntityHandle handle, std::vector<CartVect>
    // &coords_list);

  private:
    Util() {}
};

template < typename T >
inline bool Util::is_finite( T value )
{
    return moab_isfinite( (double)value );
}

}  // namespace moab

#endif