rings_test.h

Go to the documentation of this file.

#include "common.h"
using namespace std;
 
// the following headers are private...
 
 
#include "polys/monomials/ring.h"
#include "polys/monomials/p_polys.h"
 
 
 
 
class MyGlobalPrintingFixture : public GlobalPrintingFixture
{
  public:
    virtual bool setUpWorld()
    {
 
      GlobalPrintingFixture::setUpWorld();
 
 
      //TS_ASSERT_EQUALS( nRegister( n_Zp, npInitChar), n_Zp );
      //TS_ASSERT_EQUALS( nRegister( n_GF, nfInitChar), n_GF );
      //TS_ASSERT_EQUALS( nRegister( n_R, nrInitChar), n_R );
      //TS_ASSERT_EQUALS( nRegister( n_Q, nlInitChar), n_Q );
      //TS_ASSERT_EQUALS( nRegister( n_R, nrInitChar), n_R );
 
#ifdef HAVE_RINGS
      //TS_ASSERT_EQUALS( nRegister( n_Z, nrzInitChar), n_Z ); // these are UNusable at the moment!
#endif
 
      return true;
    }
};
 
 
//
// We can rely on this file being included exactly once
// and declare this global variable in the header file.
//
static MyGlobalPrintingFixture globalPrintingFixture;
 
 
namespace
{
  void PrintRing(const ring r)
  {
    rWrite(r); PrintLn();
#ifdef RDEBUG
    rDebugPrint(r); PrintLn();
#endif
  }
static inline void Test(const ring)
{
  // does nothing yet...
}
 
}
 
class PolysTestSuite : public CxxTest::TestSuite
{
public:
  void test_Z13_t()
  {
    clog << "Creating  Z/13[t]: " << endl;
 
    char* n[] = { (char*)"t"};
    ring r = rDefault( 13, 1, n);
    TS_ASSERT_DIFFERS( r, NULLp );
 
    PrintRing(r);
 
    TS_ASSERT( rField_is_Domain(r) );
    TS_ASSERT( !rField_is_Q(r) );
 
    TS_ASSERT( rField_is_Zp(r) );
    TS_ASSERT( !rField_is_Zp(r, 11) );
    TS_ASSERT( rField_is_Zp(r, 13) );
 
    TS_ASSERT_EQUALS( rVar(r), 1);
 
    Test(r);
 
    rDelete(r);
  }
 
  void test_QQ_t()
  {
    clog << "Creating  Q[s]: " << endl;
 
    char* n[] = {(char*)"s"};
    ring r = rDefault( 0, 1, n);
    TS_ASSERT_DIFFERS( r, NULLp );
 
    PrintRing(r);
 
    TS_ASSERT( rField_is_Domain(r) );
    TS_ASSERT( rField_is_Q(r) );
 
    TS_ASSERT( !rField_is_Zp(r) );
    TS_ASSERT( !rField_is_Zp(r, 11) );
 
    TS_ASSERT_EQUALS( rVar(r), 1);
 
    Test(r);
 
    rDelete(r);
  }
 
  void test_Z11_x_y_z()
  {
     clog << "Creating  Z/11[x, y, z]: " << endl;
 
     char* n[] = {(char*)"x", (char*)"y", (char*)"z"};
     ring r = rDefault( 11, 3, n);
     TS_ASSERT_DIFFERS( r, NULLp );
 
     PrintRing(r);
 
     TS_ASSERT( rField_is_Domain(r) );
     TS_ASSERT( !rField_is_Q(r) );
 
     TS_ASSERT( rField_is_Zp(r) );
     TS_ASSERT( rField_is_Zp(r, 11) );
     TS_ASSERT( !rField_is_Zp(r, 13) );
 
     TS_ASSERT_EQUALS( rVar(r), 3);
 
     Test(r);
 
     rDelete(r);
  }
 
  void test_QQ_x_y_z()
  {
     clog << "Creating  QQ[x, y, z, u]: " << endl;
 
     char* n[] = {(char*)"x", (char*)"y", (char*)"z", (char*)"u"};
     ring r = rDefault( 0, 4, n);
     TS_ASSERT_DIFFERS( r, NULLp );
 
     PrintRing(r);
 
     TS_ASSERT( rField_is_Domain(r) );
     TS_ASSERT( rField_is_Q(r) );
 
     TS_ASSERT( !rField_is_Zp(r) );
     TS_ASSERT( !rField_is_Zp(r, 11) );
 
     TS_ASSERT_EQUALS( rVar(r), 4);
 
     Test(r);
 
     rDelete(r);
   }
 
 
   void test_Z13_t_GF()
   {
     clog << "Creating  GF[t]: " << endl;
 
     char* n[] = {(char*)"t"};
 
     GFInfo param;
 
     param.GFChar= 5;
     param.GFDegree= 2;
     param.GFPar_name= (const char*)"Q";
 
     const coeffs cf = nInitChar( n_GF, &param );
 
     TS_ASSERT_DIFFERS( cf, NULLp );
 
     ring r = rDefault( cf, 1, n);  // now cf belongs to r!
     TS_ASSERT_DIFFERS( r, NULLp );
 
     PrintRing(r);
 
     TS_ASSERT( rField_is_Domain(r) );
     TS_ASSERT( !rField_is_Q(r) );
 
     TS_ASSERT( !rField_is_Zp(r) );
     TS_ASSERT( !rField_is_Zp(r, 11) );
     TS_ASSERT( !rField_is_Zp(r, 13) );
     TS_ASSERT( rField_is_GF(r) );
 
     TS_ASSERT( rField_is_GF(r, 5) );
     TS_ASSERT( !rField_is_GF(r, 25) );
 
     TS_ASSERT_EQUALS( rVar(r), 1);
 
     Test(r);
 
     rDelete(r); // kills 'cf' as well!
   }
 
 
 
 
 
  void test_QQ_x_y_z_ord_2_lp()
  {
     clog << "Creating  QQ[x, y, z, u]: " << endl;
 
     const int N = 4;
     char* n[N] = {(char*)"x", (char*)"y", (char*)"z", (char*)"u"};
 
 
     const int D = 2;
     /*order: lp,0*/
     rRingOrder_t *order = (rRingOrder_t *) omAlloc0(D* sizeof(rRingOrder_t));
     int *block0 = (int *)omAlloc0(D * sizeof(int));
     int *block1 = (int *)omAlloc0(D * sizeof(int));
     /* ringorder dp for the first block: var 1..N */
     order[0]  = ringorder_dp;
     block0[0] = 1;
     block1[0] = N;
     /* the last block: everything is 0 */
//     order[1]  = 0;
 
     ring r = rDefault(0, N, n, D, order, block0, block1);
     TS_ASSERT_DIFFERS( r, NULLp );
 
     PrintRing(r);
 
     TS_ASSERT( rField_is_Domain(r) );
     TS_ASSERT( rField_is_Q(r) );
 
     TS_ASSERT( !rField_is_Zp(r) );
     TS_ASSERT( !rField_is_Zp(r, 11) );
 
     TS_ASSERT_EQUALS( rVar(r), N);
 
     Test(r);
 
     rDelete(r);
   }
};