Merge pull request #8387 from OPM/qtcharts-summary-plots

Closes #8228 

Major refactoring of summary plotting. Now possible to create plots both with Qwt and QtChart as plotting tool.

ApplicationLibCode/ProjectDataModel/RimPlotAxisLogRangeCalculator.cpp

@@ -0,0 +1,116 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  Copyright (C) 2022-     Equinor ASA
+//
+//  ResInsight is free software: you can redistribute it and/or modify
+//  it under the terms of the GNU General Public License as published by
+//  the Free Software Foundation, either version 3 of the License, or
+//  (at your option) any later version.
+//
+//  ResInsight is distributed in the hope that it will be useful, but WITHOUT ANY
+//  WARRANTY; without even the implied warranty of MERCHANTABILITY or
+//  FITNESS FOR A PARTICULAR PURPOSE.
+//
+//  See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
+//  for more details.
+//
+/////////////////////////////////////////////////////////////////////////////////
+
+#include "RimPlotAxisLogRangeCalculator.h"
+
+#include "RiaPlotDefines.h"
+
+#include "RimPlotCurve.h"
+
+#include "cvfVector2.h"
+
+#include <cmath>
+
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+RimPlotAxisLogRangeCalculator::RimPlotAxisLogRangeCalculator( RiaDefines::PlotAxis                    axis,
+                                                              const std::vector<const RimPlotCurve*>& curves )
+    : m_axis( axis )
+    , m_curves( curves )
+{
+}
+
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+void RimPlotAxisLogRangeCalculator::computeAxisRange( double* minPositive, double* max ) const
+{
+    double minPosValue = HUGE_VAL;
+    double maxValue    = -HUGE_VAL;
+
+    for ( const RimPlotCurve* curve : m_curves )
+    {
+        double minPosCurveValue = HUGE_VAL;
+        double maxCurveValue    = -HUGE_VAL;
+
+        if ( curveValueRange( curve, &minPosCurveValue, &maxCurveValue ) )
+        {
+            if ( minPosCurveValue < minPosValue )
+            {
+                CVF_ASSERT( minPosCurveValue > 0.0 );
+                minPosValue = minPosCurveValue;
+            }
+
+            if ( maxCurveValue > maxValue )
+            {
+                maxValue = maxCurveValue;
+            }
+        }
+    }
+
+    if ( minPosValue == HUGE_VAL )
+    {
+        minPosValue = RiaDefines::minimumDefaultLogValuePlot();
+        maxValue    = RiaDefines::maximumDefaultValuePlot();
+    }
+
+    *minPositive = minPosValue;
+    *max         = maxValue;
+}
+
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+bool RimPlotAxisLogRangeCalculator::curveValueRange( const RimPlotCurve* curve, double* minPositive, double* max ) const
+{
+    if ( !curve ) return false;
+
+    if ( curve->dataSize() < 1 )
+    {
+        return false;
+    }
+
+    float minPosF = std::numeric_limits<float>::infinity();
+    float maxF    = -std::numeric_limits<float>::infinity();
+
+    int axisValueIndex = 0;
+    if ( RiaDefines::isVertical( m_axis ) )
+    {
+        axisValueIndex = 1;
+    }
+
+    for ( int i = 0; i < curve->dataSize(); ++i )
+    {
+        auto [x, y] = curve->sample( i );
+        cvf::Vec2f vec( x, y );
+        float      value = vec[axisValueIndex];
+        if ( value == HUGE_VALF ) continue;
+
+        maxF = std::max( maxF, value );
+        if ( value > 0.0f && value < minPosF )
+        {
+            minPosF = value;
+        }
+    }
+
+    *minPositive = minPosF;
+    *max         = maxF;
+
+    return true;
+}