Add random function.

src/Function.cpp

@@ -10,6 +10,10 @@
 #include <string.h>
 #include <algorithm>
 
+#include <Eigen/Core>
+#include <Eigen/Eigen>
+#include <Eigen/Dense>
+
 namespace Aurora {
     double immse(double *dataA, double *dataB, int size) {
         auto temp = new double[size];
@@ -153,5 +157,12 @@ namespace Aurora {
         return output;
     }
 
+    double *random( int size) {
+        double * data = malloc(size);
+        Eigen::Map<Eigen::VectorXd> srcV(data,size);
+        srcV.setRandom();
+        return data;
+    }
+
 
 }

src/Function.h

@@ -15,7 +15,7 @@ namespace Aurora{
     double* mulz( std::complex<double> *inputA, std::complex<double> * inputB, int size);
     double immse(double * dataA, double * dataB, int size);
     double* std(int rows, int cols, double * input);
-
+    double* random(int size);
     double* inv(int cols,double *pMatrix);
     double* real(int size, std::complex<double> * input);
     std::complex<double> * complex(int size, double * input);