#12832: Bug in cvxopt on power7
----------------------------+-----------------------------------------------
Reporter: fbissey | Owner: tbd
Type: defect | Status: new
Priority: major | Milestone: sage-5.7
Component: packages | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
----------------------------+-----------------------------------------------
Description changed by fbissey:
Old description:
> There are two doctests related to scipy failing on power7
> {{{
> ./sage -t -long "devel/sage/sage/numerical/optimize.py"
> sage -t -long "devel/sage/sage/numerical/optimize.py"
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/numerical/optimize.py",
> line 469:
> sage: sol['x']
> Expected:
> (0.999..., 1.000...)
> Got nothing
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/numerical/optimize.py",
> line 479:
> sage: sol['x']
> Expected:
> (45.000000..., 6.2499999...3, 1.00000000...)
> Got nothing
> GLPK Simplex Optimizer, v4.44
> 6 rows, 3 columns, 8 non-zeros
> Preprocessing...
> 2 rows, 2 columns, 4 non-zeros
> Scaling...
> A: min|aij| = 2.400e+01 max|aij| = 5.000e+01 ratio = 2.083e+00
> GM: min|aij| = 8.128e-01 max|aij| = 1.230e+00 ratio = 1.514e+00
> EQ: min|aij| = 6.606e-01 max|aij| = 1.000e+00 ratio = 1.514e+00
> Constructing initial basis...
> Size of triangular part = 2
> * 0: obj = -5.100000000e+01 infeas = 0.000e+00 (0)
> * 1: obj = -5.225000000e+01 infeas = 0.000e+00 (0)
> OPTIMAL SOLUTION FOUND
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/numerical/optimize.py",
> line 564:
> sage: find_fit(data, model)
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_7[5]>", line 1, in <module>
> find_fit(data, model)###line 564:
> sage: find_fit(data, model)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/numerical/optimize.py", line 655, in find_fit
> estimated_params, d = leastsq(error_function, initial_guess, args
> = (x_data, y_data))
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/scipy/optimize/minpack.py", line 324, in leastsq
> raise errors[info][1](errors[info][0])
> TypeError: Improper input parameters.
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/numerical/optimize.py",
> line 570:
> sage: fit = find_fit(data, f, parameters = [a, b, c], variables =
> [x], solution_dict = True)
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_7[7]>", line 1, in <module>
> fit = find_fit(data, f, parameters = [a, b, c], variables = [x],
> solution_dict = True)###line 570:
> sage: fit = find_fit(data, f, parameters = [a, b, c], variables =
> [x], solution_dict = True)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/numerical/optimize.py", line 655, in find_fit
> estimated_params, d = leastsq(error_function, initial_guess, args
> = (x_data, y_data))
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/scipy/optimize/minpack.py", line 324, in leastsq
> raise errors[info][1](errors[info][0])
> TypeError: Improper input parameters.
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/numerical/optimize.py",
> line 571:
> sage: fit[a], fit[b], fit[c]
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_7[8]>", line 1, in <module>
> fit[a], fit[b], fit[c]###line 571:
> sage: fit[a], fit[b], fit[c]
> NameError: name 'fit' is not defined
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/numerical/optimize.py",
> line 577:
> sage: find_fit(dataprime, a * x * log(b * x), parameters = [a, b],
> variables = [x])
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_7[10]>", line 1, in <module>
> find_fit(dataprime, a * x * log(b * x), parameters = [a, b],
> variables = [x])###line 577:
> sage: find_fit(dataprime, a * x * log(b * x), parameters = [a, b],
> variables = [x])
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/numerical/optimize.py", line 655, in find_fit
> estimated_params, d = leastsq(error_function, initial_guess, args
> = (x_data, y_data))
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/scipy/optimize/minpack.py", line 324, in leastsq
> raise errors[info][1](errors[info][0])
> TypeError: Improper input parameters.
> **********************************************************************
> }}}
> and
> {{{
> ./sage -t --long -force_lib devel/sage/sage/plot/plot3d/list_plot3d.py
> sage -t --long -force_lib "devel/sage/sage/plot/plot3d/list_plot3d.py"
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/plot/plot3d/list_plot3d.py",
> line 99:
> sage: list_plot3d(m, texture='yellow',
> interpolation_type='spline',frame_aspect_ratio=[1,1,1/3])
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_1[9]>", line 1, in <module>
> list_plot3d(m, texture='yellow',
> interpolation_type='spline',frame_aspect_ratio=[Integer(1),Integer(1),Integer(1)/Integer(3)])###line
> 99:
> sage: list_plot3d(m, texture='yellow',
> interpolation_type='spline',frame_aspect_ratio=[1,1,1/3])
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 174, in list_plot3d
> return list_plot3d_tuples(l,interpolation_type,texture,**kwds)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 458, in
> list_plot3d_tuples
> s=interpolate.bisplrep(x,y,z,[int(1)]*len(x),xmin,xmax,ymin,ymax,kx=kx,ky=ky,s=s)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/scipy/interpolate/fitpack.py", line 873, in bisplrep
> tx,ty,nxest,nyest,wrk,lwrk1,lwrk2)
> ValueError: negative dimensions are not allowed
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/plot/plot3d/list_plot3d.py",
> line 106:
> sage: list_plot3d(m, texture='yellow', interpolation_type='spline',
> degree=5, frame_aspect_ratio=[1,1,1/3])
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_1[10]>", line 1, in <module>
> list_plot3d(m, texture='yellow', interpolation_type='spline',
> degree=Integer(5),
> frame_aspect_ratio=[Integer(1),Integer(1),Integer(1)/Integer(3)])###line
> 106:
> sage: list_plot3d(m, texture='yellow', interpolation_type='spline',
> degree=5, frame_aspect_ratio=[1,1,1/3])
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 174, in list_plot3d
> return list_plot3d_tuples(l,interpolation_type,texture,**kwds)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 458, in
> list_plot3d_tuples
> s=interpolate.bisplrep(x,y,z,[int(1)]*len(x),xmin,xmax,ymin,ymax,kx=kx,ky=ky,s=s)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/scipy/interpolate/fitpack.py", line 873, in bisplrep
> tx,ty,nxest,nyest,wrk,lwrk1,lwrk2)
> ValueError: negative dimensions are not allowed
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/plot/plot3d/list_plot3d.py",
> line 271:
> sage: show(list_plot3d([[1, 1, 1, 1], [1, 2, 1, 2], [1, 1, 3, 1], [1,
> 2, 1, 4]],interpolation_type='spline'))
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_3[4]>", line 1, in <module>
> show(list_plot3d([[Integer(1), Integer(1), Integer(1),
> Integer(1)], [Integer(1), Integer(2), Integer(1), Integer(2)],
> [Integer(1), Integer(1), Integer(3), Integer(1)], [Integer(1),
> Integer(2), Integer(1), Integer(4)]],interpolation_type='spline'))###line
> 271:
> sage: show(list_plot3d([[1, 1, 1, 1], [1, 2, 1, 2], [1, 1, 3, 1], [1,
> 2, 1, 4]],interpolation_type='spline'))
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 195, in list_plot3d
> return list_plot3d_array_of_arrays(v, interpolation_type,texture,
> **kwds)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 274, in
> list_plot3d_array_of_arrays
> G = list_plot3d(m,interpolation_type,texture, **kwds)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 174, in list_plot3d
> return list_plot3d_tuples(l,interpolation_type,texture,**kwds)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 458, in
> list_plot3d_tuples
> s=interpolate.bisplrep(x,y,z,[int(1)]*len(x),xmin,xmax,ymin,ymax,kx=kx,ky=ky,s=s)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/scipy/interpolate/fitpack.py", line 873, in bisplrep
> tx,ty,nxest,nyest,wrk,lwrk1,lwrk2)
> ValueError: negative dimensions are not allowed
> **********************************************************************
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/devel/sage/sage/plot/plot3d/list_plot3d.py",
> line 342:
> sage: list_plot3d(m, texture='yellow',
> interpolation_type='spline',frame_aspect_ratio=[1,1,1/3])
> Exception raised:
> Traceback (most recent call last):
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1231, in run_one_test
> self.run_one_example(test, example, filename, compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/sagedoctest.py",
> line 38, in run_one_example
> OrigDocTestRunner.run_one_example(self, test, example, filename,
> compileflags)
> File
> "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/bin/ncadoctest.py",
> line 1172, in run_one_example
> compileflags, 1) in test.globs
> File "<doctest __main__.example_4[5]>", line 1, in <module>
> list_plot3d(m, texture='yellow',
> interpolation_type='spline',frame_aspect_ratio=[Integer(1),Integer(1),Integer(1)/Integer(3)])###line
> 342:
> sage: list_plot3d(m, texture='yellow',
> interpolation_type='spline',frame_aspect_ratio=[1,1,1/3])
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 174, in list_plot3d
> return list_plot3d_tuples(l,interpolation_type,texture,**kwds)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/sage/plot/plot3d/list_plot3d.py", line 458, in
> list_plot3d_tuples
> s=interpolate.bisplrep(x,y,z,[int(1)]*len(x),xmin,xmax,ymin,ymax,kx=kx,ky=ky,s=s)
> File "/hpc/scratch/frb15/sandbox/sage-5.0.beta10/local/lib/python
> /site-packages/scipy/interpolate/fitpack.py", line 873, in bisplrep
> tx,ty,nxest,nyest,wrk,lwrk1,lwrk2)
> ValueError: negative dimensions are not allowed
> **********************************************************************
> }}}
New description:
There are two doctests related to cvxopt failing on power7
{{{
sage -t --long -force_lib "devel/sage/doc/en/numerical_sage/cvxopt.rst"
**********************************************************************
File
"/hpc/scratch/frb15/sandbox/sage-5.7.beta4/devel/sage/doc/en/numerical_sage/cvxopt.rst",
line 129:
sage: print sol['x'] # ... below since can get -00 or +00
depending on architecture
Expected:
[ 1.00e...00]
[ 1.00e+00]
Got:
[ 9.87e-01]
[ 9.92e-01]
<BLANKLINE>
**********************************************************************
}}}
and
{{{
sage -t -long "devel/sage-main/sage/numerical/optimize.py"
**********************************************************************
File "/hpc/scratch/frb15/sandbox/sage-5.7.beta4/devel/sage-
main/sage/numerical/optimize.py", line 515:
sage: sol['x']
Expected:
(0.999..., 1.000...)
Got nothing
**********************************************************************
File "/hpc/scratch/frb15/sandbox/sage-5.7.beta4/devel/sage-
main/sage/numerical/optimize.py", line 525:
sage: sol['x']
Expected:
(45.000000..., 6.2499999...3, 1.00000000...)
Got nothing
GLPK Simplex Optimizer, v4.44
6 rows, 3 columns, 8 non-zeros
Preprocessing...
2 rows, 2 columns, 4 non-zeros
Scaling...
A: min|aij| = 2.400e+01 max|aij| = 5.000e+01 ratio = 2.083e+00
GM: min|aij| = 8.128e-01 max|aij| = 1.230e+00 ratio = 1.514e+00
EQ: min|aij| = 6.606e-01 max|aij| = 1.000e+00 ratio = 1.514e+00
Constructing initial basis...
Size of triangular part = 2
* 0: obj = -5.100000000e+01 infeas = 0.000e+00 (0)
* 1: obj = -5.225000000e+01 infeas = 0.000e+00 (0)
OPTIMAL SOLUTION FOUND
**********************************************************************
}}}
--
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/12832#comment:14>
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.