My Project
Loading...
Searching...
No Matches
Functions
fac_sqrfree.h File Reference

squarefree part and factorization over Q, Q(a) More...

#include "canonicalform.h"

Go to the source code of this file.

Functions

CFFList sortCFFList (CFFList &F)
 
CFFList sqrFreeZ (const CanonicalForm &f)
 
CanonicalForm sqrfPart (const CanonicalForm &F)
 squarefree part of a poly More...
 
CFFList sqrFreeFp (const CanonicalForm &f)
 

Detailed Description

squarefree part and factorization over Q, Q(a)

Definition in file fac_sqrfree.h.

Function Documentation

◆ sortCFFList()

CFFList sortCFFList ( CFFList F)

Definition at line 25 of file fac_sqrfree.cc.

26{
28
29 int exp;
31 CFFListIterator I = F;
33
34 // join elements with the same degree
35 while ( I.hasItem() ) {
36 f = I.getItem().factor();
37 exp = I.getItem().exp();
38 I++;
39 while ( I.hasItem() && I.getItem().exp() == exp ) {
40 f *= I.getItem().factor();
41 I++;
42 }
43 result.append( CFFactor( f, exp ) );
44 }
45
46 return result;
47}
Factor< CanonicalForm > CFFactor
FILE * f
Definition: checklibs.c:9
factory's main class
Definition: canonicalform.h:86
T & getItem() const
Definition: ftmpl_list.cc:431
void sort(int(*)(const T &, const T &))
Definition: ftmpl_list.cc:339
return result
Definition: facAbsBiFact.cc:75
static int compareFactors(const CFFactor &f, const CFFactor &g)
Definition: fac_sqrfree.cc:19
gmp_float exp(const gmp_float &a)
Definition: mpr_complex.cc:357

◆ sqrfPart()

CanonicalForm sqrfPart ( const CanonicalForm F)

squarefree part of a poly

Parameters
[in]Fsome poly

Definition at line 115 of file fac_sqrfree.cc.

116{
117 if (F.inCoeffDomain())
118 return F;
119 CFMap M;
121 CanonicalForm w, v, b;
123 int i= 1;
124 for (; i <= A.level(); i++)
125 {
126 if (!deriv (A, Variable (i)).isZero())
127 break;
128 }
129
130 w= gcd (A, deriv (A, Variable (i)));
131 b= A/w;
132 result= b;
133 if (degree (w) < 1)
134 return M (result);
135 i++;
136 for (; i <= A.level(); i++)
137 {
138 if (!deriv (w, Variable (i)).isZero())
139 {
140 b= w;
141 w= gcd (w, deriv (w, Variable (i)));
142 b /= w;
143 if (degree (b) < 1)
144 break;
146 if (degree (g) > 0)
147 result *= b/g;
148 if (degree (g) <= 0)
149 result *= b;
150 }
151 }
152 result= M (result);
153 return result;
154}
int degree(const CanonicalForm &f)
CanonicalForm deriv(const CanonicalForm &f, const Variable &x)
int i
Definition: cfEzgcd.cc:132
g
Definition: cfModGcd.cc:4090
CanonicalForm b
Definition: cfModGcd.cc:4103
CanonicalForm compress(const CanonicalForm &f, CFMap &m)
CanonicalForm compress ( const CanonicalForm & f, CFMap & m )
Definition: cf_map.cc:210
class CFMap
Definition: cf_map.h:85
bool inCoeffDomain() const
factory's class for variables
Definition: variable.h:33
const CanonicalForm & w
Definition: facAbsFact.cc:51
const Variable & v
< [in] a sqrfree bivariate poly
Definition: facBivar.h:39
bool isZero(const CFArray &A)
checks if entries of A are zero
#define A
Definition: sirandom.c:24
#define M
Definition: sirandom.c:25
int gcd(int a, int b)
Definition: walkSupport.cc:836

◆ sqrFreeFp()

CFFList sqrFreeFp ( const CanonicalForm f)

◆ sqrFreeZ()

CFFList sqrFreeZ ( const CanonicalForm f)

Definition at line 49 of file fac_sqrfree.cc.

50{
51 if ( a.inCoeffDomain() )
52 return CFFactor( a, 1 );
53 CanonicalForm aa, LcA;
54 if (isOn (SW_RATIONAL))
55 {
56 LcA= bCommonDen (a);
57 aa= a*LcA;
58 }
59 else
60 {
61 LcA= icontent (a);
62 if (lc (a).sign() < 0)
63 LcA= -LcA;
64 aa= a/LcA;
65 }
66 CanonicalForm cont = content( aa );
67 aa /= cont;
68 CanonicalForm b = aa.deriv(), c = gcd( aa, b );
69 CanonicalForm y, z, w = aa / c;
70 int i = 1;
71 CFFList F;
72 Variable v = aa.mvar();
73 CanonicalForm lcinv;
74 while ( c.degree(v) != 0 )
75 {
76 y = gcd( w, c ); z = w / y;
77 if ( degree( z, v ) > 0 )
78 {
79 if (isOn (SW_RATIONAL))
80 {
81 lcinv= 1/Lc (z);
82 z *= lcinv;
83 z *= bCommonDen (z);
84 }
85 if (lc (z).sign() < 0)
86 z= -z;
87 F.append( CFFactor( z, i ) );
88 }
89 i++;
90 w = y; c = c / y;
91 }
92 if ( degree( w,v ) > 0 )
93 {
94 if (isOn (SW_RATIONAL))
95 {
96 lcinv= 1/Lc (w);
97 w *= lcinv;
98 w *= bCommonDen (w);
99 }
100 if (lc (w).sign() < 0)
101 w= -w;
102 F.append( CFFactor( w, i ) );
103 }
104 if ( ! cont.isOne() )
105 {
106 CFFList buf= sqrFreeZ (cont);
107 buf.removeFirst();
108 F = Union( F, buf );
109 }
110 F.insert (CFFactor (LcA, 1));
111 return F;
112}
bool isOn(int sw)
switches
CanonicalForm FACTORY_PUBLIC content(const CanonicalForm &)
CanonicalForm content ( const CanonicalForm & f )
Definition: cf_gcd.cc:603
CanonicalForm lc(const CanonicalForm &f)
CanonicalForm FACTORY_PUBLIC icontent(const CanonicalForm &f)
CanonicalForm icontent ( const CanonicalForm & f )
Definition: cf_gcd.cc:74
CanonicalForm Lc(const CanonicalForm &f)
CanonicalForm bCommonDen(const CanonicalForm &f)
CanonicalForm bCommonDen ( const CanonicalForm & f )
static const int SW_RATIONAL
set to 1 for computations over Q
Definition: cf_defs.h:31
Variable mvar() const
mvar() returns the main variable of CO or Variable() if CO is in a base domain.
CF_NO_INLINE bool isOne() const
CanonicalForm deriv() const
deriv() - return the formal derivation of CO.
void append(const T &)
Definition: ftmpl_list.cc:256
void insert(const T &)
Definition: ftmpl_list.cc:193
const CanonicalForm int const CFList const Variable & y
Definition: facAbsFact.cc:53
CFFList sqrFreeZ(const CanonicalForm &a)
Definition: fac_sqrfree.cc:49
template List< Variable > Union(const List< Variable > &, const List< Variable > &)
static int sign(int x)
Definition: ring.cc:3427
int status int void * buf
Definition: si_signals.h:59