Files
rdkit/Code/RDBoost/Wrap.cpp
Paolo Tosco 19c9a3905c Enhanced generateDepictionMatching2DStructure functionality (#3811)
* - generateDepictionMatching2DStructure can be used with referencePattern smaller than reference
  to only use part of a scaffold
- adds generateDepictionMatching2DStructure overload to pass a matchVect instead of doing a substructure match
- adds allowRGroups parameter to enable using a scaffold bearing R groups as reference

* changes in response to review

* added comments

* fixes failing doctest

* - reverted change committed accidentally
- fixed get_sss_json for the case where R groups are not included in the match
- added tests for the return value of generate_aligned_coords

* Documented the value returned by GenerateDepictionMatching2DStructure (Python)

* changes in response to review

* - changes in response to review
- fixed sortMatchesByDegreeOfCoreSubstitution that was not working
- added Python wrappers for sortMatchesByDegreeOfCoreSubstitution and getMostSubstitutedCoreMatch
- added C++ and Pyhon unit tests for the above

* added missing variable initialization

Co-authored-by: Tosco, Paolo <paolo.tosco@novartis.com>
2021-02-24 05:37:31 +01:00

141 lines
4.0 KiB
C++

// $Id$
//
// Copyright (c) 2003-2006 greg Landrum and Rational Discovery LLC
//
// @@ All Rights Reserved @@
// This file is part of the RDKit.
// The contents are covered by the terms of the BSD license
// which is included in the file license.txt, found at the root
// of the RDKit source tree.
//
//
// Generic Wrapper utility functionality
//
#include "Wrap.h"
#include "pyint_api.h"
#include <RDBoost/PySequenceHolder.h>
#include <sstream>
#include <iostream>
// A helper function for dealing with errors. Throw a Python IndexError
void throw_index_error(int key) {
PyErr_SetObject(PyExc_IndexError, PyInt_FromLong(key));
python::throw_error_already_set();
}
// A helper function for dealing with errors. Throw a Python ValueError
void throw_value_error(const std::string err) {
PyErr_SetString(PyExc_ValueError, err.c_str());
python::throw_error_already_set();
}
// A helper function for dealing with errors. Throw a Python KeyError
void throw_key_error(const std::string key) {
PyErr_SetString(PyExc_KeyError, key.c_str());
python::throw_error_already_set();
}
void translate_index_error(IndexErrorException const& e) {
throw_index_error(e.index());
}
void translate_value_error(ValueErrorException const& e) {
throw_value_error(e.what());
}
void translate_key_error(KeyErrorException const& e) {
throw_key_error(e.key());
}
#ifdef INVARIANT_EXCEPTION_METHOD
// A helper function for dealing with errors. Throw a Python RuntimeError
void throw_runtime_error(const std::string err) {
PyErr_SetString(PyExc_RuntimeError, err.c_str());
python::throw_error_already_set();
}
void translate_invariant_error(Invar::Invariant const& e) {
throw_runtime_error(e.toUserString());
}
boost::dynamic_bitset<> pythonObjectToDynBitset(const python::object &obj,
boost::dynamic_bitset<>::size_type maxV) {
boost::dynamic_bitset<> res(maxV);
if (obj) {
python::stl_input_iterator<boost::dynamic_bitset<>::size_type> beg(obj), end;
while (beg != end) {
auto v = *beg;
if (v >= maxV) {
throw_value_error("list element larger than allowed value");
}
res.set(v);
++beg;
}
}
return res;
}
std::vector<std::pair<int, int>> *translateAtomMap(
const python::object &atomMap) {
PySequenceHolder<python::object> pyAtomMap(atomMap);
std::vector<std::pair<int, int>>* res;
res = nullptr;
unsigned int i;
unsigned int n = pyAtomMap.size();
if (n > 0) {
res = new std::vector<std::pair<int, int>>;
for (i = 0; i < n; ++i) {
PySequenceHolder<int> item(pyAtomMap[i]);
if (item.size() != 2) {
delete res;
res = nullptr;
throw_value_error("Incorrect format for an atomMap");
}
res->push_back(std::pair<int, int>(item[0], item[1]));
}
}
return res;
}
std::vector<std::vector<std::pair<int, int>>> translateAtomMapSeq(
const python::object &atomMapSeq) {
std::vector<std::vector<std::pair<int, int>>> aMapVec;
PySequenceHolder<python::object> pyAtomMapSeq(atomMapSeq);
for (size_t i = 0; i < pyAtomMapSeq.size(); ++i) {
std::vector<std::pair<int, int>>* res = translateAtomMap(pyAtomMapSeq[i]);
aMapVec.push_back(*res);
delete res;
}
return aMapVec;
}
RDNumeric::DoubleVector *translateDoubleSeq(const python::object &doubleSeq) {
PySequenceHolder<double> doubles(doubleSeq);
unsigned int nDoubles = doubles.size();
RDNumeric::DoubleVector *doubleVec;
doubleVec = nullptr;
unsigned int i;
if (nDoubles > 0) {
doubleVec = new RDNumeric::DoubleVector(nDoubles);
for (i = 0; i < nDoubles; ++i) {
doubleVec->setVal(i, doubles[i]);
}
}
return doubleVec;
}
std::vector<unsigned int> *translateIntSeq(const python::object &intSeq) {
PySequenceHolder<unsigned int> ints(intSeq);
std::vector<unsigned int> *intVec = nullptr;
if (ints.size() > 0) {
intVec = new std::vector<unsigned int>;
for (unsigned int i = 0; i < ints.size(); ++i) {
intVec->push_back(ints[i]);
}
}
return intVec;
}
#endif