|
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.8562492735377351919600869059664108051841-i*1.1249221758539982174957025\
45609735327269)
==> [2]:
==> (-0.3361013656343670777714001809349914480059+i*0.993421096862434190679193\
5550838062315518)
==> [3]:
==> (-0.3361013656343670777714001809349914480059+i*0.993421096862434190679193\
5550838062315518)
==> [2]:
==> [1]:
==> (-0.0945774957829047994509151973995708440987+i*0.411695549822856774098046\
8649115788679681)
==> [2]:
==> (1.353598559419160990961739112360575708109-i*0.24748564301879171662432819\
11979594330299)
==> [3]:
==> (-0.0945774957829047994509151973995708440987+i*0.411695549822856774098046\
8649115788679681)
==> [3]:
==> [1]:
==> (0.5102396697678335186091995880490086306541-i*0.0322180255640324185830203\
7742114636069658)
==> [2]:
==> (0.5102396697678335186091995880490086306541-i*0.0322180255640324185830203\
7742114636069658)
==> [3]:
==> (-0.2779459066587645175862754745812237767077+i*0.598150303028424892723585\
2887276312054417)
|
|