add constructor to QuadraticTimeMMD

examples/meta/src/statistical_testing/quadratic_time_mmd.sg

@@ -7,11 +7,9 @@ RealFeatures features_q(f_features_q)
 #![create_features]
 
 #![create_instance]
-QuadraticTimeMMD mmd()
+QuadraticTimeMMD mmd(features_p, features_q)
 GaussianKernel kernel()
 mmd.set_kernel(kernel)
-mmd.set_p(features_p)
-mmd.set_q(features_q)
 Real alpha = 0.05
 #![create_instance]
 
@@ -36,6 +34,7 @@ Real h0_rejected_permutation = mmd.perform_test(alpha)
 #![perform_test_spectrum]
 mmd.set_null_approximation_method(enum ENullApproximationMethod.MMD2_SPECTRUM)
 mmd.set_num_null_samples(200)
+mmd.spectrum_set_num_eigenvalues(5)
 Real threshold_spectrum = mmd.compute_threshold(alpha)
 Real p_value_spectrum = mmd.compute_p_value(statistic_biased)
 Real h0_rejected_spectrum = mmd.perform_test(alpha)

src/shogun/statistical_testing/QuadraticTimeMMD.cpp

@@ -46,6 +46,14 @@ CQuadraticTimeMMD::CQuadraticTimeMMD() : CMMD()
 	self = std::unique_ptr<Self>(new Self());
 }
 
+CQuadraticTimeMMD::CQuadraticTimeMMD(CFeatures* samples_from_p,
+		CFeatures* samples_from_q) : CMMD()
+{
+	self = std::unique_ptr<Self>(new Self());
+	set_p(samples_from_p);
+	set_q(samples_from_p);
+}
+
 CQuadraticTimeMMD::~CQuadraticTimeMMD()
 {
 }

src/shogun/statistical_testing/QuadraticTimeMMD.h

@@ -32,6 +32,8 @@ class CQuadraticTimeMMD : public CMMD
 	using operation = std::function<float64_t(SGMatrix<float64_t>)>;
 public:
 	CQuadraticTimeMMD();
+	CQuadraticTimeMMD(CFeatures* samples_from_p, CFeatures* samples_from_q);
+
 	virtual ~CQuadraticTimeMMD();
 
 	virtual SGVector<float64_t> sample_null() override;