number2.h File Reference

#include "kernel/mod2.h"
#include "omalloc/omalloc.h"
#include "coeffs/coeffs.h"
#include "kernel/structs.h"

Go to the source code of this file.

Functions
char *	crString (coeffs c)
void	crPrint (coeffs cf)
BOOLEAN	jjR2_CR (leftv res, leftv a)
BOOLEAN	jjCRING_Zp (leftv res, leftv a, leftv b)
BOOLEAN	jjCRING_Zm (leftv res, leftv a, leftv b)
BOOLEAN	jjEQUAL_CR (leftv res, leftv a, leftv b)

Function Documentation

crPrint()

void crPrint

(

coeffs

cf

)

Definition at line 25 of file number2.cc.

{
  char *s=crString(c);
  PrintS(s);
  omFree(s);
}

crString()

char* crString

(

coeffs

c

)

Definition at line 17 of file number2.cc.

{
  if (c==NULL)
  {
    return omStrDup("oo");
  }
  return omStrDup(nCoeffName(c));
}

jjCRING_Zm()

BOOLEAN jjCRING_Zm	(	leftv	res,
		leftv	a,
		leftv	b
	)

Definition at line 78 of file number2.cc.

{
  coeffs c1=(coeffs)a->Data();
  number i2=(number)b->Data();
  if (c1->type==n_Z)
  {
    ZnmInfo info;
    mpz_t modBase;
    mpz_init(modBase);
    nlGMP(i2,modBase,coeffs_BIGINT); // FIXME? TODO? // extern void   nlGMP(number &i, mpz_t n, const coeffs r); // to be replaced with n_MPZ(modBase,i2,coeffs_BIGINT); // ?
    info.base= (mpz_ptr)modBase;
    info.exp= 1;
    if (mpz_popcount(modBase)==1) // is a power of 2
    {
      // is exponent <=2^(8*sizeof(unsigned long))
      mp_bitcnt_t l=mpz_scan1 (modBase,0);
      if ((l>0) && (l<=8*sizeof(unsigned long)))
      {
        res->data=(void *) nInitChar(n_Z2m,(void*)(long)l);
      }
      else
      {
        mpz_set_ui(modBase,2);
        info.exp= l;
        res->data=(void *) nInitChar(n_Znm,&info);
      }
      mpz_clear(modBase);
      return FALSE;
    }
    res->data=(void *)nInitChar(n_Zn,&info);
    mpz_clear(modBase);
    return FALSE;
  }
  return TRUE;
}

jjCRING_Zp()

BOOLEAN jjCRING_Zp	(	leftv	res,
		leftv	a,
		leftv	b
	)

Definition at line 35 of file number2.cc.

{
  coeffs c1=(coeffs)a->Data();
  int    i2=(int)(long)b->Data();
  if (c1->type==n_Z)
  {
    if (i2==IsPrime(i2))
    {
      res->data=(void *)nInitChar(n_Zp,(void*)(long)i2);
    }
    else
    {
      ZnmInfo info;
      mpz_t modBase;
      mpz_init_set_ui(modBase,i2);
      info.base= modBase;
      info.exp= 1;
      if (mpz_popcount((mpz_ptr)modBase)==1) // is a power of 2
      {
        i2=SI_LOG2(i2);
        // is exponent <=2^(8*sizeof(unsigned long))
        if (i2<(8*sizeof(unsigned long)))
        {
          mpz_clear(modBase);
          res->data=(void *) nInitChar(n_Z2m,(void*)(long)i2);
          return FALSE;
        }
        else
        {
          mpz_set_ui(modBase,2);
          info.exp=i2;
          res->data=(void *) nInitChar(n_Znm,&info);
          mpz_clear(modBase);
          return FALSE;
        }
      }
      res->data=(void *)nInitChar(n_Zn,&info);
      mpz_clear(modBase);
    }
    return FALSE;
  }
  return TRUE;
}

jjEQUAL_CR()

BOOLEAN jjEQUAL_CR	(	leftv	res,
		leftv	a,
		leftv	b
	)

Definition at line 114 of file number2.cc.

{
  coeffs a2=(coeffs)a->Data();
  coeffs b2=(coeffs)b->Data();
  res->data=(void*)(long)(a2==b2);
  return FALSE;
}

jjR2_CR()

BOOLEAN jjR2_CR	(	leftv	res,
		leftv	a
	)

Definition at line 122 of file number2.cc.

{
  ring r=(ring)a->Data();
  AlgExtInfo extParam;
  extParam.r = r;
  coeffs cf=nInitChar(n_transExt,&extParam);
  res->data=(void*)cf;
  return FALSE;
}