diff --git a/CImg.h b/CImg.h
index b831c173..02570d6f 100644
--- a/CImg.h
+++ b/CImg.h
@@ -32224,8 +32224,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-25;
-
+      const Ttfloat eps = (Ttfloat)1e-15;
       if (is_empty()) { U.assign(); S.assign(); V.assign(); }
       else if (_depth!=1 || _spectrum!=1)
         throw CImgInstanceException(_cimg_instance
diff --git a/examples/dtmri_view3d.cpp b/examples/dtmri_view3d.cpp
index d1636aa5..be5b9245 100644
--- a/examples/dtmri_view3d.cpp
+++ b/examples/dtmri_view3d.cpp
@@ -317,6 +317,7 @@ int main(int argc,char **argv) {
                             x,y,z);
     }
   }
+
   float voxw = 1, voxh = 1, voxd = 1;
   std::sscanf(vsize,"%f%*c%f%*c%f",&voxw,&voxh,&voxd);
 
@@ -361,6 +362,7 @@ int main(int argc,char **argv) {
     coloredFA(x,y,z,2) = (unsigned char)b;
   }
   tensors.assign();
+
   std::fprintf(stderr,"\n- Maximum diffusivity = %g, Maximum FA = %g",eigmax,eigen.get_shared_channel(12).max());
   if (normalize) {
     std::fprintf(stderr,"\n- Normalize tensors.");