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