Create a function which makes a random matrix for CPD initial guess

btas/generic/cp.h

@@ -448,6 +448,17 @@ namespace btas {
     virtual void build_random(ind_t rank, ConvClass &converge_test, bool direct, ind_t max_als, bool calculate_epsilon,
                               double &epsilon, bool &fast_pI) = 0;
 
+    /// Create a rank \c rank initial guess using
+    /// random numbers from a uniform distribution
+
+    template<typename Generator, typename Distribution>
+    Tensor make_random_factor(ind_t row, ind_t col, Generator gen, Distribution dist){
+      Tensor a(row, col);
+      for (auto iter = a.begin(); iter != a.end(); ++iter) {
+        *(iter) = dist(gen);
+      }
+      return a;
+    }
     /// Generates V by first Multiply A^T.A then Hadamard product V(i,j) *=
     /// A^T.A(i,j);
     /// \param[in] n The mode being optimized, all other modes held constant