/usr/include/rdkit/GraphMol/ShapeHelpers/ShapeUtils.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 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 | //
// Copyright (C) 2005-2006 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.
//
#ifndef _RD_SHAPE_UTILS_H_20050128_
#define _RD_SHAPE_UTILS_H_20050128_
#include <DataStructs/DiscreteValueVect.h>
#include <vector>
namespace RDGeom {
class Point3D;
class Transform3D;
}
namespace RDKit {
class ROMol;
class Conformer;
namespace MolShapes {
//! Compute the size of the box that can fit the conformation, and offset of the
// box
//! from the origin
void computeConfDimsAndOffset(const Conformer &conf, RDGeom::Point3D &dims,
RDGeom::Point3D &offSet,
const RDGeom::Transform3D *trans = 0,
double padding = 2.5);
//! Compute the a box that will fit the confomer
/*!
\param conf The conformer of interest
\param leftBottom Storage for one extremity of the box
\param rightTop Storage for other extremity of the box
\param trans Optional transformation to be applied to the atom
coordinates
\param padding Padding added on the sides around the conformer
*/
void computeConfBox(const Conformer &conf, RDGeom::Point3D &leftBottom,
RDGeom::Point3D &rightTop,
const RDGeom::Transform3D *trans = 0, double padding = 2.5);
//! Compute the union of two boxes
void computeUnionBox(const RDGeom::Point3D &leftBottom1,
const RDGeom::Point3D &rightTop1,
const RDGeom::Point3D &leftBottom2,
const RDGeom::Point3D &rightTop2,
RDGeom::Point3D &uLeftBottom, RDGeom::Point3D &uRightTop);
//! Compute dimensions of a conformer
/*!
\param conf Conformer of interest
\param padding Padding added to the atom coordinates on all sides
\param center Optionally specify the center
\param ignoreHs if true, ignore the hydrogen atoms in computing the centroid
*/
std::vector<double> getConfDimensions(const Conformer &conf,
double padding = 2.5,
const RDGeom::Point3D *center = 0,
bool ignoreHs = true);
//! Compute the shape tanimoto distance between two molecule based on a
// predefined alignment
/*!
\param mol1 The first molecule of interest
\param mol2 The second molecule of interest
\param confId1 Conformer in the first molecule (defaults to first
conformer)
\param confId2 Conformer in the second molecule (defaults to first
conformer)
\param gridSpacing resolution of the grid used to encode the molecular shapes
\param bitsPerPoint number of bit used to encode the occupancy at each grid
point
defaults to two bits per grid point
\param vdwScale Scaling factor for the radius of the atoms to determine
the base radius
used in the encoding - grid points inside this sphere
carry the maximum occupany
\param stepSize thickness of the each layer outside the base radius, the
occupancy value is decreased
from layer to layer from the maximum value
\param maxLayers the maximum number of layers - defaults to the number
allowed the number of bits
use per grid point - e.g. two bits per grid point will
allow 3 layers
\param ignoreHs if true, ignore the hydrogen atoms in the shape encoding
process
*/
double tanimotoDistance(const ROMol &mol1, const ROMol &mol2, int confId1 = -1,
int confId2 = -1, double gridSpacing = 0.5,
DiscreteValueVect::DiscreteValueType bitsPerPoint =
DiscreteValueVect::TWOBITVALUE,
double vdwScale = 0.8, double stepSize = 0.25,
int maxLayers = -1, bool ignoreHs = true);
//! Compute the shape tanimoto distance between two conformers based on a
// predefined alignment
/*!
\param conf1 The first conformer of interest
\param conf2 The second conformer of interest
\param gridSpacing resolution of the grid used to encode the molecular shapes
\param bitsPerPoint number of bit used to encode the occupancy at each grid
point
\param vdwScale Scaling factor for the radius of the atoms to determine
the base radius
used in the encoding - grid points inside this sphere
carry the maximum occupany
\param stepSize thickness of the each layer outside the base radius, the
occupancy value is decreased
from layer to layer from the maximum value
\param maxLayers the maximum number of layers - defaults to the number
allowed the number of bits
use per grid point - e.g. two bits per grid point will
allow 3 layers
\param ignoreHs if true, ignore the hydrogen atoms in the shape encoding
process
*/
double tanimotoDistance(const Conformer &conf1, const Conformer &conf2,
double gridSpacing = 0.5,
DiscreteValueVect::DiscreteValueType bitsPerPoint =
DiscreteValueVect::TWOBITVALUE,
double vdwScale = 0.8, double stepSize = 0.25,
int maxLayers = -1, bool ignoreHs = true);
//! Compute the shape protrusion distance between two molecule based on a
// predefined alignment
/*!
\param mol1 The first molecule of interest
\param mol2 The second molecule of interest
\param confId1 Conformer in the first molecule (defaults to first
conformer)
\param confId2 Conformer in the second molecule (defaults to first
conformer)
\param gridSpacing resolution of the grid used to encode the molecular shapes
\param bitsPerPoint number of bit used to encode the occupancy at each grid
point
defaults to two bits per grid point
\param vdwScale Scaling factor for the radius of the atoms to determine
the base radius
used in the encoding - grid points inside this sphere
carry the maximum occupany
\param stepSize thickness of the each layer outside the base radius, the
occupancy value is decreased
from layer to layer from the maximum value
\param maxLayers the maximum number of layers - defaults to the number
allowed the number of bits
use per grid point - e.g. two bits per grid point will
allow 3 layers
\param ignoreHs if true, ignore the hydrogen atoms in the shape encoding
process
\param allowReordering if set the order will be automatically updated so that
the value calculated
is the protrusion of the smaller shape from the larger
one.
*/
double protrudeDistance(const ROMol &mol1, const ROMol &mol2, int confId1 = -1,
int confId2 = -1, double gridSpacing = 0.5,
DiscreteValueVect::DiscreteValueType bitsPerPoint =
DiscreteValueVect::TWOBITVALUE,
double vdwScale = 0.8, double stepSize = 0.25,
int maxLayers = -1, bool ignoreHs = true,
bool allowReordering = true);
//! Compute the shape protrusion distance between two conformers based on a
// predefined alignment
/*!
\param conf1 The first conformer of interest
\param conf2 The second conformer of interest
\param gridSpacing resolution of the grid used to encode the molecular shapes
\param bitsPerPoint number of bit used to encode the occupancy at each grid
point
\param vdwScale Scaling factor for the radius of the atoms to determine
the base radius
used in the encoding - grid points inside this sphere
carry the maximum occupany
\param stepSize thickness of the each layer outside the base radius, the
occupancy value is decreased
from layer to layer from the maximum value
\param maxLayers the maximum number of layers - defaults to the number
allowed the number of bits
use per grid point - e.g. two bits per grid point will
allow 3 layers
\param ignoreHs if true, ignore the hydrogen atoms in the shape encoding
process
\param allowReordering if set the order will be automatically updated so that
the value calculated
is the protrusion of the smaller shape from the larger
one.
*/
double protrudeDistance(const Conformer &conf1, const Conformer &conf2,
double gridSpacing = 0.5,
DiscreteValueVect::DiscreteValueType bitsPerPoint =
DiscreteValueVect::TWOBITVALUE,
double vdwScale = 0.8, double stepSize = 0.25,
int maxLayers = -1, bool ignoreHs = true,
bool allowReordering = true);
}
}
#endif
|