number.h

Go to the documentation of this file.

#ifndef __NUMBER_H__
#define __NUMBER_H__

#include <string>
#include <iostream>
#include <gmpxx.h>
#include <math.h>

#ifdef __MPF_BUILD__

namespace mpf_sdp {

typedef mpf_class Number;

#else

namespace sdp {

typedef double Number;

#endif

//
// Functions
//

template <class T, class U>
inline T convert_number(const U& num) { return (T) num; }

template <>
inline double convert_number(const mpf_class& num) { return num.get_d(); }

class double_to_mpf_class {
public:
  mpf_class operator()(double x) const { return mpf_class(x); }
};

class InsideCoefficient {
  const Number& num;

public:
  explicit InsideCoefficient(const Number& n) : num(n) { }

  friend std::ostream& operator<<(std::ostream& out, const InsideCoefficient& ic);
};

class NumberOperations {
public:
  static double abs(double num) { return num >= 0 ? num : -num; }

  static mpf_class abs(const mpf_class& num)
  {
    mpf_class ret;

    mpf_abs(ret.get_mpf_t(), num.get_mpf_t());

    return ret;
  }

  static double sqrt(double num) { return ::sqrt(num); }

  static mpf_class sqrt(const mpf_class& num)
  {
    mpf_class ret;
    
    mpf_sqrt(ret.get_mpf_t(), num.get_mpf_t());

    return ret;
  }

  static void from_string(const std::string& str, double& out);

  static void from_string(const std::string& str, mpf_class& out);
};

//
// Prototypes
//
extern std::ostream& operator<<(std::ostream& out, const InsideCoefficient& ic);

}  // namespace

#endif

// Local variables:
// mode: c++
// End: