ntupel.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/ntupel.h"
#include "coeffs/modulop.h"
#include "coeffs/longrat.h"

Go to the source code of this file.

Macros
#define	SI_BUF_SIZE   1024

Typedefs
typedef coeffs *	coeffs_array
typedef number *	number_array
typedef char	si_char_SI_BUF_SIZE[SI_BUF_SIZE]

Functions
static void	nnCoeffWrite (const coeffs r, BOOLEAN b)
static BOOLEAN	nnGreaterZero (number k, const coeffs r)
static number	nnMult (number a, number b, const coeffs r)
static number	nnInit (long l, const coeffs r)
static long	nnInt (number &n, const coeffs r)
static int	nnSize (number n, const coeffs r)
static number	nnCopy (number a, const coeffs r)
static number	nnAdd (number a, number b, const coeffs r)
static number	nnSub (number a, number b, const coeffs r)
static BOOLEAN	nnIsZero (number a, const coeffs r)
static BOOLEAN	nnIsOne (number a, const coeffs r)
static BOOLEAN	nnIsMOne (number a, const coeffs r)
static number	nnDiv (number a, number b, const coeffs r)
static number	nnInvers (number a, const coeffs r)
static number	nnInpNeg (number a, const coeffs r)
static BOOLEAN	nnGreater (number a, number b, const coeffs r)
static BOOLEAN	nnEqual (number a, number b, const coeffs r)
static void	nnWriteShort (number a, const coeffs r)
static void	nnWriteLong (number a, const coeffs r)
static const char *	nnRead (const char *s, number *a, const coeffs r)
BOOLEAN	nnDBTest (number a, const char *f, const int l, const coeffs r)
static number	nnMap (number from, const coeffs aRing, const coeffs r)
nMapFunc	nnSetMap (const coeffs, const coeffs)
static char *	nnCoeffString (const coeffs r)
static char *	nnCoeffName (const coeffs r)
static void	nnDelete (number *a, const coeffs r)
static number	nnChineseRemainder (number *, number *, int, BOOLEAN, CFArray &, const coeffs)
BOOLEAN	nnInitChar (coeffs n, void *p)
	Initialize r. More...

Variables
STATIC_VAR si_char_SI_BUF_SIZE	nnCoeffName_buf

Macro Definition Documentation

SI_BUF_SIZE

#define SI_BUF_SIZE   1024

Typedef Documentation

coeffs_array

typedef coeffs* coeffs_array

Definition at line 24 of file ntupel.cc.

number_array

typedef number* number_array

Definition at line 25 of file ntupel.cc.

si_char_SI_BUF_SIZE

typedef char si_char_SI_BUF_SIZE[SI_BUF_SIZE]

Definition at line 544 of file ntupel.cc.

Function Documentation

nnAdd()

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

static

Definition at line 151 of file ntupel.cc.

{
  number_array A=(number_array)a;
  number_array B=(number_array)b;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  i=0;
  loop
  {
    C[i]=d[i]->cfAdd(A[i],B[i],d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnChineseRemainder()

static number nnChineseRemainder ( number *  , number *  , int  , BOOLEAN  , CFArray &  , const  coeffs  )

static

Definition at line 581 of file ntupel.cc.

{
  WerrorS("not yet");
  return NULL;
}

nnCoeffName()

static char * nnCoeffName ( const coeffs  r )

static

Definition at line 546 of file ntupel.cc.

{
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  nnCoeffName_buf[0]='\0';
  loop
  {
    char *s=d[i]->cfCoeffName(d[i]);
    strncat(nnCoeffName_buf,s, SI_BUF_SIZE-strlen(nnCoeffName_buf)-strlen(s)-1);
    i++;
    if (d[i]==NULL) break;
    strncat(nnCoeffName_buf," x ", SI_BUF_SIZE-strlen(nnCoeffName_buf)-4);
  }
  return nnCoeffName_buf;
}

nnCoeffString()

static char * nnCoeffString ( const coeffs  r )

static

Definition at line 526 of file ntupel.cc.

{
  #define SI_BUF_SIZE 1024
  char buf[SI_BUF_SIZE];
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  buf[0]='\0';
  loop
  {
    char *s=d[i]->cfCoeffName(d[i]);
    strncat(buf,s, SI_BUF_SIZE-strlen(buf)-strlen(s)-1);
    i++;
    if (d[i]==NULL) break;
    strncat(buf," x ", SI_BUF_SIZE-strlen(buf)-4);
  }
  return omStrDup(buf);
}

nnCoeffWrite()

static void nnCoeffWrite ( const coeffs  r, BOOLEAN  b  )

static

Definition at line 27 of file ntupel.cc.

{
  int i=0;
  coeffs_array d=(coeffs_array)r->data;
  loop
  {
    d[i]->cfCoeffWrite(d[i],b);
    i++;
    if (d[i]==NULL) break;
    PrintS(" x ");
  }
}

nnCopy()

static number nnCopy ( number  a, const coeffs  r  )

static

Definition at line 131 of file ntupel.cc.

{
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  i=0;
  loop
  {
    C[i]=d[i]->cfCopy(A[i],d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnDBTest()

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

Definition at line 484 of file ntupel.cc.

{
  BOOLEAN bo=FALSE;
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    bo&=d[i]->cfDBTest(A[i],f,l,d[i]);
    i++;
    if (d[i]==NULL) return bo;
  }
}

nnDelete()

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

static

Definition at line 562 of file ntupel.cc.

{
  if (*a!=NULL)
  {
    coeffs_array d=(coeffs_array)r->data;
    number_array A=(number_array)(*a);
    int i=0;
    loop
    {
      d[i]->cfDelete(&(A[i]),d[i]);
      i++;
      if (d[i]==NULL) break;
    }
    omFreeSize(d,i*sizeof(number));
    *a=NULL;
  }
}

nnDiv()

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

static

Definition at line 232 of file ntupel.cc.

{
  if (nnIsZero(b,r))
  {
    WerrorS(nDivBy0);
    return nnInit(0,r);
  }
  number_array A=(number_array)a;
  number_array B=(number_array)b;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  i=0;
  loop
  {
    C[i]=d[i]->cfDiv(A[i],B[i],d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnEqual()

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

static

Definition at line 320 of file ntupel.cc.

{
  int i=0;
  number_array A=(number_array)a;
  number_array B=(number_array)b;
  coeffs_array d=(coeffs_array)r->data;
  loop
  {
    if (!d[i]->cfEqual(A[i],B[i],d[i])) return FALSE;
    i++;
    if (d[i]==NULL) return TRUE;
  }
}

nnGreater()

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

static

Definition at line 302 of file ntupel.cc.

{
  number_array A=(number_array)a;
  number_array B=(number_array)b;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    if (d[i]->type==n_R) return d[i]->cfGreater(A[i],B[i],d[i]);
    else if (d[i]->type==n_long_R) return d[i]->cfGreater(A[i],B[i],d[i]);
    i++;
    if (d[i]==NULL)
    {
      return d[i-1]->cfGreater(A[i],B[i],d[i-1]);
    }
  }
}

nnGreaterZero()

static BOOLEAN nnGreaterZero ( number  k, const coeffs  r  )

static

Definition at line 41 of file ntupel.cc.

{
  int i=0;
  coeffs_array d=(coeffs_array)r->data;
  loop
  {
    if (d[i]->type==n_R) return d[i]->cfGreaterZero(k,d[i]);
    if (d[i]->type==n_long_R) return d[i]->cfGreaterZero(k,d[i]);
    i++;
    if (d[i]==NULL) return d[i-1]->cfGreaterZero(k,d[i-1]);
  }
}

nnInit()

static number nnInit ( long  l, const coeffs  r  )

static

Definition at line 78 of file ntupel.cc.

{
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  i=0;
  loop
  {
    C[i]=d[i]->cfInit(l,d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnInitChar()

BOOLEAN nnInitChar	(	coeffs	n,
		void *	p
	)

Initialize r.

Definition at line 588 of file ntupel.cc.

{
  n->data=p;
  n->type=n_nTupel;
  n->is_field=TRUE;
  n->is_domain=TRUE;
 
  n->ch = 0;
  n->cfCoeffString = nnCoeffString;
  n->cfCoeffName = nnCoeffName;
 
  n->cfInit = nnInit;
  n->cfDelete = nnDelete;
  n->cfInt  = nnInt;
  n->cfAdd  = nnAdd;
  n->cfSub  = nnSub;
  n->cfMult = nnMult;
  n->cfDiv  = nnDiv;
  n->cfExactDiv= nnDiv;
  n->cfInpNeg   = nnInpNeg;
  n->cfInvers= nnInvers;
  n->cfCopy  = nnCopy;
  n->cfGreater = nnGreater;
  n->cfEqual = nnEqual;
  n->cfIsZero = nnIsZero;
  n->cfIsOne = nnIsOne;
  n->cfIsMOne = nnIsMOne;
  n->cfGreaterZero = nnGreaterZero;
  n->cfWriteShort = nnWriteShort;
  n->cfWriteLong = nnWriteLong;
  n->cfRead = nnRead;
  //n->cfPower = nnPower;
  n->cfSetMap = nnSetMap;
  n->cfCoeffWrite  = nnCoeffWrite;
 
  n->cfChineseRemainder=nnChineseRemainder;
  n->cfSize  = nnSize;
#ifdef LDEBUG
  n->cfDBTest=nnDBTest;
#endif
  //n->nCoeffIsEqual = ndCoeffIsEqual;
  return FALSE;
}

nnInpNeg()

static number nnInpNeg ( number  a, const coeffs  r  )

static

Definition at line 283 of file ntupel.cc.

{
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  i=0;
  loop
  {
    A[i]=d[i]->cfInpNeg(A[i],d[i]);
    i++;
    if (d[i]==NULL) return (number)A;
  }
}

nnInt()

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

static

Definition at line 100 of file ntupel.cc.

{
  int i=0;
  coeffs_array d=(coeffs_array)r->data;
  loop
  {
    if (d[i]->type==n_R) return d[i]->cfInt(n,d[i]);
    if (d[i]->type==n_long_R) return d[i]->cfInt(n,d[i]);
    i++;
    if (d[i]==NULL) return d[i-1]->cfInt(n,d[i-1]);
  }
}

nnInvers()

static number nnInvers ( number  a, const coeffs  r  )

static

Definition at line 258 of file ntupel.cc.

{
  if (nnIsZero(a,r))
  {
    WerrorS(nDivBy0);
    return nnInit(0,r);
  }
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  i=0;
  loop
  {
    C[i]=d[i]->cfInvers(A[i],d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnIsMOne()

static BOOLEAN nnIsMOne ( number  a, const coeffs  r  )

static

Definition at line 219 of file ntupel.cc.

{
  int i=0;
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  loop
  {
    if (!d[i]->cfIsMOne(A[i],d[i])) return FALSE;
    i++;
    if (d[i]==NULL) return TRUE;
  }
}

nnIsOne()

static BOOLEAN nnIsOne ( number  a, const coeffs  r  )

static

Definition at line 206 of file ntupel.cc.

{
  int i=0;
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  loop
  {
    if (!d[i]->cfIsOne(A[i],d[i])) return FALSE;
    i++;
    if (d[i]==NULL) return TRUE;
  }
}

nnIsZero()

static BOOLEAN nnIsZero ( number  a, const coeffs  r  )

static

Definition at line 193 of file ntupel.cc.

{
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    if (!d[i]->cfIsZero(A[i],d[i])) return FALSE;
    i++;
    if (d[i]==NULL) return TRUE;
  }
}

nnMap()

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

static

Definition at line 499 of file ntupel.cc.

{
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  nMapFunc nMap;
  i=0;
  loop
  {
    nMap=d[i]->cfSetMap(aRing,d[i]);
    if (nMap==NULL) Werror("no map for compoment %d",i);
    else C[i]=nMap(from,aRing,d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnMult()

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

static

Definition at line 54 of file ntupel.cc.

{
  number_array A=(number_array)a;
  number_array B=(number_array)b;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  i=0;
  loop
  {
    C[i]=d[i]->cfMult(A[i],B[i],d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnRead()

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

static

Definition at line 383 of file ntupel.cc.

{
  // reading R/longR part, setting everything else to random
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc0(i*sizeof(number));
  *a=(number)C;
  int is_zero=2; /*1: read 0, 2: read <> 0 */
  int reader=-1;
  i=0;
  loop
  {
    if ((d[i]->type==n_R)
    || (d[i]->type==n_long_R))
    {
      s=d[i]->cfRead(s,&(C[i]),d[i]);
      is_zero=d[i]->cfIsZero(C[i],d[i]);
      reader=i;
      break;
    }
    i++;
    if (d[i]==NULL)
    {
      s=d[i-1]->cfRead(s,&(C[i-1]),d[i-1]);
      is_zero=d[i-1]->cfIsZero(C[i-1],d[i-1]);
      reader=i-1;
      break;
    }
  }
  if (is_zero==0)
  { // seting random values!=0 for Z/p
    i=0;
    loop
    {
      if (i!=reader)
      {
        if (d[i]->type==n_Zp)
        {
          int l=0;
          do
          {
            l=siRand();
          } while(l%(d[i]->ch)==0);
          C[i]=d[i]->cfInit(l,d[i]);
        }
        else if ((d[i]->type==n_R)
        || (d[i]->type==n_long_R))
        {
          nMapFunc nMap=d[i]->cfSetMap(d[reader],d[i]);
          C[i]=nMap(C[reader],d[reader],d[i]);
        }
        else
        {
          WerrorS("reading is not suppiorted for such compinations of coeffs");
          break;
        }
      }
      i++;
      if (d[i]==NULL) break;
    }
  }
  else if (is_zero==1)
  { // read a zero
    i=0;
    loop
    {
      C[i]=d[i]->cfInit(0,d[i]);
      i++;
      if (d[i]==NULL) break;
    }
  }
  else if (reader==-1)/* not R or long_R*/
  {
    s=d[0]->cfRead(s,&(C[0]),d[0]);
    i=1;
    loop
    {
      if (d[i]==NULL) break;
      i++;
      {
        nMapFunc nMap=d[i]->cfSetMap(d[0],d[i]);
        C[i]=nMap(C[0],d[0],d[i]);
      }
    }
  }
  else
  {
    WerrorS("nnRead: should not happen");
  }
  return s;
}

nnSetMap()

nMapFunc nnSetMap	(	const	coeffs,
		const	coeffs
	)

Definition at line 521 of file ntupel.cc.

{
  return nnMap;
}

nnSize()

static int nnSize ( number  n, const coeffs  r  )

static

Definition at line 113 of file ntupel.cc.

{
  int s=0;
  int i=0;
  coeffs_array d=(coeffs_array)r->data;
  loop
  {
    if (d[i]->type==n_R) s+=d[i]->cfSize(n,d[i]);
    else if (d[i]->type==n_long_R) s+=d[i]->cfSize(n,d[i]);
    i++;
    if (d[i]==NULL)
    {
      if (s==0) return d[i-1]->cfSize(n,d[i-1]);
      else return s;
    }
  }
}

nnSub()

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

static

Definition at line 172 of file ntupel.cc.

{
  number_array A=(number_array)a;
  number_array B=(number_array)b;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    i++;
    if (d[i]==NULL) break;
  }
  number *C=(number*)omAlloc(i*sizeof(number));
  i=0;
  loop
  {
    C[i]=d[i]->cfSub(A[i],B[i],d[i]);
    i++;
    if (d[i]==NULL) return (number)C;
  }
}

nnWriteLong()

static void nnWriteLong ( number  a, const coeffs  r  )

static

Definition at line 348 of file ntupel.cc.

{
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  StringSetS("(");
  loop
  {
    d[i]->cfWriteLong(A[i],d[i]);
    i++;
    if (d[i]==NULL) {StringAppendS(")"); return; }
    else StringAppendS(",");
  }
}

nnWriteShort()

static void nnWriteShort ( number  a, const coeffs  r  )

static

Definition at line 334 of file ntupel.cc.

{
  number_array A=(number_array)a;
  coeffs_array d=(coeffs_array)r->data;
  int i=0;
  loop
  {
    if (d[i]->type==n_R) { d[i]->cfWriteLong(A[i],d[i]); return; }
    else if (d[i]->type==n_long_R) {d[i]->cfWriteLong(A[i],d[i]); return; }
    i++;
    if (d[i]==NULL) { d[i-1]->cfWriteLong(A[i-1],d[i-1]); return; }
  }
}

Variable Documentation

nnCoeffName_buf

STATIC_VAR si_char_SI_BUF_SIZE nnCoeffName_buf

Definition at line 545 of file ntupel.cc.