[SciPy-Dev] linalg.solveh_banded => Mac OS BLAS parameter error in DPBSV

[denis]

Folks,
   has anyone run solveh_banded on a mac ppc ?
My versions: scipy 0.7.1 numpy 1.3.0 py 2.6.4 macosx 10.4.11

If this is the wrong place or wrong form, please correct me.


# scipy.linalg.solveh_banded =>
# Parameter 8 to routine DPBSV  was incorrect
# Mac OS BLAS parameter error in DPBSV , parameter #0, (unavailable), is 0
# 144208 26 Jun  2009 _flinalg.so
# (solve_banded is in scipy/linalg/tests, solveh_banded not ...)

import sys
import numpy as np
from scipy.linalg import solve_banded, solveh_banded

N = 10
exec "\n".join( sys.argv[1:] )  # N= ...
np.random.seed(1)
np.set_printoptions( 2, threshold=100, suppress=True )  # .2f

ones = np.ones(N)
ab = np.array(( ones, 2*ones )) / 4
# ab = np.array(( ones, 2*ones, ones )) / 4
b = np.random.normal( 0, 1, N )
print "ab:", ab
print "b:", b

c, x = solveh_banded( ab, b )  # , lower=1 )
     # x = solve_banded( (1,1), ab, b )  # ok
print "x:", x

Thanks, cheers
   -- denis