Implement the optional normalization of Hadamard transform in the C++ backend...

Implement the optional normalization of Hadamard transform in the C++ backend (yet to use in wrappers).

Issue #112 related.

src/algorithm/factorization/faust_WHT.h

@@ -8,7 +8,7 @@ namespace Faust {
 	 * \brief Fast Walsh-Hadamard Transform.
 	 */
 	template<typename FPP>
-	void wht_factors(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool cloning_fact=true);
+	void wht_factors(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool cloning_fact=true, const bool norma=false);
 
 }
 #include "faust_WHT.hpp"

src/algorithm/factorization/faust_WHT.hpp

@@ -3,7 +3,7 @@
 
 namespace Faust {
 	template<typename FPP>
-		void wht_factors(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool cloning_fact)
+		void wht_factors(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool cloning_fact, const bool norma)
 		{
 			if(n == 0)
 			{
@@ -93,6 +93,12 @@ namespace Faust {
 					else
 						factors[i] = factor;
 
+				if(norma)
+				{
+					factors[0] = factor->Clone();
+					*factors[0] *= static_cast<FPP>(1.0/sqrt((float)order));
+				}
+
 			}
 		}
 }

src/faust_linear_operator/CPU/faust_MatSparse.h

@@ -86,7 +86,7 @@ template<typename FPP>
 void Faust::spgemm(const Faust::MatSparse<FPP,Cpu> & A,const Faust::MatDense<FPP,Cpu> & B, Faust::MatDense<FPP,Cpu> & C,const FPP & alpha, const FPP & beta, char  typeA, char  typeB);
 
 template<typename FPP>
-void Faust::wht_factors(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool);
+void Faust::wht_factors(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool, const bool);
 
 //! \namespace Faust
 //! \brief Faust namespace contains the principal class of the project.
@@ -106,7 +106,7 @@ namespace Faust
 		{
 
 			friend Faust::TransformHelper<FPP,Cpu>; // TODO: limit to needed member functions only
-			friend void Faust::wht_factors<>(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool);
+			friend void Faust::wht_factors<>(unsigned int n, vector<MatGeneric<FPP,Cpu>*>&  factors, const bool, const bool);
 			friend class MatDense<FPP,Cpu>;
 			//friend void MatDense<FPP,Cpu>::operator+=(const MatSparse<FPP,Cpu>& S);