Browse Source

Revert "Avoid taking root of negative number in symv_thread.c"

tags/v0.3.8^2
Martin Kroeker GitHub 6 years ago
parent
commit
8617d75548
No known key found for this signature in database GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 3 additions and 7 deletions
  1. +3
    -7
      driver/level2/symv_thread.c

+ 3
- 7
driver/level2/symv_thread.c View File

@@ -166,11 +166,7 @@ int CNAME(BLASLONG m, FLOAT *alpha, FLOAT *a, BLASLONG lda, FLOAT *x, BLASLONG i
if (nthreads - num_cpu > 1) {

double di = (double)i;
if (di * di - dnum > 0) {
width = ((BLASLONG)(-sqrt(di * di - dnum) + di) + mask) & ~mask;
} else {
width = m - i;
}
width = ((BLASLONG)(sqrt(di * di + dnum) - di) + mask) & ~mask;

if (width < 4) width = 4;
if (width > m - i) width = m - i;
@@ -216,9 +212,9 @@ int CNAME(BLASLONG m, FLOAT *alpha, FLOAT *a, BLASLONG lda, FLOAT *x, BLASLONG i

double di = (double)(m - i);
if (di * di - dnum > 0) {
width = ((BLASLONG)(-sqrt(di * di - dnum) + di) + mask) & ~mask;
width = ((BLASLONG)(-sqrt(di * di - dnum) + di) + mask) & ~mask;
} else {
width = m - i;
width = m - i;
}

if (width < 4) width = 4;


Loading…
Cancel
Save