| 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 | 
 | 
|---|
| 71 | #endif /* EMPIRICALPOTENTIAL_HPP_ */
 | 
|---|