aaron/ispc
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

use correct abs function for doubles

Browse Source
This commit is contained in:
Aaron Gutierrez
2017-05-11 16:28:11 -04:00
parent 5c0911c2a8
commit 34d26554bf
1 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
tests_ispcpp/sqrt.cpp
View File

@@ -2,6 +2,7 @@
#include <algorithm> #include <algorithm>
#include <pthread.h> #include <pthread.h>
#include <math.h> #include <math.h>
#include <cmath>
#include "CycleTimer.h" #include "CycleTimer.h"
#include "sqrt.h" #include "sqrt.h"
@@ -45,11 +46,11 @@ void sqrtSerial(int N,
double x = values[i]; double x = values[i];
double guess = initialGuess; double guess = initialGuess;
double error = abs(guess * guess * x - 1.); double error = std::abs(guess * guess * x - 1.);
while (error > kThreshold) { while (error > kThreshold) {
guess = (3. * guess - x * guess * guess * guess) * 0.5; guess = (3. * guess - x * guess * guess * guess) * 0.5;
error = fabs(guess * guess * x - 1.); error = std::abs(guess * guess * x - 1.);
} }
output[i] = x * guess; output[i] = x * guess;
@@ -66,7 +67,7 @@ static void verifyResult(int N, float* result, float* gold) {
static void verifyResult(int N, double* result, double* gold) { static void verifyResult(int N, double* result, double* gold) {
for (int i=0; i<N; i++) { for (int i=0; i<N; i++) {
if (abs(result[i] - gold[i]) > 1e-4) { if (std::abs(result[i] - gold[i]) > 1e-4) {
printf("Error: [%d] Got %f expected %f\n", i, result[i], gold[i]); printf("Error: [%d] Got %f expected %f\n", i, result[i], gold[i]);
} }
} }