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

From ApCoCoAWiki
Line 1: Line 1:
 
   <command>
 
   <command>
     <title>Numerical.BBasisOfPointsAVI</title>
+
     <title>Num.BBasisOfPointsAVI</title>
     <short_description>border basis of almost vanishing ideal for a set of points</short_description>
+
     <short_description>Border basis of almost vanishing ideal for a set of points.</short_description>
 
<syntax>
 
<syntax>
$numerical.BBasisOfPointsAVI(Points, Tau, GetO, Delta,  NormalizeType, RREFNormalizeType, RREFUseEps, RREFType):Object
+
Num.BBasisOfPointsAVI(Points:MAT, Tau:RAT, GetO:BOOL):Object
 +
Num.BBasisOfPointsAVI(Points:MAT, Tau:RAT, GetO:BOOL, Delta:RAT,  NormalizeType:INT, RREFNormalizeType:INT, RREFUseEps:BOOL, RREFType:INT):Object
 
</syntax>
 
</syntax>
    <description>
 
{{ApCoCoAServer}} Please also note that you will have to use an ApCoCoAServer with enabled BLAS/LAPACK support.
 
  
This <tt>command</tt> computes an approximate border basis of an almost vanishing ideal for a set of points using the AVI algorithm, which addresses some problems found in the ABM Algorithm.
+
<description>
<par></par>
+
<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.
The current ring has to be a ring over the rationals with a standard-degree
+
 
compatible term-ordering. The matrix Points contains the points: each
+
This function computes an approximate border basis of an almost vanishing ideal for a set of points using the AVI algorithm, which addresses some problems found in the ABM Algorithm.
point is a row in the matrix, so the number of columns must equal the
+
 
number of indeterminates in the current ring.
+
<itemize>
<par></par>
+
<item>@param <em>Points</em> The points for which a border basis is computed.</item>
Tau is a rational <formula>&gt;0</formula>
+
 
describing which singular values should be treated as 0 (smaller values for
+
<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>
tau lead to bigger errors of the polynomials evaluated at the point
+
 
set). Tau should be in the interval <formula>(0,1)</formula>. As a rule of thumb,  
+
<item>@param <em>GetO</em> A Boolean to choose the output. If it is true, the command returns a list of two values: the first contains the border basis, the second one a vector space basis of P/I comprising those power products lying outside the leading term ideal of I. If GetO is false, the function returns only the border basis (not in a list).</item>
Tau is the expected percentage of error on the input points.  
 
  
 +
<item>@return The return value depends on the parameter GetO: GetO=FALSE: The border basis of the given points as a list of polynomials. GetO=TRUE: 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>
  
GetO must be either True or False.  If it is true, the command returns a list
+
</itemize>
of two values: the first contains the approximate border basis, the
 
second one a vector space basis of <formula>P/I</formula> comprising those power products
 
lying outside the leading term ideal of I.  If GetO is false, the function
 
returns only the border basis (not in a list).
 
  
 +
The following parameters are optional:
 +
<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>
  
The last 5 parameters are optional.  
+
<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>RREFNormalizeType</em> Describes, how in each RREF steps the columns are normalized. The options correspond to the ones for NormalizeType and the default is 1 again.</item>
  
Delta must be a positiv rational. 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>@param <em>RREFUseEps</em> must be either true or false! If RREFUseEps is true, the given Delta is used within the RREF to decide if a value equals 0 or not. If this parameter is false, a replacement value for Delta is used, which is based on the matrix's norm. </item>
  
 +
<item>@param <em>RREFType</em> This must be 1 or 2. If RREFType=1, the rref operates column-wise. Otherwise it works row-wise. The default is 1.</item>
 +
</itemize>
  
NormalizeType must be one of 1,2,3, 4. The default value is 2. This parameter describes, how the input points will be normalized.
+
This command computes a border basis of an almost vanishing ideal for a set of points.  
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].
+
<par/>
With NormalizeType=2 no normalization is done at all.
+
The current ring has to be a ring over the rationals with a standard-degree
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.
+
compatible term-ordering. The matrix Points contains the points: each
The last option is NormalizeType=4. In this case, each coordinate is normalized, using the column's euclidian norm.
+
point is a row in the matrix, so the number of columns must equal the
 +
number of indeterminates in the current ring.  
  
RREFNormalizeType describes, how in each RREF steps the columns are normalized. The options correspond to the ones for NormalizeType and the default is 1 again.
 
  
RREFUseEps must be either true or false! If RREFUseEps is true, the given Delta is used within the RREF to decide if a value equals 0 or not. If this parameter is false,
 
a replacement value for Delta is used, which is based on the matrix's norm.
 
  
The last parameter is RREFType. This must be 1 or 2. If RREFType=1, the rref operates column-wise. Otherwise it works row-wise. The default is 1.
 
 
<example>
 
<example>
 
Points := Mat([[1,0,0],[0,0,1],[0,1,0]]);
 
Points := Mat([[1,0,0],[0,0,1],[0,1,0]]);
$numerical.BBasisOfPointsAVI(Points,0.001,True);
+
Num.BBasisOfPointsAVI(Points,0.001,True);
 
-- CoCoAServer: computing Cpu Time = 0
 
-- CoCoAServer: computing Cpu Time = 0
 
-------------------------------
 
-------------------------------
[[1/2x + 1/2y + 4503599627370495/9007199254740992z - 1/2, xy, 1592262918131443/2251799813685248y^2 - 1592262918131443/2251799813685248y, xz, yz, 1592262918131443/2251799813685248z^2 - 1592262918131443/2251799813685248z], [1, z, y]]
+
[[1/2x + 1/2y + 4503599627370495/9007199254740992z - 1/2, xy, 1592262918131443/2251799813685248y^2 - 1592262918131443/2251799813685248y,  
 +
xz, yz, 1592262918131443/2251799813685248z^2 - 1592262918131443/2251799813685248z], [1, z, y]]
 
-------------------------------
 
-------------------------------
 
 
</example>
 
</example>
 
     </description>
 
     </description>

Revision as of 14:02, 20 April 2009

Num.BBasisOfPointsAVI

Border basis of almost vanishing ideal for a set of points.

Syntax

Num.BBasisOfPointsAVI(Points:MAT, Tau:RAT, GetO:BOOL):Object
Num.BBasisOfPointsAVI(Points:MAT, Tau:RAT, GetO:BOOL, Delta:RAT,  NormalizeType:INT, RREFNormalizeType:INT, RREFUseEps:BOOL, RREFType: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 function computes an approximate border basis of an almost vanishing ideal for a set of points using the AVI algorithm, which addresses some problems found in the ABM Algorithm.

  • @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 GetO A Boolean to choose the output. If it is true, the command returns a list of two values: the first contains the border basis, the second one a vector space basis of P/I comprising those power products lying outside the leading term ideal of I. If GetO is false, the function returns only the border basis (not in a list).

  • @return The return value depends on the parameter GetO: GetO=FALSE: The border basis of the given points as a list of polynomials. GetO=TRUE: 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 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.

  • @param NormalizeType 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.

  • @param RREFNormalizeType Describes, how in each RREF steps the columns are normalized. The options correspond to the ones for NormalizeType and the default is 1 again.

  • @param RREFUseEps must be either true or false! If RREFUseEps is true, the given Delta is used within the RREF to decide if a value equals 0 or not. If this parameter is false, a replacement value for Delta is used, which is based on the matrix's norm.

  • @param RREFType This must be 1 or 2. If RREFType=1, the rref operates column-wise. Otherwise it works row-wise. The default is 1.

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

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

compatible term-ordering. The matrix Points contains the points: each

point is a row in the matrix, so the number of columns must equal the number of indeterminates in the current ring.


Example

Points := Mat([[1,0,0],[0,0,1],[0,1,0]]);
Num.BBasisOfPointsAVI(Points,0.001,True);
-- CoCoAServer: computing Cpu Time = 0
-------------------------------
[[1/2x + 1/2y + 4503599627370495/9007199254740992z - 1/2, xy, 1592262918131443/2251799813685248y^2 - 1592262918131443/2251799813685248y, 
xz, yz, 1592262918131443/2251799813685248z^2 - 1592262918131443/2251799813685248z], [1, z, y]]
-------------------------------

See also

Introduction to CoCoAServer

Numerical.GBasisOfPoints

Numerical.HBasisOfPoints

Numerical.GBasisOfPointsInIdeal

Numerical.BBasisOfPointsInIdeal

Numerical.HBasisOfPointsInIdeal

Numerical.FirstVanishingRelations

Numerical.FirstVanishingRelationsInIdeal