src/option.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

#include "option.hpp"
#include <iostream>

double pos (double x){
    return x>0?x:0;
}

double frac_part(double x){
    return x - floor(x);
}

int main(){
    init_alea(1);
    asian_option A(0.05, 1.0, 50.0, 0.1, 16, 45);
    int N= 10000;
    
    int d =16; 

    
    std::vector<double> result(3);
    result = monte_carlo(100, quasi_option<asian_option, sobol> (N, d, A));
    for(int i =0; i<3; i++){
        std::cout<<result[i]<<std::endl;
    }
    
    std::vector<double> result2(3);
    result2 = monte_carlo(100, quasi_option<asian_option, halton> (N, d, A));
    for(int i =0; i<3; i++){
        std::cout<<result2[i]<<std::endl;
    }


    return 0;
}