HBondGeometric::HBondGeometric(SimInfo* info,
const std::string& filename,
const std::string& sele1,
const std::string& sele2,
double rCut, double thetaCut, int nbins) :
StaticAnalyser(info, filename, nbins),
selectionScript1_(sele1), seleMan1_(info), evaluator1_(info),
selectionScript2_(sele2), seleMan2_(info), evaluator2_(info) {
setOutputName(getPrefix(filename) + ".hbg");
ff_ = info_->getForceField();
evaluator1_.loadScriptString(sele1);
if (!evaluator1_.isDynamic()) {
seleMan1_.setSelectionSet(evaluator1_.evaluate());
}
evaluator2_.loadScriptString(sele2);
if (!evaluator2_.isDynamic()) {
seleMan2_.setSelectionSet(evaluator2_.evaluate());
}
// Set up cutoff values:
rCut_ = rCut;
thetaCut_ = thetaCut;
nHBonds_.resize(nBins_);
nDonor_.resize(nBins_);
nAcceptor_.resize(nBins_);
initializeHistogram();
}
TPZNetworkSquareMidimew :: TPZNetworkSquareMidimew( const TPZComponentId& id,
const TPZString& routerId,
unsigned x,
unsigned y,
unsigned z )
: TPZNetwork(id,routerId,x,y,z)
{
unsigned diametro=getDiameter();
initializeHistogram(diametro);
}
template <typename PointFeature> void
pcl::PyramidFeatureHistogram<PointFeature>::compute ()
{
if (!initializeHistogram ())
return;
for (size_t feature_i = 0; feature_i < input_->points.size (); ++feature_i)
{
std::vector<float> feature_vector;
convertFeatureToVector (input_->points[feature_i], feature_vector);
addFeature (feature_vector);
}
is_computed_ = true;
}
/**
* Initialize a group that corresponds to an image
*/
void REM::init(int width, int height, int dimension, int percentage) {
entropy = 0;
minDist = 0;
maxDist = sqrt(pow(width, 2) + pow(height, 2));
minAngle = 0;
maxAngle = M_PI;
binDimension = dimension;
distanceBinWidth = (double)(maxDist / dimension);
angleBinWidth = (double)(maxAngle / dimension);
samplePercentage = (double) percentage / 100;
if (binDimension > 1) {
initializeHistogram();
} else {
fprintf(stderr, "ERROR: Invalid histogram dimension provided\n");
exit(-1);
}
}
