2013-08-06 02:41:10 +00:00
|
|
|
/*
|
|
|
|
|
Copyright (c) 2013 Microsoft Corporation. All rights reserved.
|
|
|
|
|
Released under Apache 2.0 license as described in the file LICENSE.
|
|
|
|
|
|
|
|
|
|
Author: Soonho Kong
|
|
|
|
|
*/
|
|
|
|
|
#pragma once
|
|
|
|
|
#include <stdint.h>
|
|
|
|
|
#include <mpfr.h>
|
2013-09-13 03:04:10 +00:00
|
|
|
#include "util/numerics/mpz.h"
|
|
|
|
|
#include "util/numerics/mpq.h"
|
|
|
|
|
#include "util/numerics/mpbq.h"
|
2013-08-06 02:41:10 +00:00
|
|
|
|
|
|
|
|
namespace lean {
|
2013-08-13 07:30:18 +00:00
|
|
|
|
|
|
|
|
void set_mpfp_rnd(bool plus_inf);
|
|
|
|
|
mpfr_rnd_t get_mpfp_rnd();
|
|
|
|
|
|
2013-08-06 02:41:10 +00:00
|
|
|
/**
|
|
|
|
|
\brief Wrapper for MPFR
|
|
|
|
|
*/
|
|
|
|
|
class mpfp {
|
2013-08-08 02:31:12 +00:00
|
|
|
friend numeric_traits<mpfp>;
|
2013-08-06 02:41:10 +00:00
|
|
|
mpfr_t m_val;
|
|
|
|
|
|
2013-08-09 03:20:53 +00:00
|
|
|
static mpz_t const & zval(mpz const & v) { return v.m_val; }
|
|
|
|
|
static mpz_t & zval(mpz & v) { return v.m_val; }
|
|
|
|
|
static mpq_t const & qval(mpq const & v) { return v.m_val; }
|
|
|
|
|
static mpq_t & qval(mpq & v) { return v.m_val; }
|
|
|
|
|
|
2013-08-06 02:41:10 +00:00
|
|
|
public:
|
2013-08-09 03:20:53 +00:00
|
|
|
friend void swap(mpfp & a, mpfp & b) { mpfr_swap(a.m_val, b.m_val); }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
|
|
|
|
// Setter functions
|
|
|
|
|
mpfp & set(mpfp const & v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set(m_val, v.m_val, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(unsigned long int const v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_ui(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(long int const v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_si(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(float const v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_flt(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(double const v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_d(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(long double const v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_ld(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(mpz_t const & v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_z(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(mpq_t const & v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_q(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(mpf_t const & v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_f(m_val, v, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(mpz const & v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_z(m_val, v.m_val, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(mpq const & v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_q(m_val, v.m_val, MPFR_RNDN); return *this;
|
|
|
|
|
}
|
|
|
|
|
mpfp & set(mpbq const & v, mpfr_rnd_t rnd = MPFR_RNDN) {
|
|
|
|
|
mpfr_set_z(m_val, v.m_num.m_val, rnd); // this = m_num
|
|
|
|
|
mpfr_div_2ui(m_val, m_val, v.m_k, rnd); // this = m_num / (2^k)
|
|
|
|
|
return *this;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
// Assignment operators
|
|
|
|
|
mpfp & operator=(mpfp const & v) { return set(v); }
|
|
|
|
|
mpfp & operator=(unsigned long int const v) { return set(v); }
|
|
|
|
|
mpfp & operator=(long int const v) { return set(v); }
|
|
|
|
|
mpfp & operator=(float const v) { return set(v); }
|
|
|
|
|
mpfp & operator=(double const v) { return set(v); }
|
|
|
|
|
mpfp & operator=(long double const v) { return set(v); }
|
|
|
|
|
mpfp & operator=(mpz_t const & v) { return set(v); }
|
|
|
|
|
mpfp & operator=(mpq_t const & v) { return set(v); }
|
|
|
|
|
mpfp & operator=(mpf_t const & v) { return set(v); }
|
|
|
|
|
mpfp & operator=(mpz const & v) { return set(v); }
|
|
|
|
|
mpfp & operator=(mpq const & v) { return set(v); }
|
|
|
|
|
mpfp & operator=(mpbq const & v) { return set(v); }
|
|
|
|
|
|
|
|
|
|
// Basic Constructors
|
|
|
|
|
mpfp() { mpfr_init(m_val); } // with default precision
|
2013-08-09 03:20:53 +00:00
|
|
|
explicit mpfp(int prec) { mpfr_init2(m_val, prec); }
|
|
|
|
|
explicit mpfp(unsigned prec) { mpfr_init2(m_val, prec); }
|
|
|
|
|
explicit mpfp(mpfr_prec_t prec) { mpfr_init2(m_val, prec); }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
|
|
|
|
// Constructors using the default precision
|
2013-08-09 03:20:53 +00:00
|
|
|
explicit mpfp(float const v ):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(double const v ):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(long double const v):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(mpz_t const & v ):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(mpq_t const & v ):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(mpf_t const & v ):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(mpz const & v ):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(mpq const & v ):mpfp() { set(v); }
|
|
|
|
|
explicit mpfp(mpbq const & v ):mpfp() { set(v); }
|
2013-08-06 02:41:10 +00:00
|
|
|
mpfp(mpfp const & v ):mpfp(mpfr_get_prec(v.m_val)) { set(v); }
|
|
|
|
|
|
|
|
|
|
// Constructors using the provided precision
|
|
|
|
|
mpfp(float const v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(double const v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(long double const v, mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(mpz_t const & v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(mpq_t const & v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(mpf_t const & v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(mpz const & v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(mpq const & v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(mpbq const & v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
mpfp(mpfp const & v , mpfr_prec_t p):mpfp(p) { set(v); }
|
|
|
|
|
|
|
|
|
|
// Constructors using the provided precision and rounding mode
|
|
|
|
|
mpfp(float const v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(double const v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(long double const v, mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(mpz_t const & v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(mpq_t const & v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(mpf_t const & v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(mpz const & v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(mpq const & v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(mpbq const & v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
mpfp(mpfp const & v , mpfr_prec_t p, mpfr_rnd_t rnd):mpfp(p) { set(v, rnd); }
|
|
|
|
|
|
2013-08-09 03:20:53 +00:00
|
|
|
mpfp(mpfp && s):mpfp(mpfr_get_prec(s.m_val)) { mpfr_swap(m_val, s.m_val); }
|
2013-09-10 20:26:14 +00:00
|
|
|
~mpfp() { mpfr_clear(m_val); mpfr_free_cache(); }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
|
|
|
|
unsigned hash() const { return static_cast<unsigned>(mpfr_get_si(m_val, MPFR_RNDN)); }
|
|
|
|
|
|
|
|
|
|
int sgn() const { return mpfr_sgn(m_val); }
|
|
|
|
|
friend int sgn(mpfp const & a) { return a.sgn(); }
|
|
|
|
|
bool is_pos() const { return sgn() > 0; }
|
|
|
|
|
bool is_neg() const { return sgn() < 0; }
|
|
|
|
|
bool is_zero() const { return sgn() == 0; }
|
|
|
|
|
bool is_nonpos() const { return !is_pos(); }
|
|
|
|
|
bool is_nonneg() const { return !is_neg(); }
|
|
|
|
|
|
|
|
|
|
void neg() { mpfr_neg(m_val, m_val, MPFR_RNDN); }
|
|
|
|
|
friend mpfp neg(mpfp a) { a.neg(); return a; }
|
|
|
|
|
|
|
|
|
|
// Math functions and theirs friend functions
|
2013-08-13 07:30:18 +00:00
|
|
|
void abs (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_abs(m_val, m_val, rnd); }
|
|
|
|
|
void exp (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_exp(m_val, m_val, rnd); }
|
|
|
|
|
void exp2 (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_exp2(m_val, m_val, rnd); }
|
|
|
|
|
void exp10(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_exp10(m_val, m_val, rnd); }
|
|
|
|
|
void log (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_log(m_val, m_val, rnd); }
|
|
|
|
|
void log2 (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_log2(m_val, m_val, rnd); }
|
|
|
|
|
void log10(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_log10(m_val, m_val, rnd); }
|
|
|
|
|
void sin (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sin(m_val, m_val, rnd); }
|
|
|
|
|
void cos (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_cos(m_val, m_val, rnd); }
|
|
|
|
|
void tan (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_tan(m_val, m_val, rnd); }
|
|
|
|
|
void sec (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sec(m_val, m_val, rnd); }
|
|
|
|
|
void csc (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_csc(m_val, m_val, rnd); }
|
|
|
|
|
void cot (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_cot(m_val, m_val, rnd); }
|
|
|
|
|
void asin (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_asin(m_val, m_val, rnd); }
|
|
|
|
|
void acos (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_acos(m_val, m_val, rnd); }
|
|
|
|
|
void atan (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_atan(m_val, m_val, rnd); }
|
|
|
|
|
void sinh (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sinh(m_val, m_val, rnd); }
|
|
|
|
|
void cosh (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_cosh(m_val, m_val, rnd); }
|
|
|
|
|
void tanh (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_tanh(m_val, m_val, rnd); }
|
|
|
|
|
void asinh(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_asinh(m_val, m_val, rnd); }
|
|
|
|
|
void acosh(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_acosh(m_val, m_val, rnd); }
|
|
|
|
|
void atanh(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_atanh(m_val, m_val, rnd); }
|
|
|
|
|
friend mpfp abs (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.abs(rnd); return a; }
|
|
|
|
|
friend mpfp exp (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.exp(rnd); return a; }
|
|
|
|
|
friend mpfp exp2 (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.exp2(rnd); return a; }
|
|
|
|
|
friend mpfp exp10(mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.exp10(rnd); return a; }
|
|
|
|
|
friend mpfp log (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.log(rnd); return a; }
|
|
|
|
|
friend mpfp log2 (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.log2(rnd); return a; }
|
|
|
|
|
friend mpfp log10(mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.log10(rnd); return a; }
|
|
|
|
|
friend mpfp sin (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.sin(rnd); return a; }
|
|
|
|
|
friend mpfp cos (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.cos(rnd); return a; }
|
|
|
|
|
friend mpfp tan (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.tan(rnd); return a; }
|
|
|
|
|
friend mpfp sec (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.sec(rnd); return a; }
|
|
|
|
|
friend mpfp csc (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.csc(rnd); return a; }
|
|
|
|
|
friend mpfp cot (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.cot(rnd); return a; }
|
|
|
|
|
friend mpfp asin (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.asin(rnd); return a; }
|
|
|
|
|
friend mpfp acos (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.acos(rnd); return a; }
|
|
|
|
|
friend mpfp atan (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.atan(rnd); return a; }
|
|
|
|
|
friend mpfp sinh (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.sinh(rnd); return a; }
|
|
|
|
|
friend mpfp cosh (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.cosh(rnd); return a; }
|
|
|
|
|
friend mpfp tanh (mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.tanh(rnd); return a; }
|
|
|
|
|
friend mpfp asinh(mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.asinh(rnd); return a; }
|
|
|
|
|
friend mpfp acosh(mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.acosh(rnd); return a; }
|
|
|
|
|
friend mpfp atanh(mpfp a, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.atanh(rnd); return a; }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-09 03:20:53 +00:00
|
|
|
void inv() { mpfr_d_div(m_val, 1.0, m_val, MPFR_RNDN); }
|
|
|
|
|
friend mpfp inv(mpfp a) { a.inv(); return a; }
|
2013-08-13 07:30:18 +00:00
|
|
|
double get_double(mpfr_rnd_t rnd = get_mpfp_rnd()) const { return mpfr_get_d(m_val, rnd); }
|
|
|
|
|
float get_float (mpfr_rnd_t rnd = get_mpfp_rnd()) const { return mpfr_get_flt(m_val, rnd); }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
|
|
|
|
friend int cmp(mpfp const & a, mpfp const & b ) { return mpfr_cmp(a.m_val, b.m_val); }
|
|
|
|
|
friend int cmp(mpfp const & a, unsigned long int const b) { return mpfr_cmp_ui(a.m_val, b); }
|
|
|
|
|
friend int cmp(mpfp const & a, long int const b ) { return mpfr_cmp_si(a.m_val, b); }
|
|
|
|
|
friend int cmp(mpfp const & a, double const b ) { return mpfr_cmp_d (a.m_val, b); }
|
|
|
|
|
friend int cmp(mpfp const & a, long double const b ) { return mpfr_cmp_ld(a.m_val, b); }
|
|
|
|
|
friend int cmp(mpfp const & a, mpz_t const & b ) { return mpfr_cmp_z (a.m_val, b); }
|
|
|
|
|
friend int cmp(mpfp const & a, mpq_t const & b ) { return mpfr_cmp_q (a.m_val, b); }
|
|
|
|
|
friend int cmp(mpfp const & a, mpf_t const & b ) { return mpfr_cmp_f (a.m_val, b); }
|
2013-08-09 03:20:53 +00:00
|
|
|
friend int cmp(mpfp const & a, mpz const & b ) { return mpfr_cmp_z (a.m_val, zval(b)); }
|
|
|
|
|
friend int cmp(mpfp const & a, mpq const & b ) { return mpfr_cmp_q (a.m_val, qval(b)); }
|
|
|
|
|
//friend int cmp(mpfp const & a, mpbq const & b) { return mpfr_cmp_(a.m_val, b); }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
|
|
|
|
friend bool operator<(mpfp const & a, mpfp const & b ) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, mpz const & b ) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, unsigned long int const b) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, long int const b ) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, double const b ) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, long double const b ) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, mpz_t const & b ) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, mpq_t const & b ) { return cmp(a, b) < 0; }
|
|
|
|
|
friend bool operator<(mpfp const & a, mpf_t const & b ) { return cmp(a, b) < 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator<(mpz const & a , mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
friend bool operator<(unsigned long int const a, mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
friend bool operator<(long int const a , mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
friend bool operator<(double const a , mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
friend bool operator<(long double const a , mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
friend bool operator<(mpz_t const & a , mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
friend bool operator<(mpq_t const & a , mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
friend bool operator<(mpf_t const & a , mpfp const & b) { return cmp(b, a) > 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator>(mpfp const & a, mpfp const & b ) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, mpz const & b ) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, unsigned long int const b) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, long int const b ) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, double const b ) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, long double const b ) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, mpz_t const & b ) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, mpq_t const & b ) { return cmp(a, b) > 0; }
|
|
|
|
|
friend bool operator>(mpfp const & a, mpf_t const & b ) { return cmp(a, b) > 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator>(mpz const & a , mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
friend bool operator>(unsigned long int const a, mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
friend bool operator>(long int const a , mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
friend bool operator>(double const a , mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
friend bool operator>(long double const a , mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
friend bool operator>(mpz_t const & a , mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
friend bool operator>(mpq_t const & a , mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
friend bool operator>(mpf_t const & a , mpfp const & b) { return cmp(b, a) < 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator<=(mpfp const & a, mpfp const & b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, mpz const & b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, unsigned long int const b) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, long int const b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, double const b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, long double const b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, mpz_t const & b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, mpq_t const & b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
friend bool operator<=(mpfp const & a, mpf_t const & b ) { return cmp(a, b) <= 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator<=(mpz const & a , mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
friend bool operator<=(unsigned long int const a, mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
friend bool operator<=(long int const a , mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
friend bool operator<=(double const a , mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
friend bool operator<=(long double const a , mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
friend bool operator<=(mpz_t const & a , mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
friend bool operator<=(mpq_t const & a , mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
friend bool operator<=(mpf_t const & a , mpfp const & b) { return cmp(b, a) >= 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator>=(mpfp const & a, mpfp const & b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, mpz const & b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, unsigned long int const b) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, long int const b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, double const b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, long double const b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, mpz_t const & b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, mpq_t const & b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
friend bool operator>=(mpfp const & a, mpf_t const & b ) { return cmp(a, b) >= 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator>=(mpz const & a , mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
friend bool operator>=(unsigned long int const a, mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
friend bool operator>=(long int const a , mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
friend bool operator>=(double const a , mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
friend bool operator>=(long double const a , mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
friend bool operator>=(mpz_t const & a , mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
friend bool operator>=(mpq_t const & a , mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
friend bool operator>=(mpf_t const & a , mpfp const & b) { return cmp(b, a) <= 0; }
|
|
|
|
|
|
|
|
|
|
friend bool operator==(mpfp const & a, mpfp const & b) { return mpfr_equal_p(a.m_val, b.m_val) != 0; }
|
2013-08-09 03:20:53 +00:00
|
|
|
friend bool operator==(unsigned long int const a, mpfp const & b) { return cmp(b, a) == 0; }
|
|
|
|
|
friend bool operator==(long int const a , mpfp const & b) { return cmp(b, a) == 0; }
|
|
|
|
|
friend bool operator==(double const a , mpfp const & b) { return cmp(b, a) == 0; }
|
|
|
|
|
friend bool operator==(long double const a , mpfp const & b) { return cmp(b, a) == 0; }
|
|
|
|
|
friend bool operator==(mpz_t const & a , mpfp const & b) { return cmp(b, a) == 0; }
|
|
|
|
|
friend bool operator==(mpq_t const & a , mpfp const & b) { return cmp(b, a) == 0; }
|
|
|
|
|
friend bool operator==(mpf_t const & a , mpfp const & b) { return cmp(b, a) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, mpz const & b ) { return cmp(a, b) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, unsigned long int const b) { return cmp(a, b) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, long int const b ) { return cmp(a, b) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, double const b ) { return cmp(a, b) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, long double const b ) { return cmp(a, b) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, mpz_t const & b ) { return cmp(a, b) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, mpq_t const & b ) { return cmp(a, b) == 0; }
|
|
|
|
|
friend bool operator==(mpfp const & a, mpf_t const & b ) { return cmp(a, b) == 0; }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-09-13 19:49:03 +00:00
|
|
|
friend bool operator!=(mpfp const & a, mpfp const & b) { return !operator==(a, b); }
|
2013-08-09 03:20:53 +00:00
|
|
|
friend bool operator!=(unsigned long int const a, mpfp const & b) { return cmp(b, a) != 0; }
|
|
|
|
|
friend bool operator!=(long int const a , mpfp const & b) { return cmp(b, a) != 0; }
|
|
|
|
|
friend bool operator!=(double const a , mpfp const & b) { return cmp(b, a) != 0; }
|
|
|
|
|
friend bool operator!=(long double const a , mpfp const & b) { return cmp(b, a) != 0; }
|
|
|
|
|
friend bool operator!=(mpz_t const & a , mpfp const & b) { return cmp(b, a) != 0; }
|
|
|
|
|
friend bool operator!=(mpq_t const & a , mpfp const & b) { return cmp(b, a) != 0; }
|
|
|
|
|
friend bool operator!=(mpf_t const & a , mpfp const & b) { return cmp(b, a) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, mpz const & b ) { return cmp(a, b) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, unsigned long int const b) { return cmp(a, b) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, long int const b ) { return cmp(a, b) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, double const b ) { return cmp(a, b) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, long double const b ) { return cmp(a, b) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, mpz_t const & b ) { return cmp(a, b) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, mpq_t const & b ) { return cmp(a, b) != 0; }
|
|
|
|
|
friend bool operator!=(mpfp const & a, mpf_t const & b ) { return cmp(a, b) != 0; }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-13 07:30:18 +00:00
|
|
|
mpfp & add(mpfp const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & add(unsigned long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add_ui(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & add(long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add_si(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & add(double const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add_d(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & add(mpz_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add_z(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & add(mpq_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add_q(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & add(mpz const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add_z(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & add(mpq const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_add_q(m_val, m_val, o.m_val, rnd); return *this; }
|
2013-08-06 02:41:10 +00:00
|
|
|
mpfp & operator+=(mpfp const & o) { return add(o); }
|
2013-08-08 02:31:12 +00:00
|
|
|
mpfp & operator+=(unsigned long int o) { return add(o); }
|
|
|
|
|
mpfp & operator+=(long int const o) { return add(o); }
|
|
|
|
|
mpfp & operator+=(double const o) { return add(o); }
|
|
|
|
|
mpfp & operator+=(mpz_t const & o) { return add(o); }
|
|
|
|
|
mpfp & operator+=(mpq_t const & o) { return add(o); }
|
|
|
|
|
mpfp & operator+=(mpz const & o) { return add(o); }
|
|
|
|
|
mpfp & operator+=(mpq const & o) { return add(o); }
|
|
|
|
|
friend mpfp operator+(mpfp a, mpfp const & b) { return a += b; }
|
|
|
|
|
friend mpfp operator+(mpfp a, unsigned long const b) { return a += b; }
|
|
|
|
|
friend mpfp operator+(mpfp a, long int const b) { return a += b; }
|
|
|
|
|
friend mpfp operator+(mpfp a, double const b) { return a += b; }
|
|
|
|
|
friend mpfp operator+(mpfp a, mpz_t const & b) { return a += b; }
|
|
|
|
|
friend mpfp operator+(mpfp a, mpq_t const & b) { return a += b; }
|
|
|
|
|
friend mpfp operator+(mpfp a, mpz const & b) { return a += b; }
|
|
|
|
|
friend mpfp operator+(mpfp a, mpq const & b) { return a += b; }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-13 07:30:18 +00:00
|
|
|
mpfp & sub(mpfp const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & sub(unsigned long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub_ui(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & sub(long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub_si(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & sub(double const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub_d(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & sub(mpz_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub_z(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & sub(mpq_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub_q(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & sub(mpz const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub_z(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & sub(mpq const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_sub_q(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rsub(unsigned long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_ui_sub(m_val, o, m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rsub(long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_si_sub(m_val, o, m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rsub(double const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_d_sub(m_val, o, m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rsub(mpz_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_z_sub(m_val, o, m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rsub(mpz const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_z_sub(m_val, o.m_val, m_val, rnd); return *this; }
|
2013-08-06 02:41:10 +00:00
|
|
|
mpfp & operator-=(mpfp const & o) { return sub(o); }
|
2013-08-08 02:31:12 +00:00
|
|
|
mpfp & operator-=(unsigned long int o) { return sub(o); }
|
|
|
|
|
mpfp & operator-=(long int const o) { return sub(o); }
|
|
|
|
|
mpfp & operator-=(double const o) { return sub(o); }
|
|
|
|
|
mpfp & operator-=(mpz_t const & o) { return sub(o); }
|
|
|
|
|
mpfp & operator-=(mpq_t const & o) { return sub(o); }
|
|
|
|
|
mpfp & operator-=(mpz const & o) { return sub(o); }
|
|
|
|
|
mpfp & operator-=(mpq const & o) { return sub(o); }
|
|
|
|
|
friend mpfp operator-(mpfp a, mpfp const & b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(mpfp a, unsigned long const b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(mpfp a, long int const b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(mpfp a, double const b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(mpfp a, mpz_t const & b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(mpfp a, mpq_t const & b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(mpfp a, mpz const & b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(mpfp a, mpq const & b) { return a -= b; }
|
|
|
|
|
friend mpfp operator-(unsigned long const a, mpfp b) { return b.rsub(a); }
|
|
|
|
|
friend mpfp operator-(long int const a, mpfp b) { return b.rsub(a); }
|
|
|
|
|
friend mpfp operator-(double const a, mpfp b) { return b.rsub(a); }
|
|
|
|
|
friend mpfp operator-(mpz_t const & a, mpfp b) { return b.rsub(a); }
|
|
|
|
|
// friend mpfp operator-(mpq_t const & a, mpfp b) { return b.rsub(a); }
|
|
|
|
|
friend mpfp operator-(mpz const & a, mpfp b) { return b.rsub(a); }
|
|
|
|
|
// friend mpfp operator-(mpq const & a, mpfp b) { return b.rsub(a); }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-13 07:30:18 +00:00
|
|
|
mpfp & mul(mpfp const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & mul(unsigned long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul_ui(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & mul(long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul_si(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & mul(double const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul_d(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & mul(mpz_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul_z(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & mul(mpq_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul_q(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & mul(mpz const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul_z(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & mul(mpq const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_mul_q(m_val, m_val, o.m_val, rnd); return *this; }
|
2013-08-06 02:41:10 +00:00
|
|
|
mpfp & operator*=(mpfp const & o) { return mul(o); }
|
2013-08-08 02:31:12 +00:00
|
|
|
mpfp & operator*=(unsigned long int o) { return mul(o); }
|
|
|
|
|
mpfp & operator*=(long int const o) { return mul(o); }
|
|
|
|
|
mpfp & operator*=(double const o) { return mul(o); }
|
|
|
|
|
mpfp & operator*=(mpz_t const & o) { return mul(o); }
|
|
|
|
|
mpfp & operator*=(mpq_t const & o) { return mul(o); }
|
|
|
|
|
mpfp & operator*=(mpz const & o) { return mul(o); }
|
|
|
|
|
mpfp & operator*=(mpq const & o) { return mul(o); }
|
|
|
|
|
friend mpfp operator*(mpfp a, mpfp const & b) { return a *= b; }
|
|
|
|
|
friend mpfp operator*(mpfp a, unsigned long const b) { return a *= b; }
|
|
|
|
|
friend mpfp operator*(mpfp a, long int const b) { return a *= b; }
|
|
|
|
|
friend mpfp operator*(mpfp a, double const b) { return a *= b; }
|
|
|
|
|
friend mpfp operator*(mpfp a, mpz_t const & b) { return a *= b; }
|
|
|
|
|
friend mpfp operator*(mpfp a, mpq_t const & b) { return a *= b; }
|
|
|
|
|
friend mpfp operator*(mpfp a, mpz const & b) { return a *= b; }
|
|
|
|
|
friend mpfp operator*(mpfp a, mpq const & b) { return a *= b; }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-13 07:30:18 +00:00
|
|
|
mpfp & div(mpfp const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & div(unsigned long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div_ui(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & div(long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div_si(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & div(double const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div_d(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & div(mpz_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div_z(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & div(mpq_t const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div_q(m_val, m_val, o, rnd); return *this; }
|
|
|
|
|
mpfp & div(mpz const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div_z(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & div(mpq const & o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_div_q(m_val, m_val, o.m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rdiv(unsigned long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_ui_div(m_val, o, m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rdiv(long int const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_si_div(m_val, o, m_val, rnd); return *this; }
|
|
|
|
|
mpfp & rdiv(double const o, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_d_div(m_val, o, m_val, rnd); return *this; }
|
2013-08-06 02:41:10 +00:00
|
|
|
mpfp & operator/=(mpfp const & o) { return div(o); }
|
2013-08-08 02:31:12 +00:00
|
|
|
mpfp & operator/=(unsigned long int o) { return div(o); }
|
|
|
|
|
mpfp & operator/=(long int const o) { return div(o); }
|
|
|
|
|
mpfp & operator/=(double const o) { return div(o); }
|
|
|
|
|
mpfp & operator/=(mpz_t const & o) { return div(o); }
|
|
|
|
|
mpfp & operator/=(mpq_t const & o) { return div(o); }
|
|
|
|
|
mpfp & operator/=(mpz const & o) { return div(o); }
|
|
|
|
|
mpfp & operator/=(mpq const & o) { return div(o); }
|
2013-08-06 02:41:10 +00:00
|
|
|
friend mpfp operator/(mpfp a, mpfp const & b) { return a /= b; }
|
2013-08-08 02:31:12 +00:00
|
|
|
friend mpfp operator/(mpfp a, unsigned long const b) { return a /= b; }
|
|
|
|
|
friend mpfp operator/(mpfp a, long int const b) { return a /= b; }
|
|
|
|
|
friend mpfp operator/(mpfp a, double const b) { return a /= b; }
|
|
|
|
|
friend mpfp operator/(mpfp a, mpz_t const & b) { return a /= b; }
|
|
|
|
|
friend mpfp operator/(mpfp a, mpq_t const & b) { return a /= b; }
|
|
|
|
|
friend mpfp operator/(mpfp a, mpz const & b) { return a /= b; }
|
|
|
|
|
friend mpfp operator/(mpfp a, mpq const & b) { return a /= b; }
|
|
|
|
|
friend mpfp operator/(unsigned long const a, mpfp b) { return b.rdiv(a); }
|
|
|
|
|
friend mpfp operator/(long int const a, mpfp b) { return b.rdiv(a); }
|
|
|
|
|
friend mpfp operator/(double const a, mpfp b) { return b.rdiv(a); }
|
|
|
|
|
// friend mpfp operator/(mpz_t const & a, mpfp b) { return a.rdiv(b); }
|
|
|
|
|
// friend mpfp operator/(mpq_t const & a, mpfp b) { return a.rdiv(b); }
|
|
|
|
|
// friend mpfp operator/(mpz const & a, mpfp b) { return a.rdiv(b); }
|
|
|
|
|
// friend mpfp operator/(mpq const & a, mpfp b) { return a.rdiv(b); }
|
|
|
|
|
|
2013-08-15 06:47:36 +00:00
|
|
|
mpfp & operator^=(unsigned long int k) { power(k); return *this; }
|
2013-08-08 02:31:12 +00:00
|
|
|
friend mpfp operator^(mpfp a, unsigned long int k) { return a ^= k; }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-09 03:20:53 +00:00
|
|
|
mpfp & operator++() { return operator+=(1lu); }
|
|
|
|
|
mpfp operator++(int) { mpfp r(*this); ++(*this); return r; }
|
|
|
|
|
|
|
|
|
|
mpfp & operator--() { return operator-=(1lu); }
|
|
|
|
|
mpfp operator--(int) { mpfp r(*this); --(*this); return r; }
|
|
|
|
|
|
2013-08-14 07:01:03 +00:00
|
|
|
mpfp operator-() const { mpfp t = *this; t.neg(); return t; }
|
|
|
|
|
|
2013-08-09 03:20:53 +00:00
|
|
|
void floor() { mpfr_floor(m_val, m_val); }
|
|
|
|
|
void ceil () { mpfr_ceil (m_val, m_val); }
|
|
|
|
|
void round() { mpfr_round(m_val, m_val); }
|
|
|
|
|
void trunc() { mpfr_trunc(m_val, m_val); }
|
2013-08-13 07:30:18 +00:00
|
|
|
void rfloor(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_rint_floor(m_val, m_val, rnd); }
|
|
|
|
|
void rceil (mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_rint_ceil (m_val, m_val, rnd); }
|
|
|
|
|
void rround(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_rint_round(m_val, m_val, rnd); }
|
|
|
|
|
void rtrunc(mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_rint_trunc(m_val, m_val, rnd); }
|
2013-08-09 03:20:53 +00:00
|
|
|
|
|
|
|
|
friend mpfp floor(mpfp const & a) { static thread_local mpfp tmp; tmp = a; tmp.floor(); return tmp; }
|
|
|
|
|
friend mpfp ceil (mpfp const & a) { static thread_local mpfp tmp; tmp = a; tmp.ceil(); return tmp; }
|
|
|
|
|
friend mpfp round(mpfp const & a) { static thread_local mpfp tmp; tmp = a; tmp.round(); return tmp; }
|
|
|
|
|
friend mpfp trunc(mpfp const & a) { static thread_local mpfp tmp; tmp = a; tmp.trunc(); return tmp; }
|
2013-08-13 07:30:18 +00:00
|
|
|
friend mpfp rfloor(mpfp const & a, mpfr_rnd_t rnd = get_mpfp_rnd()) {
|
2013-08-09 03:20:53 +00:00
|
|
|
static thread_local mpfp tmp; tmp = a; tmp.rfloor(rnd); return tmp;
|
|
|
|
|
}
|
2013-08-13 07:30:18 +00:00
|
|
|
friend mpfp rceil (mpfp const & a, mpfr_rnd_t rnd = get_mpfp_rnd()) {
|
2013-08-09 03:20:53 +00:00
|
|
|
static thread_local mpfp tmp; tmp = a; tmp.rceil(rnd); return tmp;
|
|
|
|
|
}
|
2013-08-13 07:30:18 +00:00
|
|
|
friend mpfp rround(mpfp const & a, mpfr_rnd_t rnd = get_mpfp_rnd()) {
|
2013-08-09 03:20:53 +00:00
|
|
|
static thread_local mpfp tmp; tmp = a; tmp.rround(rnd); return tmp;
|
|
|
|
|
}
|
2013-08-13 07:30:18 +00:00
|
|
|
friend mpfp rtrunc(mpfp const & a, mpfr_rnd_t rnd = get_mpfp_rnd()) {
|
2013-08-09 03:20:53 +00:00
|
|
|
static thread_local mpfp tmp; tmp = a; tmp.rtrunc(rnd); return tmp;
|
|
|
|
|
}
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-13 07:30:18 +00:00
|
|
|
void power(mpfp const & b, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_pow(m_val, m_val, b.m_val, rnd); }
|
|
|
|
|
void power(unsigned long int b, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_pow_ui(m_val, m_val, b, rnd); }
|
|
|
|
|
void power(long int b, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_pow_si(m_val, m_val, b, rnd); }
|
|
|
|
|
void power(mpz_t const & b, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_pow_z(m_val, m_val, b, rnd); }
|
|
|
|
|
void power(mpz const & b, mpfr_rnd_t rnd = get_mpfp_rnd()) { mpfr_pow_z(m_val, m_val, b.m_val, rnd); }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
2013-08-13 07:30:18 +00:00
|
|
|
friend mpfp power(mpfp a, mpfp const & b, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.power(b, rnd); return a;; }
|
|
|
|
|
friend mpfp power(mpfp a, unsigned long int b, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.power(b, rnd); return a;; }
|
|
|
|
|
friend mpfp power(mpfp a, long int b, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.power(b, rnd); return a;; }
|
|
|
|
|
friend mpfp power(mpfp a, mpz_t const & b, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.power(b, rnd); return a;; }
|
|
|
|
|
friend mpfp power(mpfp a, mpz const & b, mpfr_rnd_t rnd = get_mpfp_rnd()) { a.power(b, rnd); return a;; }
|
2013-08-06 02:41:10 +00:00
|
|
|
|
|
|
|
|
friend std::ostream & operator<<(std::ostream & out, mpfp const & v);
|
|
|
|
|
};
|
|
|
|
|
|
2013-08-08 02:31:33 +00:00
|
|
|
// Macro to implement transcendental functions
|
|
|
|
|
#define LEAN_TRANS_MPFP_FUNC(f, v, rnd) v.f(rnd);
|
|
|
|
|
|
2013-08-06 02:41:10 +00:00
|
|
|
template<>
|
|
|
|
|
class numeric_traits<mpfp> {
|
|
|
|
|
public:
|
2013-08-12 22:22:08 +00:00
|
|
|
static mpfr_rnd_t rnd() { return get_mpfp_rnd(); }
|
2013-08-06 02:41:10 +00:00
|
|
|
static bool precise() { return true; }
|
|
|
|
|
static bool is_zero(mpfp const & v) { return v.is_zero(); }
|
|
|
|
|
static bool is_pos(mpfp const & v) { return v.is_pos(); }
|
|
|
|
|
static bool is_neg(mpfp const & v) { return v.is_neg(); }
|
2013-08-12 22:22:08 +00:00
|
|
|
static void set_rounding(bool plus_inf) { set_mpfp_rnd(plus_inf); }
|
2013-08-06 02:41:10 +00:00
|
|
|
static void neg(mpfp & v) { v.neg(); }
|
2013-08-14 03:05:04 +00:00
|
|
|
static void inv(mpfp & v) { v.inv(); }
|
2013-08-09 03:20:53 +00:00
|
|
|
static void reset(mpfp & v) { v = 0.0; }
|
2013-08-06 02:41:10 +00:00
|
|
|
// v <- v^k
|
2013-08-15 06:47:36 +00:00
|
|
|
static void power(mpfp & v, unsigned k) { v.power(static_cast<unsigned long>(k)); }
|
2013-08-14 03:05:04 +00:00
|
|
|
static void abs(mpfp & v) { v.abs(); }
|
|
|
|
|
static void ceil(mpfp & v) { v.ceil(); }
|
|
|
|
|
static void floor(mpfp & v) { v.floor(); }
|
2013-08-08 02:31:33 +00:00
|
|
|
|
|
|
|
|
// constants
|
|
|
|
|
static mpfp pi_l;
|
|
|
|
|
static mpfp pi_n;
|
|
|
|
|
static mpfp pi_u;
|
|
|
|
|
static inline mpfp pi_lower() {
|
|
|
|
|
mpfr_const_pi(pi_l.m_val, MPFR_RNDD);
|
|
|
|
|
return pi_l;
|
|
|
|
|
}
|
|
|
|
|
static inline mpfp pi() {
|
|
|
|
|
mpfr_const_pi(pi_n.m_val, MPFR_RNDN);
|
|
|
|
|
return pi_n;
|
|
|
|
|
}
|
|
|
|
|
static inline mpfp pi_upper() {
|
|
|
|
|
mpfr_const_pi(pi_u.m_val, MPFR_RNDU);
|
|
|
|
|
return pi_u;
|
|
|
|
|
}
|
2013-08-14 03:06:19 +00:00
|
|
|
static inline mpfp pi_half_lower() { return pi_lower() / 2lu; }
|
|
|
|
|
static inline mpfp pi_half() { return pi() / 2lu; }
|
|
|
|
|
static inline mpfp pi_half_upper() { return pi_upper() / 2lu; }
|
|
|
|
|
static inline mpfp pi_twice_lower() { return pi_lower() * 2lu; }
|
|
|
|
|
static inline mpfp pi_twice() { return pi() * 2lu; }
|
|
|
|
|
static inline mpfp pi_twice_upper() { return pi_upper() * 2lu; }
|
2013-08-08 02:31:33 +00:00
|
|
|
|
|
|
|
|
// Transcendental functions
|
2013-08-12 22:22:08 +00:00
|
|
|
static void exp(mpfp & v) { LEAN_TRANS_MPFP_FUNC(exp, v, rnd()); }
|
|
|
|
|
static void exp2(mpfp & v) { LEAN_TRANS_MPFP_FUNC(exp2, v, rnd()); }
|
|
|
|
|
static void exp10(mpfp & v) { LEAN_TRANS_MPFP_FUNC(exp10, v, rnd()); }
|
|
|
|
|
static void log(mpfp & v) { LEAN_TRANS_MPFP_FUNC(log, v, rnd()); }
|
|
|
|
|
static void log2(mpfp & v) { LEAN_TRANS_MPFP_FUNC(log2, v, rnd()); }
|
|
|
|
|
static void log10(mpfp & v) { LEAN_TRANS_MPFP_FUNC(log10, v, rnd()); }
|
|
|
|
|
static void sin(mpfp & v) { LEAN_TRANS_MPFP_FUNC(sin, v, rnd()); }
|
|
|
|
|
static void cos(mpfp & v) { LEAN_TRANS_MPFP_FUNC(cos, v, rnd()); }
|
|
|
|
|
static void tan(mpfp & v) { LEAN_TRANS_MPFP_FUNC(tan, v, rnd()); }
|
|
|
|
|
static void sec(mpfp & v) { LEAN_TRANS_MPFP_FUNC(sec, v, rnd()); }
|
|
|
|
|
static void csc(mpfp & v) { LEAN_TRANS_MPFP_FUNC(csc, v, rnd()); }
|
|
|
|
|
static void cot(mpfp & v) { LEAN_TRANS_MPFP_FUNC(cot, v, rnd()); }
|
|
|
|
|
static void asin(mpfp & v) { LEAN_TRANS_MPFP_FUNC(asin, v, rnd()); }
|
|
|
|
|
static void acos(mpfp & v) { LEAN_TRANS_MPFP_FUNC(acos, v, rnd()); }
|
|
|
|
|
static void atan(mpfp & v) { LEAN_TRANS_MPFP_FUNC(atan, v, rnd()); }
|
|
|
|
|
static void sinh(mpfp & v) { LEAN_TRANS_MPFP_FUNC(sinh, v, rnd()); }
|
|
|
|
|
static void cosh(mpfp & v) { LEAN_TRANS_MPFP_FUNC(cosh, v, rnd()); }
|
|
|
|
|
static void tanh(mpfp & v) { LEAN_TRANS_MPFP_FUNC(tanh, v, rnd()); }
|
|
|
|
|
static void asinh(mpfp & v) { LEAN_TRANS_MPFP_FUNC(asinh, v, rnd()); }
|
|
|
|
|
static void acosh(mpfp & v) { LEAN_TRANS_MPFP_FUNC(acosh, v, rnd()); }
|
|
|
|
|
static void atanh(mpfp & v) { LEAN_TRANS_MPFP_FUNC(atanh, v, rnd()); }
|
2013-08-06 02:41:10 +00:00
|
|
|
};
|
|
|
|
|
}
|
