Impl. MatButterfly::norm() with its unit test.

misc/test/src/C++/faust_butterfly_mat.cpp.in

@@ -40,7 +40,6 @@ void testClone(MatButterfly<FPP, Cpu>& butterflyMat)
 	delete clone;
 }
 
-
 void testTranspose(MatButterfly<FPP, Cpu>& butterflyMat, MatSparse<FPP, Cpu>& spButterflyMat)
 {
 	auto size = butterflyMat.getNbRow();
@@ -155,6 +154,14 @@ void testScalMul(const MatButterfly<FPP, Cpu>& butterflyMat,  const MatSparse<FP
 	std::cout << "MatButterfly scal. mul OK" << std::endl;
 }
 
+void testNorm2(const MatButterfly<FPP, Cpu>& butterflyMat,  const MatSparse<FPP, Cpu>& spButterflyMat)
+{
+//	cout << butterflyMat.norm()  << endl;
+//	cout << spButterflyMat.norm()  << endl;
+	assert(std::abs(butterflyMat.norm() - spButterflyMat.norm()) < 1e-6);
+	std::cout << "MatButterfly fro-norm OK" << std::endl;
+}
+
 int main(int argc, char** argv)
 {
 	int log2size = 4;
@@ -220,5 +227,6 @@ int main(int argc, char** argv)
 	testType(butterflyMat);
 
 	testScalMul(butterflyMat, spButterflyMat);
+	testNorm2(butterflyMat, spButterflyMat);
 	return EXIT_SUCCESS;
 }

src/faust_linear_operator/CPU/faust_MatButterfly.hpp

@@ -292,8 +292,7 @@ namespace Faust
 	template<typename FPP>
 	matvar_t* MatButterfly<FPP, Cpu>::toMatIOVar(bool transpose, bool conjugate, const char *var_name) const
 	{
-		//TODO
-
+		throw std::runtime_error("Saving MatButterfly to a .mat file is not supported.");
 	}
 
 	template<typename FPP>
@@ -305,7 +304,8 @@ namespace Faust
 	template<typename FPP>
 	Real<FPP> MatButterfly<FPP, Cpu>::norm() const
 	{
-		//TODO
+		auto s = (D1.diagonal().array() * D1.diagonal().conjugate().array() + D2.diagonal().array() * D2.diagonal().conjugate().array()).sum();
+		return Faust::fabs(std::sqrt(s));
 	}
 
 	template<typename FPP>