AutoCorrelation3DVectorCalculator.hpp

Go to the documentation of this file.

/* 
 * AutoCorrelation3DVectorCalculator.hpp 
 *
 * This file is part of the Chemical Data Processing Toolkit
 *
 * Copyright (C) 2003 Thomas Seidel <thomas.seidel@univie.ac.at>
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public License
 * along with this library; see the file COPYING. If not, write to
 * the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
 * Boston, MA 02111-1307, USA.
 */
 
#ifndef CDPL_DESCR_AUTOCORRELATION3DVECTORCALCULATOR_HPP
#define CDPL_DESCR_AUTOCORRELATION3DVECTORCALCULATOR_HPP
 
#include <cstddef>
#include <iterator>
#include <cmath>
#include <functional>
 
#include "CDPL/Math/Matrix.hpp"
#include "CDPL/Math/VectorArray.hpp"
 
 
namespace CDPL
{
 
    namespace Descr
    {
 
        template <typename T>
        class AutoCorrelation3DVectorCalculator
        {
 
          public:
            typedef T EntityType;
 
            typedef std::function<double(const EntityType&, const EntityType&)> EntityPairWeightFunction;
 
            typedef std::function<const Math::Vector3D&(const EntityType&)> Entity3DCoordinatesFunction;
 
            AutoCorrelation3DVectorCalculator();
 
            void setStartRadius(double start_radius);
 
            double getStartRadius() const;
 
            void setRadiusIncrement(double radius_inc);
 
            double getRadiusIncrement() const;
 
            void setNumSteps(std::size_t num_steps);
 
            std::size_t getNumSteps() const;
 
            void setEntityPairWeightFunction(const EntityPairWeightFunction& func);
 
            void setEntity3DCoordinatesFunction(const Entity3DCoordinatesFunction& func);
 
            template <typename Iter, typename Vec>
            void calculate(Iter beg, Iter end, Vec& vec);
 
          private:
            template <typename Iter>
            void init(Iter beg, Iter end);
 
            double                      startRadius;
            double                      radiusIncrement;
            std::size_t                 numSteps;
            std::size_t                 numEntities;
            EntityPairWeightFunction    weightFunc;
            Entity3DCoordinatesFunction coordsFunc;
            Math::DMatrix               weightMatrix;
            Math::Vector3DArray         entityPositions;
        };
    } // namespace Descr
} // namespace CDPL
 
 
// Implementation
 
template <typename T>
CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::AutoCorrelation3DVectorCalculator():
    startRadius(0.0), radiusIncrement(0.1), numSteps(99)
{}
 
template <typename T>
void CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::setStartRadius(double start_radius)
{
    startRadius = start_radius;
}
 
template <typename T>
void CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::setRadiusIncrement(double radius_inc)
{
    radiusIncrement = radius_inc;
}
 
template <typename T>
void CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::setNumSteps(std::size_t num_steps)
{
    numSteps = num_steps;
}
 
template <typename T>
double CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::getStartRadius() const
{
    return startRadius;
}
 
template <typename T>
double CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::getRadiusIncrement() const
{
    return radiusIncrement;
}
 
template <typename T>
std::size_t CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::getNumSteps() const
{
    return numSteps;
}
 
template <typename T>
void CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::setEntityPairWeightFunction(const EntityPairWeightFunction& func)
{
    weightFunc = func;
}
 
template <typename T>
void CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::setEntity3DCoordinatesFunction(const Entity3DCoordinatesFunction& func)
{
    coordsFunc = func;
}
 
template <typename T>
template <typename Iter, typename Vec>
void CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::calculate(Iter beg, Iter end, Vec& vec)
{
    init(beg, end);
 
    for (std::size_t i = 0; i <= numSteps; i++)
        vec[i] = 0.0;
 
    Math::Vector3D tmp;
    const double   max_radius = startRadius + (numSteps + 1) * radiusIncrement;
 
    for (std::size_t i = 0; i < numEntities; i++) {
        for (std::size_t j = i; j < numEntities; j++) {
            double dist = 0.0;
 
            if (i != j) {
                tmp.assign(entityPositions[i] - entityPositions[j]);
                dist = length(tmp);
            }
 
            if (dist < startRadius || dist >= max_radius)
                continue;
 
            std::size_t idx = std::floor((dist - startRadius) / radiusIncrement);
 
            vec[idx] = weightMatrix(i, j);
        }
    }
}
 
template <typename T>
template <typename Iter>
void CDPL::Descr::AutoCorrelation3DVectorCalculator<T>::init(Iter beg, Iter end)
{
    numEntities = std::distance(beg, end);
 
    weightMatrix.resize(numEntities, numEntities, false);
    entityPositions.resize(numEntities);
 
    for (std::size_t i = 0; beg != end; i++, ++beg) {
        const EntityType& entity1 = *beg;
 
        if (coordsFunc)
            entityPositions[i] = coordsFunc(entity1);
        else
            entityPositions[i].clear(0.0);
 
        std::size_t j = i;
 
        for (Iter it = beg; it != end; ++it, j++) {
            const EntityType& entity2 = *it;
 
            if (!weightFunc)
                weightMatrix(i, j) = 1.0;
            else
                weightMatrix(i, j) = weightFunc(entity1, entity2);
        }
    }
}
 
#endif // CDPL_DESCR_AUTOCORRELATION3DVECTORCALCULATOR_HPP