Multi stat one sample

ChangeLog.md

@@ -50,6 +50,9 @@
 - *DEC*
   - Adding missing headers in some files of DEC.
     (Roland Denis, [#1349](https://github.com/DGtal-team/DGtal/pull/1349))
+
+- *Math*
+  - Fix possible division by zero in the MultiStatistics class.
 
 # DGtal 0.9.4.1
 

src/DGtal/math/MultiStatistics.ih

@@ -250,8 +250,11 @@ DGtal::MultiStatistics::terminate()
       myExp[ k ] /= mySamples[ k ];
       myExp2[ k ] /= mySamples[ k ];
       myVar[ k ] = myExp2[ k ] - myExp[ k ] * myExp[ k ];
-      myUnbiasedVar[ k ] = mySamples[ k ] * myVar[ k ] 
-	/ ( mySamples[ k ] - 1 );
+      if ( mySamples[ k ] > 1 )
+        myUnbiasedVar[ k ] = mySamples[ k ] * myVar[ k ] / ( mySamples[ k ] - 1 );
+      else
+        myUnbiasedVar[ k ] = 0;
+
     }
     myIsTerminate = true;
 }

src/DGtal/topology/doc/moduleCubicalComplex.dox

@@ -25,7 +25,7 @@ href="https://gforge.liris.cnrs.fr/projects/imagene">ImaGene</a>.
 
 
 The following programs are related to this documentation:
-cubical-complex-collapse.cpp, cubical-complex-illustrations.cpp, 
+cubical-complex-collapse.cpp, cubical-complex-illustrations.cpp,
 testCubicalComplex.cpp, cubicalComplexThinning.cpp.
 
 @section dgtal_ccomplex_sec1  Introduction to cubical complexes