#12832: Bug in scipy on power7
----------------------------+-----------------------------------------------
Reporter: fbissey | Owner: tbd
Type: defect | Status: new
Priority: major | Milestone: sage-5.0
Component: packages | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
----------------------------+-----------------------------------------------
Comment (by zimmerma):
I investigated the first error in {{{optimize.py}}}:
{{{
c=vector(RDF,[-4,-5])
G=matrix(RDF,[[2,1],[1,2],[-1,0],[0,-1]])
h=vector(RDF,[3,3,0,0])
from cvxopt.base import matrix as m
from cvxopt import solvers
solvers.options['show_progress']=True
c_=m(c.base_extend(RDF).numpy())
G_=m(G.base_extend(RDF).numpy())
h_=m(h.base_extend(RDF).numpy())
sol=solvers.lp(c_,G_,h_,solver=None)
}}}
The first numerical difference comes from the call
{{{base.gemm(spmatrix(W['di'], list(range(ml)), list(range(ml))), G,
F['Gs'], partial = True)}}} in {{{local/lib/python2.7/site-
packages/cvxopt/misc.py}}}, where the inputs are:
{{{
W['di']= [ 3.06e+00]
[ 3.32e+00]
[ 1.29e+00]
[ 1.05e+00]
ml= 4
G= [ 2.00e+00 1.00e+00]
[ 1.00e+00 2.00e+00]
[-1.00e+00 0.00e+00]
[ 0.00e+00 -1.00e+00]
}}}
and the expected output is:
{{{
F['Gs']= [ 6.11e+00 3.06e+00]
[ 3.32e+00 6.63e+00]
[-1.29e+00 0.00e+00]
[ 0.00e+00 -1.05e+00]
}}}
whereas we get:
{{{
F['Gs']= [ 8.11e+00 4.06e+00]
[ 4.32e+00 8.63e+00]
[-2.29e+00 0.00e+00]
[ 0.00e+00 -2.05e+00]
}}}
Does anybody know where the {{{base.gemm}}} function is defined?
Paul
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/12832#comment:2>
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 post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
