#11705: Port Sage to SUSE Linux Power 7 (ppc64).
---------------------------------------------------+------------------------
Reporter: was | Owner: drkirkby
Type: enhancement | Status: new
Priority: major | Milestone: sage-5.7
Component: porting | Resolution:
Keywords: sd32 sd35.5 | Work issues:
Report Upstream: N/A | Reviewers:
Authors: Paul Zimmermann, Jeroen Demeyer | Merged in:
Dependencies: #12829, #12832, #14098 | Stopgaps:
---------------------------------------------------+------------------------
Comment (by fbissey):
That's repeating from earlier comments but this is now a long thread:
{{{
./sage -t --long -force_lib devel/sage/sage/functions/special.py
sage -t --long -force_lib "devel/sage/sage/functions/special.py"
;;;
;;; Stack overflow.
;;; Jumping to the outermost toplevel prompt
;;;
**********************************************************************
File
"/hpc/scratch/frb15/sandbox/sage-5.7.beta4/devel/sage/sage/functions/special.py",
line 1614:
sage: elliptic_f(RR(pi/2), 0.5)
Expected:
1.8540746773
Got:
sage165
;;;
;;; Stack overflow.
;;; Jumping to the outermost toplevel prompt
;;;
**********************************************************************
File
"/hpc/scratch/frb15/sandbox/sage-5.7.beta4/devel/sage/sage/functions/special.py",
line 1623:
sage: elliptic_f(RR(pi/2), 0.5)
Expected:
1.8540746773
Got:
sage167
**********************************************************************
2 items had failures:
1 of 5 in __main__.example_44
1 of 5 in __main__.example_45
***Test Failed*** 2 failures.
For whitespace errors, see the file
/hpc/home/frb15/.sage//tmp/special_36988.py
[31.9 s]
----------------------------------------------------------------------
}}}
A maxima/ecl problem, I'll try to get more details so we can forward this
to the ecl people (or maxima but I suspect the fault is in ecl).
{{{
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>
**********************************************************************
}}}
could be considered noise, but it is of a worrying kind especially when
combined with the next one:
{{{
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
**********************************************************************
}}}
I have digged deeper at what is happening under the hood in the first case
comparing with 5.6 on my imac.
The imac run
{{{
sage: c=vector(RDF,[-4,-5])
sage: G=matrix(RDF,[[2,1],[1,2],[-1,0],[0,-1]])
sage: h=vector(RDF,[3,3,0,0])
sage: sol=linear_program(c,G,h)
sage: sol['x']
(0.999999939767, 1.00000001047)
sage: from cvxopt.base import matrix as m
sage: from cvxopt import solvers
sage: solvers.options['show_progress']=True
sage: c_=m(c.base_extend(RDF).numpy())
sage: G_=m(G.base_extend(RDF).numpy())
sage: h_=m(h.base_extend(RDF).numpy())
sage: solvers.lp(c_,G_,h_,solver=None)
pcost dcost gap pres dres k/t
0: -8.1000e+00 -1.8300e+01 4e+00 0e+00 8e-01 1e+00
1: -8.8055e+00 -9.4357e+00 2e-01 1e-16 4e-02 3e-02
2: -8.9981e+00 -9.0049e+00 2e-03 1e-16 5e-04 4e-04
3: -9.0000e+00 -9.0000e+00 2e-05 1e-16 5e-06 4e-06
4: -9.0000e+00 -9.0000e+00 2e-07 1e-16 5e-08 4e-08
Optimal solution found.
{'status': 'optimal', 'dual slack': 2.7986120477842287e-08, 'iterations':
4, 'residual as primal infeasibility certificate': None, 'relative gap':
2.7257517543273638e-08, 'dual objective': -9.00000049248484, 'residual as
dual infeasibility certificate': None, 'gap': 2.453176527488781e-07, 's':
<4x1 matrix, tc='d'>, 'primal infeasibility': 1.0136264998440638e-16,
'dual infeasibility': 4.812169483742963e-08, 'primal objective':
-8.99999981140672, 'primal slack': 3.929722790979621e-08, 'y': <0x1
matrix, tc='d'>, 'x': <2x1 matrix, tc='d'>, 'z': <4x1 matrix, tc='d'>}
}}}
The power7 one
{{{
sage: c=vector(RDF,[-4,-5])
sage: G=matrix(RDF,[[2,1],[1,2],[-1,0],[0,-1]])
sage: h=vector(RDF,[3,3,0,0])
sage: sol=linear_program(c,G,h)
sage: sol['x']
sage: from cvxopt.base import matrix as m
sage: from cvxopt import solvers
sage: solvers.options['show_progress']=True
sage: c_=m(c.base_extend(RDF).numpy())
sage: G_=m(G.base_extend(RDF).numpy())
sage: h_=m(h.base_extend(RDF).numpy())
sage: solvers.lp(c_,G_,h_,solver=None)
pcost dcost gap pres dres k/t
0: -8.1000e+00 -1.8300e+01 4e+00 0e+00 8e-01 1e+00
1: -8.2446e+00 -9.9047e+00 3e-01 8e-02 1e-01 6e-02
2: -7.6587e+00 -8.5347e+00 9e-02 1e+00 5e-02 2e-02
3: -7.0226e+00 -7.7110e+00 6e-02 1e+01 1e-01 1e-02
4: -6.4985e+00 -7.0882e+00 5e-02 2e+01 2e-01 9e-03
5: -5.5876e+00 -6.0350e+00 4e-02 2e+02 3e-01 5e-03
6: -4.8663e+00 -5.2180e+00 3e-02 3e+02 4e-01 3e-03
7: -4.1532e+00 -4.4125e+00 2e-02 2e+03 5e-01 2e-03
8: -3.6054e+00 -3.8041e+00 2e-02 4e+03 6e-01 1e-03
9: -3.1924e+00 -3.3415e+00 1e-02 2e+04 6e-01 8e-04
10: -2.8113e+00 -2.9258e+00 1e-02 4e+04 7e-01 5e-04
11: -2.5706e+00 -2.6588e+00 8e-03 2e+05 7e-01 4e-04
12: -2.2854e+00 -2.3535e+00 6e-03 4e+05 7e-01 3e-04
13: -2.1174e+00 -2.1709e+00 5e-03 2e+06 8e-01 3e-04
14: -1.9020e+00 -1.9438e+00 4e-03 4e+06 8e-01 2e-04
15: -1.7759e+00 -1.8093e+00 3e-03 2e+07 8e-01 2e-04
16: -1.6170e+00 -1.6435e+00 2e-03 4e+07 8e-01 1e-04
17: -1.5242e+00 -1.5457e+00 2e-03 2e+08 8e-01 1e-04
18: -1.4088e+00 -1.4261e+00 2e-03 4e+08 8e-01 1e-04
19: -1.3438e+00 -1.3581e+00 1e-03 2e+09 8e-01 1e-04
20: -1.2614e+00 -1.2731e+00 1e-03 5e+09 9e-01 8e-05
21: -1.2205e+00 -1.2303e+00 1e-03 2e+10 9e-01 8e-05
22: -1.1641e+00 -1.1723e+00 8e-04 5e+10 9e-01 7e-05
23: -1.1457e+00 -1.1526e+00 7e-04 3e+11 9e-01 6e-05
24: -1.1120e+00 -1.1178e+00 6e-04 6e+11 9e-01 5e-05
25: -1.1175e+00 -1.1225e+00 5e-04 3e+12 9e-01 5e-05
26: -1.1071e+00 -1.1114e+00 4e-04 7e+12 9e-01 5e-05
27: -1.1440e+00 -1.1477e+00 4e-04 3e+13 9e-01 5e-05
28: -1.1644e+00 -1.1676e+00 4e-04 8e+13 9e-01 4e-05
29: -1.2555e+00 -1.2584e+00 3e-04 3e+14 9e-01 4e-05
30: -1.3345e+00 -1.3370e+00 3e-04 1e+15 9e-01 4e-05
31: -1.5684e+00 -1.5705e+00 3e-04 4e+15 8e-01 5e-05
32: -1.8504e+00 -1.8522e+00 3e-04 1e+16 8e-01 5e-05
33: -3.2927e+00 -3.2940e+00 4e-04 5e+16 6e-01 6e-05
34: -7.2895e+00 -7.2897e+00 4e-04 1e+17 2e-01 6e-05
35: -8.2349e+00 -8.2350e+00 2e-04 2e+17 8e-02 3e-05
36: -8.1449e+00 -8.1450e+00 3e-04 3e+18 9e-02 2e-05
37: -8.1562e+00 -8.1563e+00 2e-04 5e+18 9e-02 3e-05
38: -8.5850e+00 -8.5850e+00 2e-04 4e+19 5e-02 1e-05
39: -8.7398e+00 -8.7398e+00 1e-04 7e+19 3e-02 7e-06
40: -8.5796e+00 -8.5796e+00 1e-04 2e+21 5e-02 9e-06
41: -8.5844e+00 -8.5844e+00 9e-05 2e+21 5e-02 1e-05
42: -8.7635e+00 -8.7635e+00 5e-05 1e+22 3e-02 4e-06
43: -8.6210e+00 -8.6210e+00 6e-05 3e+23 4e-02 5e-06
44: -8.6127e+00 -8.6127e+00 4e-05 5e+23 4e-02 6e-06
45: -8.7308e+00 -8.7308e+00 4e-05 4e+24 3e-02 3e-06
46: -8.8323e+00 -8.8323e+00 2e-05 7e+24 2e-02 2e-06
47: -8.7174e+00 -8.7174e+00 3e-05 2e+26 3e-02 2e-06
48: -8.6973e+00 -8.6973e+00 3e-05 2e+26 3e-02 3e-06
49: -8.8058e+00 -8.8058e+00 2e-05 2e+27 2e-02 1e-06
50: -8.8546e+00 -8.8546e+00 1e-05 4e+27 2e-02 9e-07
51: -8.7514e+00 -8.7514e+00 2e-05 1e+29 3e-02 1e-06
52: -8.7405e+00 -8.7405e+00 1e-05 1e+29 3e-02 1e-06
53: -8.8161e+00 -8.8161e+00 1e-05 1e+30 2e-02 7e-07
54: -8.8781e+00 -8.8781e+00 7e-06 2e+30 1e-02 5e-07
55: -8.7917e+00 -8.7917e+00 1e-05 5e+31 2e-02 6e-07
56: -8.7732e+00 -8.7732e+00 7e-06 7e+31 3e-02 8e-07
57: -8.8288e+00 -8.8288e+00 7e-06 7e+32 2e-02 5e-07
58: -8.9006e+00 -8.9006e+00 4e-06 1e+33 1e-02 3e-07
59: -8.8275e+00 -8.8275e+00 6e-06 2e+34 2e-02 3e-07
60: -8.8035e+00 -8.8035e+00 4e-06 3e+34 2e-02 4e-07
61: -8.8360e+00 -8.8360e+00 5e-06 3e+35 2e-02 3e-07
62: -8.9251e+00 -8.9251e+00 1e-06 4e+35 8e-03 1e-07
63: -8.8661e+00 -8.8661e+00 3e-06 8e+36 1e-02 2e-07
64: -8.8364e+00 -8.8364e+00 3e-06 1e+37 2e-02 2e-07
65: -8.8262e+00 -8.8262e+00 4e-06 2e+38 2e-02 2e-07
66: -8.9471e+00 -8.9471e+00 6e-07 2e+38 6e-03 8e-08
67: -8.9304e+00 -8.9304e+00 7e-07 3e+38 8e-03 8e-08
68: -8.8850e+00 -8.8850e+00 1e-06 5e+39 1e-02 9e-08
69: -8.8583e+00 -8.8583e+00 1e-06 1e+40 2e-02 1e-07
70: -8.8071e+00 -8.8071e+00 3e-06 2e+41 2e-02 1e-07
71: -8.9335e+00 -8.9335e+00 6e-07 2e+41 7e-03 1e-07
72: -8.9106e+00 -8.9106e+00 6e-07 3e+41 1e-02 1e-07
73: -8.8518e+00 -8.8518e+00 1e-06 5e+42 2e-02 9e-08
74: -8.8289e+00 -8.8289e+00 9e-07 8e+42 2e-02 1e-07
75: -8.8343e+00 -8.8343e+00 1e-06 1e+44 2e-02 8e-08
76: -8.9463e+00 -8.9463e+00 2e-07 1e+44 6e-03 3e-08
77: -8.9336e+00 -8.9336e+00 2e-07 3e+44 7e-03 3e-08
78: -8.9035e+00 -8.9035e+00 4e-07 4e+45 1e-02 3e-08
79: -8.8828e+00 -8.8828e+00 4e-07 8e+45 1e-02 4e-08
80: -8.8369e+00 -8.8369e+00 8e-07 2e+47 2e-02 4e-08
81: -8.8225e+00 -8.8225e+00 7e-07 2e+47 2e-02 7e-08
82: -8.8762e+00 -8.8762e+00 8e-07 3e+48 1e-02 5e-08
83: -8.9445e+00 -8.9445e+00 2e-07 3e+48 6e-03 2e-08
84: -8.9062e+00 -8.9062e+00 4e-07 5e+49 1e-02 2e-08
85: -8.8852e+00 -8.8852e+00 4e-07 9e+49 1e-02 3e-08
86: -8.8567e+00 -8.8567e+00 9e-07 2e+51 2e-02 4e-08
87: -8.9491e+00 -8.9491e+00 2e-07 1e+51 6e-03 2e-08
88: -8.9309e+00 -8.9309e+00 2e-07 2e+51 8e-03 2e-08
89: -8.8818e+00 -8.8818e+00 5e-07 4e+52 1e-02 3e-08
90: -8.8660e+00 -8.8660e+00 4e-07 7e+52 1e-02 4e-08
91: -8.8989e+00 -8.8989e+00 4e-07 7e+53 1e-02 3e-08
92: -8.9524e+00 -8.9524e+00 1e-07 9e+53 5e-03 1e-08
93: -8.9187e+00 -8.9187e+00 2e-07 1e+55 9e-03 1e-08
94: -8.9004e+00 -8.9004e+00 2e-07 3e+55 1e-02 2e-08
95: -8.8722e+00 -8.8722e+00 5e-07 5e+56 1e-02 2e-08
96: -8.9530e+00 -8.9530e+00 1e-07 4e+56 5e-03 1e-08
97: -8.9359e+00 -8.9359e+00 1e-07 6e+56 7e-03 2e-08
98: -8.8923e+00 -8.8923e+00 3e-07 1e+58 1e-02 2e-08
99: -8.8786e+00 -8.8786e+00 2e-07 2e+58 1e-02 2e-08
100: -8.9075e+00 -8.9075e+00 3e-07 2e+59 1e-02 2e-08
Terminated (maximum number of iterations reached).
{'dual infeasibility': 0.010220406013259782,
'dual objective': -8.907506570132597,
'dual slack': 9.311702859901095e-68,
'gap': 2.543051531583804e-07,
'iterations': 100,
'primal infeasibility': 1.9410386424538202e+59,
'primal objective': -8.907506586345107,
'primal slack': 8.257930718191955e-09,
'relative gap': 2.85495329914458e-08,
'residual as dual infeasibility certificate': 2.1791043583727054e+58,
'residual as primal infeasibility certificate': None,
's': <4x1 matrix, tc='d'>,
'status': 'unknown',
'x': <2x1 matrix, tc='d'>,
'y': <0x1 matrix, tc='d'>,
'z': <4x1 matrix, tc='d'>}
}}}
We have a method that completely fail to converge. I had done some
exploration in #12832 when I thought it was a problem in scipy and I will
probably continue to track that one over there after a good editing of the
issue description.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11705#comment:136>
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.
