2 | | > Nicolas Thiery just tested the doctests of our book with Sage 8.1beta5 and several do fail (while all tests pass with Sage 8.0). Here are the issues: |
3 | | > |
4 | | > 1) the {{{lil_matrix}}} do not work any more. Cf the example on page 299 of the english translation: |
5 | | > {{{ |
6 | | > sage: from scipy.sparse.linalg.dsolve import * |
7 | | > sage: from scipy.sparse import lil_matrix |
8 | | > sage: from numpy import array |
9 | | > sage: n = 200 |
10 | | > sage: n2 = n*n |
11 | | > sage: A = lil_matrix((n2, n2)) |
12 | | > sage: h2 = 1./float((n+1)^2) |
13 | | > sage: for i in range(0,n2): |
14 | | > ....: A[i,i]=4*h2+1. |
15 | | > ....: if i+1<n2: A[i,int(i+1)]=-h2 |
16 | | > ....: if i>0: A[i,int(i-1)]=-h2 |
17 | | > ....: if i+n<n2: A[i,int(i+n)]=-h2 |
18 | | > ....: if i-n>=0: A[i,int(i-n)]=-h2 |
19 | | > sage: Acsc = A.tocsc() |
20 | | > sage: b = array([1 for i in range(0,n2)]) |
21 | | > sage: solve = factorized(Acsc) # LU factorisation |
22 | | > sage: S = solve(b) # resolution |
23 | | > }}} |
24 | | > Now if fails at the {{{A = lil_matrix((n2, n2))}}} line with: |
25 | | > {{{ |
26 | | > TypeError: unrecognized lil_matrix constructor usage |
27 | | > }}} |