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
const char *	ngfRead (const char *s, number *a, const coeffs r)
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 number	ngfSub (number a, number b, const coeffs R)
static number	ngfMult (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)
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 153 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 389 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 407 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 96 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 82 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 186 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 283 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 256 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 246 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 236 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 39 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 499 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->cfSub     = ngfSub;
  n->cfMult    = ngfMult;
  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;
#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 458 of file gnumpfl.cc.

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

ngfInt()

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

static

Definition at line 50 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 133 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 276 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 266 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 61 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 449 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 441 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 414 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 432 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 423 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 175 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 122 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 207 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 228 of file gnumpfl.cc.

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

ngfRead()

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

Definition at line 312 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 402 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 464 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 68 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 164 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 371 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");
  }
}