diff --git a/Makefile.am b/Makefile.am
index 7907e9d9..bbec2d2c 100644
--- a/Makefile.am
+++ b/Makefile.am
@@ -128,8 +128,6 @@ UTIL_FILES   =			\
     maths.cpp			\
     maths.hpp			\
     maths.ipp			\
-    maths/fast.hpp		\
-    maths/fast.ipp		\
     matrix.cpp			\
     matrix.hpp			\
     matrix.ipp			\
@@ -369,7 +367,6 @@ TEST_BIN = 			\
     test/json_types		\
     test/ray 			\
     test/maths			\
-    test/maths_fast		\
     test/matrix			\
     test/md2			\
     test/md4			\
diff --git a/maths/fast.hpp b/maths/fast.hpp
deleted file mode 100644
index 9dfe7f50..00000000
--- a/maths/fast.hpp
+++ /dev/null
@@ -1,44 +0,0 @@
-/*
- * Licensed under the Apache License, Version 2.0 (the "License");
- * you may not use this file except in compliance with the License.
- * You may obtain a copy of the License at
- *
- *     http://www.apache.org/licenses/LICENSE-2.0
- *
- * Unless required by applicable law or agreed to in writing, software
- * distributed under the License is distributed on an "AS IS" BASIS,
- * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
- * See the License for the specific language governing permissions and
- * limitations under the License.
- *
- * Copyright 2015 Danny Robson <danny@nerdcruft.net>
- */
-
-#ifndef __UTIL_MATHS_FAST_HPP
-#define __UTIL_MATHS_FAST_HPP
-
-namespace util { namespace maths { namespace fast {
-    ///////////////////////////////////////////////////////////////////////////
-    constexpr float pow (float, float);
-
-    constexpr float pow2 (float);
-
-
-    ///////////////////////////////////////////////////////////////////////////
-    constexpr float exp (float);
-
-
-    ///////////////////////////////////////////////////////////////////////////
-    constexpr float log (float);
-
-    constexpr float log2 (float);
-
-
-    ///////////////////////////////////////////////////////////////////////////
-    constexpr float sqrt (float);
-    constexpr float invsqrt (float);
-} } }
-
-#include "fast.ipp"
-
-#endif
diff --git a/maths/fast.ipp b/maths/fast.ipp
deleted file mode 100644
index a63b10c5..00000000
--- a/maths/fast.ipp
+++ /dev/null
@@ -1,107 +0,0 @@
-/*
- * Licensed under the Apache License, Version 2.0 (the "License");
- * you may not use this file except in compliance with the License.
- * You may obtain a copy of the License at
- *
- *     http://www.apache.org/licenses/LICENSE-2.0
- *
- * Unless required by applicable law or agreed to in writing, software
- * distributed under the License is distributed on an "AS IS" BASIS,
- * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
- * See the License for the specific language governing permissions and
- * limitations under the License.
- *
- * Copyright 2015 Danny Robson <danny@nerdcruft.net>
- */
-
-#ifdef __UTIL_MATHS_FAST_IPP
-#error
-#endif
-#define __UTIL_MATHS_FAST_IPP
-    
-
-namespace util { namespace maths { namespace fast {
-    ///////////////////////////////////////////////////////////////////////////
-    constexpr float
-    pow2 (float p)
-    {
-        float offset = (p < 0) ? 1.0f : 0.0f;
-        float clipp = (p < -126) ? -126.0f : p;
-        int32_t w = static_cast<int32_t> (clipp);
-        float z = clipp - w + offset;
-        union { uint32_t i; float f; } v = { static_cast<uint32_t> ( (1 << 23) * (clipp + 121.2740575f + 27.7280233f / (4.84252568f - z) - 1.49012907f * z) ) };
-
-        return v.f;
-
-        //float clipp = (p < -126) ? -126.0f : p;
-        //union { uint32_t i; float f; } v = { static_cast<uint32_t> ( (1 << 23) * (clipp + 126.94269504f) ) };
-        //return v.f;
-    }
-
-
-    constexpr float
-    log2 (float x)
-    {
-        union { float f; uint32_t i; } vx = { x };
-        union { uint32_t i; float f; } mx = { (vx.i & 0x007FFFFF) | 0x3f000000 };
-        float y = vx.i;
-        y *= 1.1920928955078125e-7f;
-
-        return y - 124.22551499f
-            - 1.498030302f * mx.f 
-            - 1.72587999f / (0.3520887068f + mx.f);
-
-
-        //union { float f; uint32_t i; } vx = { x };
-        //float y = vx.i;
-        //y *= 1.1920928955078125e-7f;
-        //return y - 126.94269504f;
-    }
-
-
-    constexpr float
-    log (float x)
-    {
-        return 0.69314718f * log2 (x);
-    }
-
-
-    constexpr float
-    pow (float a, float b)
-    {
-        return pow2 (b * log2 (a));
-    }
-
-
-    constexpr float
-    sqrt (float x)
-    {
-        //return pow (x, 0.5f);
-        union { float f; int32_t i; } u  = { x };
-
-        int32_t v = u.i;
-        v  -= 1 << 23;
-        v >>= 1;
-        v  += 1 << 29;
-
-        u.i = v;
-        return u.f;
-    }
-
-
-    constexpr float
-    invsqrt (float x)
-    {
-        union
-        {
-            float x;
-            int32_t i;
-        } u = { x };
-
-        u.i = 0x5f3759df - (u.i >> 1);
-
-        // refine estimate. repeat as necessary.
-        u.x = u.x * (1.5f - x * 0.5f * u.x * u.x);
-        return u.x;
-    }
-} } }
diff --git a/test/maths_fast.cpp b/test/maths_fast.cpp
deleted file mode 100644
index 24f5b798..00000000
--- a/test/maths_fast.cpp
+++ /dev/null
@@ -1,39 +0,0 @@
-#include "tap.hpp"
-
-#include "maths/fast.hpp"
-
-
-constexpr float
-threshold (float a, float b)
-{
-    constexpr float PARTS = 100;
-    return std::abs (a + b) / 2 / PARTS;
-}
-
-
-int
-main (void)
-{
-    util::TAP::logger tap;
-
-
-    {
-        auto a = util::maths::fast::log2 (3.456f);
-        auto b = std::log2 (3.456f);
-        tap.expect_lt (std::abs (a - b), threshold (a, b), "fast log2");
-    }
-
-    {
-        auto a = util::maths::fast::pow2 (-100.f);
-        auto b = std::pow (2.f, -100.f);
-        tap.expect_lt (std::abs (a - b), threshold (a, b), "fast pow2");
-    }
-
-    {
-        auto a = util::maths::fast::pow (0.8f, 100.f);
-        auto b = std::pow (0.8f, 100.f);
-        tap.expect_lt (std::abs (a - b), threshold (a, b), "fast pow");
-    }
-
-    return tap.status ();
-}