| [6bb72a] | 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 |  | 
|---|
| [4f82f8] | 19 | #include "FunctionApproximation/FunctionArgument.hpp" | 
|---|
| [fdd23a] | 20 | #include "FunctionApproximation/FunctionModel.hpp" | 
|---|
|  | 21 | #include "Potentials/SerializablePotential.hpp" | 
|---|
| [4f82f8] | 22 |  | 
|---|
| [6bb72a] | 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 | */ | 
|---|
| [fdd23a] | 32 | class EmpiricalPotential : | 
|---|
|  | 33 | public FunctionModel, | 
|---|
|  | 34 | public SerializablePotential | 
|---|
| [6bb72a] | 35 | { | 
|---|
|  | 36 | public: | 
|---|
| [3ccea3] | 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 |  | 
|---|
| [6bb72a] | 46 | /** Default constructor for class EmpiricalPotential. | 
|---|
|  | 47 | * | 
|---|
|  | 48 | */ | 
|---|
| [fdd23a] | 49 | EmpiricalPotential(const ParticleTypes_t &_ParticleTypes) : | 
|---|
|  | 50 | SerializablePotential(_ParticleTypes) | 
|---|
|  | 51 | {} | 
|---|
| [6bb72a] | 52 |  | 
|---|
|  | 53 | /** Destructor for class EmpiricalPotential. | 
|---|
|  | 54 | * | 
|---|
|  | 55 | */ | 
|---|
|  | 56 | virtual ~EmpiricalPotential() {} | 
|---|
| [3ccea3] | 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; | 
|---|
| [713888] | 74 |  | 
|---|
|  | 75 | protected: | 
|---|
|  | 76 | /** Default constructor for class EmpiricalPotential. | 
|---|
|  | 77 | * | 
|---|
|  | 78 | * Callable only by derived functions. | 
|---|
|  | 79 | * | 
|---|
|  | 80 | */ | 
|---|
|  | 81 | EmpiricalPotential() | 
|---|
|  | 82 | {} | 
|---|
|  | 83 |  | 
|---|
|  | 84 |  | 
|---|
| [6bb72a] | 85 | }; | 
|---|
|  | 86 |  | 
|---|
|  | 87 | #endif /* EMPIRICALPOTENTIAL_HPP_ */ | 
|---|