Changes between Initial Version and Version 1 of Ticket #23572, comment 15


Ignore:
Timestamp:
09/15/17 10:02:31 (2 years ago)
Author:
jdemeyer
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • Ticket #23572, comment 15

    initial v1  
    11Replying to [comment:12 zimmerma]:
    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 > }}}
    282> Aren't {{{lil_matrix}}} tested in the current doctests?
    293