849bc7df93
this is meant for those nasty keywords where the dimension of an entry depends on a user-defined value of some field of a (potentially different) keyword. One example for this are the surface rates of the produced fluids for the .CON(INJ|PROD).* keywords which exhibit different units depending on whether the user choses to control for the surface gas or the liquid rate. the approach taken in this patch is to convert all numbers to NaN if the unit is queried in SI (via item->getSIDouble()). It might be more desireable to throw an exception in this case, but this approach would be more elaborate and NaNs should be quickly noticeable by the users of this code. I tried to adapt all implemented keywords, but it's quite likely that I missed some...
77 lines
2.1 KiB
C++
77 lines
2.1 KiB
C++
/*
|
|
Copyright 2013 Statoil ASA.
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM 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.
|
|
|
|
OPM 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 for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#include <opm/parser/eclipse/Units/Dimension.hpp>
|
|
|
|
#include <string>
|
|
#include <stdexcept>
|
|
#include <cmath>
|
|
|
|
namespace Opm {
|
|
|
|
Dimension::Dimension() {
|
|
|
|
}
|
|
|
|
Dimension::Dimension(const std::string& name, double SIfactor)
|
|
{
|
|
for (auto iter = name.begin(); iter != name.end(); ++iter) {
|
|
if (!isalpha(*iter) && (*iter) != '1')
|
|
throw std::invalid_argument("Invalid dimension name");
|
|
}
|
|
m_name = name;
|
|
m_SIfactor = SIfactor;
|
|
}
|
|
|
|
|
|
double Dimension::getSIScaling() const {
|
|
if (!std::isfinite(m_SIfactor))
|
|
throw std::logic_error("The DeckItem contains a field with a context dependent unit. "
|
|
"Use getRawDoubleData() and convert the returned value manually!");
|
|
return m_SIfactor;
|
|
}
|
|
|
|
const std::string& Dimension::getName() const {
|
|
return m_name;
|
|
}
|
|
|
|
|
|
Dimension * Dimension::newComposite(const std::string& dim , double SIfactor) {
|
|
Dimension * dimension = new Dimension();
|
|
dimension->m_name = dim;
|
|
dimension->m_SIfactor = SIfactor;
|
|
|
|
return dimension;
|
|
}
|
|
|
|
|
|
bool Dimension::equal(const Dimension& other) const {
|
|
if (m_name != other.m_name)
|
|
return false;
|
|
if (m_SIfactor == other.m_SIfactor)
|
|
return true;
|
|
if (std::isnan(m_SIfactor) && std::isnan(other.m_SIfactor))
|
|
return true;
|
|
return false;
|
|
}
|
|
|
|
}
|
|
|
|
|