flintcf_Q.cc File Reference

#include <ctype.h>
#include "misc/auxiliary.h"
#include <flint/flint.h>
#include <flint/fmpz.h>
#include <flint/fmpq.h>
#include <flint/fmpq_poly.h>
#include "factory/factory.h"
#include "coeffs/coeffs.h"
#include "coeffs/numbers.h"
#include "coeffs/longrat.h"

Go to the source code of this file.

Typedefs
typedef fmpq_poly_struct *	fmpq_poly_ptr
typedef fmpz *	fmpz_ptr

Functions
static char *	nlEatLong (char *s, mpz_ptr i)
static BOOLEAN	CoeffIsEqual (const coeffs r, n_coeffType n, void *)
static void	SetChar (const coeffs)
static number	Mult (number a, number b, const coeffs)
static number	Sub (number a, number b, const coeffs)
static number	Add (number a, number b, const coeffs)
static number	Div (number a, number b, const coeffs)
static number	ExactDiv (number a, number b, const coeffs)
static number	Init (long i, const coeffs)
static number	InitMPZ (mpz_t i, const coeffs)
static int	Size (number n, const coeffs)
static long	Int (number &n, const coeffs)
static void	MPZ (mpz_t result, number &n, const coeffs)
static number	Neg (number a, const coeffs)
static number	Invers (number a, const coeffs)
static number	Copy (number a, const coeffs)
static BOOLEAN	IsOne (number a, const coeffs)
static BOOLEAN	IsZero (number a, const coeffs)
static void	WriteShort (number a, const coeffs r)
static const char *	Read (const char *st, number *a, const coeffs r)
static void	Normalize (number &a, const coeffs)
static BOOLEAN	Greater (number a, number b, const coeffs)
static BOOLEAN	Equal (number a, number b, const coeffs)
static BOOLEAN	IsMOne (number k, const coeffs)
static BOOLEAN	GreaterZero (number, const coeffs)
static void	Power (number a, int i, number *result, const coeffs)
static number	GetDenom (number &n, const coeffs)
static number	GetNumerator (number &n, const coeffs)
static number	Gcd (number a, number b, const coeffs)
static number	ExtGcd (number a, number b, number *s, number *t, const coeffs)
static number	Lcm (number, number, const coeffs)
static void	Delete (number *a, const coeffs)
static nMapFunc	SetMap (const coeffs, const coeffs)
static number	Farey (number, number, const coeffs)
static number	ChineseRemainder (number *, number *, int, BOOLEAN, CFArray &, const coeffs)
static int	ParDeg (number x, const coeffs)
static number	Parameter (const int, const coeffs)
static void	WriteFd (number a, const ssiInfo *d, const coeffs)
static number	ReadFd (const ssiInfo *d, const coeffs)
static number	ConvFactoryNSingN (const CanonicalForm, const coeffs)
static CanonicalForm	ConvSingNFactoryN (number, BOOLEAN, const coeffs)
char *	CoeffName (const coeffs r)
coeffs	flintQInitCfByName (char *s, n_coeffType n)
static BOOLEAN	DBTest (number, const char *, const int, const coeffs)
static void	KillChar (coeffs cf)
BOOLEAN	flintQ_InitChar (coeffs cf, void *infoStruct)

Typedef Documentation

fmpq_poly_ptr

typedef fmpq_poly_struct* fmpq_poly_ptr

Definition at line 24 of file flintcf_Q.cc.

fmpz_ptr

typedef fmpz* fmpz_ptr

Definition at line 25 of file flintcf_Q.cc.

Function Documentation

Add()

static number Add ( number  a, number  b, const  coeffs  )

static

Definition at line 70 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_add(res,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
  return (number)res;
}

ChineseRemainder()

static number ChineseRemainder ( number *  , number *  , int  , BOOLEAN  , CFArray &  , const  coeffs  )

static

Definition at line 439 of file flintcf_Q.cc.

{
  WerrorS("not yet: ChineseRemainder");
  return NULL;
}

CoeffIsEqual()

static BOOLEAN CoeffIsEqual ( const coeffs  r, n_coeffType  n, void *    )

static

Definition at line 48 of file flintcf_Q.cc.

{
  return (r->type==n);
}

CoeffName()

char * CoeffName

(

const coeffs

r

)

Definition at line 519 of file flintcf_Q.cc.

{
  STATIC_VAR char CoeffName_flint_Q[20];
  sprintf(CoeffName_flint_Q,"flintQp[%s]",r->pParameterNames[0]);
  return (char*)CoeffName_flint_Q;
 
}

ConvFactoryNSingN()

static number ConvFactoryNSingN ( const  CanonicalForm, const  coeffs  )

static

Definition at line 509 of file flintcf_Q.cc.

{
  WerrorS("not yet: ConvFactoryNSingN");
  return NULL;
}

ConvSingNFactoryN()

static CanonicalForm ConvSingNFactoryN ( number  , BOOLEAN  , const  coeffs  )

static

Definition at line 514 of file flintcf_Q.cc.

{
  WerrorS("not yet: ConvSingNFactoryN");
  return CanonicalForm(0);
}

Copy()

static number Copy ( number  a, const  coeffs  )

static

Definition at line 200 of file flintcf_Q.cc.

{
 fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
 fmpq_poly_init(res);
 fmpq_poly_set(res,(fmpq_poly_ptr)a);
 return (number)res;
}

DBTest()

static BOOLEAN DBTest ( number  , const char *  , const int  , const  coeffs  )

static

Definition at line 544 of file flintcf_Q.cc.

{
  return TRUE;
}

Delete()

static void Delete ( number *  a, const  coeffs  )

static

Definition at line 400 of file flintcf_Q.cc.

{
  if ((*a)!=NULL)
  {
    fmpq_poly_clear((fmpq_poly_ptr)*a);
    omFree(*a);
    *a=NULL;
  }
}

Div()

static number Div ( number  a, number  b, const  coeffs  )

static

Definition at line 77 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  if(fmpq_poly_is_zero((fmpq_poly_ptr)b))
  {
     WerrorS(nDivBy0);
  }
  else
  {
    fmpq_poly_div(res,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
    fmpq_poly_t mod;
    fmpq_poly_init(mod);
    fmpq_poly_rem(mod,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
    if (!fmpq_poly_is_zero((fmpq_poly_ptr)mod))
    {
      WerrorS("cannot divide");
    }
    fmpq_poly_clear(mod);
  }
  return (number)res;
}

Equal()

static BOOLEAN Equal ( number  a, number  b, const  coeffs  )

static

Definition at line 320 of file flintcf_Q.cc.

{
  return (fmpq_poly_equal((fmpq_poly_ptr)a,(fmpq_poly_ptr)b));
}

ExactDiv()

static number ExactDiv ( number  a, number  b, const  coeffs  )

static

Definition at line 99 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  if(fmpq_poly_is_zero((fmpq_poly_ptr)b))
  {
     WerrorS(nDivBy0);
  }
  else
    fmpq_poly_div(res,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
  return (number)res;
}

ExtGcd()

static number ExtGcd ( number  a, number  b, number *  s, number *  t, const  coeffs  )

static

Definition at line 386 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_init((fmpq_poly_ptr)*s);
  fmpq_poly_init((fmpq_poly_ptr)*t);
  fmpq_poly_xgcd(res,(fmpq_poly_ptr)*s,(fmpq_poly_ptr)*t,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
  return (number)res;
}

Farey()

static number Farey ( number  , number  , const  coeffs  )

static

Definition at line 434 of file flintcf_Q.cc.

{
  WerrorS("not yet: Farey");
  return NULL;
}

flintQ_InitChar()

BOOLEAN flintQ_InitChar	(	coeffs	cf,
		void *	infoStruct
	)

Definition at line 554 of file flintcf_Q.cc.

{
  char *pp=(char*)infoStruct;
  cf->cfCoeffName    = CoeffName;
  cf->nCoeffIsEqual  = CoeffIsEqual;
  cf->cfKillChar     = KillChar;
  cf->cfSetChar      = SetChar;
  cf->ch=0; //char 0
  cf->cfMult         = Mult;
  cf->cfSub          = Sub;
  cf->cfAdd          = Add;
  cf->cfDiv          = Div;
  cf->cfExactDiv     = ExactDiv; // ???
  cf->cfInit         =Init;
  cf->cfInitMPZ      =InitMPZ;
  cf->cfSize         = Size;
  cf->cfInt          = Int;
  cf->cfMPZ          = MPZ;
  cf->cfInpNeg       = Neg;
  cf->cfInvers       = Invers;
  cf->cfCopy         = Copy;
  cf->cfRePart       = Copy;
  // default: cf->cfImPart       = ndReturn0;
  cf->cfWriteLong    = WriteShort; //WriteLong;
  cf->cfWriteShort = WriteShort;
  cf->cfRead         = Read;
  cf->cfNormalize    = Normalize;
 
  //cf->cfDivComp=
  //cf->cfIsUnit=
  //cf->cfGetUnit=
  //cf->cfDivBy=
 
  cf->cfGreater=Greater;
  cf->cfEqual  =Equal;
  cf->cfIsZero =IsZero;
  cf->cfIsOne  =IsOne;
  cf->cfIsMOne =IsMOne;
  cf->cfGreaterZero=GreaterZero;
 
  cf->cfPower        = Power;
  cf->cfGetDenom     = GetDenom;
  cf->cfGetNumerator = GetNumerator;
  cf->cfGcd          = Gcd;
  cf->cfExtGcd         = ExtGcd;
  cf->cfLcm          = Lcm;
  cf->cfDelete       = Delete;
  cf->cfSetMap       = SetMap;
  // default: cf->cfInpMult
  // default: cf->cfInpAdd
  cf->cfFarey        =Farey;
  cf->cfChineseRemainder=ChineseRemainder;
  cf->cfParDeg = ParDeg;
  cf->cfParameter = Parameter;
  //  cf->cfClearContent = ClearContent;
  //  cf->cfClearDenominators = ClearDenominators;
  cf->convFactoryNSingN=ConvFactoryNSingN;
  cf->convSingNFactoryN=ConvSingNFactoryN;
  cf->cfWriteFd  = WriteFd;
  cf->cfReadFd = ReadFd;
#ifdef LDEBUG
  cf->cfDBTest       = DBTest;
#endif
 
  cf->iNumberOfParameters = 1;
  char **pn=(char**)omAlloc0(sizeof(char*));
  pn[0]=omStrDup(pp);
  cf->pParameterNames = (const char **)pn;
  cf->has_simple_Inverse= FALSE;
  cf->has_simple_Alloc= FALSE;
  cf->is_field=FALSE;
 
  return FALSE;
}

flintQInitCfByName()

coeffs flintQInitCfByName	(	char *	s,
		n_coeffType	n
	)

Definition at line 526 of file flintcf_Q.cc.

{
  const char start[]="flintQp[";
  const int start_len=strlen(start);
  if (strncmp(s,start,start_len)==0)
  {
    s+=start_len;
    char st[10];
    int l=sscanf(s,"%s",st);
    if (l==1)
    {
      while (st[strlen(st)-1]==']') st[strlen(st)-1]='\0';
      return nInitChar(n,(void*)st);
    }
  }
  return NULL;
}

Gcd()

static number Gcd ( number  a, number  b, const  coeffs  )

static

Definition at line 379 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_gcd(res,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
  return (number)res;
}

GetDenom()

static number GetDenom ( number &  n, const  coeffs  )

static

Definition at line 362 of file flintcf_Q.cc.

{
 fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
 fmpq_poly_init(res);
 fmpz_ptr den=fmpq_poly_denref((fmpq_poly_ptr)n);
 fmpq_poly_set_fmpz(res,den);
 return (number)res;
}

GetNumerator()

static number GetNumerator ( number &  n, const  coeffs  )

static

Definition at line 370 of file flintcf_Q.cc.

{
 fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
 fmpq_poly_init(res);
 fmpq_poly_set(res,(fmpq_poly_ptr)n);
 fmpz_ptr den=fmpq_poly_denref(res);
 fmpq_poly_scalar_mul_fmpz(res,res,den);
 return (number)res;
}

Greater()

static BOOLEAN Greater ( number  a, number  b, const  coeffs  )

static

Definition at line 316 of file flintcf_Q.cc.

{
  return (fmpq_poly_cmp((fmpq_poly_ptr)a,(fmpq_poly_ptr)b)>0);
}

GreaterZero()

static BOOLEAN GreaterZero ( number  , const  coeffs  )

static

Definition at line 349 of file flintcf_Q.cc.

{
  // does it have a leading sign?
  // no: 0 and 1 do not have, everything else is in (...)
  return TRUE;
}

Init()

static number Init ( long  i, const  coeffs  )

static

Definition at line 120 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_set_si(res,i);
  return (number)res;
}

InitMPZ()

static number InitMPZ ( mpz_t  i, const  coeffs  )

static

Definition at line 127 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_fit_length(res, 1);
  fmpz_set_mpz(res->coeffs, i);
  fmpz_one(res->den);
  _fmpq_poly_set_length(res, 1);
  _fmpq_poly_normalise(res);
  return (number)res;
}

Int()

static long Int ( number &  n, const  coeffs  )

static

Definition at line 142 of file flintcf_Q.cc.

{
  if (fmpq_poly_degree((fmpq_poly_ptr)n)==0)
  {
    fmpq_t m;
    fmpq_init(m);
    fmpq_poly_get_coeff_fmpq(m,(fmpq_poly_ptr)n,0);
    long nl=fmpz_get_si(fmpq_numref(m));
    if (fmpz_cmp_si(fmpq_numref(m),nl)!=0) nl=0;
    long dl=fmpz_get_si(fmpq_denref(m));
    if ((dl!=1)||(fmpz_cmp_si(fmpq_denref(m),dl)!=0)) nl=0;
    fmpq_clear(m);
    return nl;
  }
  return 0;
}

Invers()

static number Invers ( number  a, const  coeffs  )

static

Definition at line 180 of file flintcf_Q.cc.

{
  if(fmpq_poly_is_zero((fmpq_poly_ptr)a))
  {
    WerrorS(nDivBy0);
    return NULL;
  }
  if (fmpq_poly_degree((fmpq_poly_ptr)a)==0)
  {
    fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
    fmpq_poly_init(res);
    fmpq_poly_inv(res,(fmpq_poly_ptr)a);
    return (number)res;
  }
  else
  {
    WerrorS("not invertable");
    return NULL;
  }
}

IsMOne()

static BOOLEAN IsMOne ( number  k, const  coeffs  )

static

Definition at line 332 of file flintcf_Q.cc.

{
  if (fmpq_poly_length((fmpq_poly_ptr)k)>0) return FALSE;
  fmpq_poly_canonicalise((fmpq_poly_ptr)k);
  fmpq_t m;
  fmpq_init(m);
  fmpq_poly_get_coeff_fmpq(m,(fmpq_poly_ptr)k,0);
  BOOLEAN result=TRUE;
  if (fmpz_cmp_si(fmpq_numref(m),(long)-1)!=0) result=FALSE;
  else
  {
    int dl=(int)fmpz_get_si(fmpq_denref(m));
    if ((dl!=1)||(fmpz_cmp_si(fmpq_denref(m),(long)dl)!=0)) result=FALSE;
  }
  fmpq_clear(m);
  return (result);
}

IsOne()

static BOOLEAN IsOne ( number  a, const  coeffs  )

static

Definition at line 328 of file flintcf_Q.cc.

{
  return fmpq_poly_is_one((fmpq_poly_ptr)a);
}

IsZero()

static BOOLEAN IsZero ( number  a, const  coeffs  )

static

Definition at line 324 of file flintcf_Q.cc.

{
  return fmpq_poly_is_zero((fmpq_poly_ptr)a);
}

KillChar()

static void KillChar ( coeffs  cf )

static

Definition at line 549 of file flintcf_Q.cc.

{
  omFree((ADDRESS)(cf->pParameterNames[0]));
  omFreeSize(cf->pParameterNames,sizeof(char*));
}

Lcm()

static number Lcm ( number  , number  , const  coeffs  )

static

Definition at line 395 of file flintcf_Q.cc.

{
  WerrorS("not yet: Lcm");
  return NULL;
}

MPZ()

static void MPZ ( mpz_t  result, number &  n, const  coeffs  )

static

Definition at line 158 of file flintcf_Q.cc.

{
  mpz_init(result);
  if (fmpq_poly_degree((fmpq_poly_ptr)n)==0)
  {
    fmpq_t m;
    fmpq_init(m);
    fmpq_poly_get_coeff_fmpq(m,(fmpq_poly_ptr)n,0);
    mpz_t den;
    mpz_init(den);
    fmpq_get_mpz_frac(result,den,m);
    int dl=(int)mpz_get_si(den);
    if ((dl!=1)||(mpz_cmp_si(den,(long)dl)!=0)) mpz_set_ui(result,0);
    mpz_clear(den);
    fmpq_clear(m);
  }
}

Mult()

static number Mult ( number  a, number  b, const  coeffs  )

static

Definition at line 56 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_mul(res,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
  return (number)res;
}

Neg()

static number Neg ( number  a, const  coeffs  )

static

Definition at line 175 of file flintcf_Q.cc.

{
  fmpq_poly_neg((fmpq_poly_ptr)a,(fmpq_poly_ptr)a);
  return a;
}

nlEatLong()

static char * nlEatLong ( char *  s, mpz_ptr  i  )

static

Definition at line 29 of file flintcf_Q.cc.

{
  const char * start=s;
 
  while (*s >= '0' && *s <= '9') s++;
  if (*s=='\0')
  {
    mpz_set_str(i,start,10);
  }
  else
  {
    char c=*s;
    *s='\0';
    mpz_set_str(i,start,10);
    *s=c;
  }
  return s;
}

Normalize()

static void Normalize ( number &  a, const  coeffs  )

static

Definition at line 312 of file flintcf_Q.cc.

{
  fmpq_poly_canonicalise((fmpq_poly_ptr)a);
}

Parameter()

static number Parameter ( const int  , const  coeffs  )

static

Definition at line 448 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_set_coeff_si(res,1,1);
  return (number)res;
}

ParDeg()

static int ParDeg ( number  x, const  coeffs  )

static

Definition at line 444 of file flintcf_Q.cc.

{
  return fmpq_poly_degree((fmpq_poly_ptr)x);
}

Power()

static void Power ( number  a, int  i, number *  result, const  coeffs  )

static

Definition at line 355 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  *result=(number)res;
  fmpq_poly_pow((fmpq_poly_ptr)(*result),(fmpq_poly_ptr)a,i);
}

Read()

static const char * Read ( const char *  st, number *  a, const coeffs  r  )

static

Definition at line 265 of file flintcf_Q.cc.

{
// we only read "monomials" (i.e. [-][digits][parameter]),
// everythings else (+,*,^,()) is left to the singular interpreter
  char *s=(char *)st;
  *a=(number)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init((fmpq_poly_ptr)(*a));
  BOOLEAN neg=FALSE;
  if (*s=='-') { neg=TRUE; s++;}
  if (isdigit(*s))
  {
    mpz_t z;
    mpz_init(z);
    fmpz_t z1;
    fmpz_init(z1);
    s=nlEatLong((char *)s, z);
    fmpz_set_mpz(z1,z);
    fmpq_poly_set_fmpz((fmpq_poly_ptr)(*a),z1);
    if (*s == '/')
    {
      s++;
      s=nlEatLong((char *)s, z);
      fmpz_set_mpz(z1,z);
      fmpq_poly_scalar_div_fmpz((fmpq_poly_ptr)(*a),(fmpq_poly_ptr)(*a),z1);
    }
    fmpz_clear(z1);
    mpz_clear(z);
  }
  else if(strncmp(s,r->pParameterNames[0],strlen(r->pParameterNames[0]))==0)
  {
    fmpq_poly_set_coeff_si((fmpq_poly_ptr)(*a),1,1);
    s+=strlen(r->pParameterNames[0]);
    if(isdigit(*s))
    {
      int i=1;
      s=nEati(s,&i,0);
      if (i!=1)
      {
        fmpq_poly_set_coeff_si((fmpq_poly_ptr)(*a),1,0);
        fmpq_poly_set_coeff_si((fmpq_poly_ptr)(*a),i,1);
      }
    }
  }
  if (neg)
    fmpq_poly_neg((fmpq_poly_ptr)(*a),(fmpq_poly_ptr)(*a));
  return s;
}

ReadFd()

static number ReadFd ( const ssiInfo *  d, const  coeffs  )

static

Definition at line 479 of file flintcf_Q.cc.

{
  // format: len a_len .. a_0
  fmpq_poly_ptr aa=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(aa);
  int l=s_readint(d->f_read);
  mpz_t tmp;
  mpz_init(tmp);
  fmpq_t m;
  fmpq_init(m);
  fmpz_t num,den;
  fmpz_init(num);
  fmpz_init(den);
  for (int i=l;i>=0;i--)
  {
    s_readmpz_base (d->f_read,tmp, SSI_BASE);
    fmpz_set_mpz(num, tmp);
    s_readmpz_base (d->f_read,tmp, SSI_BASE);
    fmpz_set_mpz(den, tmp);
    fmpq_set_fmpz_frac(m,num,den);
    fmpq_poly_set_coeff_fmpq(aa,i,m);
  }
  mpz_clear(tmp);
  fmpz_clear(den);
  fmpz_clear(num);
  fmpq_clear(m);
  return (number)aa;
}

SetChar()

static void SetChar ( const  coeffs )

static

Definition at line 52 of file flintcf_Q.cc.

{
  // dummy
}

SetMap()

static nMapFunc SetMap ( const  coeffs, const  coeffs  )

static

Definition at line 409 of file flintcf_Q.cc.

{
  WerrorS("not yet: SetMap");
  return NULL;
}

Size()

static int Size ( number  n, const  coeffs  )

static

Definition at line 138 of file flintcf_Q.cc.

{
  return fmpq_poly_degree((fmpq_poly_ptr)n);
}

Sub()

static number Sub ( number  a, number  b, const  coeffs  )

static

Definition at line 63 of file flintcf_Q.cc.

{
  fmpq_poly_ptr res=(fmpq_poly_ptr)omAlloc(sizeof(fmpq_poly_t));
  fmpq_poly_init(res);
  fmpq_poly_sub(res,(fmpq_poly_ptr)a,(fmpq_poly_ptr)b);
  return (number)res;
}

WriteFd()

static void WriteFd ( number  a, const ssiInfo *  d, const  coeffs  )

static

Definition at line 455 of file flintcf_Q.cc.

{
  // format: len a_len(num den) .. a_0
  fmpq_poly_ptr aa=(fmpq_poly_ptr)a;
  int l=fmpq_poly_length(aa);
  fprintf(d->f_write,"%d ",l);
  fmpq_t m;
  fmpq_init(m);
  mpz_t num,den;
  mpz_init(num);
  mpz_init(den);
  for(int i=l; i>=0; i--)
  {
    fmpq_poly_get_coeff_fmpq(m,(fmpq_poly_ptr)a,i);
    fmpq_get_mpz_frac(num,den,m);
    mpz_out_str (d->f_write,SSI_BASE, num);
    fputc(' ',d->f_write);
    mpz_out_str (d->f_write,SSI_BASE, den);
    fputc(' ',d->f_write);
  }
  mpz_clear(den);
  mpz_clear(num);
  fmpq_clear(m);
}

WriteShort()

static void WriteShort ( number  a, const coeffs  r  )

static

Definition at line 218 of file flintcf_Q.cc.

{
  //fmpq_poly_print_pretty((fmpq_poly_ptr)a,r->pParameterNames[0]);
  if (IsOne(a,r)) StringAppendS("1");
  else if (IsZero(a,r)) StringAppendS("0");
  else
  {
  StringAppendS("(");
  fmpq_t m;
  fmpq_init(m);
  BOOLEAN need_plus=FALSE;
  for(int i=fmpq_poly_length((fmpq_poly_ptr)a);i>=0;i--)
  {
    fmpq_poly_get_coeff_fmpq(m,(fmpq_poly_ptr)a,i);
    if (!fmpq_is_zero(m))
    {
      if (need_plus && (fmpq_cmp_ui(m,0)>0))
        StringAppendS("+");
      need_plus=TRUE;
      int l=fmpz_sizeinbase(fmpq_numref(m),10);
      l=si_max(l,(int)fmpz_sizeinbase(fmpq_denref(m),10));
      l+=2;
      char *s=(char*)omAlloc(l);
      char *z=fmpz_get_str(s,10,fmpq_numref(m));
      if ((i==0)
      ||(fmpz_cmp_si(fmpq_numref(m),1)!=0)
      ||(fmpz_cmp_si(fmpq_denref(m),1)!=0))
      {
        StringAppendS(z);
        if (fmpz_cmp_si(fmpq_denref(m),1)!=0)
        {
          StringAppendS("/");
          z=fmpz_get_str(s,10,fmpq_denref(m));
          StringAppendS(z);
        }
        if (i!=0) StringAppendS("*");
      }
      if (i>1)
        StringAppend("%s^%d",r->pParameterNames[0],i);
      else if (i==1)
        StringAppend("%s",r->pParameterNames[0]);
    }
  }
  fmpq_clear(m);
  StringAppendS(")");
  }
}