diff --git a/src/tests/util/numerics/CMakeLists.txt b/src/tests/util/numerics/CMakeLists.txt
index 8d34b39f7..1ccdc92a8 100644
--- a/src/tests/util/numerics/CMakeLists.txt
+++ b/src/tests/util/numerics/CMakeLists.txt
@@ -4,3 +4,6 @@ add_test(mpq ${CMAKE_CURRENT_BINARY_DIR}/mpq)
 add_executable(mpbq mpbq.cpp)
 target_link_libraries(mpbq ${EXTRA_LIBS})
 add_test(mpbq ${CMAKE_CURRENT_BINARY_DIR}/mpbq)
+add_executable(mpfp mpfp.cpp)
+target_link_libraries(mpfp ${EXTRA_LIBS})
+add_test(mpfp ${CMAKE_CURRENT_BINARY_DIR}/mpfp)
diff --git a/src/tests/util/numerics/mpfp.cpp b/src/tests/util/numerics/mpfp.cpp
new file mode 100644
index 000000000..451026b18
--- /dev/null
+++ b/src/tests/util/numerics/mpfp.cpp
@@ -0,0 +1,48 @@
+/*
+Copyright (c) 2013 Microsoft Corporation. All rights reserved.
+Released under Apache 2.0 license as described in the file LICENSE.
+
+Author: Soonho Kong
+*/
+#include <iostream>
+#include "test.h"
+#include "mpfp.h"
+using namespace lean;
+
+static void tst1() {
+    mpfp a(3.141592, 64);
+    std::cout << "a            = |" << a << "|" << std::endl;
+    std::cout << "exp(a, UP  ) = |" << exp(a, MPFR_RNDU) << "|" << std::endl;
+    std::cout << "exp(a, NEAR) = |" << exp(a, MPFR_RNDN) << "|" << std::endl;
+    std::cout << "exp(a, DOWN) = |" << exp(a, MPFR_RNDD) << "|" << std::endl;
+
+    mpfp b = mpfp(5.141592, 128);
+    std::cout << "b            = |" << b << "|" << std::endl;
+    std::cout << "exp(b, UP  ) = |" << exp(b, MPFR_RNDU) << "|" << std::endl;
+    std::cout << "exp(b, NEBR) = |" << exp(b, MPFR_RNDN) << "|" << std::endl;
+    std::cout << "exp(b, DOWN) = |" << exp(b, MPFR_RNDD) << "|" << std::endl;
+
+    mpfp c = mpfp(6.141592, 256);
+    std::cout << "c            = |" << c << "|" << std::endl;
+    std::cout << "exp(c, UP  ) = |" << exp(c, MPFR_RNDU) << "|" << std::endl;
+    std::cout << "exp(c, NEAR) = |" << exp(c, MPFR_RNDN) << "|" << std::endl;
+    std::cout << "exp(c, DOWN) = |" << exp(c, MPFR_RNDD) << "|" << std::endl;
+
+    std::cout << "a            = |" << a << "|" << std::endl;
+    std::cout << "exp(a, UP  ) = |" << exp(a, MPFR_RNDU) << "|" << std::endl;
+    std::cout << "exp(a, NEAR) = |" << exp(a, MPFR_RNDN) << "|" << std::endl;
+    std::cout << "exp(a, DOWN) = |" << exp(a, MPFR_RNDD) << "|" << std::endl;
+
+    std::cout << "b            = |" << b << "|" << std::endl;
+    std::cout << "exp(b, UP  ) = |" << exp(b, MPFR_RNDU) << "|" << std::endl;
+    std::cout << "exp(b, NEAR) = |" << exp(b, MPFR_RNDN) << "|" << std::endl;
+    std::cout << "exp(b, DOWN) = |" << exp(b, MPFR_RNDD) << "|" << std::endl;
+}
+
+
+int main() {
+    continue_on_violation(true);
+    tst1();
+    return has_violations() ? 1 : 0;
+
+}