D.8.7.13 num_radical1

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.

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));