#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 could reduce the issue to the following:
{{{
import numpy
from cvxopt.base import spmatrix
from cvxopt.base import matrix as m
from cvxopt import umfpack
RealNumber=float
W = [ 3.06e+00, 3.32e+00, 1.29e+00, 1.05e+00]
ml=4
W2 = spmatrix(W, list(range(ml)), list(range(ml)))
G= m([[2.00e+00, 1.00e+00, -1.00e+00, 0.00e+00],
[1.00e+00, 2.00e+00, 0.00e+00, -1.00e+00]])
F = m([[0.0,0.0,0.0,0.0],[0.0,0.0,0.0,0.0]])
from cvxopt.base import gemm
import cvxopt
cvxopt.base.gemm(W2,G,F,partial = True)
print F
cvxopt.base.gemm(W2,G,F,partial = True)
print F
}}}
The expected output is twice the same matrix F:
{{{
[ 6.12e+00 3.06e+00]
[ 3.32e+00 6.64e+00]
[-1.29e+00 0.00e+00]
[ 0.00e+00 -1.05e+00]
[ 6.12e+00 3.06e+00]
[ 3.32e+00 6.64e+00]
[-1.29e+00 0.00e+00]
[ 0.00e+00 -1.05e+00]
}}}
On silius with sage-5.0.rc1 I get instead:
{{{
[ 6.12e+00 3.06e+00]
[ 3.32e+00 6.64e+00]
[-1.29e+00 0.00e+00]
[ 0.00e+00 -1.05e+00]
[ 1.22e+01 6.12e+00]
[ 6.64e+00 1.33e+01]
[-2.58e+00 0.00e+00]
[ 0.00e+00 -2.10e+00]
}}}
thus it seems the product is **accumulated** to F.
Paul
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/12832#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 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.
