Difference between revisions of "ApCoCoA-1:Num.SubABM"

From ApCoCoAWiki
(Updated example)
Line 1: Line 1:
 
   <command>
 
   <command>
 
     <title>Num.SubABM</title>
 
     <title>Num.SubABM</title>
     <short_description>Computes a border basis of an almost vanishing sub-ideal for a set of points and ideal.</short_description>
+
     <short_description>Computes a border basis of an almost vanishing sub-ideal for a set of points and an ideal.</short_description>
 
<syntax>
 
<syntax>
 
Num.SubABM(Points:MAT, Tau:RAT, GBasis:LIST):Object
 
Num.SubABM(Points:MAT, Tau:RAT, GBasis:LIST):Object
Line 9: Line 9:
 
<em>Please note:</em> The function(s) explained on this page is/are using the <em>ApCoCoAServer</em>. You will have to start the ApCoCoAServer in order to use it/them.
 
<em>Please note:</em> The function(s) explained on this page is/are using the <em>ApCoCoAServer</em>. You will have to start the ApCoCoAServer in order to use it/them.
 
<par/>
 
<par/>
This command computes a border basis of an almost vanishing sub-ideal for a set of points and ideal.
+
This command computes a border basis of an almost vanishing sub-ideal for a set of points and an ideal.
 
<par/>
 
<par/>
The current ring has to be a ring over the rationals with a standard-degree
+
The current ring has to be a ring over the rational numbers with a standard-degree
compatible term-ordering. The matrix Points contains the points: each
+
compatible term-ordering. Each row in the matrix <tt>Points</tt> represents one point, so the number of columns must equal the
point is a row in the matrix, so the number of columns must equal the
 
 
number of indeterminates in the current ring.  
 
number of indeterminates in the current ring.  
  
Line 21: Line 20:
 
<item>@param <em>Tau</em> A positive rational number describing which singular values should be treated as 0 (smaller values for tau lead to bigger errors of the polynomials evaluated at the point set). Tau should be in the interval (0,1). As a rule of thumb, Tau is the expected percentage of error on the input points. </item>
 
<item>@param <em>Tau</em> A positive rational number describing which singular values should be treated as 0 (smaller values for tau lead to bigger errors of the polynomials evaluated at the point set). Tau should be in the interval (0,1). As a rule of thumb, Tau is the expected percentage of error on the input points. </item>
  
<item>@param <em>GBasis</em> A homogeneous Groebner Basis in the current ring. This basis defines the ideal we compute the approximate vanishing ideal's basis in.  Warning: for reasons of efficiency the function does not check the validity of GBasis.</item>
+
<item>@param <em>GBasis</em> A homogeneous Groebner Basis in the current ring. This basis defines the ideal in which we compute the basis of the approximate vanishing ideal.  <em>Warning:</em> For reasons of efficiency the function does not check the validity of <tt>GBasis</tt>.</item>
  
 
<item>@return A list of two results. First the border basis as a list of polynomials, second the vector space basis of P/I as a list of terms.</item>
 
<item>@return A list of two results. First the border basis as a list of polynomials, second the vector space basis of P/I as a list of terms.</item>
Line 29: Line 28:
 
The following parameters are optional:
 
The following parameters are optional:
 
<itemize>
 
<itemize>
<item>@param <em>Delta</em> A positiv rational number. Delta describes the computing precision. In different steps, it is crucial, if a value is 0 or not. The algorithm assumes every value in [-Delta, Delta] to be 0. The default value for Delta is 0.00000000001.</item>
+
<item>@param <em>Delta</em> A positive rational number which describes the computing precision. In different steps, it is crucial, if a value is 0 or not. The algorithm assumes every value in <tt>[-Delta, Delta]</tt> to be 0. The default value for <tt>Delta</tt> is 0.00000000001.</item>
  
<item>@param <em>NormalizeType</em> A integer of the range 1..4. The default value is 1. This parameter describes, if / how the input points are normalized. If NormalizeType equals 1, each coordinate is divided by the maximal absolut value of the matrix's corresponding column. This ensures that all point's coordinates are in [-1,1]. With NormalizeType=2 no normalization is done at all. NormalizeType=3 shifts each coordinate to [-1,1]. So it's minimum is mapped to -1 and the maximum to one, describing a unique affine mapping. The last option is NormalizeType=4. In this case, each coordinate is normalized, using the column's euclidian norm. Due to backward compatibility, the default is 1, although 3 is in most cases a better choice.</item>
+
<item>@param <em>NormalizeType</em> A integer of the set <tt>{1,2,3,4}</tt>. The default value is 1. This parameter describes, if and where required the input points are normalized. If <tt>NormalizeType</tt> equals 1, each coordinate of a point is divided by the maximal absolute value of all coordinates of this point. This ensures that all coordinates of the points are within <tt>[-1,1]</tt>. With <tt>NormalizeType=2</tt> no normalization is done at all. <tt>NormalizeType=3</tt> shifts each coordinate to <tt>[-1,1]</tt>, i.e. the minimal coordinate of a point is mapped to -1 and the maximal coordinate to 1, which describes a unique affine mapping. The last option is <tt>NormalizeType=4</tt>. In this case, each point is normalized by its euclidean norm. Although <tt>NormalizeType=3</tt> is in most cases a better choice, the default value is due to backward compatibility 1.</item>
  
 
</itemize>
 
</itemize>

Revision as of 13:20, 8 July 2009

Num.SubABM

Computes a border basis of an almost vanishing sub-ideal for a set of points and an ideal.

Syntax

Num.SubABM(Points:MAT, Tau:RAT, GBasis:LIST):Object
Num.SubABM(Points:MAT, Tau:RAT, GBasis:LIST, Delta:RAT, NormalizeType:INT):Object

Description

Please note: The function(s) explained on this page is/are using the ApCoCoAServer. You will have to start the ApCoCoAServer in order to use it/them.

This command computes a border basis of an almost vanishing sub-ideal for a set of points and an ideal.

The current ring has to be a ring over the rational numbers with a standard-degree

compatible term-ordering. Each row in the matrix Points represents one point, so the number of columns must equal the

number of indeterminates in the current ring.

  • @param Points The points for which a border basis is computed.

  • @param Tau A positive rational number describing which singular values should be treated as 0 (smaller values for tau lead to bigger errors of the polynomials evaluated at the point set). Tau should be in the interval (0,1). As a rule of thumb, Tau is the expected percentage of error on the input points.

  • @param GBasis A homogeneous Groebner Basis in the current ring. This basis defines the ideal in which we compute the basis of the approximate vanishing ideal. Warning: For reasons of efficiency the function does not check the validity of GBasis.

  • @return A list of two results. First the border basis as a list of polynomials, second the vector space basis of P/I as a list of terms.


The following parameters are optional:

  • @param Delta A positive rational number which describes the computing precision. In different steps, it is crucial, if a value is 0 or not. The algorithm assumes every value in [-Delta, Delta] to be 0. The default value for Delta is 0.00000000001.

  • @param NormalizeType A integer of the set {1,2,3,4}. The default value is 1. This parameter describes, if and where required the input points are normalized. If NormalizeType equals 1, each coordinate of a point is divided by the maximal absolute value of all coordinates of this point. This ensures that all coordinates of the points are within [-1,1]. With NormalizeType=2 no normalization is done at all. NormalizeType=3 shifts each coordinate to [-1,1], i.e. the minimal coordinate of a point is mapped to -1 and the maximal coordinate to 1, which describes a unique affine mapping. The last option is NormalizeType=4. In this case, each point is normalized by its euclidean norm. Although NormalizeType=3 is in most cases a better choice, the default value is due to backward compatibility 1.


Example

Use P::=QQ[x,y,z];

Points := Mat([[2/3,0,0],[0,1,0],[0,0,1/3]]);
R:=Num.SubABM(Points, 0.1, [1,x]);
Dec(R[1],2);
R[2];

-- CoCoAServer: computing Cpu Time = 0.015
-------------------------------
[<quotes>0.41 x +0.27 y +0.82 z -0.27  </quotes>, <quotes>0.94 z^2 -0.31 z -0.00  </quotes>, <quotes>1 yz </quotes>, <quotes>1 xz </quotes>, <quotes>0.70 y^2 -0.70 y -0.00  </quotes>, <quotes>1 xy </quotes>]
-------------------------------
[1, z, y]
-------------------------------

See also

Introduction to CoCoAServer

Num.ABM