/////////////////////////////////////////////////////////////////////////////////
//
//  Copyright (C) 2018-     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.
//
/////////////////////////////////////////////////////////////////////////////////

#pragma once

#include "RiaPorosityModel.h"

#include "RigEclipseToStimPlanCellTransmissibilityCalculator.h"

#include "cvfMatrix4.h"

#include <map>
#include <memory>

class QString;

class RimEclipseCase;
class RigFractureGrid;
class RigTransmissibilityCondenser;
class RimFracture;

//==================================================================================================
///
//==================================================================================================
class RigEclipseToStimPlanCalculator
{
public:
    explicit RigEclipseToStimPlanCalculator( const RimEclipseCase*  caseToApply,
                                             cvf::Mat4d             fractureTransform,
                                             double                 skinFactor,
                                             double                 cDarcy,
                                             const RigFractureGrid& fractureGrid,
                                             const RimFracture*     fracture );

    void appendDataToTransmissibilityCondenser( bool useFiniteConductivityInFracture, RigTransmissibilityCondenser* condenser ) const;

    // Returns the area intersecting eclipse cells open for flow, from both active and inactive cells
    // Truncated parts of the fracture are not included
    double totalEclipseAreaOpenForFlow() const;

    double areaWeightedMatrixPermeability() const;
    double areaWeightedWidth() const;
    double areaWeightedConductivity() const;
    double longestYSectionOpenForFlow() const;

private:
    void computeValues();

private:
    const RimEclipseCase*  m_case;
    const RimFracture*     m_fracture;
    double                 m_cDarcy;
    double                 m_fractureSkinFactor;
    cvf::Mat4d             m_fractureTransform;
    const RigFractureGrid& m_fractureGrid;

    std::map<size_t, std::unique_ptr<RigEclipseToStimPlanCellTransmissibilityCalculator>> m_singleFractureCellCalculators;
};
