scipy/scipy — pull request #3206 by bjodah (2014-01-12 15:25:21+00:00)
Title: Updated integrate.ode doc (lband, uband)

The reference to `rband` in the docstring is puzzling to me (right?).
I also added info about the dimension requirements for the Jacobian. That extra lband is used as work space
if my memory serves me correctly (I vaguely remember this is done in LAPACK too but I couldn't find any reference about it).

Example of lband=uband=4 when using a matrix dimension of (lband+uband+1, len(y)):

0-th dimension must be fixed to 13 but got 9
rv_cb_arr is NULL
Call-back cb_jac_in_dvode__user__routines failed.
Traceback (most recent call last):
File "four_species.py", line 63, in
argh.dispatch_command(main)
File "/usr/local/lib/python2.7/dist-packages/argh/dispatching.py", line 224, in dispatch_command
dispatch(parser, *args, **kwargs)
File "/usr/local/lib/python2.7/dist-packages/argh/dispatching.py", line 121, in dispatch
for line in lines:
File "/usr/local/lib/python2.7/dist-packages/argh/dispatching.py", line 197, in _execute_command
for line in result:
File "/usr/local/lib/python2.7/dist-packages/argh/dispatching.py", line 180, in _call
result = args.function(*positional, **keywords)
File "four_species.py", line 55, in main
tout, yout, info = run(sys, y0, t0, tend, nt)
File "/home/bjorn/vc/chemreac/chemreac/integrate.py", line 69, in run
runner.integrate(tout[i])
File "/usr/local/lib/python2.7/dist-packages/scipy/integrate/_ode.py", line 334, in integrate
self.f_params, self.jac_params)
File "/usr/local/lib/python2.7/dist-packages/scipy/integrate/_ode.py", line 636, in run
args[5:]))
SystemError: NULL result without error in PyObject_Call

changing the dimension to (2*lband+uband+1, len(y)) fixes this. It would be great if someone could confirm that this is the case.