Home Online Manual
Top
Back: extend
Forward: locusdg
FastBack: general_lib
FastForward: inout_lib
Up: grobcov_lib
Top: Singular Manual
Contents: Table of Contents
Index: Index
About: About this document

D.2.4.10 locus

Procedure from library grobcov.lib (see grobcov_lib).

Options:
The algorithm allows the following options as pair of arguments: "movdim", d : by default movdim is 2 but it can be set to other values, and represents the number of mever variables. they should be given as the last variables of the ring.
"version", v : There are two versions of the algorithm. ("version",1) is a full algorithm that always distinguishes correctly between 'Normal' and 'Special' components, whereas \("version",0) can decalre a component as 'Normal' being really 'Special', but is more effective. By default ("version",1) is used when the number of variables is less than 4 and 0 if not. The user can force to use one or other version, but it is not recommended. "system", ideal F: if the initial system is passed as an argument. This is actually not used. "comments", c: by default it is 0, but it can be set to 1. Usually locus problems have mover coordinates, variables and tracer coordinates. The mover coordinates are to be placed as the last variables, and by default, its number is 2. If one consider locus problems in higer dimensions, the number of mover coordinates (placed as the last variables) is to be given as an option.

Return:
The locus. The output is a list of the components ( C_1,.. C_n ) where C_i = (p_i,(p_i1,..p_is_i), type_i,l evel_i ) and type_i can be 'Normal', 'Special', Accumulation', 'Degenerate'. The 'Special' components return more information, namely the antiimage of the component, that is 0-dimensional for these kind of components.
Normal components: for each point in the component,
the number of solutions in the variables is finite, and the solutions depend on the point in the component.
Special components: for each point in the component, the number of solutions in the variables is finite. The antiimage of the component is 0-dimensional.
Accumlation points: are 0-dimensional components whose antiimage is not zero-dimansional.
Degenerate components: are components of dimension greater than 0 whose antiimage is not-zero-diemansional.
The components are given in canonical P-representation. The levels of a class of locus are 1,
because they represent locally closed. sets.

Note:
It can only be called after computing the grobcov of the parametrical ideal in generic representation ('ext',0), which is the default.
The basering R, must be of the form Q[a_1,..,a_m][x_1,..,x_n].

Example:
 
LIB "grobcov.lib";
ring R=(0,a,b),(x,y),dp;
short=0;
// Concoid
ideal S96=x^2+y^2-4,(b-2)*x-a*y+2*a,(a-x)^2+(b-y)^2-1;
// System S96=
S96;
==> S96[1]=x^2+y^2-4
==> S96[2]=(b-2)*x+(-a)*y+(2*a)
==> S96[3]=x^2+y^2+(-2*a)*x+(-2*b)*y+(a^2+b^2-1)
locus(grobcov(S96));
==> [1]:
==>    [1]:
==>       _[1]=(a^4+2*a^2*b^2-9*a^2+b^4-9*b^2+4*b+12)
==>    [2]:
==>       [1]:
==>          _[1]=1
==>    [3]:
==>       Normal
==>    [4]:
==>       1
==> [2]:
==>    [1]:
==>       _[1]=(a^2+b^2-4*b+3)
==>    [2]:
==>       [1]:
==>          _[1]=1
==>    [3]:
==>       [1]:
==>          Special
==>       [2]:
==>          y-2,x
==>    [4]:
==>       1
kill R;
// ********************************************
ring R=(0,a,b),(x4,x3,x2,x1),dp;
short=0;
ideal S=(x1-3)^2+(x2-1)^2-9,
(4-x2)*(x3-3)+(x1-3)*(x4-1),
(3-x1)*(x3-x1)+(4-x2)*(x4-x2),
(4-x4)*a+(x3-3)*b+3*x4-4*x3,
(a-x1)^2+(b-x2)^2-(x1-x3)^2-(x2-x4)^2;
// System S=
S;
==> S[1]=x2^2+x1^2-2*x2-6*x1+1
==> S[2]=-x3*x2+x4*x1-3*x4+4*x3+3*x2-x1-9
==> S[3]=-x4*x2+x2^2-x3*x1+x1^2+4*x4+3*x3-4*x2-3*x1
==> S[4]=(-a+3)*x4+(b-4)*x3+(4*a-3*b)
==> S[5]=-x4^2-x3^2+2*x4*x2+2*x3*x1+(-2*b)*x2+(-2*a)*x1+(a^2+b^2)
locus(grobcov(S));
==> locus detected that the mover must avoid point (x1-3,x2-4) in order to ob\
   tain the correct locus
==> [1]:
==>    [1]:
==>       _[1]=(a^2-6*a+b^2+b+7)
==>    [2]:
==>       [1]:
==>          _[1]=(b-2)
==>          _[2]=(a^2-6*a+13)
==>    [3]:
==>       Normal
==>    [4]:
==>       1
==> [2]:
==>    [1]:
==>       _[1]=(a^4-12*a^3+2*a^2*b^2-13*a^2*b+236*a^2-12*a*b^2+78*a*b-1200*a+\
   b^4-13*b^3+60*b^2-85*b+1495)
==>    [2]:
==>       [1]:
==>          _[1]=(b-4)
==>          _[2]=(a-3)
==>    [3]:
==>       Normal
==>    [4]:
==>       1
==> [3]:
==>    [1]:
==>       _[1]=(b-2)
==>       _[2]=(a^2-6*a+13)
==>    [2]:
==>       [1]:
==>          _[1]=1
==>    [3]:
==>       Accumulation
==>    [4]:
==>       1
==> [4]:
==>    [1]:
==>       _[1]=(b-4)
==>       _[2]=(a-3)
==>    [2]:
==>       [1]:
==>          _[1]=1
==>    [3]:
==>       Accumulation
==>    [4]:
==>       1
kill R;
"********************************************";
==> ********************************************
ring R=(0,x,y),(x1,x2),dp;
short=0;
ideal S=-(x - 5)*(x1 - 1) - (x2 - 2)*(y - 2),
(x1 - 5)^2 + (x2 - 2)^2 - 4,
-2*(x - 5)*(x2 - 2) + 2*(x1 - 5)*(y - 2);
// System S=
S;
==> S[1]=(-x+5)*x1+(-y+2)*x2+(x+2*y-9)
==> S[2]=x1^2+x2^2-10*x1-4*x2+25
==> S[3]=(2*y-4)*x1+(-2*x+10)*x2+(4*x-10*y)
locus(grobcov(S));
==> [1]:
==>    [1]:
==>       _[1]=(3*x^2-30*x-y^2+4*y+71)
==>    [2]:
==>       [1]:
==>          _[1]=(y-2)
==>          _[2]=(x-5)
==>    [3]:
==>       [1]:
==>          Special
==>       [2]:
==>          x1-4,x2^2-4*x2+1
==>    [4]:
==>       1
==> [2]:
==>    [1]:
==>       _[1]=(y-2)
==>       _[2]=(x-5)
==>    [2]:
==>       [1]:
==>          _[1]=1
==>    [3]:
==>       Accumulation
==>    [4]:
==>       1