From bbd2320c56865d08d975a5019c9de3c9a3a11975 Mon Sep 17 00:00:00 2001
From: Krad <kradchen123@163.com>
Date: Tue, 25 Apr 2023 09:41:55 +0800
Subject: [PATCH] Add random function.

---
 src/Function.cpp | 11 +++++++++++
 src/Function.h   |  2 +-
 2 files changed, 12 insertions(+), 1 deletion(-)

diff --git a/src/Function.cpp b/src/Function.cpp
index df47a4e..3716780 100644
--- a/src/Function.cpp
+++ b/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;
+    }
+
 
 }
\ No newline at end of file
diff --git a/src/Function.h b/src/Function.h
index a29f57d..225741e 100644
--- a/src/Function.h
+++ b/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);