diff --git a/compute/pzlaset.c b/compute/pzlaset.c index 519457c3adbf513c458d2ac34c208a2e4374dcf8..3536b3561e40b8e64c2130de4088af263476aca1 100644 --- a/compute/pzlaset.c +++ b/compute/pzlaset.c @@ -70,26 +70,18 @@ void chameleon_pzlaset(cham_uplo_t uplo, } } else if (uplo == ChamUpper) { - for (j = 1; j < A->nt; j++){ + for (j = 0; j < A->nt; j++){ tempjn = j == A->nt-1 ? A->n-j*A->nb : A->nb; - for (i = 0; i < chameleon_min(j, A->mt); i++){ + for (i = 0; i < chameleon_min(j+1, A->mt); i++){ tempim = i == A->mt-1 ? A->m-i*A->mb : A->mb; ldai = BLKLDD(A, i); INSERT_TASK_zlaset( &options, - ChamUpperLower, tempim, tempjn, alpha, alpha, + ChamUpperLower, tempim, tempjn, + alpha, (i == j) ? beta : alpha, A(i, j), ldai); } } - for (j = 0; j < minmn; j++){ - tempjm = j == A->mt-1 ? A->m-j*A->mb : A->mb; - tempjn = j == A->nt-1 ? A->n-j*A->nb : A->nb; - ldaj = BLKLDD(A, j); - INSERT_TASK_zlaset( - &options, - ChamUpper, tempjm, tempjn, alpha, beta, - A(j, j), ldaj); - } } else { for (i = 0; i < A->mt; i++){ @@ -99,19 +91,11 @@ void chameleon_pzlaset(cham_uplo_t uplo, tempjn = j == A->nt-1 ? A->n-j*A->nb : A->nb; INSERT_TASK_zlaset( &options, - ChamUpperLower, tempim, tempjn, alpha, alpha, + ChamUpperLower, tempim, tempjn, + alpha, (i == j) ? beta : alpha, A(i, j), ldai); } } - for (j = 0; j < minmn; j++){ - tempjm = j == A->mt-1 ? A->m-j*A->mb : A->mb; - tempjn = j == A->nt-1 ? A->n-j*A->nb : A->nb; - ldaj = BLKLDD(A, j); - INSERT_TASK_zlaset( - &options, - ChamUpperLower, tempjm, tempjn, alpha, beta, - A(j, j), ldaj); - } } RUNTIME_options_finalize(&options, chamctxt); }