ARingZZpTest.cpp

Go to the documentation of this file.

// Copyright 2012-2013 Michael E. Stillman

 
#include <cstdio>
#include <string>
#include <iostream>
#include <sstream>
#include <memory>
#include <gtest/gtest.h>
#include <mpfr.h>
 
#include "reader.hpp"
#include "ZZp.hpp"
#include "aring-zzp-ffpack.hpp"
#include "aring-zzp.hpp"
#include "ARingTest.hpp"
 
static bool maxH_initialized = false;
static mpz_t maxH;
 
gmp_ZZ getRandomInteger()
{
  if (!maxH_initialized)
    {
      maxH_initialized = true;
      mpz_init(maxH);
      mpz_set_str(maxH, "100000000000", 10);
    }
  return rawRandomInteger(maxH);
}
 
template <>
void getElement<M2::ARingZZp>(const M2::ARingZZp& R,
                              int index,
                              M2::ARingZZp::ElementType& result)
{
  if (index < 50)
    R.set_from_long(result, index - 25);
  else
    {
      gmp_ZZ a = getRandomInteger();
      R.set_from_mpz(result, a);
    }
}
 
template <typename RT>
void testCoerceToLongInteger(const RT& R)
{
  long top1 = (R.characteristic() > 10000 ? 10000 : R.characteristic());
  long bottom2 = (R.characteristic() > 10000 ? R.characteristic() - 10000 : 0);
  for (long i = 0; i < top1; i++)
    {
      typename RT::ElementType a;
      R.init(a);
      R.set_from_long(a, i);
      long b = R.coerceToLongInteger(a);
      if (b < 0) b += R.characteristic();
      EXPECT_EQ(b, i);
    }
  for (long i = bottom2; i < R.characteristic(); i++)
    {
      typename RT::ElementType a;
      R.init(a);
      R.set_from_long(a, i);
      long b = R.coerceToLongInteger(a);
      if (b < 0) b += R.characteristic();
      EXPECT_EQ(b, i);
    }
}
TEST(RingZZp, create)
{
  const Z_mod* R = Z_mod::create(101);
  EXPECT_FALSE(R == nullptr);
  buffer o;
  o << "Ring being tested: ";
  R->text_out(o);
  fprintf(stdout, "%s\n", o.str());
}
 
TEST(ARingZZp, create)
{
  M2::ARingZZp R(101);
 
  M2::ARingZZp::ElementType a;
  buffer o;
 
  ARingElementGenerator<M2::ARingZZp> gen(R);
  R.init(a);
  gen.nextElement(a);
 
  EXPECT_EQ(ringName(R), "AZZ/101");
  EXPECT_EQ(R.cardinality(), 101);
  EXPECT_EQ(R.characteristic(), 101);
  // Now check what the generator is, as an integer
  R.init(a);
  R.set_var(a, 0);
  R.elem_text_out(o, a, true, true, false);
  std::cout << "generator is " << o.str() << std::endl;
  R.clear(a);
}
 
TEST(ARingZZp, arithmetic101)
{
  M2::ARingZZp R(101);
  testFiniteField(R, ntrials);
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZp, arithmetic32749)
{
  M2::ARingZZp R(32749);
  testFiniteField(R, ntrials);
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZp, arithmetic2)
{
  M2::ARingZZp R(2);
  testFiniteField(R, ntrials);
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZp, arithmetic3)
{
  M2::ARingZZp R(3);
  testFiniteField(R, ntrials);
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZp, fromStream)
{
  std::istringstream i("+1234 +345 -235*a");
  M2::ARingZZp R(32003);
  M2::ARingZZp::ElementType a;
  R.init(a);
  while (true)
    {
      while (isspace(i.peek())) i.get();
 
      if (!isdigit(i.peek()) && i.peek() != '+' && i.peek() != '-') break;
 
      fromStream(i, R, a);
 
      buffer o;
      R.elem_text_out(o, a);
      std::cout << o.str() << " peek: " << i.peek() << std::endl;
    }
  R.clear(a);
}
 
template <>
void getElement<M2::ARingZZpFFPACK>(const M2::ARingZZpFFPACK& R,
                                    int index,
                                    M2::ARingZZpFFPACK::ElementType& result)
{
  if (index < 50)
    R.set_from_long(result, index - 25);
  else
    {
      gmp_ZZ a = getRandomInteger();
      R.set_from_mpz(result, a);
    }
}
 
TEST(ARingZZpFFPACK, create)
{
  M2::ARingZZpFFPACK R(101);
 
  EXPECT_EQ(ringName(R), "ZZpFPACK(101,1)");
  testSomeMore(R);
 
  std::cout << "max modulus for ffpack zzp: "
            << M2::ARingZZpFFPACK::getMaxModulus() << std::endl;
 
  M2::ARingZZpFFPACK::ElementType a;
  R.init(a);
  R.set_from_long(a, 99);
  R.set_from_long(a, 101);
  R.set_from_long(a, 103);
  R.clear(a);
}
 
TEST(ARingZZpFFPACK, arithmetic101)
{
  M2::ARingZZpFFPACK R(101);
  testFiniteField(R, ntrials);
}
 
// TODO: commented out because it appears wrong.  Perhaps p=2 isn't allowed
// here?
// strange: does not fail on my thinkpad...
TEST(ARingZZpFFPACK, arithmetic2)
{
  M2::ARingZZpFFPACK R(2);
  testFiniteField(R, ntrials);
 
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZpFFPACK, arithmetic3)
{
  M2::ARingZZpFFPACK R(3);
  testFiniteField(R, ntrials);
 
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZpFFPACK, arithmetic66000007)
{
  M2::ARingZZpFFPACK R(66000007);
 
  testCoerceToLongInteger(R);
 
  testCoercions(R);
  testNegate(R, ntrials);
  testAdd(R, ntrials);
  testSubtract(R, ntrials);
  testMultiply(R, ntrials);
  testDivide(R, ntrials);
  testReciprocal(R, ntrials);
  testPower(R, ntrials);
  testAxioms(R, ntrials);
}
 
// TODO: commented out because it takes too long:
// Actually: now this characteristic seems too big?!
TEST(ARingZZpFFPACK, arithmetic67108859)
{
  M2::ARingZZpFFPACK R(67108859);
 
  testCoerceToLongInteger(R);
 
  testCoercions(R);
  testNegate(R, ntrials);
  testAdd(R, ntrials);
  testSubtract(R, ntrials);
  testMultiply(R, ntrials);
  testDivide(R, ntrials);
  testReciprocal(R, ntrials);
  testPower(R, ntrials);
  testAxioms(R, ntrials);
}
 
TEST(ARingZZpFFPACK, arithmetic33500479)
{
  M2::ARingZZpFFPACK R(33500479);
 
  testCoerceToLongInteger(R);
 
  testCoercions(R);
  testNegate(R, ntrials);
  testAdd(R, ntrials);
  testSubtract(R, ntrials);
  testMultiply(R, ntrials);
  testDivide(R, ntrials);
  testReciprocal(R, ntrials);
  testPower(R, ntrials);
  testAxioms(R, ntrials);
}
 
TEST(ARingZZp, read)
{
  std::string a = "-42378489327498312749c3";
  std::istringstream i(a);
 
  M2::ARingZZp R(101);
  M2::Reader<M2::ARingZZp> reader(R);
  M2::ARingZZp::ElementType b, c;
  R.init(b);
  R.init(c);
  reader.read(i, b);
  R.set_from_long(c, 3);
 
  EXPECT_TRUE(R.is_equal(b, c));
}

// Flint ZZ/p arithmetic ///
#include "../aring-zzp-flint.hpp"
template <>
void getElement<M2::ARingZZpFlint>(const M2::ARingZZpFlint& R,
                                   int index,
                                   M2::ARingZZpFlint::ElementType& result)
{
  if (index < 50)
    R.set_from_long(result, index - 25);
  else
    {
      gmp_ZZ a = getRandomInteger();
      R.set_from_mpz(result, a);
    }
}
 
TEST(ARingZZpFlint, create)
{
  M2::ARingZZpFlint R(101);
 
  EXPECT_EQ(ringName(R), "AZZFlint/101");
  testSomeMore(R);
 
  M2::ARingZZpFlint::ElementType a;
  R.init(a);
  R.set_from_long(a, 99);
  R.set_from_long(a, 101);
  R.set_from_long(a, 103);
  R.clear(a);
}
 
TEST(ARingZZpFlint, arithmetic101)
{
  M2::ARingZZpFlint R(101);
  testFiniteField(R, ntrials);
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZpFlint, arithmetic2)
{
  M2::ARingZZpFlint R(2);
  testFiniteField(R, ntrials);
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZpFlint, arithmetic3)
{
  M2::ARingZZpFlint R(3);
  testFiniteField(R, ntrials);
  testCoerceToLongInteger(R);
}
 
TEST(ARingZZpFlint, arithmetic66000007)
{
  M2::ARingZZpFlint R(66000007);
 
  testCoerceToLongInteger(R);
 
  testCoercions(R);
  testNegate(R, ntrials);
  testAdd(R, ntrials);
  testSubtract(R, ntrials);
  testMultiply(R, ntrials);
  testDivide(R, ntrials);
  testReciprocal(R, ntrials);
  testPower(R, ntrials);
  testAxioms(R, ntrials);
}
 
TEST(ARingZZpFlint, arithmetic67108859)
{
  M2::ARingZZpFlint R(67108859);
 
  testCoerceToLongInteger(R);
 
  testCoercions(R);
  testNegate(R, ntrials);
  testAdd(R, ntrials);
  testSubtract(R, ntrials);
  testMultiply(R, ntrials);
  testDivide(R, ntrials);
  testReciprocal(R, ntrials);
  testPower(R, ntrials);
  testAxioms(R, ntrials);
}
 
TEST(ARingZZpFlint, arithmetic33500479)
{
  M2::ARingZZpFlint R(33500479);
 
  testCoerceToLongInteger(R);
 
  testCoercions(R);
  testNegate(R, ntrials);
  testAdd(R, ntrials);
  testSubtract(R, ntrials);
  testMultiply(R, ntrials);
  testDivide(R, ntrials);
  testReciprocal(R, ntrials);
  testPower(R, ntrials);
  testAxioms(R, ntrials);
}
 
TEST(ARingZZpFlint, arithmetic9223372036854775783)
{
  // largest prime < 2^63
  if (sizeof(unsigned long) <= 4)
    std::cout << "seems to be a 32bit machine: skipping the test" << std::endl;
  else
    {
      M2::ARingZZpFlint R(9223372036854775783L);
 
      testCoerceToLongInteger(R);
 
      testCoercions(R);
      testNegate(R, ntrials);
      testAdd(R, ntrials);
      testSubtract(R, ntrials);
      testMultiply(R, ntrials);
      testDivide(R, ntrials);
      testReciprocal(R, ntrials);
      //  testPower(R, ntrials);  // this test fails: as it expects the
      //  characteristic to fit into an int.
      testAxioms(R, ntrials);
    }
}
 
// Even though flint handles primes up to 2^64-1, M2 assumes that the
// characteristic is < 2^63.  If needed, we could change the type of
// the characteristic to unsigned long, but that would have further
// complications.
TEST(ARingZZpFlint, arithmetic18446744073709551557)
{
  // largest prime < 2^64
  if (sizeof(unsigned long) <= 4)
    std::cout << "seems to be a 32bit machine: skipping the test" << std::endl;
  else
    {
      M2::ARingZZpFlint R(18446744073709551557UL);
 
      //    testCoerceToLongInteger(R); // this fails for charac > 2^63
 
      testCoercions(R);
      testNegate(R, ntrials);
      testAdd(R, ntrials);
      testSubtract(R, ntrials);
      testMultiply(R, ntrials);
      testDivide(R, ntrials);
      testReciprocal(R, ntrials);
      //  testPower(R, ntrials);  // this test fails: as it expects the
      //  characteristic to fit into an int.
      testAxioms(R, ntrials);
    }
}
 
// Local Variables:
// compile-command: "make -C $M2BUILDDIR/Macaulay2/e/unit-tests check  "
// indent-tabs-mode: nil
// End: