def load_standard_laurent_system():
"""
Returns the Laurent polynomials stored in the system container
for standard double precision arithmetic.
"""
from phcpy.phcpy2c3 import py2c_syscon_number_of_standard_Laurentials
from phcpy.phcpy2c3 import py2c_syscon_load_standard_Laurential
dim = py2c_syscon_number_of_standard_Laurentials()
result = []
for ind in range(1, dim + 1):
result.append(py2c_syscon_load_standard_Laurential(ind))
return result
def load_standard_laurent_system():
"""
Returns the Laurent polynomials stored in the system container
for standard double precision arithmetic.
"""
from phcpy.phcpy2c3 import py2c_syscon_number_of_standard_Laurentials
from phcpy.phcpy2c3 import py2c_syscon_load_standard_Laurential
dim = py2c_syscon_number_of_standard_Laurentials()
result = []
for ind in range(1, dim+1):
result.append(py2c_syscon_load_standard_Laurential(ind))
return result
def initial_form(pols, normal):
r"""
Returns the initial form of the polynomials in *pols*
with respect to the inner normal with coordinates in *normal*.
"""
from phcpy.phcpy2c3 import py2c_syscon_clear_standard_Laurent_system
from phcpy.phcpy2c3 \
import py2c_syscon_initialize_number_of_standard_Laurentials
from phcpy.phcpy2c3 import py2c_syscon_store_standard_Laurential
from phcpy.phcpy2c3 import py2c_syscon_load_standard_Laurential
from phcpy.phcpy2c3 import py2c_giftwrap_initial_form
py2c_syscon_clear_standard_Laurent_system()
dim = max(len(pols), len(normal))
py2c_syscon_initialize_number_of_standard_Laurentials(dim)
for i in range(len(pols)):
nchar = len(pols[i])
fail = py2c_syscon_store_standard_Laurential(nchar, dim, i+1, pols[i])
strnrm = str(tuple(normal))
fail = py2c_giftwrap_initial_form(len(normal), len(strnrm), strnrm)
result = []
for i in range(len(pols)):
result.append(py2c_syscon_load_standard_Laurential(i+1))
return result
def initial_form(pols, normal):
r"""
Returns the initial form of the polynomials in *pols*
with respect to the inner normal with coordinates in *normal*.
"""
from phcpy.phcpy2c3 import py2c_syscon_clear_standard_Laurent_system
from phcpy.phcpy2c3 \
import py2c_syscon_initialize_number_of_standard_Laurentials
from phcpy.phcpy2c3 import py2c_syscon_store_standard_Laurential
from phcpy.phcpy2c3 import py2c_syscon_load_standard_Laurential
from phcpy.phcpy2c3 import py2c_giftwrap_initial_form
py2c_syscon_clear_standard_Laurent_system()
dim = max(len(pols), len(normal))
py2c_syscon_initialize_number_of_standard_Laurentials(dim)
for i in range(len(pols)):
nchar = len(pols[i])
fail = py2c_syscon_store_standard_Laurential(nchar, dim, i + 1,
pols[i])
strnrm = str(tuple(normal))
fail = py2c_giftwrap_initial_form(len(normal), len(strnrm), strnrm)
result = []
for i in range(len(pols)):
result.append(py2c_syscon_load_standard_Laurential(i + 1))
return result
