ZZ.hpp

Go to the documentation of this file.

// Copyright 1995 Michael E. Stillman.
 
#ifndef _ZZ_hh_
#define _ZZ_hh_

 
#include "error.h"
#include "ring.hpp"
#include <cstddef>
 
namespace M2 {
class ARingZZGMP;
};
 
// The following lines are here only to remove complaints about old style casts
// from gmp
extern "C" inline int mask_mpz_cmp_si(mpz_srcptr x, long int i)
{
  return mpz_cmp_si(x, i);
}
extern "C" inline int mask_mpq_cmp_si(mpq_srcptr x, long int i, long int j)
{
  return mpq_cmp_si(x, i, j);
}

class RingZZ : public Ring
{
  friend class M2::ARingZZGMP;
 
  mpz_ptr new_elem() const;
 
  M2::ARingZZGMP *coeffR;
 
 protected:
  virtual ~RingZZ() {}
 public:
  typedef mpz_ptr element_type;

  // Creation of globalZZ is done in PolyRing::make_trivial_ZZ_poly_ring
  // These two routines should not be called from elsewhere
  RingZZ() {}
  bool initialize_ZZ(const PolynomialRing *deg_ring);
 
  RingZZ *cast_to_RingZZ() { return this; }
  const RingZZ *cast_to_RingZZ() const { return this; }
  M2::ARingZZGMP *get_ARing() const
  {
    return coeffR;
  }
 
  virtual MutableMatrix *makeMutableMatrix(size_t nrows,
                                           size_t ncols,
                                           bool dense) const;
 
  // The following are all the routines required by 'ring'
  virtual bool is_ZZ() const { return true; }
  virtual CoefficientType coefficient_type() const { return COEFF_ZZ; }
  virtual void text_out(buffer &o) const;
 
  static unsigned int mod_ui(mpz_srcptr n, unsigned int p);
  static std::pair<bool, int> get_si(mpz_srcptr n);
 
  // If the base ring of a is ZZ:
  // To get a bignum from a RingElement a, use: a.get_value().get_mpz()
  // To get a bignum from an ring_elem  a, use: a.get_mpz()
 
  virtual unsigned int computeHashValue(const ring_elem a) const;
 
  virtual std::pair<bool, long> coerceToLongInteger(ring_elem a) const;
 
  virtual ring_elem from_long(long n) const;
  virtual ring_elem from_int(mpz_srcptr n) const;
  virtual bool from_rational(mpq_srcptr q, ring_elem &result) const;
  virtual bool promote(const Ring *R,
                       const ring_elem f,
                       ring_elem &result) const;
  virtual bool lift(const Ring *R, const ring_elem f, ring_elem &result) const;
 
  virtual ring_elem preferred_associate(ring_elem f) const;
  virtual bool lower_associate_divisor(ring_elem &f, ring_elem g) const;
 
  virtual void lower_content(ring_elem &c, ring_elem g)
      const;  // c is a content elem, g is in ring
 
  int is_positive(const ring_elem a) const;
 
  virtual bool is_unit(const ring_elem f) const;
  virtual bool is_zero(const ring_elem f) const;
  virtual bool is_equal(const ring_elem f, const ring_elem g) const;
  virtual int compare_elems(const ring_elem f, const ring_elem g) const;
 
  virtual ring_elem copy(const ring_elem f) const;
  virtual void remove(ring_elem &f) const;
 
#if 0  
  void internal_negate_to(ring_elem &f) const;
  void internal_add_to(ring_elem &f, ring_elem &g) const;
  void internal_subtract_to(ring_elem &f, ring_elem &g) const;
#endif
  
  virtual ring_elem negate(const ring_elem f) const;
  virtual ring_elem add(const ring_elem f, const ring_elem g) const;
  virtual ring_elem subtract(const ring_elem f, const ring_elem g) const;
  virtual ring_elem mult(const ring_elem f, const ring_elem g) const;
  virtual ring_elem power(const ring_elem f, mpz_srcptr n) const;
  virtual ring_elem power(const ring_elem f, int n) const;
 
  virtual ring_elem invert(const ring_elem f) const;
 
  virtual ring_elem divide(const ring_elem f, const ring_elem g) const;
 
  virtual ring_elem remainder(const ring_elem f, const ring_elem g) const;
  virtual ring_elem quotient(const ring_elem f, const ring_elem g) const;
  virtual ring_elem remainderAndQuotient(const ring_elem f,
                                         const ring_elem g,
                                         ring_elem &quot) const;
 
  ring_elem gcd(const ring_elem f, const ring_elem g) const;
  ring_elem gcd_extended(const ring_elem f,
                         const ring_elem g,
                         ring_elem &u,
                         ring_elem &v) const;
 
  virtual void syzygy(const ring_elem a,
                      const ring_elem b,
                      ring_elem &x,
                      ring_elem &y) const;
 
  virtual ring_elem random() const;
 
  virtual void elem_text_out(buffer &o,
                             const ring_elem f,
                             bool p_one = true,
                             bool p_plus = false,
                             bool p_parens = false) const;
 
  virtual ring_elem eval(const RingMap *map,
                         const ring_elem f,
                         int first_var) const;
};
 
unsigned int computeHashValue_mpz(mpz_srcptr a);
#endif
 
// Local Variables:
// compile-command: "make -C $M2BUILDDIR/Macaulay2/e "
// indent-tabs-mode: nil
// End: