Use BDCSVD instead of JacobiSVD in MatDense::best_low_rank() (the former is...

Use BDCSVD instead of JacobiSVD in MatDense::best_low_rank() (the former is faster than the latter on large matrices).

src/faust_linear_operator/CPU/faust_MatDense.hpp

@@ -1260,7 +1260,7 @@ bool MatDense<FPP,Cpu>::eq_rows(const MatDense<FPP, Cpu> & other, faust_unsigned
 template<typename FPP>
 void MatDense<FPP, Cpu>::best_low_rank(const int &r, MatDense<FPP,Cpu> &bestX, MatDense<FPP, Cpu> &bestY) const
 {
-	Eigen::JacobiSVD<Eigen::Matrix<FPP, Eigen::Dynamic, Eigen::Dynamic>> svd(this->mat, Eigen::ComputeThinU | Eigen::ComputeThinV);
+	Eigen::BDCSVD<Eigen::Matrix<FPP, Eigen::Dynamic, Eigen::Dynamic>> svd(this->mat, Eigen::ComputeThinU | Eigen::ComputeThinV);
 	if(bestX.getNbRow() != this->getNbRow() || r != bestX.getNbCol())
 		bestX.resize(this->getNbRow(), r);
 	if(bestY.getNbRow() != this->getNbRow() || r != bestY.getNbCol())