mirror of
https://github.com/OPM/ResInsight.git
synced 2025-02-10 23:46:00 -06:00
8defab9bc1
#8401 Flow Vector Result: Handle empty results for first time step Use OpenMP for vector calculations Janitor: Use field of double instead of float to avoid rounding errors. If float is used in field, the current editor is has rounding effects that makes the text representation of 0.1 turn into 0.09999999....
83 lines
3.0 KiB
C++
83 lines
3.0 KiB
C++
/////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// Copyright (C) Statoil ASA
|
|
// Copyright (C) Ceetron Solutions AS
|
|
//
|
|
// 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 "RigResultModifierFactory.h"
|
|
|
|
#include "RigCaseCellResultsData.h"
|
|
#include "RigEclipseCaseData.h"
|
|
#include "RigResultModifier.h"
|
|
|
|
#include <cmath>
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
cvf::ref<RigResultModifier> RigResultModifierFactory::createResultModifier( RigEclipseCaseData* eclipseCase,
|
|
size_t gridIndex,
|
|
RiaDefines::PorosityModelType porosityModel,
|
|
size_t timeStepIndex,
|
|
const RigEclipseResultAddress& resVarAddr )
|
|
{
|
|
if ( !eclipseCase ) return nullptr;
|
|
|
|
if ( !eclipseCase->results( porosityModel ) || !eclipseCase->activeCellInfo( porosityModel ) )
|
|
{
|
|
return nullptr;
|
|
}
|
|
|
|
if ( !resVarAddr.isValid() )
|
|
{
|
|
return nullptr;
|
|
}
|
|
|
|
RigGridBase* grid = eclipseCase->grid( gridIndex );
|
|
if ( !grid )
|
|
{
|
|
return nullptr;
|
|
}
|
|
|
|
auto scalarSetResults = eclipseCase->results( porosityModel )->modifiableCellScalarResultTimesteps( resVarAddr );
|
|
|
|
if ( timeStepIndex >= scalarSetResults->size() )
|
|
{
|
|
return nullptr;
|
|
}
|
|
|
|
std::vector<double>* resultValues = nullptr;
|
|
if ( timeStepIndex < scalarSetResults->size() )
|
|
{
|
|
resultValues = &( scalarSetResults->at( timeStepIndex ) );
|
|
}
|
|
|
|
if ( resultValues->empty() ) return nullptr;
|
|
|
|
bool useGlobalActiveIndex = eclipseCase->results( porosityModel )->isUsingGlobalActiveIndex( resVarAddr );
|
|
if ( useGlobalActiveIndex )
|
|
{
|
|
cvf::ref<RigResultModifier> object =
|
|
new RigActiveCellsResultModifier( grid, eclipseCase->activeCellInfo( porosityModel ), resultValues );
|
|
return object;
|
|
}
|
|
else
|
|
{
|
|
cvf::ref<RigResultModifier> object = new RigAllGridCellsResultModifier( grid, resultValues );
|
|
return object;
|
|
}
|
|
}
|