| 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 |  | 
|---|
| 21 | /** An EmpiricalPotential is a function that is given a vector of objects as | 
|---|
| 22 | * arguments which it uses to evaluate an internal function and returns a | 
|---|
| 23 | * value representing the energy of this configuration indicated by the | 
|---|
| 24 | * arguments. | 
|---|
| 25 | * | 
|---|
| 26 | * It is to be used inside an std::accumulate function after a vector of | 
|---|
| 27 | * arguments (i.e. a vector of a vector) has been prepared initially. | 
|---|
| 28 | * | 
|---|
| 29 | */ | 
|---|
| 30 | class EmpiricalPotential | 
|---|
| 31 | { | 
|---|
| 32 | public: | 
|---|
| 33 | //!> typedef for the argument vector as input to the function | 
|---|
| 34 | typedef std::vector<argument_t> arguments_t; | 
|---|
| 35 | //!> typedef for a single result degree of freedom | 
|---|
| 36 | typedef double result_t; | 
|---|
| 37 | //!> typedef for the result vector as returned by the function | 
|---|
| 38 | typedef std::vector<result_t> results_t; | 
|---|
| 39 | //!> typedef for the components of the derivative | 
|---|
| 40 | typedef std::vector<result_t> derivative_components_t; | 
|---|
| 41 |  | 
|---|
| 42 | /** Default constructor for class EmpiricalPotential. | 
|---|
| 43 | * | 
|---|
| 44 | */ | 
|---|
| 45 | EmpiricalPotential() {} | 
|---|
| 46 |  | 
|---|
| 47 | /** Destructor for class EmpiricalPotential. | 
|---|
| 48 | * | 
|---|
| 49 | */ | 
|---|
| 50 | virtual ~EmpiricalPotential() {} | 
|---|
| 51 |  | 
|---|
| 52 | /** Evaluates the function with the given \a arguments and the current set of | 
|---|
| 53 | * parameters. | 
|---|
| 54 | * | 
|---|
| 55 | * \param arguments set of arguments as input variables to the function | 
|---|
| 56 | * \return result of the function | 
|---|
| 57 | */ | 
|---|
| 58 | virtual results_t operator()(const arguments_t &arguments) const=0; | 
|---|
| 59 |  | 
|---|
| 60 | /** Evaluates the derivative of the function with the given \a arguments | 
|---|
| 61 | * for each component. | 
|---|
| 62 | * | 
|---|
| 63 | * \param arguments set of arguments as input variables to the function | 
|---|
| 64 | * \return result vector containing the derivative with respect to each | 
|---|
| 65 | *         input comonent of the function | 
|---|
| 66 | */ | 
|---|
| 67 | virtual derivative_components_t derivative(const arguments_t &arguments) const=0; | 
|---|
| 68 | }; | 
|---|
| 69 |  | 
|---|
| 70 | #endif /* EMPIRICALPOTENTIAL_HPP_ */ | 
|---|