Replace nearbyint with roundl in CL

Originally, CL used FORTRANs NINT function which always rounds up for .5, while nearbyint rounds towards the even number in this case.
iraf-community · Jan 17, 2024 · 844c79f · 844c79f
1 parent bebbccb
commit 844c79f
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/pkg/cl/unop.c b/pkg/cl/unop.c
@@ -260,7 +260,7 @@ unop (int opcode)
 	    break;
 	case OP_NINT:
 	    if (in_type == OT_REAL)
-		iresult = nearbyint (rval);
+		iresult = roundl (rval);
 	    else
 		iresult = ival;
 	    break;

diff --git a/pkg/ecl/unop.c b/pkg/ecl/unop.c
@@ -294,7 +294,7 @@ unop (int opcode)
 	    break;
 	case OP_NINT:
 	    if (in_type == OT_REAL)
-		iresult = nearbyint (rval);
+		iresult = roundl (rval);
 	    else
 		iresult = ival;
 	    break;