clapconv.h File Reference

#include "polys/monomials/ring.h"
#include "factory/factory.h"

Go to the source code of this file.

Functions
poly	convFactoryPSingP (const CanonicalForm &f, const ring r)
CanonicalForm	convSingPFactoryP (poly p, const ring r)
int	convFactoryISingI (const CanonicalForm &f)
CanonicalForm	convSingAPFactoryAP (poly p, const Variable &a, const ring r)
poly	convFactoryAPSingAP (const CanonicalForm &f, const ring r)
poly	convFactoryAPSingAP_R (const CanonicalForm &f, int par_start, int var_start)
CanonicalForm	convSingAFactoryA (poly p, const Variable &a, const ring r)
poly	convFactoryASingA (const CanonicalForm &f, const ring r)
CanonicalForm	convSingTrPFactoryP (poly p, const ring r)
BOOLEAN	convSingTrP (poly p, const ring r)
poly	convFactoryPSingTrP (const CanonicalForm &f, const ring r)

Function Documentation

convFactoryAPSingAP()

poly convFactoryAPSingAP	(	const CanonicalForm &	f,
		const ring	r
	)

Definition at line 184 of file clapconv.cc.

{
  return convFactoryAPSingAP_R(f,0,rPar(r),r);
}

convFactoryAPSingAP_R()

poly convFactoryAPSingAP_R	(	const CanonicalForm &	f,
		int	par_start,
		int	var_start
	)

convFactoryASingA()

poly convFactoryASingA	(	const CanonicalForm &	f,
		const ring	r
	)

Definition at line 284 of file clapconv.cc.

{
  poly a=NULL;
  for( CFIterator i=f; i.hasTerms(); i++)
  {
    number n= convFactoryNSingAN( i.coeff(), r );
    if (n_IsZero(n,r->cf->extRing->cf))
    {
      n_Delete(&n,r->cf->extRing->cf);
    }
    else
    {
      poly t= p_Init (r->cf->extRing);
      pGetCoeff(t)=n;
      p_SetExp(t,1,i.exp(),r->cf->extRing);
      //p_Setm(t,r->cf->extRing);// not needed for rings with 1 variable
      a=p_Add_q(a,t,r->cf->extRing);
    }
  }
  if (a!=NULL)
  {
    if( r->cf->extRing != NULL )
      if (r->cf->extRing->qideal->m[0]!=NULL)
      {
        poly l=r->cf->extRing->qideal->m[0];
        if (p_GetExp(a,1,r->cf->extRing) >= p_GetExp(l,1,r->cf->extRing))
          a = p_PolyDiv (a, l, FALSE, r->cf->extRing); // ???
      }
  }
  return a;
}

convFactoryISingI()

int convFactoryISingI

(

const CanonicalForm &

f

)

Definition at line 142 of file clapconv.cc.

{
  if (!f.isImm()) WerrorS("int overflow in det");
  return f.intval();
}

convFactoryPSingP()

poly convFactoryPSingP	(	const CanonicalForm &	f,
		const ring	r
	)

Definition at line 40 of file clapconv.cc.

{
  if (f.isZero()) return NULL;
  int n = rVar(r)+1;
  /* ASSERT( level( f ) <= pVariables, "illegal number of variables" ); */
  int * exp = (int*)omAlloc0(n*sizeof(int));
  sBucket_pt result_bucket=sBucketCreate(r);
  conv_RecPP( f, exp, result_bucket, r );
  poly result; int dummy;
  sBucketDestroyMerge(result_bucket,&result,&dummy);
  omFreeSize((ADDRESS)exp,n*sizeof(int));
  return result;
}

convFactoryPSingTrP()

poly convFactoryPSingTrP	(	const CanonicalForm &	f,
		const ring	r
	)

Definition at line 366 of file clapconv.cc.

{
  if (f.isZero()) return NULL;
  int n = rVar(r)+1;
  int * exp = (int*)omAlloc0(n*sizeof(int));
  poly result = NULL;
  convRecTrP( f, exp, result , rPar(r), r );
  omFreeSize((ADDRESS)exp,n*sizeof(int));
  return result;
}

convSingAFactoryA()

CanonicalForm convSingAFactoryA	(	poly	p,
		const Variable &	a,
		const ring	r
	)

Definition at line 232 of file clapconv.cc.

{
  CanonicalForm result = 0;
  int e;
 
  while ( p!=NULL )
  {
    CanonicalForm term;
    if ( rField_is_Zp_a(r) )
    {
      term = n_Int( p_GetCoeff( p, r->cf->extRing ), r->cf->extRing->cf );
    }
    else
    {
      if ( SR_HDL(p_GetCoeff( p, r->cf->extRing )) & SR_INT )
        term = SR_TO_INT(p_GetCoeff( p, r->cf->extRing )) ;
      else
      {
        if ( p_GetCoeff( p, r->cf->extRing )->s == 3 )
        {
          mpz_t dummy;
          mpz_init_set( dummy, (p_GetCoeff( p,r->cf->extRing )->z) );
          term = make_cf( dummy );
        }
        else
        {
          // assume s==0 or s==1
          mpz_t num, den;
          On(SW_RATIONAL);
          mpz_init_set( num, (p_GetCoeff( p, r->cf->extRing )->z) );
          mpz_init_set( den, (p_GetCoeff( p, r->cf->extRing )->n) );
          term = make_cf( num, den, ( p_GetCoeff( p, r->cf->extRing )->s != 1 ));
        }
      }
    }
    if ( (e = p_GetExp( p, 1, r->cf->extRing )) != 0 )
      term *= power( a , e );
    result += term;
    p = pNext( p );
  }
  return result;
}

convSingAPFactoryAP()

CanonicalForm convSingAPFactoryAP	(	poly	p,
		const Variable &	a,
		const ring	r
	)

Definition at line 148 of file clapconv.cc.

{
  CanonicalForm result = 0;
  int e, n = r-> N;
  int off=rPar(r);
 
  if (!rField_is_Zp_a(r))
    On(SW_RATIONAL);
  while ( p!=NULL)
  {
    CanonicalForm term=convSingAFactoryA(((poly)p_GetCoeff(p, r->cf->extRing)),a, r);
    for ( int i = 1; i <= n; i++ )
    {
      if ( (e = p_GetExp( p, i, r )) != 0 )
        term *= CanonicalForm( Variable( i + off), e );
    }
    result += term;
    pIter( p );
  }
  return result;
}

convSingPFactoryP()

CanonicalForm convSingPFactoryP	(	poly	p,
		const ring	r
	)

Definition at line 136 of file clapconv.cc.

{
  BOOLEAN setChar=TRUE;
  return convSingPFactoryP_intern(p,pLength(p),setChar,r);
}

convSingTrP()

BOOLEAN convSingTrP	(	poly	p,
		const ring	r
	)

Definition at line 352 of file clapconv.cc.

{
  while ( p!=NULL )
  {
    n_Normalize(p_GetCoeff(p, r), r->cf);
 
    // test if denominator is constant
    if (!p_IsConstant(DEN ((fraction)p_GetCoeff (p,r)),r->cf->extRing))
      return FALSE;
    pIter(p);
  }
  return TRUE;
}

convSingTrPFactoryP()

CanonicalForm convSingTrPFactoryP	(	poly	p,
		const ring	r
	)

Definition at line 316 of file clapconv.cc.

{
  CanonicalForm result = 0;
  int e, n = rVar(r);
  int offs = rPar(r);
 
  while ( p!=NULL )
  {
    //n_Normalize(p_GetCoeff(p, r), r->cf);
 
    // test if denominator is constant
    if (!errorreported && !p_IsConstant(DEN ((fraction)p_GetCoeff (p,r)),r->cf->extRing))
      WerrorS("conversion error: denominator!= 1");
 
    CanonicalForm term=convSingPFactoryP(NUM ((fraction)p_GetCoeff(p, r)),r->cf->extRing);
 
    // if denominator is not NULL it should be a constant at this point
    if (DEN ((fraction)p_GetCoeff(p,r)) != NULL)
    {
      CanonicalForm den= convSingPFactoryP(DEN ((fraction)p_GetCoeff(p, r)),r->cf->extRing);
      if (rChar (r) == 0)
        On (SW_RATIONAL);
      term /= den;
    }
 
    for ( int i = n; i > 0; i-- )
    {
      if ( (e = p_GetExp( p, i,r )) != 0 )
        term = term * power( Variable( i + offs ), e );
    }
    result += term;
    p = pNext( p );
  }
  return result;
}