gnumpfl.cc File Reference

#include "misc/auxiliary.h"
#include "omalloc/omalloc.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 152 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 384 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 409 of file gnumpfl.cc.

{
  static 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 402 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 502 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 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 185 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 278 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 251 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 241 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 231 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 507 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 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);
  }
  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 271 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 261 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 462 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 454 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 416 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 445 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 424 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 174 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 202 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 223 of file gnumpfl.cc.

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

ngfRead()

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

Definition at line 307 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 397 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 471 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 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 163 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 366 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");
  }