opm-common/python/cxx/deck_keyword.cpp

#include <boost/python/suite/indexing/vector_indexing_suite.hpp>

#include <opm/parser/eclipse/Units/UnitSystem.hpp>

#include <opm/parser/eclipse/Parser/ParserKeyword.hpp>

#include <opm/parser/eclipse/Deck/DeckValue.hpp>
#include <opm/parser/eclipse/Deck/DeckItem.hpp>
#include <opm/parser/eclipse/Deck/DeckKeyword.hpp>
#include <opm/parser/eclipse/Deck/DeckRecord.hpp>
#include <opm/parser/eclipse/Utility/Typetools.hpp>

#include "export.hpp"
#include "converters.hpp"

#include <iostream>


namespace {

/* DeckKeyword */
const DeckRecord& (DeckKeyword::*getRecord)(size_t index) const = &DeckKeyword::getRecord;


const DeckItem& getItem(const DeckRecord& record, size_t index) {
  return record.getItem(index);
}

py::list item_to_pylist( const DeckItem& item )
{
    switch (item.getType())
    {
    case type_tag::integer:
        return iterable_to_pylist( item.getData< int >() );
        break;
    case type_tag::fdouble:
        throw py::type_error("Double list access must be specified by either 'get_raw_data_list' or 'get_SI_data_list'.");
        break;
    case type_tag::string:
        return iterable_to_pylist( item.getData< std::string >() );
        break;
    default:
        throw std::logic_error( "Type not set." );
        break;
    }
}


py::list raw_data_to_pylist( const DeckItem& item) {
    return iterable_to_pylist( item.getData<double>() );
}


py::list SI_data_to_pylist( const DeckItem& item) {
    return iterable_to_pylist( item.getSIDoubleData() );
}


bool is_int(const std::string& s)
{
    return !s.empty() && std::find_if(s.begin(),
        s.end(), [](char c) { return !std::isdigit(c); }) == s.end();
}


void push_string_as_deck_value(std::vector<DeckValue>& record, const std::string str) {

    std::size_t star_pos = str.find('*');
    if (star_pos != std::string::npos) {
        int multiplier = 1;

        std::string mult_str = str.substr(0, star_pos);

        if (mult_str.length() > 0) {
            if (is_int(mult_str))
                multiplier = std::stoi( mult_str );
            else
                throw py::type_error();
        }

        std::string value_str = str.substr(star_pos + 1, str.length());
        DeckValue value;

        if (value_str.length() > 0) {
            if (is_int(value_str))
                value = DeckValue( stoi(value_str) );
            else
                value = DeckValue( stod(value_str) );
        }

        for (int i = 0; i < multiplier; i++)
            record.push_back( value );

    }
    else
        record.push_back( DeckValue(str) );

}


py::array_t<int> get_int_array(const DeckKeyword& kw) {
    return convert::numpy_array( kw.getIntData() );
}

py::array_t<double> get_raw_array(const DeckKeyword& kw) {
    return convert::numpy_array( kw.getRawDoubleData() );
}

py::array_t<double> get_SI_array(const DeckKeyword& kw) {
    return convert::numpy_array( kw.getSIDoubleData() );
}

}

void python::common::export_DeckKeyword(py::module& module) {
    py::class_< DeckKeyword >( module, "DeckKeyword")
        .def(py::init<const ParserKeyword& >())

        .def(py::init([](const ParserKeyword& parser_keyword, py::list record_list, UnitSystem& active_system, UnitSystem& default_system) {

            std::vector< std::vector<DeckValue> > value_record_list;

            for (py::handle record_obj : record_list) {
                 py::list record = record_obj.cast<py::list>();
                 std::vector<DeckValue> value_record;

                 for (const py::handle& value_obj : record) {

                      try {
                          int val_int = value_obj.cast<int>();
                          value_record.push_back( DeckValue(val_int) );
                          continue;
                      }
                      catch (const std::exception& e_int) {}

                      try {
                           double val_double = value_obj.cast<double>();
                           value_record.push_back( DeckValue(val_double) );
                           continue;
                      }
                      catch (const std::exception& e_double) {}

                      try {
                           std::string val_string = value_obj.cast<std::string>();
                           push_string_as_deck_value(value_record, val_string);
                           continue;
                      }
                      catch (const std::exception& e_string) {}

                      throw py::type_error("DeckKeyword: tried to add unkown type to record.");

                 }
                 value_record_list.push_back( value_record );
             }
             return DeckKeyword(parser_keyword, value_record_list, active_system, default_system);
         }  )  )

        .def( "__repr__", &DeckKeyword::name )
        .def( "__str__", &str<DeckKeyword> )
        .def("__iter__",  [] (const DeckKeyword &keyword) { return py::make_iterator(keyword.begin(), keyword.end()); }, py::keep_alive<0,1>())
        .def( "__getitem__", getRecord, ref_internal)
        .def( "__len__", &DeckKeyword::size )
        .def_property_readonly("name", &DeckKeyword::name )

    .def(py::init([](const ParserKeyword& parser_keyword, py::array_t<int> py_data) {
            return DeckKeyword(parser_keyword, convert::vector(py_data));
        } ) )

    .def(py::init([](const ParserKeyword& parser_keyword, py::array_t<double> py_data, UnitSystem& active_system, UnitSystem& default_system) {
            return DeckKeyword(parser_keyword, convert::vector(py_data), active_system, default_system);
        } ) )

    .def("get_int_array", &get_int_array)
    .def("get_raw_array", &get_raw_array)
    .def("get_SI_array", &get_SI_array)
         ;


    py::class_< DeckRecord >( module, "DeckRecord")
        .def( "__repr__", &str<DeckRecord> )
        .def( "__iter__", +[] (const DeckRecord& record) { return py::make_iterator(record.begin(), record.end()); }, py::keep_alive<0,1>())
        .def( "__getitem__", &getItem, ref_internal)
        .def( "__len__", &DeckRecord::size )
        ;


    py::class_< DeckItem >(module, "DeckItem")
        .def( "__len__", &DeckItem::data_size )
        .def("get_str", &DeckItem::get<std::string>)
        .def("get_int", &DeckItem::get<int>)
        .def("get_raw", &DeckItem::get<double>)
        .def("get_SI", &DeckItem::getSIDouble)
        .def("get_data_list", &item_to_pylist)
        .def("get_raw_data_list", &raw_data_to_pylist)
        .def("get_SI_data_list", &SI_data_to_pylist)
        .def("has_value", &DeckItem::hasValue)
        .def("defaulted", &DeckItem::defaultApplied)
        ;


}