.

Author: dtschump

CImg.h

@@ -21532,13 +21532,13 @@ namespace cimg_library {
               _cimg_mp_return_nan();
             }
 
-            if (!std::strncmp(ss,"eig(",4)) { // Matrix eigenvalues/eigenvector
-              _cimg_mp_op("Function 'eig()'");
-              arg1 = compile(ss4,se1,depth1,0,block_flags);
+            if (!std::strncmp(ss,"eigen(",6)) { // Matrix eigenvalues/eigenvector
+              _cimg_mp_op("Function 'eigen()'");
+              arg1 = compile(ss6,se1,depth1,0,block_flags);
               _cimg_mp_check_matrix_square(arg1,1);
               p1 = (unsigned int)cimg::round(std::sqrt((float)size(arg1)));
               pos = vector((p1 + 1)*p1);
-              CImg<ulongT>::vector((ulongT)mp_matrix_eig,pos,arg1,p1).move_to(code);
+              CImg<ulongT>::vector((ulongT)mp_matrix_eigen,pos,arg1,p1).move_to(code);
               return_comp = true;
               _cimg_mp_return(pos);
             }
@@ -28751,7 +28751,7 @@ namespace cimg_library {
         return cimg::type<double>::nan();
       }
 
-      static double mp_matrix_eig(_cimg_math_parser& mp) {
+      static double mp_matrix_eigen(_cimg_math_parser& mp) {
         double *ptrd = &_mp_arg(1) + 1;
         const double *ptr1 = &_mp_arg(2) + 1;
         const unsigned int k = (unsigned int)mp.opcode[3];
@@ -33440,7 +33440,7 @@ namespace cimg_library {
     const CImg<T>& SVD(CImg<t>& U, CImg<t>& S, CImg<t>& V, const bool sorting=true,
                        const unsigned int max_iteration=40, const float lambda=0) const {
       typedef _cimg_Ttfloat Ttfloat;
-      const Ttfloat eps = (Ttfloat)1e-15;
+      const Ttfloat eps = (Ttfloat)1e-8f;
       if (is_empty()) { U.assign(); S.assign(); V.assign(); }
       else if (_depth!=1 || _spectrum!=1)
         throw CImgInstanceException(_cimg_instance
@@ -33475,7 +33475,7 @@ namespace cimg_library {
               U(i,i) = f - g;
               for (int j = l; j<width(); ++j) {
                 s = 0;
-                for (int k=i; k<height(); ++k) s+=U(i,k)*U(j,k);
+                for (int k = i; k<height(); ++k) s+=U(i,k)*U(j,k);
                 f = s/h;
                 for (int k = i; k<height(); ++k) U(j,k)+=f*U(i,k);
               }
@@ -33548,8 +33548,8 @@ namespace cimg_library {
             bool flag = true;
             for (l = k; l>=1; --l) {
               nm = l - 1;
-              if ((cimg::abs(rv1[l]) + anorm)==anorm) { flag = false; break; }
-              if ((cimg::abs(S[nm]) + anorm)==anorm) break;
+              if (l==1 || cimg::abs(rv1[l])<=eps*anorm) { flag = false; break; }
+              if (cimg::abs(S[nm])<=eps*anorm) break;
             }
             if (flag) {
               c = 0;

html/header.html

@@ -23,7 +23,7 @@
     <div class="header">
 <a href="index.html"><img alt="Logo" src="img/logo_header.jpg" class="center_image" style="margin-top:1em;"/></a>
       <h2  style="padding-bottom: 1em">
-        Latest stable version: <b><a href="http://cimg.eu/files/CImg_.zip">3.5.5</a></b> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Current pre-release: <b><a href="http://cimg.eu/files/CImg_latest.zip">3.6.0</a></b> (2025/07/18)
+        Latest stable version: <b><a href="http://cimg.eu/files/CImg_.zip">3.5.5</a></b> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Current pre-release: <b><a href="http://cimg.eu/files/CImg_latest.zip">3.6.0</a></b> (2025/07/19)
       </h2>
 
       <hr/>

html/header_doxygen.html

@@ -26,7 +26,7 @@
     <div class="header">
 <a href="../index.html"><img alt="Logo" src="../img/logo_header.jpg" class="center_image" style="margin-top:1em;"/></a>
       <h2  style="padding-bottom: 1em">
-        Latest stable version: <b><a href="http://cimg.eu/files/CImg_.zip">3.5.5</a></b> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Current pre-release: <b><a href="http://cimg.eu/files/CImg_latest.zip">3.6.0</a></b> (2025/07/18)
+        Latest stable version: <b><a href="http://cimg.eu/files/CImg_.zip">3.5.5</a></b> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Current pre-release: <b><a href="http://cimg.eu/files/CImg_latest.zip">3.6.0</a></b> (2025/07/19)
       </h2>
 
       <hr/>