The user can specify the number of decimal places desired to refine solutions.
Let us compute a square root with a working precision of 64 decimal places:
i1 : R = CC[x,y]; S = {x^2 - 1/3, x*y - 1}; roots = solveSystem(S);
using temporary files /var/folders/40/dy88l5qd361391m_3v2m51bm0000gn/T/M2-44686-0/0PHCinput and /var/folders/40/dy88l5qd361391m_3v2m51bm0000gn/T/M2-44686-0/0PHCoutput
|
i4 : r0 = roots#0#Coordinates#1 o4 = 1.73205080756888 o4 : CC (of precision 53) |
i5 : newRoots = refineSolutions(S,roots,64)
writing input system to /var/folders/40/dy88l5qd361391m_3v2m51bm0000gn/T/M2-44686-0/1PHCinput
appending solutions to /var/folders/40/dy88l5qd361391m_3v2m51bm0000gn/T/M2-44686-0/1PHCinput
preparing input data for phc -v in /var/folders/40/dy88l5qd361391m_3v2m51bm0000gn/T/M2-44686-0/3PHCbatch
using temporary file /var/folders/40/dy88l5qd361391m_3v2m51bm0000gn/T/M2-44686-0/2PHCoutput for storing refined solutions
solutions in Maple format in /var/folders/40/dy88l5qd361391m_3v2m51bm0000gn/T/M2-44686-0/5PHCsolutions
o5 = {{.57735, 1.73205}, {-.57735, -1.73205}}
o5 : List
|
i6 : newRoots#0 -- recall that solutions are of type Point
o6 = {.57735, 1.73205}
o6 : Point
|
i7 : r1 = newRoots#0#Coordinates#1 o7 = 1.73205080756887735414922460641658082304678611432812598193257967 o7 : CC (of precision 213) |