I have heavily altered the PHCpack interface for sage, and posted the relevant files as sage-trac 399:
http://www.sagemath.org:9002/sage_trac/ticket/399 Comments are very welcome. In order to get it working, I have temporarily removed some potential functionality - the ability to pass in arbitrary command lists to phc. On the positive side, what is there works, and I have added some parsing of the raw phc output to get at the final complex solutions of the system. Here is part of the Examples section to give an idea of how it works: EXAMPLES: sage: from sage.interfaces.phc import phc sage: R.<x,y> = PolynomialRing(QQ,2) sage: testsys = [x^2 + 1, x*y - 1] sage: v = phc.blackbox(testsys, R) sage: v.solutions() [[1.00000000000000*I, -1.00000000000000*I], [-1.00000000000000*I, 1.00000000000000*I]] sage: v.solution_dicts() [{y: 1.00000000000000*I, x: -1.00000000000000*I}, {y: -1.00000000000000*I, x: 1.00000000000000*I}] sage: residuals = [[test_equation.subs(sol) for test_equation in testsys] for sol in v.solution_dicts()] sage: residuals [[0, 0], [0, 0]] -Marshall Hampton --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~----------~----~----~----~------~----~------~--~---
