| 1 | /* | 
|---|
| 2 | * EmpiricalPotential.hpp | 
|---|
| 3 | * | 
|---|
| 4 | *  Created on: Sep 26, 2012 | 
|---|
| 5 | *      Author: heber | 
|---|
| 6 | */ | 
|---|
| 7 |  | 
|---|
| 8 | #ifndef EMPIRICALPOTENTIAL_HPP_ | 
|---|
| 9 | #define EMPIRICALPOTENTIAL_HPP_ | 
|---|
| 10 |  | 
|---|
| 11 |  | 
|---|
| 12 | // include config.h | 
|---|
| 13 | #ifdef HAVE_CONFIG_H | 
|---|
| 14 | #include <config.h> | 
|---|
| 15 | #endif | 
|---|
| 16 |  | 
|---|
| 17 | #include <vector> | 
|---|
| 18 |  | 
|---|
| 19 | #include "FunctionApproximation/FunctionArgument.hpp" | 
|---|
| 20 | #include "FunctionApproximation/FunctionModel.hpp" | 
|---|
| 21 | #include "Potentials/SerializablePotential.hpp" | 
|---|
| 22 |  | 
|---|
| 23 | /** An EmpiricalPotential is a function that is given a vector of objects as | 
|---|
| 24 | * arguments which it uses to evaluate an internal function and returns a | 
|---|
| 25 | * value representing the energy of this configuration indicated by the | 
|---|
| 26 | * arguments. | 
|---|
| 27 | * | 
|---|
| 28 | * It is to be used inside an std::accumulate function after a vector of | 
|---|
| 29 | * arguments (i.e. a vector of a vector) has been prepared initially. | 
|---|
| 30 | * | 
|---|
| 31 | */ | 
|---|
| 32 | class EmpiricalPotential : | 
|---|
| 33 | public FunctionModel, | 
|---|
| 34 | public SerializablePotential | 
|---|
| 35 | { | 
|---|
| 36 | public: | 
|---|
| 37 | //!> typedef for the argument vector as input to the function | 
|---|
| 38 | typedef std::vector<argument_t> arguments_t; | 
|---|
| 39 | //!> typedef for a single result degree of freedom | 
|---|
| 40 | typedef double result_t; | 
|---|
| 41 | //!> typedef for the result vector as returned by the function | 
|---|
| 42 | typedef std::vector<result_t> results_t; | 
|---|
| 43 | //!> typedef for the components of the derivative | 
|---|
| 44 | typedef std::vector<result_t> derivative_components_t; | 
|---|
| 45 |  | 
|---|
| 46 | /** Default constructor for class EmpiricalPotential. | 
|---|
| 47 | * | 
|---|
| 48 | */ | 
|---|
| 49 | EmpiricalPotential(const ParticleTypes_t &_ParticleTypes) : | 
|---|
| 50 | SerializablePotential(_ParticleTypes) | 
|---|
| 51 | {} | 
|---|
| 52 |  | 
|---|
| 53 | /** Destructor for class EmpiricalPotential. | 
|---|
| 54 | * | 
|---|
| 55 | */ | 
|---|
| 56 | virtual ~EmpiricalPotential() {} | 
|---|
| 57 |  | 
|---|
| 58 | /** Evaluates the function with the given \a arguments and the current set of | 
|---|
| 59 | * parameters. | 
|---|
| 60 | * | 
|---|
| 61 | * \param arguments set of arguments as input variables to the function | 
|---|
| 62 | * \return result of the function | 
|---|
| 63 | */ | 
|---|
| 64 | virtual results_t operator()(const arguments_t &arguments) const=0; | 
|---|
| 65 |  | 
|---|
| 66 | /** Evaluates the derivative of the function with the given \a arguments | 
|---|
| 67 | * for each component. | 
|---|
| 68 | * | 
|---|
| 69 | * \param arguments set of arguments as input variables to the function | 
|---|
| 70 | * \return result vector containing the derivative with respect to each | 
|---|
| 71 | *         input comonent of the function | 
|---|
| 72 | */ | 
|---|
| 73 | virtual derivative_components_t derivative(const arguments_t &arguments) const=0; | 
|---|
| 74 |  | 
|---|
| 75 | protected: | 
|---|
| 76 | /** Default constructor for class EmpiricalPotential. | 
|---|
| 77 | * | 
|---|
| 78 | * Callable only by derived functions. | 
|---|
| 79 | * | 
|---|
| 80 | */ | 
|---|
| 81 | EmpiricalPotential() | 
|---|
| 82 | {} | 
|---|
| 83 |  | 
|---|
| 84 |  | 
|---|
| 85 | }; | 
|---|
| 86 |  | 
|---|
| 87 | #endif /* EMPIRICALPOTENTIAL_HPP_ */ | 
|---|