CPowerMultiplier Class Reference

#include <ncSAMult.h>

Public Types
typedef CPower	CExponent

Public Member Functions
	CPowerMultiplier (ring r)
virtual	~CPowerMultiplier ()
CSpecialPairMultiplier *	GetPair (int i, int j) const
CSpecialPairMultiplier *&	GetPair (int i, int j)
virtual poly	MultiplyEE (const CExponent expLeft, const CExponent expRight)
virtual poly	MultiplyME (const poly pMonom, const CExponent expRight)
virtual poly	MultiplyEM (const CExponent expLeft, const poly pMonom)
poly	MultiplyPE (const poly pPoly, const CExponent expRight)
poly	MultiplyEP (const CExponent expLeft, const poly pPoly)
poly	MultiplyPEDestroy (poly pPoly, const CExponent expRight)
poly	MultiplyEPDestroy (const CExponent expLeft, poly pPoly)
Public Member Functions inherited from CMultiplier< CPower >
	CMultiplier (ring rBaseRing)
virtual	~CMultiplier ()
ring	GetBasering () const
int	NVars () const
poly	LM (const poly pTerm, const ring r, int i=1) const
poly	MultiplyTE (const poly pTerm, const CPower expRight)
poly	MultiplyET (const CPower expLeft, const poly pTerm)

Private Attributes
CSpecialPairMultiplier **	m_specialpairs

Additional Inherited Members
Protected Attributes inherited from CMultiplier< CPower >
const ring	m_basering
const int	m_NVars

Detailed Description

Definition at line 163 of file ncSAMult.h.

Member Typedef Documentation

CExponent

typedef CPower CPowerMultiplier::CExponent

Definition at line 194 of file ncSAMult.h.

Constructor & Destructor Documentation

CPowerMultiplier()

CPowerMultiplier::CPowerMultiplier

(

ring

r

)

Definition at line 832 of file ncSAMult.cc.

                                        : CMultiplier<CPower>(r)
{
#if OUTPUT
  PrintS("CPowerMultiplier::CPowerMultiplier(ring)!");
  PrintLn();
#endif
 
  m_specialpairs = (CSpecialPairMultiplier**)omAlloc0( ((NVars() * (NVars()-1)) / 2) * sizeof(CSpecialPairMultiplier*) );
 
  for( int i = 1; i < NVars(); i++ )
    for( int j = i + 1; j <= NVars(); j++ )
      GetPair(i, j) = AnalyzePair(GetBasering(), i, j); // factory method!
}

~CPowerMultiplier()

CPowerMultiplier::~CPowerMultiplier ( )

virtual

Definition at line 847 of file ncSAMult.cc.

{
#if OUTPUT
  PrintS("CPowerMultiplier::~CPowerMultiplier()!");
  PrintLn();
#endif
 
  omFreeSize((ADDRESS)m_specialpairs, ((NVars() * (NVars()-1)) / 2) * sizeof(CSpecialPairMultiplier*) );
}

Member Function Documentation

GetPair() [1/2]

CSpecialPairMultiplier *& CPowerMultiplier::GetPair ( int  i, int  j  )

inline

Definition at line 183 of file ncSAMult.h.

    {
      assume( m_specialpairs != NULL );
      assume( i > 0 );
      assume( i < j );
      assume( j <= NVars() );
 
      return m_specialpairs[( (NVars() * ((i)-1) - ((i) * ((i)-1))/2 + (j)-1) - (i) )];
    }

GetPair() [2/2]

CSpecialPairMultiplier * CPowerMultiplier::GetPair ( int  i, int  j  ) const

inline

Definition at line 173 of file ncSAMult.h.

    {
      assume( m_specialpairs != NULL );
      assume( i > 0 );
      assume( i < j );
      assume( j <= NVars() );
 
      return m_specialpairs[( (NVars() * ((i)-1) - ((i) * ((i)-1))/2 + (j)-1) - (i) )];
    }

MultiplyEE()

poly CPowerMultiplier::MultiplyEE ( const CExponent  expLeft, const CExponent  expRight  )

virtual

Implements CMultiplier< CPower >.

Definition at line 991 of file ncSAMult.cc.

{
#if OUTPUT
  PrintS("CPowerMultiplier::MultiplyEE)!");
  PrintLn();
#endif
 
  const int i = expRight.Var, j = expLeft.Var;
  const int ei = expRight.Power, ej = expLeft.Power;
 
#if OUTPUT
  Print("Input: var(%d)^{%d} * var(%d)^{%d}", j, ej, i, ei);
  PrintLn();
#endif
 
  assume(1 <= i);
  assume(j <= NVars());
  assume(1 <= j);
  assume(i <= NVars());
  assume(ei > 0);
  assume(ej > 0);
 
  if( i >= j )
  {
    const ring r = GetBasering();
 
    poly product = p_One(r);
    p_SetExp(product, j, ej, r);
    p_SetExp(product, i, ei, r);
    p_Setm(product, r);
 
    return product;
 
  } else
  {
    assume(i <  j);
 
    // No Cache Lookup!? :(
 
    CSpecialPairMultiplier* pSpecialMultiplier = GetPair(i, j);
 
    // Special case?
    if( pSpecialMultiplier != NULL )
    {
      assume( pSpecialMultiplier->GetI() == i );
      assume( pSpecialMultiplier->GetJ() == j );
      assume( pSpecialMultiplier->GetBasering() == GetBasering() );
 
      return pSpecialMultiplier->MultiplyEE(ej, ei);
    } else
    {
      // Perform general NC Multiplication:
      // TODO
 
      WerrorS("Sorry the general case is not implemented this way yet!!!");
      assume(0);
 
      // poly product = NULL;
    }
  }
 
  return NULL;
}

MultiplyEM()

poly CPowerMultiplier::MultiplyEM ( const CExponent  expLeft, const poly  pMonom  )

virtual

Implements CMultiplier< CPower >.

Definition at line 924 of file ncSAMult.cc.

{
  const ring r = GetBasering();
 
  // TODO: as above! (difference due to Left/Right semmantics!)
  const int j = expLeft.Var;
  const int n = expLeft.Power;
 
#if OUTPUT
  Print("CPowerMultiplier::MultiplyEM(var(%d)^{%d} * monom)!", j, n);
  PrintLn();
  PrintS("Monom: "); p_Write(pMonom, r);
#endif
 
  assume( (j > 0) && (j <= NVars()));
 
  if( n == 0 )
    return p_Head(pMonom, r); // Copy?!?
 
 
  int v = 1; // NVars();
  int e = p_GetExp(pMonom, v, r);
 
  while((v < j) && (e == 0))
    e = p_GetExp(pMonom, ++v, r);
 
  if( v == j )
  {
    poly p = p_Head(pMonom, r);
    p_SetExp(p, j, e + n, r);
    p_Setm(p, r);
 
    return p;
  }
 
  assume( v < j );
  assume( e > 0 );
 
 
  // And now the General Case: v > j!
 
  poly p = MultiplyEE( expLeft, CPower(v, e) ); // Easy way!
 
  ++v;
 
  while(v <= NVars())
  {
    e = p_GetExp(pMonom, v, r);
 
    if( e > 0 )
      p = MultiplyPEDestroy(p, CPower(v, e));
 
    ++v;
  }
 
#if OUTPUT
  PrintS("CPowerMultiplier::MultiplyEM() ===> ");
  p_Write(p, r);
#endif
 
  return p;
 
}

MultiplyEP()

poly CPowerMultiplier::MultiplyEP ( const CExponent  expLeft, const poly  pPoly  )

inline

Definition at line 223 of file ncSAMult.h.

    {
      bool bUsePolynomial = TEST_OPT_NOT_BUCKETS || (pLength(pPoly) < MIN_LENGTH_BUCKET);
      CPolynomialSummator sum(GetBasering(), bUsePolynomial);
 
      for( poly q = pPoly; q !=NULL; q = pNext(q) )
        sum += MultiplyET(expLeft, q);
 
      return sum;
    }

MultiplyEPDestroy()

poly CPowerMultiplier::MultiplyEPDestroy ( const CExponent  expLeft, poly  pPoly  )

inline

Definition at line 247 of file ncSAMult.h.

    {
      bool bUsePolynomial = TEST_OPT_NOT_BUCKETS || (pLength(pPoly) < MIN_LENGTH_BUCKET);
      CPolynomialSummator sum(GetBasering(), bUsePolynomial);
 
      for( ; pPoly!=NULL; pPoly  = p_LmDeleteAndNext(pPoly, GetBasering()) )
        sum += MultiplyET(expLeft, pPoly);
 
      return sum;
    }

MultiplyME()

poly CPowerMultiplier::MultiplyME ( const poly  pMonom, const CExponent  expRight  )

virtual

Implements CMultiplier< CPower >.

Definition at line 860 of file ncSAMult.cc.

{
  const int j = expRight.Var;
  const int n = expRight.Power;
 
  const ring r = GetBasering();
 
#if OUTPUT
  Print("CPowerMultiplier::MultiplyME(monom * var(%d)^{%d})!", j, n);
  PrintLn();
  PrintS("Monom: "); p_Write(pMonom, r);
#endif
 
  assume( (j > 0) && (j <= NVars()));
 
  if( n == 0 )
    return p_Head(pMonom, r); // Copy?!?
 
 
  int v = NVars();
  int e = p_GetExp(pMonom, v, r);
 
  while((v > j) && (e == 0))
    e = p_GetExp(pMonom, --v, r);
 
  // TODO: review this!
  if( v == j )
  {
    poly p = p_Head(pMonom, r);
    p_SetExp(p, v, e + n, r);
    p_Setm(p, r);
 
    return p;
  }
 
  assume( v > j );
  assume( e > 0 );
 
  // And now the General Case: v > j!
 
  poly p = MultiplyEE( CPower(v, e), expRight ); // Easy way!
 
  --v;
 
  while(v > 0)
  {
    e = p_GetExp(pMonom, v, GetBasering());
 
    if( e > 0 )
      p = MultiplyEPDestroy(CPower(v, e), p);
 
    --v;
  }
 
#if OUTPUT
  PrintS("CPowerMultiplier::MultiplyME() ===> ");
  p_Write(p, GetBasering());
#endif
 
  return p;
}

MultiplyPE()

poly CPowerMultiplier::MultiplyPE ( const poly  pPoly, const CExponent  expRight  )

inline

Definition at line 211 of file ncSAMult.h.

    {
      bool bUsePolynomial = TEST_OPT_NOT_BUCKETS || (pLength(pPoly) < MIN_LENGTH_BUCKET);
      CPolynomialSummator sum(GetBasering(), bUsePolynomial);
 
      for( poly q = pPoly; q !=NULL; q = pNext(q) )
        sum += MultiplyTE(q, expRight);
 
      return sum;
    }

MultiplyPEDestroy()

poly CPowerMultiplier::MultiplyPEDestroy ( poly  pPoly, const CExponent  expRight  )

inline

Definition at line 235 of file ncSAMult.h.

    {
      bool bUsePolynomial = TEST_OPT_NOT_BUCKETS || (pLength(pPoly) < MIN_LENGTH_BUCKET);
      CPolynomialSummator sum(GetBasering(), bUsePolynomial);
 
      for( ; pPoly!=NULL; pPoly  = p_LmDeleteAndNext(pPoly, GetBasering()) )
        sum += MultiplyTE(pPoly, expRight);
 
      return sum;
    }

Field Documentation

m_specialpairs

CSpecialPairMultiplier** CPowerMultiplier::m_specialpairs

private

Definition at line 166 of file ncSAMult.h.

---

The documentation for this class was generated from the following files:

• libpolys/polys/nc/ncSAMult.h
• libpolys/polys/nc/ncSAMult.cc