Fixed harmless macro bug in level-1m operations.

Details:
- Fixed some inconsistent usage of n_iter_max and n_iter in the two
  bl2_set_dims_incs_uplo_[12]m macros. The right thing ended up happening
  despite the bug, which is why I had not discovered it until now.
This commit is contained in:
Field G. Van Zee
2012-12-20 11:35:14 -06:00
parent 8945db6ec9
commit 2fecc88ca2
8 changed files with 33 additions and 31 deletions

View File

@@ -133,8 +133,8 @@ void PASTEMAC3(cha,chx,chy,varname)( \
ctype_y* y1; \
uplo_t uplox_eff; \
conj_t conjx; \
dim_t n_iter_max, n_iter; \
dim_t n_elem_max, n_elem; \
dim_t n_iter; \
dim_t n_elem, n_elem_max; \
inc_t ldx, incx; \
inc_t ldy, incy; \
dim_t j, i; \
@@ -152,7 +152,7 @@ void PASTEMAC3(cha,chx,chy,varname)( \
/* Set various loop parameters. */ \
bl2_set_dims_incs_uplo_2m( diagoffx, diagx, transx, \
uplox, m, n, rs_x, cs_x, rs_y, cs_y, \
uplox_eff, n_elem_max, n_iter_max, incx, ldx, incy, ldy, \
uplox_eff, n_elem_max, n_iter, incx, ldx, incy, ldy, \
ij0, n_shift ); \
\
if ( bl2_is_zeros( uplox_eff ) ) return; \

View File

@@ -119,8 +119,8 @@ void PASTEMAC2(chx,chy,varname)( \
ctype_y* y1; \
uplo_t uplox_eff; \
conj_t conjx; \
dim_t n_iter_max, n_iter; \
dim_t n_elem_max, n_elem; \
dim_t n_iter; \
dim_t n_elem, n_elem_max; \
inc_t ldx, incx; \
inc_t ldy, incy; \
dim_t j, i; \
@@ -135,7 +135,7 @@ void PASTEMAC2(chx,chy,varname)( \
/* Set various loop parameters. */ \
bl2_set_dims_incs_uplo_2m( diagoffx, diagx, transx, \
uplox, m, n, rs_x, cs_x, rs_y, cs_y, \
uplox_eff, n_elem_max, n_iter_max, incx, ldx, incy, ldy, \
uplox_eff, n_elem_max, n_iter, incx, ldx, incy, ldy, \
ij0, n_shift ); \
\
if ( bl2_is_zeros( uplox_eff ) ) return; \
@@ -154,9 +154,9 @@ void PASTEMAC2(chx,chy,varname)( \
y1 = y_cast + (j )*ldy + (0 )*incy; \
\
PASTEMAC2(chx,chy,kername)( conjx, \
n_elem, \
x1, incx, \
y1, incy ); \
n_elem, \
x1, incx, \
y1, incy ); \
} \
} \
else \
@@ -172,9 +172,9 @@ void PASTEMAC2(chx,chy,varname)( \
y1 = y_cast + (ij0+j )*ldy + (0 )*incy; \
\
PASTEMAC2(chx,chy,kername)( conjx, \
n_elem, \
x1, incx, \
y1, incy ); \
n_elem, \
x1, incx, \
y1, incy ); \
} \
} \
else if ( bl2_is_lower( uplox_eff ) ) \
@@ -189,9 +189,9 @@ void PASTEMAC2(chx,chy,varname)( \
y1 = y_cast + (j )*ldy + (ij0+i )*incy; \
\
PASTEMAC2(chx,chy,kername)( conjx, \
n_elem, \
x1, incx, \
y1, incy ); \
n_elem, \
x1, incx, \
y1, incy ); \
} \
} \
\

View File

@@ -115,8 +115,8 @@ void PASTEMAC2(chx,chy,varname)( \
ctype_y* y1; \
uplo_t uplox_eff; \
conj_t conjx; \
dim_t n_iter_max, n_iter; \
dim_t n_elem_max, n_elem; \
dim_t n_iter; \
dim_t n_elem, n_elem_max; \
inc_t ldx, incx; \
inc_t ldy, incy; \
dim_t j, i; \
@@ -129,7 +129,7 @@ void PASTEMAC2(chx,chy,varname)( \
/* Set various loop parameters. */ \
bl2_set_dims_incs_uplo_2m( diagoffx, BLIS_NONUNIT_DIAG, transx, \
uplox, m, n, rs_x, cs_x, rs_y, cs_y, \
uplox_eff, n_elem_max, n_iter_max, incx, ldx, incy, ldy, \
uplox_eff, n_elem_max, n_iter, incx, ldx, incy, ldy, \
ij0, n_shift ); \
\
if ( bl2_is_zeros( uplox_eff ) ) return; \

View File

@@ -133,8 +133,8 @@ void PASTEMAC3(cha,chx,chy,varname)( \
ctype_y* y1; \
uplo_t uplox_eff; \
conj_t conjx; \
dim_t n_iter_max, n_iter; \
dim_t n_elem_max, n_elem; \
dim_t n_iter; \
dim_t n_elem, n_elem_max; \
inc_t ldx, incx; \
inc_t ldy, incy; \
dim_t j, i; \
@@ -163,7 +163,7 @@ void PASTEMAC3(cha,chx,chy,varname)( \
/* Set various loop parameters. */ \
bl2_set_dims_incs_uplo_2m( diagoffx, diagx, transx, \
uplox, m, n, rs_x, cs_x, rs_y, cs_y, \
uplox_eff, n_elem_max, n_iter_max, incx, ldx, incy, ldy, \
uplox_eff, n_elem_max, n_iter, incx, ldx, incy, ldy, \
ij0, n_shift ); \
\
if ( bl2_is_zeros( uplox_eff ) ) return; \

View File

@@ -121,8 +121,8 @@ void PASTEMAC2(chb,chx,varname)( \
ctype_x* x_cast = x; \
ctype_x* x1; \
uplo_t uplox_eff; \
dim_t n_iter_max, n_iter; \
dim_t n_elem_max, n_elem; \
dim_t n_iter; \
dim_t n_elem, n_elem_max; \
inc_t ldx, incx; \
dim_t j, i; \
dim_t ij0, n_shift; \

View File

@@ -117,8 +117,8 @@ void PASTEMAC2(chb,chx,varname)( \
ctype_x* x_cast = x; \
ctype_x* x1; \
uplo_t uplox_eff; \
dim_t n_iter_max, n_iter; \
dim_t n_elem_max, n_elem; \
dim_t n_iter; \
dim_t n_elem, n_elem_max; \
inc_t ldx, incx; \
dim_t j, i; \
dim_t ij0, n_shift; \

View File

@@ -513,8 +513,8 @@
// argument.
#define bl2_set_dims_incs_uplo_1m( diagoffa, diaga, \
uploa, m, n, rs_a, cs_a, \
uplo_eff, n_elem, n_iter, inca, lda, \
uploa, m, n, rs_a, cs_a, \
uplo_eff, n_elem_max, n_iter, inca, lda, \
ij0, n_shift ) \
{ \
/* If matrix A is entirely "unstored", that is, if either:
@@ -529,6 +529,7 @@
{ \
doff_t diagoffa_use = diagoffa; \
doff_t diagoff_eff; \
dim_t n_iter_max; \
\
if ( bl2_is_unit_diag( diaga ) ) \
bl2_shift_diag_offset_to_shrink_uplo( uploa, diagoffa_use ); \
@@ -599,8 +600,8 @@
#define bl2_set_dims_incs_uplo_2m( \
diagoffa, diaga, transa, \
uploa, m, n, rs_a, cs_a, rs_b, cs_b, \
uplo_eff, n_elem_max, n_iter_max, inca, lda, incb, ldb, \
uploa, m, n, rs_a, cs_a, rs_b, cs_b, \
uplo_eff, n_elem_max, n_iter, inca, lda, incb, ldb, \
ij0, n_shift \
) \
{ \
@@ -616,6 +617,7 @@
{ \
doff_t diagoffa_use = diagoffa; \
doff_t diagoff_eff; \
dim_t n_iter_max; \
\
if ( bl2_is_unit_diag( diaga ) ) \
bl2_shift_diag_offset_to_shrink_uplo( uploa, diagoffa_use ); \

View File

@@ -97,8 +97,8 @@ void PASTEMAC(ch,varname)( \
ctype omega; \
double max_m_n; \
uplo_t uplox_eff; \
dim_t n_iter_max, n_iter; \
dim_t n_elem_max, n_elem; \
dim_t n_iter; \
dim_t n_elem, n_elem_max; \
inc_t ldx, incx; \
dim_t j, i; \
dim_t ij0, n_shift; \