OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-01-05 21:53:27 -06:00
Code Issues Packages Projects Releases Wiki Activity

#8900 Summary Data: fix error in h5-data for time scale above 250 years

Browse Source 
The usual implementation of chrono::system_clock::to_time_t() uses nanoseconds
which will overflow on data with long time spans. Fixed by lowering resolution
to seconds, and adding custom conversion method.

Fixes #8900.
This commit is contained in:
Kristian Bendiksen 2022-05-27 10:20:09 +02:00 committed by Magne Sjaastad
parent 5c224241e1
commit c9aa9285d2
1 changed files with 10 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
ApplicationLibCode/FileInterface/RifHdf5SummaryReader.cpp
View File

@ -90,24 +90,24 @@ std::vector<time_t> RifHdf5SummaryReader::timeSteps() const
double time_unit = 24 * 3600;
using namespace std::chrono;
using TP = time_point<system_clock>;
using DoubSec = duration<double, seconds::period>;
using TP = time_point<system_clock, DoubSec>;
auto timeDeltasInDays = values( "TIME" );
std::vector<std::chrono::system_clock::time_point> timePoints;
// Add custom method to convert from time_point to time_t. The usual implementation of
// chrono::system_clock::to_time_t() uses nanoseconds which will overflow on data with
// long time spans.
auto convertTimePointToTimeT = []( const TP& value ) {
return std::chrono::duration_cast<std::chrono::seconds>( value.time_since_epoch() ).count();
};
TP startDat = std::chrono::system_clock::from_time_t( startDate() );
auto startDat = std::chrono::system_clock::from_time_t( startDate() );
timePoints.reserve( timeDeltasInDays.size() );
for ( const auto& t : timeDeltasInDays )
{
timePoints.push_back( startDat + duration_cast<TP::duration>( DoubSec( t * time_unit ) ) );
}
for ( const auto& d : timePoints )
{
auto timeAsTimeT = std::chrono::system_clock::to_time_t( d );
TP timePoint = startDat + duration_cast<TP::duration>( DoubSec( t * time_unit ) );
std::time_t timeAsTimeT = convertTimePointToTimeT( timePoint );
times.push_back( timeAsTimeT );
}