Top
Back: num_radical_via_randlincom
Forward: num_radical2
FastBack:
FastForward:
Up: recover_lib
Top: Singular Manual
Contents: Table of Contents
Index: Index
About: About this document

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.


Top Back: num_radical_via_randlincom Forward: num_radical2 FastBack: FastForward: Up: recover_lib Top: Singular Manual Contents: Table of Contents Index: Index About: About this document
            User manual for Singular version 4.3.2, 2023, generated by texi2html.