gmpfield.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GMPFIELD_HH
#define DUNE_GMPFIELD_HH
 
#include <iostream>
#include <string>
#include <type_traits>
 
#if HAVE_GMP || DOXYGEN
 
#include <gmpxx.h>
 
#include <dune/common/promotiontraits.hh>
#include <dune/common/typetraits.hh>
 
namespace Dune
{
 
  template< unsigned int precision >
  class GMPField
    : public mpf_class
  {
    typedef mpf_class Base;
 
  public:
    GMPField ()
      : Base(0,precision)
    {}
 
    GMPField ( const char* str )
      : Base(str,precision)
    {}
 
    GMPField ( const std::string& str )
      : Base(str,precision)
    {}
 
    template< class T,
              typename EnableIf = typename std::enable_if<
                std::is_convertible<T, mpf_class>::value>::type
              >
    GMPField ( const T &v )
      : Base( v,precision )
    {}
 
    // type conversion operators
    operator double () const
    {
      return this->get_d();
    }
 
  };
 
  template <unsigned int precision>
  struct IsNumber<GMPField<precision>>
    : public std::integral_constant<bool, true> {
  };
 
  template< unsigned int precision1, unsigned int precision2 >
  struct PromotionTraits<GMPField<precision1>, GMPField<precision2>>
  {
    typedef GMPField<(precision1 > precision2 ? precision1 : precision2)> PromotedType;
  };
 
  template< unsigned int precision >
  struct PromotionTraits<GMPField<precision>,GMPField<precision>>
  {
    typedef GMPField<precision> PromotedType;
  };
 
  template< unsigned int precision, class T >
  struct PromotionTraits<GMPField<precision>, T>
  {
    typedef GMPField<precision> PromotedType;
  };
 
  template< class T, unsigned int precision >
  struct PromotionTraits<T, GMPField<precision>>
  {
    typedef GMPField<precision> PromotedType;
  };
}
 
#endif // HAVE_GMP
 
#endif // #ifndef DUNE_GMPFIELD_HH