mirror of
https://github.com/rdkit/rdkit.git
synced 2026-06-04 21:54:27 +08:00
* add ROMol::atomBonds() and ROMol::atomNeighbors() methods * remove some warnings * start using the new code * add default for those template params * some more applications * get the SWIG builds working * get rid of extraneous ref * remove extraneous comments
60 lines
1.5 KiB
C++
60 lines
1.5 KiB
C++
//
|
|
// Copyright (C) 2021 Greg Landrum and other RDKit contributors
|
|
//
|
|
// @@ 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.
|
|
//
|
|
|
|
#include "catch.hpp"
|
|
|
|
#include <GraphMol/RDKitBase.h>
|
|
#include <GraphMol/SmilesParse/SmilesParse.h>
|
|
#include <GraphMol/SmilesParse/SmilesWrite.h>
|
|
#include <algorithm>
|
|
|
|
using namespace RDKit;
|
|
|
|
TEST_CASE("mol.atoms()") {
|
|
const auto m = "CC(C)CO"_smiles;
|
|
REQUIRE(m);
|
|
unsigned int ccount = 0;
|
|
for (const auto atom : m->atoms()) {
|
|
if (atom->getAtomicNum() == 6) {
|
|
++ccount;
|
|
}
|
|
}
|
|
CHECK(ccount == 4);
|
|
}
|
|
|
|
TEST_CASE("mol.atomNeighbors()") {
|
|
const auto m = "CC(C)CO"_smiles;
|
|
REQUIRE(m);
|
|
unsigned int count = 0;
|
|
for (const auto atom : m->atomNeighbors(m->getAtomWithIdx(1))) {
|
|
count += atom->getDegree();
|
|
}
|
|
CHECK(count == 4);
|
|
for (auto atom : m->atomNeighbors(m->getAtomWithIdx(1))) {
|
|
atom->setAtomicNum(7);
|
|
}
|
|
MolOps::sanitizeMol(*m);
|
|
CHECK(MolToSmiles(*m) == "NC(N)NO");
|
|
}
|
|
|
|
TEST_CASE("mol.atomBonds()") {
|
|
const auto m = "CC(=C)CO"_smiles;
|
|
REQUIRE(m);
|
|
double count = 0;
|
|
for (const auto bond : m->atomBonds(m->getAtomWithIdx(1))) {
|
|
count += bond->getBondTypeAsDouble();
|
|
}
|
|
CHECK(count == 4);
|
|
for (auto bond : m->atomBonds(m->getAtomWithIdx(1))) {
|
|
bond->setBondType(Bond::BondType::SINGLE);
|
|
}
|
|
MolOps::sanitizeMol(*m);
|
|
CHECK(MolToSmiles(*m) == "CC(C)CO");
|
|
} |