[SciPy-dev] sparse solvers in scipy rock

Gabriel Gellner ggellner@uoguelph...
Sun Jan 20 16:01:21 CST 2008


Yes! Please, if you have time, make a cookbook entry . . . I would love to see
how cython is used in this way.

Gabriel

On Tue, Jan 15, 2008 at 11:16:10AM -0700, Fernando Perez wrote:
> On Jan 15, 2008 6:20 AM, Ondrej Certik <ondrej@certik.cz> wrote:
> > Hi,
> >
> > yesterday I wrote some solver for 1D time dependent schrodinger
> > equation to my school,
> > using explicit and implicit methods and using scipy sparse solvers,
> > because I was lazy to use
> > some fortran routines for inverting the tridiagonal matrix, and I was
> > surprised how nice it works in scipy.
> > You did a great job. I am convinced that is the way to call all sparse solvers.
> >
> > Here is the code:
> >
> > http://hg.certik.cz/schrod1D/
> >
> > I based it on the example:
> >
> > http://www.scipy.org/Cookbook/SchrodingerFDTD
> >
> > which uses central differences explicit method, but I wrote it my way,
> > using complex numbers directly and Cython to speed it up.
> >
> > I also implemented euler differences explicit method,
> > which totally blows up after a few iterations, and also an implicit
> > method, which works really well.
> >
> > Where do you think I could put it together with some documentation how
> > to play with it? to the above wiki, or should I create a new wiki for
> > that?
> > It's an example how to use Cython+Numpy+scipy sparse solvers (superlu
> > currently).
> 
> I'd suggest making another Cookbook entry for it, and linking both
> Schrodinger ones to each other, so that an interested party can read
> both the simpler, pure python FDTD approach and your more
> sophisticated one as well.
> 
> The better the cookbook gets, the easier it will be for newcomers to
> find examples that are similart to their needs, from which they can
> learn and get started.
> 
> Thanks!
> 
> Cheers,
> 
> f
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-dev


More information about the Scipy-dev mailing list