29 for(
int i=0;
i<
n;
i++){
34 CPPL_INT size =(
kl+
ku+1)*n;
60 const CPPL_INT& _kl,
const CPPL_INT& _ku)
63 if( _m<0 || _n<0 || _kl<0 || _ku<0 || _m<_kl || _n<_ku ){
65 std::cerr <<
"It is impossible to make a matrix you ordered. " << std::endl
66 <<
"Your input was (" << _m <<
"," << _n <<
","<< _ku <<
"," << _kl <<
")." << std::endl;
78 for(
int i=0;
i<
n;
i++){
101 {CPPL_VERBOSE_REPORT;
CPPL_INT ku
upper band width
friend _zgematrix i(const zgbmatrix &)
CPPL_INT kl
lower band width
comple ** darray
array of pointers of column head addresses
CPPL_INT n
matrix column size
CPPL_INT ku
upper band width
CPPL_INT m
matrix row size
Complex Double-precision General Band Matrix Class.
(DO NOT USE) Smart-temporary Complex Double-precision General Band Matrix Class
CPPL_INT kl
lower band width
CPPL_INT m
matrix row size
comple * array
1D array to store matrix data
comple ** darray
array of pointers of column head addresses
CPPL_INT n
matrix column size
comple * array
1D array to store matrix data