1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
// ************************************************************************************************
//
// BornAgain: simulate and fit reflection and scattering
//
//! @file Tests/Functional/Mumufit/ClassicalTestFunctions.cpp
//! @brief Implements set of ObjectiveTestFunctions.
//!
//! @homepage http://www.bornagainproject.org
//! @license GNU General Public License v3 or higher (see COPYING)
//! @copyright Forschungszentrum Jülich GmbH 2018
//! @authors Scientific Computing Group at MLZ (see CITATION, AUTHORS)
//
// ************************************************************************************************
#include "Tests/Functional/Mumufit/ClassicalTestFunctions.h"
#include <cassert>
#include <cmath>
//! RosenBrock function: F(x,y) = 100 (y-x^2)^2 + (1-x)^2
//! start point: F(-1.2,1.0) = 24.20
//! minimum : F(1.0,1.0) = 0.
double TestFunctions::RosenBrock(const std::vector<double>& par)
{
assert(par.size() == 2);
const double x = par[0];
const double y = par[1];
const double tmp1 = y - x * x;
const double tmp2 = 1 - x;
return 100 * tmp1 * tmp1 + tmp2 * tmp2;
}
//! Fourth-degree polynomial
//! F(w,x,y,z) = 100 (y-w^2)^2 + (w-1)^2 + 90 (z-y^2)^2
//! + (1-y)^2 + 10.1 [(x-1)^2 + (z-1)^2]
//! + 19.8 (x-1)(z-1)
//!
//! start point: F(-3,-1,-3,-1) = 19192
//! minimum : F(1,1,1,1) = 0.
double TestFunctions::WoodFour(const std::vector<double>& par)
{
assert(par.size() == 4);
const double w = par[0];
const double x = par[1];
const double y = par[2];
const double z = par[3];
const double w1 = w - 1;
const double x1 = x - 1;
const double y1 = y - 1;
const double z1 = z - 1;
const double tmp1 = x - w * w;
const double tmp2 = z - y * y;
return 100 * tmp1 * tmp1 + w1 * w1 + 90 * tmp2 * tmp2 + y1 * y1 + 10.1 * (x1 * x1 + z1 * z1)
+ 19.8 * x1 * z1;
}
//! Decaying sinus from lmfit tutorial.
double TestFunctions::DecayingSin(double x, const std::vector<double>& par)
{
assert(par.size() == 4);
const double amp = par[0];
const double freq = par[1];
const double phase = par[2];
const double decay = par[3];
return amp * std::sin(x * freq + phase) * std::exp(-x * x * decay);
}
|
