Add two unit tests for matfaust and pyfaust DFT factorization (factparams.ParamsHierarchicalDFT).

misc/test/src/Matlab/FaustFactoryTest.m

@@ -558,6 +558,20 @@ classdef FaustFactoryTest < matlab.unittest.TestCase
 			this.verifyLessThanOrEqual(err, 0.20)
 		end
 
+		function test_hierarchical_dft(this)
+			disp('Test hierarchical dft')
+			import matfaust.fact.hierarchical
+			import matfaust.dft
+			DFT = full(dft(32))
+			F = hierarchical(DFT, 'dft', 'backend', 2020)
+			err = norm(full(F)-DFT)/norm(DFT)
+			this.verifyLessThanOrEqual(err, 1e-3)
+			DFT = full(dft(32))
+			F = hierarchical(DFT, 'dft', 'backend', 2016)
+			err = norm(full(F)-DFT)/norm(DFT)
+			this.verifyLessThanOrEqual(err, 1e-3)
+		end
+
 	end
 
 	methods

misc/test/src/Python/test_FaustPy.py

@@ -1605,6 +1605,17 @@ class TestFaustFactory(unittest.TestCase):
        F = hierarchical_mhtp(M, param, mhtp_param)
        self.assertLessEqual(norm((F.toarray()-M), "fro")/norm(M,"fro"), 0.5)
 
+    def test_hierarchical_dft(self):
+        print("Test hierarchical dft")
+        from pyfaust import dft
+        from pyfaust.fact import hierarchical
+        DFT = dft(32).toarray()
+        F = hierarchical(DFT, 'dft', backend=2020)
+        err = norm(F.toarray()-DFT)/norm(DFT)
+        self.assertLessEqual(err, 1e-6)
+        F = hierarchical(DFT, 'dft', backend=2016)
+        err = norm(F.toarray()-DFT)/norm(DFT)
+        self.assertLessEqual(err, 1e-6)
 
 
 if __name__ == "__main__":