|
D.8.7.13 num_radical1
Procedure from library recover.lib (see recover_lib).
- Usage:
- num_radical1(P,D,C); list P, int D, bigint C
P a list of lists representing a witness point set representing an ideal I
D should be a bound to the degree of the elements of the components
C the number with which the images of the Veronese embeddings are multiplied
- Return:
- list of ideals: each of the ideals a prime component of the radical of I
- Remarks:
- This procedure merely calls num_prime_decom1 with the same input and then
intersects the returned components.
- Note:
- Should only be called from a ring over the complex numbers.
Example:
| LIB "recover.lib";
//First, we write the input file for bertini and compute the radical symbolically.
ring r=0,(x,y,z),dp;
ideal I=4xy2-4z3,-2x2y+5xz2;
ideal Rad=radical(I);
writeBertiniInput(I,100);
//Then we attempt to compute the radical via the hybrid algorithm.
ring R=(complex,100,i),(x,y,z),dp;
system("sh","bertini input");
list P=getWitnessSet();
int D=2;
bigint C=bigint(10)**30;
ideal Rad1=num_radical1(P,D,C);
//Lastly, we compare the results.
Rad1;
ideal Rad=fetch(r,Rad);
Rad;
reduce(Rad,std(Rad1));
reduce(Rad1,std(Rad));
| See also:
num_prime_decom1;
num_radical2.
|