| 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 | | > }}} |
