#14628: Zero solution does not result in zero
-----------------------------+----------------------------------------------
Reporter: gagern | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-5.10
Component: symbolics | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
-----------------------------+----------------------------------------------
Comment (by kcrisman):
Unfortunately, I can't debug this more now because of the annoying thing
that Maxima doesn't allow copying and pasting of more than a certain
number of characters, and I can't remember the workaround right now. What
I *can* say is that it is something about our form of the imaginary part
of z:
{{{
(%i1) display2d:false;
(%o1) false
(%i2) declare(phi,real);
(%o2) done
(%i3)
f:(%i+1)*sqrt(-%i)/((%e^(%i*phi)-1)*((%i+1)/(%e^(%i*phi)-1)-%i*%e^(-%i*phi)));
(%o3) sqrt(-%i)*(%i+1)/(((%i+1)/(%e^(%i*phi)-1)-%i*%e^-(%i*phi))
*(%e^(%i*phi)-1))
(%i4) solve(imagpart(f),phi);
(%o4) [sin(phi) = -sqrt(-cos(phi)^2+cos(phi)+1),
sin(phi) = sqrt(-cos(phi)^2+cos(phi)+1)]
}}}
on which `to_poly_solve` also works correctly.
{{{
(%i5) load(to_poly_solve);
Loading maxima-grobner $Revision: 1.6 $ $Date: 2009-06-02 07:49:49 $
(%o5)
"/Users/.../sage-5.9.rc1/local/share/maxima/5.29.1/share/to_poly_solve/to_poly_solve.mac"
(%i6) %solve(imagpart(f),phi);
(%o6) %union([phi = 2*%pi*%z8+%pi/2])
}}}
It ''is'' strange that the error is so dramatically wrong.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/14628#comment:3>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.
