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
101
102
103
104
105
106
107
108
109
110
111
112
113
114
|
/*
* $Id: ForceField.i 1779 2011-07-01 04:43:39Z glandrum $
*
* Copyright (c) 2010, Novartis Institutes for BioMedical Research Inc.
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions are
* met:
*
* * Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
* * Redistributions in binary form must reproduce the above
* copyright notice, this list of conditions and the following
* disclaimer in the documentation and/or other materials provided
* with the distribution.
* * Neither the name of Novartis Institutes for BioMedical Research Inc.
* nor the names of its contributors may be used to endorse or promote
* products derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
* OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
* DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
%{
#include <GraphMol/ForceFieldHelpers/UFF/AtomTyper.h>
#include <GraphMol/ForceFieldHelpers/UFF/Builder.h>
#include <ForceField/ForceField.h>
#include <ForceField/Contrib.h>
#include <ForceField/UFF/AngleBend.h>
#include <ForceField/UFF/BondStretch.h>
#include <ForceField/UFF/Contribs.h>
#include <ForceField/UFF/DistanceConstraint.h>
#include <ForceField/UFF/Nonbonded.h>
#include <ForceField/UFF/Params.h>
#include <ForceField/UFF/TorsionAngle.h>
#include <boost/tuple/tuple.hpp>
#include <GraphMol/ROMol.h>
%}
// We have trouble with the definition as it is -- make it a Point3D_Vect instead
%ignore ForceFields::ForceField::positions();
%ignore ForceFields::ForceField::positions() const;
%rename(get) ForceFields::UFF::ParamCollection::operator();
%include <ForceField/ForceField.h>
%include <ForceField/Contrib.h>
%include <ForceField/UFF/AngleBend.h>
%include <ForceField/UFF/BondStretch.h>
%include <ForceField/UFF/Contribs.h>
%include <ForceField/UFF/DistanceConstraint.h>
%include <ForceField/UFF/Nonbonded.h>
%include <ForceField/UFF/Params.h>
%include <ForceField/UFF/TorsionAngle.h>
%include <GraphMol/ForceFieldHelpers/UFF/AtomTyper.h>
%include <GraphMol/ForceFieldHelpers/UFF/Builder.h>
%template(FF_Contrib_Vect) std::vector<ForceFields::ContribPtr>;
%extend ForceFields::ForceField {
std::vector<RDGeom::Point3D *> &positions3D () {
return ((std::vector<RDGeom::Point3D *> &) ($self)->positions());
}
static int UFFOptimizeMolecule(RDKit::ROMol &mol, int maxIters=200,
double vdwThresh=10.0, int confId=-1,
bool ignoreInterfragInteractions=true ) {
ForceFields::ForceField *ff=RDKit::UFF::constructForceField(mol,vdwThresh, confId, ignoreInterfragInteractions);
ff->initialize();
int res=ff->minimize(maxIters);
delete ff;
return res;
}
static ForceFields::ForceField *UFFGetMoleculeForceField(RDKit::ROMol &mol,
double vdwThresh=10.0,
int confId=-1,
bool ignoreInterfragInteractions=true ) {
ForceFields::ForceField *ff=RDKit::UFF::constructForceField(mol,vdwThresh, confId, ignoreInterfragInteractions);
ff->initialize();
return ff;
}
static bool UFFHasAllMoleculeParams(const RDKit::ROMol &mol){
RDKit::UFF::AtomicParamVect types;
bool foundAll;
boost::tie(types,foundAll)=RDKit::UFF::getAtomTypes(mol);
return foundAll;
}
/* From GraphMol/ForceFieldHelpers/UFF/AtomTyper.h */
static void UFFAddAtomChargeFlags(const RDKit::Atom *atom, std::string &atomKey, bool tolerateChargeMismatch=true) {
RDKit::UFF::Tools::addAtomChargeFlags(atom, atomKey, tolerateChargeMismatch);
}
static std::string UFFGetAtomLabel(const RDKit::Atom *atom) {
return RDKit::UFF::Tools::getAtomLabel(atom);
}
/* From GraphMol/ForceFieldHelpers/UFF/AtomTyper.h */
static std::pair<RDKit::UFF::AtomicParamVect,bool> UFFGetAtomTypes(const RDKit::ROMol &mol, const std::string ¶mData="") {
return RDKit::UFF::getAtomTypes(mol, paramData);
}
}
|
