1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
|
/*
TO DO:
Find a good way to input the name of the new variable and ajust this to helpText.
*/
#include "parser.h"
#include "printer.h"
#include "polynomial.h"
#include "division.h"
#include "buchberger.h"
#include "wallideal.h"
#include "lp.h"
#include "reversesearch.h"
#include "termorder.h"
#include "gfanapplication.h"
#include "wallideal.h"
class HomogenizeApplication : public GFanApplication
{
SimpleOption optionIdeal;
SimpleOption optionHomogenizationVector;
public:
const char *helpText()
{
return "This program homogenises a list of polynomials by introducing an extra variable. The name of the variable to be introduced is read from the input after the list of polynomials. Without the -w option the homogenisation is done with respect to total degree.\n"
"Example:\n"
"Input:\n"
"Q[x,y]{y-1}\n"
"Output:\n"
"Q[x,y,z]{y-z}\n";
;
}
HomogenizeApplication():
optionIdeal("-i","Treat input as an ideal. This will make the program compute the homogenisation of the input ideal. This is done by computing a degree Groebner basis and homogenising it."),
optionHomogenizationVector("-w","Specify a homogenisation vector. The length of the vector must be the same as the number of variables in the ring. The vector is read from the input after the list of polynomials.\n")
{
registerOptions();
}
char *name()
{
return "_homogenize";
}
int main()
{
FileParser P(Stdin);
PolynomialSet g=P.parsePolynomialSetWithRing();
string homogenizingVariableName=P.parseVariableName();
IntegerVector grading;
if(optionHomogenizationVector.getValue())
grading=P.parseIntegerVector();
else
grading=IntegerVector::allOnes(g.numberOfVariablesInRing());
if(optionIdeal.getValue())
{
WeightReverseLexicographicTermOrder t(grading);
buchberger(&g,t);
}
PolynomialSet h(g.getRing().withVariablesAppended(homogenizingVariableName));
AsciiPrinter(Stdout).printPolynomialRing(h.getRing());
AsciiPrinter(Stdout).printNewLine();
for(PolynomialSet::const_iterator i=g.begin();i!=g.end();i++)
{
h.push_back(i->homogenization(h.getRing(),&grading));
}
AsciiPrinter(Stdout).printPolynomialSet(h);
return 0;
}
};
static HomogenizeApplication theApplication;
|
