diff --git a/py/formatfloat.c b/py/formatfloat.c
index 35cd5d51af705d329d0f4062681212b32697b2be..b61a958a25972753d3d6fd3af8a363709821c834 100644
--- a/py/formatfloat.c
+++ b/py/formatfloat.c
@@ -30,6 +30,7 @@
 #include <assert.h>
 #include <stdlib.h>
 #include <stdint.h>
+#include <math.h>
 #include "py/formatfloat.h"
 
 /***********************************************************************
@@ -82,7 +83,6 @@ static inline int fp_isless1(float x) { union floatbits fb = {x}; return fb.u <
 #define FPROUND_TO_ONE 0.999999999995
 #define FPDECEXP 256
 #define FPMIN_BUF_SIZE 7 // +9e+199
-#include <math.h>
 #define fp_signbit(x) signbit(x)
 #define fp_isspecial(x) 1
 #define fp_isnan(x) isnan(x)
@@ -122,7 +122,7 @@ int mp_format_float(FPTYPE f, char *buf, size_t buf_size, char fmt, int prec, ch
         }
         return buf_size >= 2;
     }
-    if (fp_signbit(f)) {
+    if (fp_signbit(f) && !isnan(f)) {
         *s++ = '-';
         f = -f;
     } else {
diff --git a/tests/float/complex1.py b/tests/float/complex1.py
index 85441054532454745e177b56abadf226b5ba4caf..479b4b3485b269fc5a7ce85420b1782039a407f9 100644
--- a/tests/float/complex1.py
+++ b/tests/float/complex1.py
@@ -59,6 +59,7 @@ ans = (-1.2) ** -3.4; print("%.5g %.5g" % (ans.real, ans.imag))
 
 # check printing of inf/nan
 print(float('nan') * 1j)
+print(float('-nan') * 1j)
 print(float('inf') * (1 + 1j))
 print(float('-inf') * (1 + 1j))
 
diff --git a/tests/float/float1.py b/tests/float/float1.py
index 137dacc233bf63aeee27537e257b74ed2d3a3a19..c64f965a7d7cd8ae41d712f47634abee07861473 100644
--- a/tests/float/float1.py
+++ b/tests/float/float1.py
@@ -21,6 +21,7 @@ print(float("INF"))
 print(float("infinity"))
 print(float("INFINITY"))
 print(float("nan"))
+print(float("-nan"))
 print(float("NaN"))
 try:
     float("")