Display input configuration if verbosity level > 1 in GivensFGFT(Parallel)Gen C++ code.

src/algorithm/factorization/faust_GivensFGFTGen.hpp

@@ -66,6 +66,8 @@ const vector<int>& GivensFGFTGen<FPP,DEVICE,FPP2,FPP4>::get_ord_indices()
 template<typename FPP, FDevice DEVICE, typename FPP2, typename FPP4>
 void GivensFGFTGen<FPP,DEVICE,FPP2,FPP4>::compute_facts()
 {
+
+
 	is_D_ordered = false; // facts (re)computed then D must be reordered
 	ite = 0;
 	bool stopping = false;
@@ -86,6 +88,16 @@ void GivensFGFTGen<FPP,DEVICE,FPP2,FPP4>::compute_facts()
 			break;
 		}
 	}
+	if(verbosity > 1)
+	{
+		std::cout << "GivensFGFTGen::compute_facts() end" << std::endl;
+		std::cout << "J: " << J << std::endl;
+		std::cout << "tol: " << stoppingError << std::endl;
+		std::cout << "stopcrit is error: " << stoppingCritIsError << std::endl;
+		std::cout << "relErr: " << errIsRel << std::endl;
+		std::cout << "order: " << is_D_ordered << std::endl;
+		std::cout << "enable_large_Faust: " << enable_large_Faust << std::endl;
+	}
 }
 
 template<typename FPP, FDevice DEVICE, typename FPP2, typename FPP4>

src/algorithm/factorization/faust_GivensFGFTParallelGen.hpp

@@ -5,8 +5,28 @@
 #endif
 
 template<typename FPP, FDevice DEVICE, typename FPP2, typename FPP4>
-GivensFGFTParallelGen<FPP,DEVICE,FPP2,FPP4>::GivensFGFTParallelGen(int t, Faust::GivensFGFTGen<FPP, DEVICE, FPP2, FPP4> & alg) : alg(alg), t(t), fact_nrots(0) 
+GivensFGFTParallelGen<FPP,DEVICE,FPP2,FPP4>::GivensFGFTParallelGen(int t, Faust::GivensFGFTGen<FPP, DEVICE, FPP2, FPP4> & alg) : alg(alg), t(t), fact_nrots(0)
 {
+
+	if(alg.verbosity > 1)
+	{
+		std::cout << "GivensFGFTGenParallelGen ctor:" << std::endl;
+		std::cout << "J: " << alg.J << std::endl;
+		std::cout << "tol: " << alg.stoppingError << std::endl;
+		std::cout << "stopcrit is error: " << alg.stoppingCritIsError << std::endl;
+		std::cout << "relErr: " << alg.errIsRel << std::endl;
+		std::cout << "order: " << alg.D_order_dir << std::endl;
+		std::cout << "enable_large_Faust: " << alg.enable_large_Faust << std::endl;
+		auto dLap = dynamic_cast<MatDense<FPP, DEVICE>*>(&alg.Lap);
+		if(dLap)
+			std::cout << "matrix norm: " << dLap->norm() << std::endl;
+		else
+		{
+			auto sLap = dynamic_cast<MatSparse<FPP, DEVICE>*>(&alg.Lap);
+			if(sLap)
+				std::cout << "matrix norm: " << sLap->norm() << std::endl;
+		}
+	}
 }
 
 template<typename FPP, FDevice DEVICE, typename FPP2, typename FPP4>