aring-qq-flint.hpp

Go to the documentation of this file.

// Copyright 2013 Michael E. Stillman
 
#ifndef _aring_QQ_flint_hpp_
#define _aring_QQ_flint_hpp_

 
#include "interface/gmp-util.h"  // for mpz_reallocate_limbs
 
// The following needs to be included before any flint files are included.
#include <M2/gc-include.h>
 
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wconversion"
#include <flint/flint.h>  // for fmpq_denref, fmpq_numref, fmpq, flin...
#include <flint/fmpq.h>   // for fmpq_init, fmpq_set, fmpq_set_mpq
#include <flint/fmpz.h>   // for fmpz_get_ui, fmpz_cmp_si, fmpz_is_one
#pragma GCC diagnostic pop
 
#include "aring.hpp"
#include "buffer.hpp"
#include "ringelem.hpp"
#include "exceptions.hpp"
 
// promote needs ring.hpp.  After moving promote out, remove it here!
#include "ring.hpp"
 
namespace M2 {
 
class ARingQQFlint : public SimpleARing<ARingQQFlint>
{
 public:
  static const RingID ringID = ring_QQFlint;
 
  typedef fmpq ElementType;
  typedef ElementType elem;
 
  ARingQQFlint();
  ~ARingQQFlint();
 
 public:
  // ring informational
  size_t characteristic() const { return 0; }
  size_t cardinality() const { return static_cast<size_t>(-1); }
  unsigned int computeHashValue(const ElementType& a) const
  {
    unsigned long numhash = fmpz_get_ui(fmpq_numref(&a));
    unsigned long denhash = fmpz_get_ui(fmpq_denref(&a));
    return static_cast<unsigned int>(13253 * numhash + 7647 * denhash);
  }

 
  bool is_pm_one(const ElementType& f) const
  {
    return fmpq_is_one(&f) or ((fmpz_cmp_si(fmpq_numref(&f), -1) == 0) and
                               fmpz_is_one(fmpq_denref(&f)));
  }
  bool is_unit(const ElementType& f) const { return not is_zero(f); }
  bool is_zero(const ElementType& f) const { return fmpq_is_zero(&f); }

 
  bool is_equal(const ElementType& f, const ElementType& g) const
  {
    return fmpq_equal(&f, &g);
  }
  int compare_elems(const ElementType& f, const ElementType& g) const
  {
    int cmp = fmpq_cmp(&f, &g);
    if (cmp > 0) return 1;
    if (cmp < 0) return -1;
    return 0;
  }


 
  void init_set(ElementType& result, const ElementType& a) const
  {
    fmpq_init(&result);
    fmpq_set(&result, &a);
  }
 
  void init(ElementType& result) const { fmpq_init(&result); }
  static void clear(ElementType& result) { fmpq_clear(&result); }
  void set(ElementType& result, const ElementType& a) const
  {
    fmpq_set(&result, &a);
  }
 
  void set_zero(ElementType& result) const { fmpq_zero(&result); }
  void set_from_long(ElementType& result, long a) const
  {
    fmpq_set_si(&result, a, 1);
  }
 
  void set_from_mpz(ElementType& result, mpz_srcptr a) const
  {
    // printf("ARingQQFlint::calling set_from_mpz\n");
    fmpz_set_mpz(fmpq_numref(&result), a);
    fmpz_one(fmpq_denref(&result));
  }
 
  bool set_from_mpq(ElementType& result, mpq_srcptr a) const
  {
    fmpq_set_mpq(&result, a);
    return true;
  }
 
  bool set_from_BigReal(ElementType& result, gmp_RR a) const
  {
    (void) result;
    (void) a;
    return false;
  }
 
  void set_var(ElementType& result, int v) const
  {
    (void) v;
    fmpq_set_si(&result, 1, 1);
  }


  void negate(ElementType& result, const ElementType& a) const
  {
    fmpq_neg(&result, &a);
  }
 
  bool invert(ElementType& result, const ElementType& a) const
  {
    if (is_unit(a))
      {
        fmpq_inv(&result, &a);
        return true;
      }
    set_zero(result);
    return false;
  }
 
  void add(ElementType& result,
           const ElementType& a,
           const ElementType& b) const
  {
    fmpq_add(&result, &a, &b);
  }
 
  void subtract(ElementType& result,
                const ElementType& a,
                const ElementType& b) const
  {
    fmpq_sub(&result, &a, &b);
  }
 
  void subtract_multiple(ElementType& result,
                         const ElementType& a,
                         const ElementType& b) const
  {
    fmpq_submul(&result, &a, &b);
  }
 
  void mult(ElementType& result,
            const ElementType& a,
            const ElementType& b) const
  {
    fmpq_mul(&result, &a, &b);
  }

  void divide(ElementType& result,
              const ElementType& a,
              const ElementType& b) const
  {
    if (is_zero(b)) throw exc::division_by_zero_error();
    fmpq_div(&result, &a, &b);
  }
 
  void power(ElementType& result, const ElementType& a, long n) const
  {
    return fmpq_pow_si(&result, &a, n);
  }
 
  void power_mpz(ElementType& result,
                 const ElementType& a,
                 mpz_srcptr n) const
  {
    std::pair<bool, int> n1 = RingZZ::get_si(n);
    if (n1.first)
      power(result, a, n1.second);
    else
      throw exc::engine_error("exponent too large");
  }
 
  void syzygy(const ElementType& a,
              const ElementType& b,
              ElementType& x,
              ElementType& y) const;

  void swap(ElementType& a, ElementType& b) const { fmpq_swap(&a, &b); }
  void random(ElementType& result) const
  {
    fmpq_randtest(&result, mRandomState, mMaxHeight);
  }


  void text_out(buffer& o) const { o << "QQFlint"; }
  void elem_text_out(buffer& o,
                     const ElementType& a,
                     bool p_one = true,
                     bool p_plus = false,
                     bool p_parens = false) const;

 
  void to_ring_elem(ring_elem& result, const ElementType& a) const
  {
    mpq_ptr b = getmemstructtype(mpq_ptr);
    mpq_init(b);
    fmpq_get_mpq(b, &a);
    mpz_reallocate_limbs(mpq_numref(b));
    mpz_reallocate_limbs(mpq_denref(b));
    result = ring_elem(b);
  }
 
  void from_ring_elem(ElementType& result, const ring_elem& a) const
  {
    fmpq_set_mpq(&result, a.get_mpq());
  }

  const ElementType from_ring_elem_const(const ring_elem& a) const
  {
    mpq_srcptr a1 = a.get_mpq();
    fmpq result;
    result.num = PTR_TO_COEFF(mpq_numref(a1));
    result.den = PTR_TO_COEFF(mpq_denref(a1));
    return result;
  }

 
  bool promote(const Ring* Rf, const ring_elem f, ElementType& result) const
  {
    // printf("ARingQQFlint::calling promote\n");
    // Rf = ZZ ---> QQ
    if (Rf->is_ZZ())
      {
        set_from_mpz(result, f.get_mpz());
        return true;
      }
    return false;
  }
 
  bool lift(const Ring* Rg, const ElementType& f, ring_elem& result) const
  {
    (void) Rg;
    (void) f;
    (void) result;
    return false;
  }
 
  // map : this --> target(map)
  //       primelem --> map->elem(first_var)
  // evaluate map(f)
  void eval(const RingMap* map,
            const ElementType& f,
            int first_var,
            ring_elem& result) const;
 
 private:
  mutable flint_rand_t mRandomState;
  long int mMaxHeight;
};
};
 
#endif
 
// Local Variables:
// compile-command: "make -C $M2BUILDDIR/Macaulay2/e  "
// indent-tabs-mode: nil
// End: