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

D.8.7.7 getWitnessSet

Procedure from library recover.lib (see recover_lib).

Usage:
getWitnessSet();

Assume:
There is a text-document "main_data" in the current directory which was produced by Bertini.
The basefield is the field of real numbers or the field of complex numbers.

Return:
list; a list P of lists p_i of numbers: P a set of witness points

Note:
Reads the file "main_data", searches the strings containing the witness points, and converts them into floating point numbers.

Example:
 
LIB "recover.lib";
//First, we write the input file for bertini, then run bertini
ring r=0,(x,y,z),dp;
ideal I=(x-y)*(y-z)*(x-z);
writeBertiniInput(I,40);
system("sh","bertini input");
==> 
==>    Bertini(TM) v1.6
==>     (May 22, 2018)
==> 
==>  D.J. Bates, J.D. Hauenstein,
==>  A.J. Sommese, C.W. Wampler
==> 
==> (using GMP v6.0.0, MPFR v3.1.2)
==> 
==> 
==> 
==> NOTE: You have requested to use adaptive path tracking.  Please make sure\
    that you have
==> setup the following tolerances appropriately:
==> CoeffBound: 6.000000000000e+00, DegreeBound: 3.000000000000e+00
==> AMPSafetyDigits1: 1, AMPSafetyDigits2: 1, AMPMaxPrec: 160
==> 
==> 
==> Tracking regeneration codim 1 of 1: 3 paths to track.
==> Tracking path 0 of 3
==> Tracking path 1 of 3
==> Tracking path 2 of 3
==> 
==> Sorting codimension 1 of 1: 3 paths to sort.
==> Sorting 0 of 3
==> Sorting 1 of 3
==> Sorting 2 of 3
==> 
==> 
==> ************ Regenerative Cascade Summary ************
==> 
==> NOTE: nonsingular vs singular is based on rank deficiency and identical e\
   ndpoints
==> 
==> |codim|   paths   |witness superset| nonsingular | singular |nonsolutions\
   | inf endpoints | other bad endpoints
==> ----------------------------------------------------------------------------------------------------------------
==> | 1   |   3       |   3            |  3          |  0       |  0         \
   |   0           |  0
==> ----------------------------------------------------------------------------------------------------------------
==> |total|   3
==> 
==> ****************************************************
==> 
==> 
==> 
==> *************** Witness Set Summary ****************
==> 
==> NOTE: nonsingular vs singular is based on rank deficiency and identical e\
   ndpoints
==> 
==> |codim| witness points | nonsingular | singular 
==> -------------------------------------------------
==> | 1   |   3            |  3          |  0       
==> -------------------------------------------------
==> 
==> ****************************************************
==> 
==> 
==> Calculating traces for codimension 1.
==> Calculating 0 of 3
==> Calculating 1 of 3
==> Calculating 2 of 3
==> 
==> Using combinatorial trace test to decompose codimension 1.
==> 
==> 
==> ************* Witness Set Decomposition *************
==> 
==> | dimension | components | classified | unclassified
==> -----------------------------------------------------
==> |   1       |   0        |   0        |  3
==> -----------------------------------------------------
==> 
==> ************** Decomposition by Degree **************
==> 
==> *****************************************************
==> 
==> 0
//Then we change the ring and extract the witness set from main_data
ring R=(complex,40,i),(x,y,z),dp;
list P=getWitnessSet();
P;
==> [1]:
==>    [1]:
==> (-0.3133771220271943131647016495574044329905+i*0.149319306273712443760775\
   4574813232743642)
==>    [2]:
==> (-0.3133771220271943131647016495574044329905+i*0.149319306273712443760775\
   4574813232743642)
==>    [3]:
==> (-0.3236835696859500265582033518080813162881+i*0.120979768405647429138007\
   6875729169467213)
==> [2]:
==>    [1]:
==> (-0.3212684582929880880341864797248958078041+i*0.136898103093073877150095\
   9867882795547697)
==>    [2]:
==> (-0.3086806001515909537016402824032076176217+i*0.164519057109862805358972\
   1143137415991986)
==>    [3]:
==> (-0.3212684582929880880341864797248958078041+i*0.136898103093073877150095\
   9867882795547697)
==> [3]:
==>    [1]:
==> (-0.4654825244434110956527195748061499827706+i*0.255549561350552455480763\
   9314453979096486)
==>    [2]:
==> (-0.1245166868485723086630707453638664230182+i*0.081799393365234509591583\
   34524040056787389)
==>    [3]:
==> (-0.1245166868485723086630707453638664230182+i*0.081799393365234509591583\
   34524040056787389)


Top Back: gaussColWithoutPerm Forward: writeBertiniInput FastBack: FastForward: Up: recover_lib Top: Singular Manual Contents: Table of Contents Index: Index About: About this document
            User manual for Singular version 4-2-1p3, 2021, generated by texi2html.