fix: preserve sign in math/base/special/{ceil2,floor2}

lib/node_modules/@stdlib/math/base/special/ceil2/lib/main.js

@@ -72,7 +72,7 @@ function ceil2( x ) {
 
 	// If provided the smallest subnormal, no rounding possible:
 	if ( p === MIN_EXP_SUBNORMAL ) {
-		return x;
+		return sign * x;
 	}
 	// Determine a power of two which rounds the input value toward positive infinity:
 	if ( sign === -1.0 ) {

lib/node_modules/@stdlib/math/base/special/ceil2/src/main.c

@@ -58,7 +58,7 @@ double stdlib_base_ceil2( const double x ) {
 
 	// If provided the smallest subnormal, no rounding possible:
 	if ( p == STDLIB_CONSTANT_FLOAT64_MIN_BASE2_EXPONENT_SUBNORMAL ) {
-		return xc;
+		return x;
 	}
 
 	// Determine a power of two which rounds the input value toward positive infinity:

lib/node_modules/@stdlib/math/base/special/ceil2/test/test.js

@@ -82,6 +82,12 @@ tape( 'the function returns the minimum double-precision floating-point value if
 	t.end();
 });
 
+tape( 'the function returns the negative minimum double-precision floating-point value if provided the negative minimum double-precision floating-point value', function test( t ) {
+	var v = ceil2( -SMALLEST_SUBNORMAL );
+	t.strictEqual( v, -SMALLEST_SUBNORMAL, 'returns negative smallest subnormal' );
+	t.end();
+});
+
 tape( 'the function overflows if provided a sufficiently large positive value', function test( t ) {
 	var x;
 	var v;

lib/node_modules/@stdlib/math/base/special/ceil2/test/test.native.js

@@ -91,6 +91,12 @@ tape( 'the function returns the minimum double-precision floating-point value if
 	t.end();
 });
 
+tape( 'the function returns the negative minimum double-precision floating-point value if provided the negative minimum double-precision floating-point value', opts, function test( t ) {
+	var v = ceil2( -SMALLEST_SUBNORMAL );
+	t.strictEqual( v, -SMALLEST_SUBNORMAL, 'returns negative smallest subnormal' );
+	t.end();
+});
+
 tape( 'the function overflows if provided a sufficiently large positive value', opts, function test( t ) {
 	var x;
 	var v;

lib/node_modules/@stdlib/math/base/special/floor2/lib/main.js

@@ -72,7 +72,7 @@ function floor2( x ) {
 
 	// If provided the smallest subnormal, no rounding possible:
 	if ( p === MIN_EXP_SUBNORMAL ) {
-		return x;
+		return sign * x;
 	}
 	// Determine a power of two which rounds the input value toward negative infinity:
 	if ( sign === 1.0 ) {

lib/node_modules/@stdlib/math/base/special/floor2/src/main.c

@@ -57,7 +57,7 @@ double stdlib_base_floor2( const double x ) {
 
 	// If provided the smallest subnormal, no rounding possible:
 	if ( p == STDLIB_CONSTANT_FLOAT64_MIN_BASE2_EXPONENT_SUBNORMAL ) {
-		return xc;
+		return x;
 	}
 	// Determine a power of two which rounds the input value toward negative infinity:
 	if ( sign == 1.0 ) {

lib/node_modules/@stdlib/math/base/special/floor2/test/test.js

@@ -82,6 +82,12 @@ tape( 'the function returns the minimum double-precision floating-point value if
 	t.end();
 });
 
+tape( 'the function returns the negative minimum double-precision floating-point value if provided the negative minimum double-precision floating-point value', function test( t ) {
+	var v = floor2( -SMALLEST_SUBNORMAL );
+	t.strictEqual( v, -SMALLEST_SUBNORMAL, 'returns negative smallest subnormal' );
+	t.end();
+});
+
 tape( 'the function overflows if provided a sufficiently large negative value', function test( t ) {
 	var x;
 	var v;

lib/node_modules/@stdlib/math/base/special/floor2/test/test.native.js

@@ -91,6 +91,12 @@ tape( 'the function returns the minimum double-precision floating-point value if
 	t.end();
 });
 
+tape( 'the function returns the negative minimum double-precision floating-point value if provided the negative minimum double-precision floating-point value', opts, function test( t ) {
+	var v = floor2( -SMALLEST_SUBNORMAL );
+	t.strictEqual( v, -SMALLEST_SUBNORMAL, 'returns negative smallest subnormal' );
+	t.end();
+});
+
 tape( 'the function overflows if provided a sufficiently large negative value', opts, function test( t ) {
 	var x;
 	var v;