In supress V1 for root poly, Algo OK

micmacIGN · Dec 21, 2024 · 81c6749 · 81c6749
1 parent 88628cf
commit 81c6749
Show file tree

Hide file tree

Showing 4 changed files with 336 additions and 18 deletions.
diff --git a/MMVII/ExternalInclude/eigen-3.4.0/unsupported/Eigen/src/Polynomials/PolynomialUtils.h b/MMVII/ExternalInclude/eigen-3.4.0/unsupported/Eigen/src/Polynomials/PolynomialUtils.h
@@ -53,7 +53,8 @@ T poly_eval( const Polynomials& poly, const T& x )
   {
     T val=poly[0];
     T inv_x = T(1)/x;
-    for( DenseIndex i=1; i<poly.size(); ++i ){
+    // for( DenseIndex i=1; i<poly.size(); ++i ){  MPD::MMVII
+    for( DenseIndex i=1; i<DenseIndex(poly.size()); ++i ){
       val = val*inv_x + poly[i]; }
 
     return numext::pow(x,(T)(poly.size()-1)) * val;

diff --git a/MMVII/include/MMVII_nums.h b/MMVII/include/MMVII_nums.h
@@ -1008,12 +1008,16 @@ template <class Type> class  cPolynom
 
 
            Type  Value(const Type & aVal) const;
+           /// return som(|a_k x^k|) , used for some bounding stuffs
+           Type  AbsValue(const Type & aVal) const;
 
            cPolynom<Type> operator * (const cPolynom<Type> & aP2) const;
            cPolynom<Type> operator + (const cPolynom<Type> & aP2) const;
            cPolynom<Type> operator - (const cPolynom<Type> & aP2) const;
            cPolynom<Type> operator * (const  Type & aVal) const;
-           std::vector<Type> RealRoots(const Type & aTol,int ItMax);
+           cPolynom<Type> Deriv() const;
+
+           std::vector<Type> RealRoots(const Type & aTol,int ItMax) const;
 
 
            Type&   operator [] (size_t aK) {return mVCoeffs[aK];}

diff --git a/MMVII/src/MMV1/Numerics.cpp b/MMVII/src/MMV1/Numerics.cpp
@@ -54,16 +54,6 @@ INST_V1ROOTS(tREAL8)
 INST_V1ROOTS(tREAL16)
 
 
-/*
-template <class Type> void
-     RealRootsOfRealPolynome
-     (
-         ElSTDNS vector<Type> &  Sols,
-         const ElPolynome<Type>  &aPol,
-         Type                    tol,
-         INT                     ItMax
-     )
-     */
 
 
 };