diff --git a/src/Function.cpp b/src/Function.cpp
index 5578324..9b2fe1d 100644
--- a/src/Function.cpp
+++ b/src/Function.cpp
@@ -149,6 +149,10 @@ namespace Aurora {
         return (double *) mkl_malloc(size * sizeof(std::complex<double>), 64);
     }
 
+    void free(void* ptr){
+        mkl_free(ptr);
+    }
+
     double * mul(double scalar, double *input, int size) {
         double* output = malloc(size);
         vdMulI(size,input,1,&scalar,0,output,1);
diff --git a/src/Function.h b/src/Function.h
index d92b455..ced7354 100644
--- a/src/Function.h
+++ b/src/Function.h
@@ -9,6 +9,7 @@
 
 namespace Aurora{
     double* malloc(int size,bool complex = false);
+    void free(void* ptr);
     double* mul( double scalar, double * input, int size);
     double* mul( double* inputA, double * inputB, int size);
     double* mulz( std::complex<double> *inputA, std::complex<double> * inputB, int size);