check of input in dsbtrd and zhbtrd

Post here if you want to report a bug to the LAPACK team

check of input in dsbtrd and zhbtrd

Postby Yana Gurieva » Tue Sep 18, 2012 11:30 pm

I've noticed that the input parameters N and KD of dsbtrd abd zhbtrd are checked to be positive. Should'n it also be checked that KD<=N as KD is the number of matrix super/sub-diagonals?
Yana Gurieva
Posts: 7
Joined: Mon Apr 11, 2011 1:26 am

Re: check of input in dsbtrd and zhbtrd

Postby rodney » Wed Sep 19, 2012 10:53 pm


You are right, it would be a reasonable thing to check that KD were not too large for the specified N. The convention that LAPACK and BLAS uses for all of the banded matrix routines is to instead check that the leading dimension of the band matrix, LDAB in this case, is large enough for the specified number of bands. Internally, the routines then typically cut the band depth to be compatible with N if necessary. These routines, for example, internally use KDN = MIN( N-1, KD ) as the band size to ensure that N and KD are compatible.

Posts: 48
Joined: Thu Feb 10, 2011 8:20 pm
Location: Colorado College

Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest