flintcf_Q.cc File Reference

#include <ctype.h>
#include "misc/auxiliary.h"
#include <flint/flint.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 22 of file flintcf_Q.cc.

fmpz_ptr

typedef fmpz* fmpz_ptr

Definition at line 23 of file flintcf_Q.cc.

Function Documentation

Add()

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

static

Definition at line 68 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 450 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 46 of file flintcf_Q.cc.

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

CoeffName()

char * CoeffName

(

const coeffs

r

)

Definition at line 526 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 516 of file flintcf_Q.cc.

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

ConvSingNFactoryN()

static CanonicalForm ConvSingNFactoryN ( number  , BOOLEAN  , const  coeffs  )

static

Definition at line 521 of file flintcf_Q.cc.

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

Copy()

static number Copy ( number  a, const  coeffs  )

static

Definition at line 202 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 551 of file flintcf_Q.cc.

{
  return TRUE;
}

Delete()

static void Delete ( number *  a, const  coeffs  )

static

Definition at line 411 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 75 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 324 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 97 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 397 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 445 of file flintcf_Q.cc.

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

flintQ_InitChar()

BOOLEAN flintQ_InitChar	(	coeffs	cf,
		void *	infoStruct
	)

Definition at line 561 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 533 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 390 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 373 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 381 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 320 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 360 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 118 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 125 of file flintcf_Q.cc.

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

Int()

static long Int ( number &  n, const  coeffs  )

static

Definition at line 136 of file flintcf_Q.cc.

{
  if (fmpq_poly_degree((fmpq_poly_ptr)n)==0)
  {
    mpq_t m;
    mpq_init(m);
    fmpq_poly_get_coeff_mpq(m,(fmpq_poly_ptr)n,0);
    mpz_t num,den;
    mpz_init(num);
    mpz_init(den);
    mpq_get_num(num,m);
    mpq_get_den(den,m);
    long nl=mpz_get_si(num);
    if (mpz_cmp_si(num,nl)!=0) nl=0;
    long dl=mpz_get_si(den);
    if ((dl!=1)||(mpz_cmp_si(den,dl)!=0)) nl=0;
    mpz_clear(num);
    mpz_clear(den);
    mpq_clear(m);
    return nl;
  }
  return 0;
}

Invers()

static number Invers ( number  a, const  coeffs  )

static

Definition at line 182 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 336 of file flintcf_Q.cc.

{
  if (fmpq_poly_length((fmpq_poly_ptr)k)>0) return FALSE;
  fmpq_poly_canonicalise((fmpq_poly_ptr)k);
  mpq_t m;
  mpq_init(m);
  fmpq_poly_get_coeff_mpq(m,(fmpq_poly_ptr)k,0);
  mpz_t num,den;
  mpz_init(num);
  mpq_get_num(num,m);
  BOOLEAN result=TRUE;
  if (mpz_cmp_si(num,(long)-1)!=0) result=FALSE;
  else
  {
    mpz_init(den);
    mpq_get_den(den,m);
    int dl=(int)mpz_get_si(den);
    if ((dl!=1)||(mpz_cmp_si(den,(long)dl)!=0)) result=FALSE;
    mpz_clear(den);
  }
  mpz_clear(num);
  mpq_clear(m);
  return (result);
}

IsOne()

static BOOLEAN IsOne ( number  a, const  coeffs  )

static

Definition at line 332 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 328 of file flintcf_Q.cc.

{
  return fmpq_poly_is_zero((fmpq_poly_ptr)a);
}

KillChar()

static void KillChar ( coeffs  cf )

static

Definition at line 556 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 406 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 159 of file flintcf_Q.cc.

{
  mpz_init(result);
  if (fmpq_poly_degree((fmpq_poly_ptr)n)==0)
  {
    mpq_t m;
    mpq_init(m);
    fmpq_poly_get_coeff_mpq(m,(fmpq_poly_ptr)n,0);
    mpz_t den;
    mpz_init(den);
    mpq_get_num(result,m);
    mpq_get_den(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);
    mpq_clear(m);
  }
}

Mult()

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

static

Definition at line 54 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 177 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 27 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 316 of file flintcf_Q.cc.

{
  fmpq_poly_canonicalise((fmpq_poly_ptr)a);
}

Parameter()

static number Parameter ( const int  , const  coeffs  )

static

Definition at line 459 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 455 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 366 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 274 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);
    s=nlEatLong((char *)s, z);
    fmpq_poly_set_mpz((fmpq_poly_ptr)(*a),z);
    if (*s == '/')
    {
      s++;
      s=nlEatLong((char *)s, z);
      fmpq_poly_scalar_div_mpz((fmpq_poly_ptr)(*a),(fmpq_poly_ptr)(*a),z);
    }
    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 491 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 nm;
  mpz_init(nm);
  mpq_t m;
  mpq_init(m);
  for (int i=l;i>=0;i--)
  {
 
    s_readmpz_base (d->f_read,nm, SSI_BASE);
    mpq_set_num(m,nm);
    s_readmpz_base (d->f_read,nm, SSI_BASE);
    mpq_set_den(m,nm);
    fmpq_poly_set_coeff_mpq(aa,i,m);
  }
  mpz_clear(nm);
  mpq_clear(m);
  return (number)aa;
}

SetChar()

static void SetChar ( const  coeffs )

static

Definition at line 50 of file flintcf_Q.cc.

{
  // dummy
}

SetMap()

static nMapFunc SetMap ( const  coeffs, const  coeffs  )

static

Definition at line 420 of file flintcf_Q.cc.

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

Size()

static int Size ( number  n, const  coeffs  )

static

Definition at line 132 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 61 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 466 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);
  mpq_t m;
  mpq_init(m);
  mpz_t num,den;
  mpz_init(num);
  mpz_init(den);
  for(int i=l; i>=0; i--)
  {
    fmpq_poly_get_coeff_mpq(m,(fmpq_poly_ptr)a,i);
    mpq_get_num(num,m);
    mpq_get_den(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);
  mpq_clear(m);
}

WriteShort()

static void WriteShort ( number  a, const coeffs  r  )

static

Definition at line 220 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("(");
  mpq_t m;
  mpq_init(m);
  mpz_t num,den;
  mpz_init(num);
  mpz_init(den);
  BOOLEAN need_plus=FALSE;
  for(int i=fmpq_poly_length((fmpq_poly_ptr)a);i>=0;i--)
  {
    fmpq_poly_get_coeff_mpq(m,(fmpq_poly_ptr)a,i);
    mpq_get_num(num,m);
    mpq_get_den(den,m);
    if (mpz_sgn1(num)!=0)
    {
      if (need_plus && (mpz_sgn1(num)>0))
        StringAppendS("+");
      need_plus=TRUE;
      int l=mpz_sizeinbase(num,10);
      l=si_max(l,(int)mpz_sizeinbase(den,10));
      l+=2;
      char *s=(char*)omAlloc(l);
      char *z=mpz_get_str(s,10,num);
      if ((i==0)
      ||(mpz_cmp_si(num,1)!=0)
      ||(mpz_cmp_si(den,1)!=0))
      {
        StringAppendS(z);
        if (mpz_cmp_si(den,1)!=0)
        {
          StringAppendS("/");
          z=mpz_get_str(s,10,den);
          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]);
    }
  }
  mpz_clear(den);
  mpz_clear(num);
  mpq_clear(m);
  StringAppendS(")");
  }
}