tarf_monte_carlo 3.44 → 3.50

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 3983fc474126c8a34f0fce99d3f127b0057f5bf92e35949278a7018bc4c6eba0
-  data.tar.gz: 3071b3646c65f592b9b2dfa75a8cfee0dac798b20a4b6d9d21db8758068c2191
+  metadata.gz: d378d28787178efe7bf21a742d34f732738ca1dd0d6cc88aa82726676fb1c767
+  data.tar.gz: 51ccc1a2352a018d9db80dd1b95acb51739c8ff51e6e80f73a8799444f381b21
 SHA512:
-  metadata.gz: a457cd81fbd653af84c587af0de58495646567a9e9e0616f95226592db687e6122c109c62559867f4589bd81ae2a3bada9f79768f6ee9d1c418230868c5579bc
-  data.tar.gz: eada7cec1fa54069644d41817988789080d127fe003918e5a7fd9d40f24b9001bece1c9090dd0defd006a168716158c731f9e47fc00b2cedc31f92c26cf68522
+  metadata.gz: 6cac6595d2997a12b7d39835fbe6e072e6ab052ad90aa9f0579fb4d3e3355bb0ef947ed818b36308ccde5239aef682f299844fda9fdd4d1f6925dc2238f642fb
+  data.tar.gz: 6df19f6a90f00d43c43b76a86e287ba2358ea63b47e68b75aa5fbc104ed562408dfbc0de249609b4446dafdda6a2fc54ab29cadc00f0d85d6d4025be02bd1b0a

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -43,8 +43,11 @@
 
 #define FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE 13
 #define FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE 14
+#define FX_AMERICAN_BARRIER_WINDOW_KNOCKIN 15
+#define FX_AMERICAN_BARRIER_WINDOW_KNOCKOUT 16
 
 #define DATAPOINTS 200 // data for plotting
+#define BARRIER_DP 4 // data for plotting
 #define INTERVAL 50
 #define SIM_LIMIT 196 // 196 + 4 = 200 simulations nedded
 
@@ -142,6 +145,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   double *USts_array  = ( double* ) malloc( NL * sizeof(double) );
   double *TempNs_array = ( double* ) malloc( NL * sizeof(double) );
   double *Barrier_array = ( double* ) malloc( NL * sizeof(double) );
+  double *Barrier2_array = ( double* ) malloc( NL * sizeof(double) );
   double *Rebate_array = ( double* ) malloc( NL * sizeof(double) );
 
   VALUE Ls = rb_hash_aref(MCInputs, rb_str_new2("leverage_ratios") );
@@ -156,6 +160,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   VALUE LSts = rb_hash_aref(MCInputs, rb_str_new2("lower_strikes") );
   VALUE USts = rb_hash_aref(MCInputs, rb_str_new2("upper_strikes") );
   VALUE Brs = rb_hash_aref(MCInputs, rb_str_new2("barriers") );
+  VALUE Brs2 = rb_hash_aref(MCInputs, rb_str_new2("barriers2") );
   VALUE Rbts = rb_hash_aref(MCInputs, rb_str_new2("rebates") );
 
   for (leg = 0; leg < NL; ++leg) {
@@ -187,6 +192,9 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     for (leg = 0; leg < NL; ++leg) {
       Barrier_array[leg]  = NUM2DBL( rb_ary_entry(Brs, leg) );
       Rebate_array[leg]   = NUM2DBL( rb_ary_entry(Rbts, leg) );
+      if(KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKIN || KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKOUT){
+        Barrier2_array[leg]  = NUM2DBL( rb_ary_entry(Brs2, leg) );
+      }
     }
   }
 
@@ -195,15 +203,15 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   //
   double **metrics;
   if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
-    metrics = ( double** ) malloc( 3 * sizeof(double*) );
-    for( metric = 0; metric < 4; metric++ ) {
+    metrics = ( double** ) malloc( BARRIER_DP * sizeof(double*) );
+    for( metric = 0; metric < BARRIER_DP; metric++ ) {
       metrics[metric] = ( double* ) malloc( DATAPOINTS * sizeof(double) );
       for(leg = 0; leg < DATAPOINTS; leg++) {
         metrics[metric][leg] = 0.0;
       }
     }
   }
-  else{  
+  else{
     metrics = ( double** ) malloc( DATAPOINTS * sizeof(double*) );
     for( metric = 0; metric < DATAPOINTS; metric++ ) {
       metrics[metric] = ( double* ) malloc( NL * sizeof(double) );
@@ -261,6 +269,16 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
         if(knockedLeg_dash == -1 && knockin_dash >= 0){
           knockedLeg_dash = leg;
         }
+      } else if(KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKIN || KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKOUT){
+        profit_loss = Spot;
+        profit_loss_dash = Spot_dash;
+
+        if(knockedLeg == -1 && (Spot <= *(Barrier_array + leg) || Spot >= *(Barrier2_array + leg))){
+          knockedLeg = leg;
+        }
+        if(knockedLeg_dash == -1 && (Spot_dash <= *(Barrier_array + leg) || Spot_dash >= *(Barrier2_array + leg))){
+          knockedLeg_dash = leg;
+        }
       } else if ( KType == DOUBLE_STRIKE_POINTS || KType == DOUBLE_STRIKE_ABSOLUTE || KType == DOUBLE_STRIKE_LEGS ) {
         if ( Spot < *( LSts_array + leg ) ) {
           profit_loss = Spot - (*( LSts_array + leg ));
@@ -418,14 +436,11 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
         sim_dash_pos[NL-1] = get_equivalent_rebate(rebate_conversion_sign, ( *( Rebate_array + 0 ) ), sim_dash[NL-1], dir_sign);
       }
       if ((point_pos < DATAPOINTS) && (sim_count % INTERVAL) == 0){
-        metrics[0][point_pos] = sim[NL-1];
-        metrics[1][point_pos] = (knockedLeg >= 0) ? ( *( Xs_array + (knockedLeg) ) ) : 0;
+        metrics[0][point_pos]   = sim[NL-1];
+        metrics[1][point_pos]   = (knockedLeg >= 0) ? ( *( Xs_array + (knockedLeg) ) ) : 0;
+        metrics[3][point_pos]   = knockedLeg;
         metrics[0][point_pos+1] = sim_dash[NL-1];
         metrics[1][point_pos+1] = (knockedLeg_dash >= 0) ? ( *( Xs_array + (knockedLeg_dash) ) ) : 0;
-
-        metrics[2][point_pos]   = sim_pos[NL-1];
-        metrics[2][point_pos+1] = sim_dash_pos[NL-1];
-        metrics[3][point_pos] = knockedLeg;
         metrics[3][point_pos+1] = knockedLeg_dash;
       }
     } else if(KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
@@ -445,10 +460,12 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       }
 
       if ((point_pos < DATAPOINTS) && (sim_count % INTERVAL) == 0){
-        metrics[0][point_pos] = (knockedLeg >= 0) ? sim[knockedLeg] : sim[NL-1];
-        metrics[1][point_pos] = ( *( Xs_array + 0 ) );
+        metrics[0][point_pos]   = (knockedLeg >= 0) ? sim[knockedLeg] : sim[NL-1];
+        metrics[1][point_pos]   = ( *( Xs_array + 0 ) );
+        metrics[3][point_pos]   = knockedLeg;
         metrics[0][point_pos+1] = (knockedLeg_dash >= 0) ? sim[knockedLeg_dash] : sim[NL-1];
         metrics[1][point_pos+1] = ( *( Xs_array + 0 ) );
+        metrics[3][point_pos+1] = knockedLeg_dash;
       }
 
     } else if( KType == ABSOLUTE || KType == PIVOT_ABSOLUTE || KType == COLLAR_ABSOLUTE || KType == DOUBLE_STRIKE_ABSOLUTE ) {
@@ -581,8 +598,8 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     
     if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
       if ((point_pos < DATAPOINTS) && (sim_count % INTERVAL) == 0){
-        // metrics[2][point_pos]   = sim_pos_sum;
-        // metrics[2][point_pos+1] = sim_dash_pos_sum;
+        metrics[2][point_pos]   = sim_pos_sum;
+        metrics[2][point_pos+1] = sim_dash_pos_sum;
         point_pos += 2;
       }
     }
@@ -631,7 +648,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   VALUE final_metrics = rb_ary_new();
   
   if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
-    for(metric = 0; metric < 4; metric++) {
+    for(metric = 0; metric < BARRIER_DP; metric++) {
       VALUE leg_metrics = rb_ary_new();
 
       for(leg = 0; leg < DATAPOINTS; leg++) {
@@ -663,7 +680,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   //
   // free your arrays from memory once you're done using them
   //
-  int dp = (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE) ? 4 : DATAPOINTS; 
+  int dp = (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE) ? BARRIER_DP : DATAPOINTS; 
   for(metric = 0; metric < dp; metric++) {
     free( metrics[metric] );
   }

data/lib/tarf_monte_carlo/version.rb

@@ -3,5 +3,5 @@
 # gem yank tarf_monte_carlo -v 2.3
 
 module TarfMonteCarlo
-  VERSION = "3.44"
+  VERSION = "3.50"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.44'
+  version: '3.50'
 platform: ruby
 authors:
 - Vivek Routh
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2020-10-06 00:00:00.000000000 Z
+date: 2020-10-07 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -140,7 +140,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.0.8
+rubygems_version: 3.0.6
 signing_key: 
 specification_version: 4
 summary: Monte Carlo Simulation.