#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.


Reply via email to