diff --git a/relapack/src/chetrf.c b/relapack/src/chetrf.c index 8cd3c07742..4351065a19 100644 --- a/relapack/src/chetrf.c +++ b/relapack/src/chetrf.c @@ -19,8 +19,8 @@ void RELAPACK_chetrf( float *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (CREC_SPLIT rounds up) + const blasint cleanlWork = *n * CREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/chetrf_rook.c b/relapack/src/chetrf_rook.c index 9ed1261cff..0738e6245b 100644 --- a/relapack/src/chetrf_rook.c +++ b/relapack/src/chetrf_rook.c @@ -19,8 +19,8 @@ void RELAPACK_chetrf_rook( float *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (CREC_SPLIT rounds up) + const blasint cleanlWork = *n * CREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/csytrf.c b/relapack/src/csytrf.c index 807c91eced..75cb85c37f 100644 --- a/relapack/src/csytrf.c +++ b/relapack/src/csytrf.c @@ -19,8 +19,8 @@ void RELAPACK_csytrf( float *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (CREC_SPLIT rounds up) + const blasint cleanlWork = *n * CREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/csytrf_rook.c b/relapack/src/csytrf_rook.c index 105c6b8b69..ec50170043 100644 --- a/relapack/src/csytrf_rook.c +++ b/relapack/src/csytrf_rook.c @@ -19,8 +19,8 @@ void RELAPACK_csytrf_rook( float *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (CREC_SPLIT rounds up) + const blasint cleanlWork = *n * CREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/dsytrf.c b/relapack/src/dsytrf.c index ba869ad118..63027e980b 100644 --- a/relapack/src/dsytrf.c +++ b/relapack/src/dsytrf.c @@ -19,8 +19,8 @@ void RELAPACK_dsytrf( double *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (DREC_SPLIT rounds up) + const blasint cleanlWork = *n * DREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/dsytrf_rook.c b/relapack/src/dsytrf_rook.c index fcdc2809ff..8c045b69c3 100644 --- a/relapack/src/dsytrf_rook.c +++ b/relapack/src/dsytrf_rook.c @@ -19,8 +19,8 @@ void RELAPACK_dsytrf_rook( double *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (DREC_SPLIT rounds up) + const blasint cleanlWork = *n * DREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/ssytrf.c b/relapack/src/ssytrf.c index 5f8e033913..9bd3997328 100644 --- a/relapack/src/ssytrf.c +++ b/relapack/src/ssytrf.c @@ -18,8 +18,8 @@ void RELAPACK_ssytrf( float *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (SREC_SPLIT rounds up) + const blasint cleanlWork = *n * SREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/ssytrf_rook.c b/relapack/src/ssytrf_rook.c index b40f12271a..512a84b3a7 100644 --- a/relapack/src/ssytrf_rook.c +++ b/relapack/src/ssytrf_rook.c @@ -19,8 +19,8 @@ void RELAPACK_ssytrf_rook( float *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (SREC_SPLIT rounds up) + const blasint cleanlWork = *n * SREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/zhetrf.c b/relapack/src/zhetrf.c index 3d458fecf8..abe581da76 100644 --- a/relapack/src/zhetrf.c +++ b/relapack/src/zhetrf.c @@ -19,8 +19,8 @@ void RELAPACK_zhetrf( double *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (ZREC_SPLIT rounds up) + const blasint cleanlWork = *n * ZREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/zhetrf_rook.c b/relapack/src/zhetrf_rook.c index 605e3a77f7..e90a49de9b 100644 --- a/relapack/src/zhetrf_rook.c +++ b/relapack/src/zhetrf_rook.c @@ -19,8 +19,8 @@ void RELAPACK_zhetrf_rook( double *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (ZREC_SPLIT rounds up) + const blasint cleanlWork = *n * ZREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/zsytrf.c b/relapack/src/zsytrf.c index 59daba02f7..c0b8378642 100644 --- a/relapack/src/zsytrf.c +++ b/relapack/src/zsytrf.c @@ -19,8 +19,8 @@ void RELAPACK_zsytrf( double *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (ZREC_SPLIT rounds up) + const blasint cleanlWork = *n * ZREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1; diff --git a/relapack/src/zsytrf_rook.c b/relapack/src/zsytrf_rook.c index 0fd8e70335..5256771f64 100644 --- a/relapack/src/zsytrf_rook.c +++ b/relapack/src/zsytrf_rook.c @@ -19,8 +19,8 @@ void RELAPACK_zsytrf_rook( double *Work, const blasint *lWork, blasint *info ) { - // Required work size - const blasint cleanlWork = *n * (*n / 2); + // Required work size (ZREC_SPLIT rounds up) + const blasint cleanlWork = *n * ZREC_SPLIT(*n); blasint minlWork = cleanlWork; #if XSYTRF_ALLOW_MALLOC minlWork = 1;