diff --git a/compute/pzgered.c b/compute/pzgered.c index 1051ee91f45370bba6cefbdb21d3b17a6862c99f..c43e9a0e4704b958bbe6380032f8c32ee3b060fe 100644 --- a/compute/pzgered.c +++ b/compute/pzgered.c @@ -162,7 +162,7 @@ void chameleon_pzgered( cham_uplo_t uplo, double prec, CHAM_desc_t *A, RUNTIME_option_t options; CHAM_desc_t Wcol; CHAM_desc_t Welt; - double gnorm, lnorm, threshold, eps; + double gnorm, threshold, eps; int workmt, worknt; int m, n; @@ -233,16 +233,14 @@ void chameleon_pzgered( cham_uplo_t uplo, double prec, CHAM_desc_t *A, int nmax = ( uplo == ChamLower ) ? chameleon_min(m+1, A->nt) : A->nt; for(n = nmin; n < nmax; n++) { - CHAM_tile_t *tile = A->get_blktile( A, m, n ); - int tempnn = ( n == (A->nt-1) ) ? A->n - n * A->nb : A->nb; /* - * u_{high} = 1e-16 (later should be application accuracy) - * u_{low} = 1e-8 - * ||A_{i,j}||_F < u_{high} * || A ||_F / (nt * u_{low}) - * ||A_{i,j}||_F < threshold / u_{low} - */ + * u_{high} = 1e-16 (later should be application accuracy) + * u_{low} = 1e-8 + * ||A_{i,j}||_F < u_{high} * || A ||_F / (nt * u_{low}) + * ||A_{i,j}||_F < threshold / u_{low} + */ INSERT_TASK_zgered( &options, threshold, tempmm, tempnn, A( m, n ), W( &Wcol, m, n ) ); } diff --git a/compute/pzhered.c b/compute/pzhered.c index 97e171b3f11dd6a40ccb44e7d2f85ee41324b62f..cc32f4243c81dd2e63656cf9972850e53dddbd83 100644 --- a/compute/pzhered.c +++ b/compute/pzhered.c @@ -187,7 +187,7 @@ void chameleon_pzhered( cham_trans_t trans, cham_uplo_t uplo, double prec, CHAM_ RUNTIME_option_t options; CHAM_desc_t Wcol; CHAM_desc_t Welt; - double gnorm, lnorm, threshold, eps; + double gnorm, threshold, eps; int workmt, worknt; int m, n; @@ -263,17 +263,14 @@ void chameleon_pzhered( cham_trans_t trans, cham_uplo_t uplo, double prec, CHAM_ for (n = nmin; n < nmax; n++) { - CHAM_tile_t *tile = A->get_blktile(A, m, n); - int tempnn = (n == (A->nt - 1)) ? A->n - n * A->nb : A->nb; /* - * u_{high} = 1e-16 (later should be application accuracy) - * u_{low} = 1e-8 - * ||A_{i,j}||_F < u_{high} * || A ||_F / (nt * u_{low}) - * ||A_{i,j}||_F < threshold / u_{low} - */ - + * u_{high} = 1e-16 (later should be application accuracy) + * u_{low} = 1e-8 + * ||A_{i,j}||_F < u_{high} * || A ||_F / (nt * u_{low}) + * ||A_{i,j}||_F < threshold / u_{low} + */ INSERT_TASK_zgered( &options, threshold, tempmm, tempnn, A( m, n ), W( &Wcol, m, n ) ); }