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 *	ngfCoeffString (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 nMapFunc	ngfSetMap (const coeffs src, const coeffs dst)
static void	ngfCoeffWrite (const coeffs r, BOOLEAN)
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 383 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 408 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;

ngfCoeffString()

static char* ngfCoeffString ( const coeffs  r )

static

Definition at line 401 of file gnumpfl.cc.

{
  char *s=(char*)omAlloc(30);
  snprintf(s,30,"Float(%d,%d)",r->float_len,r->float_len2);
  return s;

ngfCoeffWrite()

static void ngfCoeffWrite ( const coeffs  r, BOOLEAN    )

static

Definition at line 501 of file gnumpfl.cc.

{
  Print("Float(%d,%d)", r->float_len,r->float_len2);  /* long R */

ngfCopy()

static number ngfCopy ( number  a, const coeffs  r  )

static

Definition at line 95 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 81 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 184 of file gnumpfl.cc.

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

ngfEatFloatNExp()

static char* ngfEatFloatNExp ( char *  s )

static

Definition at line 277 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 250 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 240 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 230 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 38 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 506 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->cfCoeffString=ngfCoeffString;
  n->cfCoeffName=ngfCoeffName;
 
  n->cfDelete  = ngfDelete;
  //n->cfNormalize=ndNormalize;
  n->cfInit   = ngfInit;
  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;
  n->cfCoeffWrite = ngfCoeffWrite;
#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;

ngfInt()

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

static

Definition at line 49 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 132 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
 
  gmp_float* f= NULL;
  if (((gmp_float*)a)->isZero() )
  {
    WerrorS(nDivBy0);
  }
  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 270 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 260 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 60 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 415 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 423 of file gnumpfl.cc.

{
  assume( getCoeffType(r) == n_long_R );
  assume( aRing->rep == n_rep_gap_gmp);
 
  if ( from != NULL )
  {
    if (SR_HDL(from) & SR_INT)
    {
      gmp_float f_i= gmp_float(SR_TO_INT(from));
      gmp_float *res=new gmp_float(f_i);
      return (number)res;
    }
    gmp_float f_i=(mpz_ptr)from;
    gmp_float *res=new gmp_float(f_i);
    return (number)res;
  }
  else
    return NULL;

ngfMult()

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

static

Definition at line 173 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 121 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 201 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 222 of file gnumpfl.cc.

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

ngfRead()

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

Definition at line 306 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 396 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 470 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, 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 67 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 162 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 365 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");
  }