gnumpfl.cc File Reference

#include "misc/auxiliary.h"
#include "reporter/reporter.h"
#include "coeffs/coeffs.h"
#include "coeffs/numbers.h"
#include "coeffs/mpr_complex.h"
#include "coeffs/longrat.h"
#include "coeffs/shortfl.h"
#include "coeffs/gnumpfl.h"
#include "coeffs/gnumpc.h"
#include "coeffs/modulop.h"

Go to the source code of this file.

Data Structures
union	nf

Functions
static number	ngfInit (long i, const coeffs r)
static long	ngfInt (number &i, const coeffs r)
static BOOLEAN	ngfIsZero (number a, const coeffs r)
static int	ngfSize (number n, const coeffs r)
static void	ngfDelete (number *a, const coeffs r)
static number	ngfCopy (number a, const coeffs r)
static number	ngfNeg (number a, const coeffs r)
static number	ngfInvers (number a, const coeffs r)
static number	ngfAdd (number a, number b, const coeffs R)
static void	ngfInpAdd (number &a, number b, const coeffs R)
static number	ngfSub (number a, number b, const coeffs R)
static number	ngfMult (number a, number b, const coeffs R)
static void	ngfInpMult (number &a, number b, const coeffs R)
static number	ngfDiv (number a, number b, const coeffs r)
static number	ngfPower (number x, int exp, const coeffs r)
static void	ngfPower (number x, int exp, number *u, const coeffs r)
static BOOLEAN	ngfGreaterZero (number a, const coeffs r)
static BOOLEAN	ngfGreater (number a, number b, const coeffs r)
static BOOLEAN	ngfEqual (number a, number b, const coeffs r)
static BOOLEAN	ngfIsOne (number a, const coeffs r)
static BOOLEAN	ngfIsMOne (number a, const coeffs r)
static char *	ngfEatFloatNExp (char *s)
const char *	ngfRead (const char *start, number *a, const coeffs r)
static void	ngfWrite (number a, const coeffs r)
static BOOLEAN	ngfCoeffIsEqual (const coeffs r, n_coeffType n, void *parameter)
static void	ngfSetChar (const coeffs r)
static char *	ngfCoeffName (const coeffs r)
static number	ngfMapQ (number from, const coeffs src, const coeffs dst)
static number	ngfMapZ (number from, const coeffs aRing, const coeffs r)
static number	ngfMapR (number from, const coeffs src, const coeffs dst)
static number	ngfMapP (number from, const coeffs src, const coeffs dst)
static number	ngfMapC (number from, const coeffs src, const coeffs dst)
static number	ngfInitMPZ (mpz_t m, const coeffs)
static nMapFunc	ngfSetMap (const coeffs src, const coeffs dst)
BOOLEAN	ngfInitChar (coeffs n, void *parameter)
	Initialize r. More...

Function Documentation

ngfAdd()

static number ngfAdd ( number  a, number  b, const coeffs  R  )

static

Definition at line 151 of file gnumpfl.cc.

{
  assume( getCoeffType(R) == n_long_R );
 
  gmp_float* r= new gmp_float( (*(gmp_float*)a) + (*(gmp_float*)b) );
  return (number)r;
}

ngfCoeffIsEqual()

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

static

Definition at line 401 of file gnumpfl.cc.

{
  if (n==n_long_R)
  {
    LongComplexInfo* p = (LongComplexInfo *)(parameter);
    if ((p!=NULL)
    && (p->float_len == r->float_len)
    && (p->float_len2 == r->float_len2))
      return TRUE;
  }
  return FALSE;
}

ngfCoeffName()

static char * ngfCoeffName ( const coeffs  r )

static

Definition at line 419 of file gnumpfl.cc.

{
  STATIC_VAR char ngfCoeffName_buf[30];
  snprintf(ngfCoeffName_buf,30,"Float(%d,%d)",r->float_len,r->float_len2);
  return ngfCoeffName_buf;
}

ngfCopy()

static number ngfCopy ( number  a, const coeffs  r  )

static

Definition at line 94 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  gmp_float* b= new gmp_float( *(gmp_float*)a );
  return (number)b;
}

ngfDelete()

static void ngfDelete ( number *  a, const coeffs  r  )

static

Definition at line 80 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  if ( *a != NULL )
  {
    delete *(gmp_float**)a;
    *a=NULL;
  }
}

ngfDiv()

static number ngfDiv ( number  a, number  b, const coeffs  r  )

static

Definition at line 198 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
  
  gmp_float* f;
  if ( ((gmp_float*)b)->isZero() )
  {
    // a/0 = error
    WerrorS(nDivBy0);
    f= new gmp_float( 0 );
  }
  else
  {
    f= new gmp_float( (*(gmp_float*)a) / (*(gmp_float*)b) );
  }
  return (number)f;
}

ngfEatFloatNExp()

static char * ngfEatFloatNExp ( char *  s )

static

Definition at line 295 of file gnumpfl.cc.

{
  char *start= s;
 
  // eat floats (mantissa) like:
  //   0.394394993, 102.203003008,  .300303032, pssibly starting with -
  if (*s == '-') s++;
  while ((*s >= '0' && *s <= '9')||(*s == '.')) s++;
 
  // eat the exponent, starts with 'e' followed by '+', '-'
  // and digits, like:
  //   e-202, e+393, accept also E7
  if ( (s != start) && ((*s == 'e')||(*s=='E')))
  {
    if (*s=='E') *s='e';
    s++; // skip 'e'/'E'
    if ((*s == '+') || (*s == '-')) s++;
    while ((*s >= '0' && *s <= '9')) s++;
  }
 
  return s;
}

ngfEqual()

static BOOLEAN ngfEqual ( number  a, number  b, const coeffs  r  )

static

Definition at line 268 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  return ( (*(gmp_float*)a) == (*(gmp_float*)b) );
}

ngfGreater()

static BOOLEAN ngfGreater ( number  a, number  b, const coeffs  r  )

static

Definition at line 258 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  return ( (*(gmp_float*)a) > (*(gmp_float*)b) );
}

ngfGreaterZero()

static BOOLEAN ngfGreaterZero ( number  a, const coeffs  r  )

static

Definition at line 248 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  return (((gmp_float*)a)->sign() > 0);
}

ngfInit()

static number ngfInit ( long  i, const coeffs  r  )

static

Definition at line 37 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  gmp_float* n= new gmp_float( (double)i );
  return (number)n;
}

ngfInitChar()

BOOLEAN ngfInitChar	(	coeffs	n,
		void *	parameter
	)

Initialize r.

Definition at line 511 of file gnumpfl.cc.

{
  assume( getCoeffType(n) == n_long_R );
 
  n->is_field=TRUE;
  n->is_domain=TRUE;
  n->rep=n_rep_gmp_float;
 
  //n->cfKillChar = ndKillChar; /* dummy */
 
  n->cfSetChar = ngfSetChar;
  n->ch = 0;
  n->cfCoeffName=ngfCoeffName;
 
  n->cfDelete  = ngfDelete;
  //n->cfNormalize=ndNormalize;
  n->cfInit    = ngfInit;
  n->cfInitMPZ = ngfInitMPZ;
  n->cfInt     = ngfInt;
  n->cfAdd     = ngfAdd;
  n->cfInpAdd  = ngfInpAdd;
  n->cfSub     = ngfSub;
  n->cfMult    = ngfMult;
  n->cfInpMult = ngfInpMult;
  n->cfDiv     = ngfDiv;
  n->cfExactDiv= ngfDiv;
  n->cfInpNeg  = ngfNeg;
  n->cfInvers  = ngfInvers;
  n->cfCopy    = ngfCopy;
  n->cfGreater = ngfGreater;
  n->cfEqual   = ngfEqual;
  n->cfIsZero  = ngfIsZero;
  n->cfIsOne   = ngfIsOne;
  n->cfIsMOne  = ngfIsMOne;
  n->cfGreaterZero = ngfGreaterZero;
  n->cfWriteLong  = ngfWrite;
  n->cfRead    = ngfRead;
  n->cfPower   = ngfPower;
  n->cfSetMap = ngfSetMap;
  n->cfSize   = ngfSize;
#ifdef LDEBUG
  //n->cfDBTest  = ndDBTest; // not yet implemented: ngfDBTest
#endif
 
  n->nCoeffIsEqual = ngfCoeffIsEqual;
 
  if( parameter != NULL)
  {
    LongComplexInfo* p = (LongComplexInfo*)parameter;
 
    n->float_len = p->float_len;
    n->float_len2 = p->float_len2;
  } else // default values, just for testing!
  {
    n->float_len = SHORT_REAL_LENGTH;
    n->float_len2 = SHORT_REAL_LENGTH;
  }
 
  assume( n->float_len2 >= SHORT_REAL_LENGTH );
 
  assume( n_NumberOfParameters(n) == 0 );
  assume( n_ParameterNames(n) == NULL );
 
  return FALSE;
}

ngfInitMPZ()

static number ngfInitMPZ ( mpz_t  m, const  coeffs  )

static

Definition at line 470 of file gnumpfl.cc.

{
  gmp_float *res=new gmp_float(m);
  return (number)res;
}

ngfInpAdd()

static void ngfInpAdd ( number &  a, number  b, const coeffs  R  )

static

Definition at line 159 of file gnumpfl.cc.

{
  assume( getCoeffType(R) == n_long_R );
 
  (*(gmp_float*)a) += (*(gmp_float*)b);
}

ngfInpMult()

static void ngfInpMult ( number &  a, number  b, const coeffs  R  )

static

Definition at line 188 of file gnumpfl.cc.

{
  assume( getCoeffType(R) == n_long_R );
 
  (*(gmp_float*)a) *= (*(gmp_float*)b);
}

ngfInt()

static long ngfInt ( number &  i, const coeffs  r  )

static

Definition at line 48 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  double d=(double)*(gmp_float*)i;
  if (d<0.0)
    return (long)(d-0.5);
  else
    return (long)(d+0.5);
}

ngfInvers()

static number ngfInvers ( number  a, const coeffs  r  )

static

Definition at line 131 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  gmp_float* f= NULL;
  if (((gmp_float*)a)->isZero() )
  {
    WerrorS(nDivBy0);
    f= new gmp_float( 0 );
  }
  else
  {
    f= new gmp_float( gmp_float(1) / (*(gmp_float*)a) );
  }
  return (number)f;
}

ngfIsMOne()

static BOOLEAN ngfIsMOne ( number  a, const coeffs  r  )

static

Definition at line 288 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  return ((gmp_float*)a)->isMOne();
}

ngfIsOne()

static BOOLEAN ngfIsOne ( number  a, const coeffs  r  )

static

Definition at line 278 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  return ((gmp_float*)a)->isOne();
}

ngfIsZero()

static BOOLEAN ngfIsZero ( number  a, const coeffs  r  )

static

Definition at line 59 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  return ( ((gmp_float*)a)->isZero() );
}

ngfMapC()

static number ngfMapC ( number  from, const coeffs  src, const coeffs  dst  )

static

Definition at line 461 of file gnumpfl.cc.

{
  assume( getCoeffType(dst) == n_long_R );
  assume( getCoeffType(src) ==  n_long_C );
 
  gmp_float *res=new gmp_float(((gmp_complex*)from)->real());
  return (number)res;
}

ngfMapP()

static number ngfMapP ( number  from, const coeffs  src, const coeffs  dst  )

static

Definition at line 453 of file gnumpfl.cc.

{
  assume( getCoeffType(dst) == n_long_R );
  assume( getCoeffType(src) ==  n_Zp );
 
  return ngfInit(npInt(from,src), dst); // FIXME? TODO? // extern int     npInt         (number &n, const coeffs r);
}

ngfMapQ()

static number ngfMapQ ( number  from, const coeffs  src, const coeffs  dst  )

static

Definition at line 426 of file gnumpfl.cc.

{
  assume( getCoeffType(dst) == n_long_R );
  assume( src->rep == n_rep_gap_rat );
 
  gmp_float *res=new gmp_float(numberFieldToFloat(from,QTOF));
  return (number)res;
}

ngfMapR()

static number ngfMapR ( number  from, const coeffs  src, const coeffs  dst  )

static

Definition at line 444 of file gnumpfl.cc.

{
  assume( getCoeffType(dst) == n_long_R );
  assume( getCoeffType(src) == n_R );
 
  gmp_float *res=new gmp_float((double)nf(from).F());
  return (number)res;
}

ngfMapZ()

static number ngfMapZ ( number  from, const coeffs  aRing, const coeffs  r  )

static

Definition at line 435 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
  assume( aRing->rep == n_rep_gmp);
 
  gmp_float *res=new gmp_float((mpz_ptr)from);
  return (number)res;
}

ngfMult()

static number ngfMult ( number  a, number  b, const coeffs  R  )

static

Definition at line 180 of file gnumpfl.cc.

{
  assume( getCoeffType(R) == n_long_R );
 
  gmp_float* r= new gmp_float( (*(gmp_float*)a) * (*(gmp_float*)b) );
  return (number)r;
}

ngfNeg()

static number ngfNeg ( number  a, const coeffs  r  )

static

Definition at line 120 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  *(gmp_float*)a= -(*(gmp_float*)a);
  return (number)a;
}

ngfPower() [1/2]

static number ngfPower ( number  x, int  exp, const coeffs  r  )

static

Definition at line 219 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  if ( exp == 0 )
  {
    gmp_float* n = new gmp_float(1);
    return (number)n;
  }
  else if ( ngfIsZero(x, r) ) // 0^e, e>0
  {
    return ngfInit(0, r);
  }
  else if ( exp == 1 )
  {
    return ngfCopy(x,r);
  }
  return (number) ( new gmp_float( (*(gmp_float*)x)^exp ) );
}

ngfPower() [2/2]

static void ngfPower ( number  x, int  exp, number *  u, const coeffs  r  )

static

Definition at line 240 of file gnumpfl.cc.

{
  *u = ngfPower(x, exp, r);
}

ngfRead()

const char * ngfRead	(	const char *	start,
		number *	a,
		const coeffs	r
	)

Definition at line 324 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R or getCoeffType(r) == n_long_C);
 
  char *s= (char *)start;
 
  //Print("%s\n",s);
 
  s= ngfEatFloatNExp( s );
 
  if (*s=='\0')  // 0
  {
    if ( *(gmp_float**)a == NULL ) (*(gmp_float**)a)= new gmp_float();
    (*(gmp_float**)a)->setFromStr(start);
  }
  else if (s==start)  // 1
  {
    if ( *(gmp_float**)a != NULL )  delete (*(gmp_float**)a);
    (*(gmp_float**)a)= new gmp_float(1);
  }
  else
  {
    gmp_float divisor(1.0);
    char *start2=s;
    if ( *s == '/' )
    {
      s++;
      s= ngfEatFloatNExp( (char *)s );
      if (s!= start2+1)
      {
        char tmp_c=*s;
        *s='\0';
        divisor.setFromStr(start2+1);
        *s=tmp_c;
      }
      else
      {
        Werror("wrong long real format: %s",start2);
      }
    }
    char c=*start2;
    *start2='\0';
    if ( *(gmp_float**)a == NULL ) (*(gmp_float**)a)= new gmp_float();
    (*(gmp_float**)a)->setFromStr(start);
    *start2=c;
    if (divisor.isZero())
    {
      WerrorS(nDivBy0);
    }
    else
      (**(gmp_float**)a) /= divisor;
  }
 
  return s;
}

ngfSetChar()

static void ngfSetChar ( const coeffs  r )

static

Definition at line 414 of file gnumpfl.cc.

{
  setGMPFloatDigits(r->float_len, r->float_len2);
}

ngfSetMap()

static nMapFunc ngfSetMap ( const coeffs  src, const coeffs  dst  )

static

Definition at line 476 of file gnumpfl.cc.

{
  assume( getCoeffType(dst) == n_long_R );
 
  if (src->rep==n_rep_gap_rat) /*Q, Z*/
  {
    return ngfMapQ;
  }
  if (src->rep==n_rep_gap_gmp) /*Q, bigint*/
  {
    return ngfMapQ;
  }
  if (src->rep==n_rep_gmp) /* Z*/
  {
    return ngfMapZ;
  }
  if ((src->rep==n_rep_gmp_float) && nCoeff_is_long_R(src))
  {
    return ndCopyMap; //ngfCopyMap;
  }
  if ((src->rep==n_rep_float) && nCoeff_is_R(src))
  {
    return ngfMapR;
  }
  if ((src->rep==n_rep_gmp_complex) && nCoeff_is_long_C(src))
  {
    return ngfMapC;
  }
  if ((src->rep==n_rep_int) && nCoeff_is_Zp(src))
  {
    return ngfMapP;
  }
  return NULL;
}

ngfSize()

static int ngfSize ( number  n, const coeffs  r  )

static

Definition at line 66 of file gnumpfl.cc.

{
  long i = ngfInt(n, r);
  /* basically return the largest integer in n;
     only if this happens to be zero although n != 0,
     return 1;
     (this code ensures that zero has the size zero) */
  if ((i == 0) && (ngfIsZero(n,r) == FALSE)) i = 1;
  return ABS(i);
}

ngfSub()

static number ngfSub ( number  a, number  b, const coeffs  R  )

static

Definition at line 169 of file gnumpfl.cc.

{
  assume( getCoeffType(R) == n_long_R );
 
  gmp_float* r= new gmp_float( (*(gmp_float*)a) - (*(gmp_float*)b) );
  return (number)r;
}

ngfWrite()

static void ngfWrite ( number  a, const coeffs  r  )

static

Definition at line 383 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  char *out;
  if ( a != NULL )
  {
    out= floatToStr(*(gmp_float*)a, r->float_len);
    StringAppendS(out);
    //omFreeSize((void *)out, (strlen(out)+1)* sizeof(char) );
    omFree( (void *)out );
  }
  else
  {
    StringAppendS("0");
  }
}