|
|
Ublas : |
Subject: Re: [ublas] Fwd: boost/numeric/bindings/lapack
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2008-11-12 18:16:53
> I've received the following patch suggestion from Julius Ziegler. What he
> points out does look fishy. Since I've never worked with banded matrices in
> LAPACK, I hope someone on the list is better equipped to give a definitive
> answer.
The fix from Julius Ziegler is correct (according to the documentation of SGBTRF and ublas/banded_matrix), but it doesn't go far enough. It should also change "assert(ku > 0)" into "assert(ku >= 0)", because that is all that is required by SGBTRF.