From a4e321400b534d7f547cad8127ff6da45e1ac872 Mon Sep 17 00:00:00 2001 From: "Steven G. Johnson" Date: Fri, 3 Aug 2018 13:00:10 -0400 Subject: [PATCH 1/2] fabs -> fabsl Fixes two calls that were using `fabs` on a `long double` argument rather than `fabsl`, which looks like it is doing an unintentional truncation to `double` precision. --- interface/rotg.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/interface/rotg.c b/interface/rotg.c index 092554299..69443a5a0 100644 --- a/interface/rotg.c +++ b/interface/rotg.c @@ -22,8 +22,8 @@ void CNAME(FLOAT *DA, FLOAT *DB, FLOAT *C, FLOAT *S){ long double s; long double r, roe, z; - long double ada = fabs(da); - long double adb = fabs(db); + long double ada = fabsl(da); + long double adb = fabsl(db); long double scale = ada + adb; #ifndef CBLAS From 165f00c159cf0c4e7e6eef8f656fa68e7cda4ea2 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sat, 4 Aug 2018 20:14:51 +0200 Subject: [PATCH 2/2] fabs -> fabsl --- interface/zrotg.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/interface/zrotg.c b/interface/zrotg.c index 187343d41..8caa411fc 100644 --- a/interface/zrotg.c +++ b/interface/zrotg.c @@ -14,7 +14,7 @@ void NAME(FLOAT *DA, FLOAT *DB, FLOAT *C, FLOAT *S){ long double db_i = *(DB + 1); long double r; - long double ada = fabs(da_r) + fabs(da_i); + long double ada = fabsl(da_r) + fabsl(da_i); PRINT_DEBUG_NAME;