/usr/include/rdkit/GraphMol/FMCS/MaximumCommonSubgraph.h is in librdkit-dev 201603.5-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 | //
// Copyright (C) 2014 Novartis Institutes for BioMedical Research
//
// @@ 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.
//
#pragma once
#include <vector>
#include <string>
#include <stdexcept>
#include "../RDKitBase.h"
#include "FMCS.h"
#include "DebugTrace.h" // algorithm filter definitions
#include "SeedSet.h"
#include "Target.h"
#include "SubstructureCache.h"
#include "DuplicatedSeedCache.h"
#include "MatchTable.h"
#include "TargetMatch.h"
#include "RingMatchTableSet.h"
namespace RDKit {
bool FinalChiralityCheckFunction(const short unsigned c1[],
const short unsigned c2[], const ROMol& mol1,
const FMCS::Graph& query, const ROMol& mol2,
const FMCS::Graph& target,
const MCSParameters* p);
namespace FMCS {
class MaximumCommonSubgraph {
struct MCS { // current result. Reference to a fragment of source molecule
std::vector<const Atom*> Atoms;
std::vector<const Bond*> Bonds;
std::vector<unsigned> AtomsIdx;
std::vector<unsigned> BondsIdx; // need for results and size() only !
const ROMol* QueryMolecule;
std::vector<Target> Targets;
};
unsigned long long To;
MCSProgressData Stat;
MCSParameters Parameters;
unsigned ThresholdCount; // min number of matching
std::vector<const ROMol*> Molecules;
#ifdef FAST_SUBSTRUCT_CACHE
std::vector<unsigned> QueryAtomLabels; // for code Morgan. Value based on
// current functor and parameters
std::vector<unsigned> QueryBondLabels; // for code Morgan. Value based on
// current functor and parameters
SubstructureCache HashCache;
MatchTable QueryAtomMatchTable;
MatchTable QueryBondMatchTable;
RingMatchTableSet RingMatchTables;
#endif
#ifdef DUP_SUBSTRUCT_CACHE
DuplicatedSeedCache DuplicateCache;
#endif
const ROMol* QueryMolecule;
unsigned QueryMoleculeMatchedBonds;
unsigned QueryMoleculeMatchedAtoms;
std::vector<Target> Targets;
SeedSet Seeds;
MCS McsIdx;
public:
#ifdef VERBOSE_STATISTICS_ON
ExecStatistics VerboseStatistics;
#endif
MaximumCommonSubgraph(const MCSParameters* params);
~MaximumCommonSubgraph() { clear(); }
MCSResult find(const std::vector<ROMOL_SPTR>& mols);
const ROMol& getQueryMolecule() const { return *QueryMolecule; }
unsigned getMaxNumberBonds() const { return McsIdx.BondsIdx.size(); }
unsigned getMaxNumberAtoms() const { return McsIdx.AtomsIdx.size(); }
// internal:
bool checkIfMatchAndAppend(Seed& seed);
private:
void clear() {
Targets.clear();
Molecules.clear();
To = nanoClock();
}
void init();
void makeInitialSeeds();
bool createSeedFromMCS(size_t newQueryTarget, Seed& seed);
bool growSeeds(); // returns false if canceled
std::string generateResultSMARTS(const MCS& McsIdx) const;
bool match(Seed& seed);
bool matchIncrementalFast(Seed& seed, unsigned itarget);
};
}
} // namespace RDKit
|