Add Faust::Vect::Rand/setRand().

src/faust_linear_operator/CPU/faust_Vect.h

@@ -193,6 +193,8 @@ namespace Faust
         bool equality(Vect<FPP,Cpu> const &x, FPP precision) const;
 		FPP mean();
 
+		static Faust::Vect<FPP, Cpu>* rand(faust_unsigned_int size);
+		void setRand();
 
         // friend algebra
         friend void Faust::gemv<>(const Faust::MatDense<FPP,Cpu> & A,const Faust::Vect<FPP,Cpu> & x,Faust::Vect<FPP,Cpu> & y,const FPP & alpha, const FPP & beta, char typeA);

src/faust_linear_operator/CPU/faust_Vect.hpp

@@ -309,4 +309,27 @@ FPP Faust::Vect<FPP, Cpu>::normL1() const
 		sum += std::abs(vec.data()[j]);
 	return sum;
 }
+
+	template<typename FPP>
+Faust::Vect<FPP, Cpu>* Faust::Vect<FPP, Cpu>::rand(faust_unsigned_int size)
+{
+	Faust::Vect<FPP, Cpu>* v = nullptr;
+	try
+	{
+		v = new Faust::Vect<FPP, Cpu>(size);
+		v->setRand();
+	}
+	catch(std::bad_alloc e)
+	{
+		//		std::cerr << "Out of memory." << std::endl;
+	}
+	return v;
+}
+
+	template<typename FPP>
+void Faust::Vect<FPP, Cpu>::setRand()
+{
+	vec = Eigen::Matrix<FPP, Eigen::Dynamic, 1>::Random(this->size(), 1);
+}
+
 #endif