shortfl.cc File Reference

#include "misc/auxiliary.h"
#include "misc/mylimits.h"
#include "reporter/reporter.h"
#include "coeffs/numbers.h"
#include "coeffs/coeffs.h"
#include "coeffs/mpr_complex.h"
#include "coeffs/shortfl.h"
#include "coeffs/longrat.h"
#include <cmath>

Go to the source code of this file.

Data Structures
union	nf

Macros
#define	SR_HDL(A)   ((long)(A))
#define	IS_INT(A)   ((A)->s==3)
#define	IS_IMM(A)   (SR_HDL(A) & SR_INT)
#define	GET_NOM(A)   ((A)->z)
#define	GET_DENOM(A)   ((A)->n)

Functions
static nMapFunc	nrSetMap (const coeffs src, const coeffs dst)
	Get a mapping function from src into the domain of this type: n_R. More...
static number	nrMapQ (number from, const coeffs r, const coeffs aRing)
SI_FLOAT	nrFloat (number n)
	Converts a n_R number into a float. Needed by Maps. More...
static BOOLEAN	nrGreaterZero (number k, const coeffs r)
static number	nrMult (number a, number b, const coeffs r)
static number	nrInit (long i, const coeffs r)
static long	nrInt (number &n, const coeffs r)
static number	nrAdd (number a, number b, const coeffs r)
static number	nrSub (number a, number b, const coeffs r)
static BOOLEAN	nrIsZero (number a, const coeffs r)
static BOOLEAN	nrIsOne (number a, const coeffs r)
static BOOLEAN	nrIsMOne (number a, const coeffs r)
static number	nrDiv (number a, number b, const coeffs r)
static number	nrInvers (number c, const coeffs r)
static number	nrNeg (number c, const coeffs r)
static BOOLEAN	nrGreater (number a, number b, const coeffs r)
static BOOLEAN	nrEqual (number a, number b, const coeffs r)
static void	nrWrite (number a, const coeffs r)
static const char *	nrRead (const char *s, number *a, const coeffs r)
static BOOLEAN	nrDBTest (number a, const char *f, const int l, const coeffs r)
static number	nrMapP (number from, const coeffs aRing, const coeffs r)
static number	nrMapLongR (number from, const coeffs aRing, const coeffs r)
static number	nrMapC (number from, const coeffs aRing, const coeffs r)
static number	nrMapZ (number from, const coeffs aRing, const coeffs r)
static char *	nrCoeffString (const coeffs r)
static char *	nrCoeffName (const coeffs r)
BOOLEAN	nrInitChar (coeffs n, void *p)
	Initialize r. More...

Variables
static const SI_FLOAT	nrEps = 1.0e-3

Macro Definition Documentation

GET_DENOM

#define GET_DENOM

(

A

)

((A)->n)

GET_NOM

#define GET_NOM

(

A

)

((A)->z)

IS_IMM

#define IS_IMM

(

A

)

(SR_HDL(A) & SR_INT)

IS_INT

#define IS_INT

(

A

)

((A)->s==3)

SR_HDL

#define SR_HDL

(

A

)

((long)(A))

Function Documentation

nrAdd()

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

static

Definition at line 94 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  SI_FLOAT x = nf(a).F();
  SI_FLOAT y = nf(b).F();
  SI_FLOAT f = x + y;
  if (x > 0.0)
  {
    if (y < 0.0)
    {
      x = f / (x - y);
      if (x < 0.0)
        x = -x;
      if (x < nrEps)
        f = 0.0;
    }
  }
  else
  {
    if (y > 0.0)
    {
      x = f / (y - x);
      if (x < 0.0)
        x = -x;
      if (x < nrEps)
        f = 0.0;
    }
  }
  return nf(f).N();
}

nrCoeffName()

static char * nrCoeffName ( const coeffs  r )

static

Definition at line 653 of file shortfl.cc.

{
  return (char*)"Float()";
}

nrCoeffString()

static char * nrCoeffString ( const coeffs  r )

static

Definition at line 648 of file shortfl.cc.

{
  return omStrDup("Float()");
}

nrDBTest()

static BOOLEAN nrDBTest ( number  a, const char *  f, const int  l, const coeffs  r  )

static

Definition at line 370 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  return TRUE;
}

nrDiv()

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

static

Definition at line 183 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  SI_FLOAT n = nf(b).F();
  if (n == 0.0)
  {
    WerrorS(nDivBy0);
    return nf((SI_FLOAT)0.0).N();
  }
  else
    return nf(nf(a).F() / n).N();
}

nrEqual()

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

static

Definition at line 224 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  number x = nrSub(a,b,r);
  return nf(x).F() == nf((SI_FLOAT)0.0).F();
}

nrFloat()

SI_FLOAT nrFloat

(

number

n

)

Converts a n_R number into a float. Needed by Maps.

Definition at line 48 of file shortfl.cc.

{
  return nf(n).F();
}

nrGreater()

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

static

Definition at line 217 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  return nf(a).F() > nf(b).F();
}

nrGreaterZero()

static BOOLEAN nrGreaterZero ( number  k, const coeffs  r  )

static

Definition at line 53 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  return nf(k).F() >= 0.0;
}

nrInit()

static number nrInit ( long  i, const coeffs  r  )

static

Definition at line 70 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  SI_FLOAT f = (SI_FLOAT)i;
  return nf(nf(f).F()).N();
}

nrInitChar()

BOOLEAN nrInitChar	(	coeffs	n,
		void *	p
	)

Initialize r.

Definition at line 658 of file shortfl.cc.

{
  assume( getCoeffType(n) == n_R );
 
  assume( p == NULL );
 
  n->is_field=TRUE;
  n->is_domain=TRUE;
  n->rep=n_rep_float;
 
  //n->cfKillChar = ndKillChar; /* dummy */
  n->ch = 0;
  n->cfCoeffName = nrCoeffName;
  n->cfCoeffString = nrCoeffString;
 
  n->cfInit = nrInit;
  n->cfInt  = nrInt;
  n->cfAdd   = nrAdd;
  n->cfSub   = nrSub;
  n->cfMult  = nrMult;
  n->cfDiv   = nrDiv;
  n->cfExactDiv= nrDiv;
  n->cfInpNeg   = nrNeg;
  n->cfInvers= nrInvers;
  //n->cfCopy  = ndCopy;
  n->cfGreater = nrGreater;
  n->cfEqual = nrEqual;
  n->cfIsZero = nrIsZero;
  n->cfIsOne = nrIsOne;
  n->cfIsMOne = nrIsMOne;
  n->cfGreaterZero = nrGreaterZero;
  n->cfWriteLong = nrWrite;
  n->cfRead = nrRead;
  //n->cfPower = nrPower;
  n->cfSetMap = nrSetMap;
 
    /* nName= ndName; */
    /*nSize  = ndSize;*/
#ifdef LDEBUG
  n->cfDBTest=nrDBTest; // not yet implemented: nrDBTest;
#endif
 
  //n->nCoeffIsEqual = ndCoeffIsEqual;
 
  n->float_len = SHORT_REAL_LENGTH;
  n->float_len2 = SHORT_REAL_LENGTH;
 
  // TODO: Any variables?
  return FALSE;
}

nrInt()

static long nrInt ( number &  n, const coeffs  r  )

static

Definition at line 81 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  long i;
  SI_FLOAT f = nf(n).F();
  if (((SI_FLOAT)(-MAX_INT_VAL-1) <= f) || ((SI_FLOAT)MAX_INT_VAL >= f))
    i = (long)f;
  else
    i = 0;
  return i;
}

nrInvers()

static number nrInvers ( number  c, const coeffs  r  )

static

Definition at line 197 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  SI_FLOAT n = nf(c).F();
  if (n == 0.0)
  {
    WerrorS(nDivBy0);
    return nf((SI_FLOAT)0.0).N();
  }
  return nf(1.0 / n).N();
}

nrIsMOne()

static BOOLEAN nrIsMOne ( number  a, const coeffs  r  )

static

Definition at line 174 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  SI_FLOAT aa=nf(a).F()+1.0;
  if (aa<0.0) aa=-aa;
  return (aa<nrEps);
}

nrIsOne()

static BOOLEAN nrIsOne ( number  a, const coeffs  r  )

static

Definition at line 165 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  SI_FLOAT aa=nf(a).F()-1.0;
  if (aa<0.0) aa=-aa;
  return (aa<nrEps);
}

nrIsZero()

static BOOLEAN nrIsZero ( number  a, const coeffs  r  )

static

Definition at line 158 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  return (0.0 == nf(a).F());
}

nrMapC()

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

static

Definition at line 397 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
  assume( getCoeffType(aRing) == n_long_C );
 
  gmp_float h = ((gmp_complex*)from)->real();
  SI_FLOAT t =(SI_FLOAT)mpf_get_d((mpf_srcptr)&h);
  return nf(t).N();
}

nrMapLongR()

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

static

Definition at line 388 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
  assume( getCoeffType(aRing) == n_long_R );
 
  SI_FLOAT t =(SI_FLOAT)mpf_get_d((mpf_srcptr)from);
  return nf(t).N();
}

nrMapP()

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

static

Definition at line 378 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
  assume( getCoeffType(aRing) ==  n_Zp );
 
  int i = (int)((long)from);
  SI_FLOAT f = (SI_FLOAT)i;
  return nf(f).N();
}

nrMapQ()

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

static

Definition at line 408 of file shortfl.cc.

{
/* in longrat.h
#define SR_INT    1
#define mpz_size1(A) (ABS((A)->_mp_size))
*/
#define SR_HDL(A) ((long)(A))
#define IS_INT(A) ((A)->s==3)
#define IS_IMM(A) (SR_HDL(A) & SR_INT)
#define GET_NOM(A) ((A)->z)
#define GET_DENOM(A) ((A)->n)
 
  assume( getCoeffType(r) == n_R );
  assume( aRing->rep == n_rep_gap_rat );
 
  if (IS_IMM(from))
  {
    SI_FLOAT f = (SI_FLOAT)SR_TO_INT(from);
    return nf(nf(f).F()).N();
  }
  else
  {
    /* read out the enumerator */
    if (IS_INT(from))
    {
      mpf_t e;
      mpf_init(e);
      mpf_set_z(e,GET_NOM(from));
      SI_FLOAT f = mpf_get_d(e);
      mpf_clear(e);
      return nf(nf(f).F()).N();
    }
    else /*quotient*/
    {
      mpf_t z,n,q;
      mpf_init(z);
      mpf_init(n);
      mpf_init(q);
      mpf_set_z(z,GET_NOM(from));
      mpf_set_z(n,GET_DENOM(from));
      mpf_div(q,z,n);
      mpf_clear(z);
      mpf_clear(n);
      SI_FLOAT f = mpf_get_d(q);
      mpf_clear(q);
      return nf(nf(f).F()).N();
    }
  }
}

nrMapZ()

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

static

Definition at line 458 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
  assume( aRing->rep == n_rep_gap_gmp );
 
  mpz_ptr z;
  mpz_ptr zz=NULL;
  if (IS_IMM(from))
  {
     zz=(mpz_ptr)omAlloc(sizeof(mpz_t));
     mpz_init_set_si(zz,SR_TO_INT(from));
     z=zz;
  }
  else
  {
    /* read out the enumerator */
    z=(mpz_ptr)from;
  }
 
  int i = mpz_size1(z);
  mpf_t e;
  mpf_init(e);
  mpf_set_z(e,z);
  int sign= mpf_sgn(e);
  mpf_abs (e, e);
 
  if (zz!=NULL)
  {
    mpz_clear(zz);
    omFreeSize(zz,sizeof(mpz_t));
  }
  if(i>4)
  {
    WerrorS("float overflow");
    return nf(0.0).N();
  }
  double basis;
  signed long int exp;
  basis = mpf_get_d_2exp(&exp, e);
  SI_FLOAT f= sign*ldexp(basis,exp);
  mpf_clear(e);
  return nf(f).N();
}

nrMult()

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

static

Definition at line 60 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  return nf(nf(a).F() * nf(b).F()).N();
}

nrNeg()

static number nrNeg ( number  c, const coeffs  r  )

static

Definition at line 210 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  return nf(-nf(c).F()).N();
}

nrRead()

static const char * nrRead ( const char *  s, number *  a, const coeffs  r  )

static

Definition at line 299 of file shortfl.cc.

{
 
  assume( getCoeffType(r) == n_R );
 
  static const char *nIllegalChar="illegal character in number";
 
  const char *t;
  const char *start=s;
  SI_FLOAT z1,z2;
  SI_FLOAT n=1.0;
 
  s = nrEatr(s, &z1);
  if (*s == '/')
  {
    if (s==start) { WerrorS(nIllegalChar);return s; }
    s++;
    s = nrEatr(s, &z2);
    if (z2==0.0)
      WerrorS(nDivBy0);
    else
      z1 /= z2;
  }
  else if (*s =='.')
  {
    if (s==start) { WerrorS(nIllegalChar);return s; }
    s++;
    t = s;
    while (*t >= '0' && *t <= '9')
    {
      t++;
      n *= 10.0;
    }
    s = nrEatr(s, &z2);
    z1 = (z1*n + z2) / n;
    if (*s=='e')
    {
      int e=0; /* exponent */
      int si=1;/* sign of exponent */
      s++;
      if (*s=='+') s++;
      else if (*s=='-') {s++; si=-1; }
      while (*s >= '0' && *s <= '9')
      {
        e=e*10+(*s)-'0';
        s++;
      }
      if (si==1)
      {
        while (e>0) {z1*=10.0; e--; }
      }
      else
      {
        while (e>0) {z1/=10.0; e--; }
      }
    }
  }
  *a = nf(z1).N();
  return s;
}

nrSetMap()

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

static

Get a mapping function from src into the domain of this type: n_R.

Definition at line 617 of file shortfl.cc.

{
  assume( getCoeffType(dst) == n_R );
 
  if (src->rep==n_rep_gap_rat) /*Q, Z */
  {
    return nrMapQ;
  }
  if (src->rep==n_rep_gap_gmp) /*Q, Z */
  {
    return nrMapZ;
  }
  if ((src->rep==n_rep_gmp_float) && nCoeff_is_long_R(src))
  {
    return nrMapLongR;
  }
  if ((src->rep==n_rep_float) && nCoeff_is_R(src))
  {
    return ndCopyMap;
  }
  if ((src->rep==n_rep_int) && nCoeff_is_Zp(src))
  {
    return nrMapP;
  }
  if ((src->rep==n_rep_gmp_complex) && nCoeff_is_long_C(src))
  {
    return nrMapC;
  }
  return NULL;
}

nrSub()

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

static

Definition at line 126 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  SI_FLOAT x = nf(a).F();
  SI_FLOAT y = nf(b).F();
  SI_FLOAT f = x - y;
  if (x > 0.0)
  {
    if (y > 0.0)
    {
      x = f / (x + y);
      if (x < 0.0)
        x = -x;
      if (x < nrEps)
        f = 0.0;
    }
  }
  else
  {
    if (y < 0.0)
    {
      x = f / (x + y);
      if (x < 0.0)
        x = -x;
      if (x < nrEps)
        f = 0.0;
    }
  }
  return nf(f).N();
}

nrWrite()

static void nrWrite ( number  a, const coeffs  r  )

static

Definition at line 232 of file shortfl.cc.

{
  assume( getCoeffType(r) == n_R );
 
  //#if SIZEOF_DOUBLE == SIZEOF_LONG
  //char ch[16];
  //int n = sprintf(ch,"%12.6e", nf(a).F());
  //#else
  char ch[11];
  int n = sprintf(ch,"%9.3e", nf(a).F());
  //#endif
  if (ch[0] == '-')
  {
    char* chbr = new char[n+3];
    memcpy(&chbr[2],&ch[1],n-1);
    chbr[0] = '-';
    chbr[1] = '(';
    chbr[n+1] = ')';
    chbr[n+2] = '\0';
    StringAppendS(chbr);
    delete[] chbr;
  }
  else
    StringAppend("(%s)",ch);
}

Variable Documentation

nrEps

const SI_FLOAT nrEps = 1.0e-3

static

Definition at line 30 of file shortfl.cc.