src/opti.cpp


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

#include "opti.hpp"
#include "option.hpp"


double f (const std::vector<double> &X, std::vector<double> &grad, void *params) {
    int d = X.size();
    double *p =(double *)params;
    double r=p[0] ; double T=p[1]; double S0=p[2];
    double V=p[3]; double K=p[4];
    double norm = 0;
    for(int i=0; i<d; i++){
        norm+=X[i]*X[i];
    }
    asian_option A(r, T, S0, V, d, K);
    return log(A(X)) - 0.5*norm;
};


std::vector<double> argmax(double r, double T, double S0, double V, double K, int d){
    
    double params[5] = {r, T, S0, V, K};
    
    nlopt::opt opt(nlopt::LN_COBYLA, d);

    opt.set_max_objective(f, &params);

    opt.set_xtol_rel(1e-4);

    std::vector<double> x(d,0);
    std::vector<double> g(0);

    double maxf;
    nlopt::result result = opt.optimize(x, maxf);

    return x;
}