From b64727816c06689a1269b5e9ffc4a897dfea02d9 Mon Sep 17 00:00:00 2001
From: Magne Sjaastad <magne.sjaastad@ceetronsolutions.com>
Date: Fri, 15 Mar 2019 11:50:21 +0100
Subject: [PATCH] Add reference to documentation

---
 .../ResultStatisticsCache/RigStatisticsMath.cpp     | 13 ++++++++++---
 1 file changed, 10 insertions(+), 3 deletions(-)

diff --git a/ApplicationCode/ResultStatisticsCache/RigStatisticsMath.cpp b/ApplicationCode/ResultStatisticsCache/RigStatisticsMath.cpp
index 08c16ad958..4e2505d6bd 100644
--- a/ApplicationCode/ResultStatisticsCache/RigStatisticsMath.cpp
+++ b/ApplicationCode/ResultStatisticsCache/RigStatisticsMath.cpp
@@ -17,11 +17,13 @@
 /////////////////////////////////////////////////////////////////////////////////
 
 #include "RigStatisticsMath.h"
+
+#include "cvfBase.h"
+#include "cvfMath.h"
+
 #include <algorithm>
 #include <assert.h>
 #include <math.h>
-#include "cvfBase.h"
-#include "cvfMath.h"
 
 //--------------------------------------------------------------------------------------------------
 /// A function to do basic statistical calculations 
@@ -254,7 +256,12 @@ double RigHistogramCalculator::calculatePercentil(double pVal)
         {
             double domainValueAtEndOfBin = m_min + (binIdx+1) * binWidth;
             double unusedFractionOfLastBin = (double)(accObsCount - pValObservationCount)/binObsCount;
-            return domainValueAtEndOfBin - unusedFractionOfLastBin*binWidth;
+
+            double histogramBasedEstimate = domainValueAtEndOfBin - unusedFractionOfLastBin*binWidth;
+
+            // See https://resinsight.org/docs/casegroupsandstatistics/#percentile-methods for details
+
+            return histogramBasedEstimate;
         }
     }
     assert(false);