diff --git a/graph_framework/special_functions.hpp b/graph_framework/special_functions.hpp
index e250b0a2308c9d53fadd70afc16bf7ce2c10f2d1..7eadcfa3bd76e15a3716bf9c56cd619d20c3e722 100644
--- a/graph_framework/special_functions.hpp
+++ b/graph_framework/special_functions.hpp
@@ -18,7 +18,6 @@ template<typename T>
 using complex_type = complex<T>;
 
 #else
-#include <iostream>
 #include <complex>
 #include <cfloat>
 #include <cmath>
@@ -1583,7 +1582,7 @@ namespace special {
 #endif
     complex_type<T> erfi(const complex_type<T> z) {
 //  Avoids NaN instead of doing i<T>*z and -i*temp;
-        const complex_type<T> temp = erf_complex<T> (complex_type<T> (-imag<T> (z), real<T> (z)));
+        const complex_type<T> temp = erf_complex<T> (complex_type<T> (-imag(z), real(z)));
         return complex_type<T> (imag(temp), -real(temp));
     }
 }