tarf_monte_carlo 3.46 → 3.52

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 189ce30c759cdc40444b5ffc5443e8d93680adb392be409beb2f8d6445e1f667
-  data.tar.gz: 93b441a9d9481e0595b9b029e28792d2a4bf08299b0bd0ca3396a680bb85184c
+  metadata.gz: 981ac06c2a4fa081cd670f85bc08c2efc714ed6d057618b3347b0e2d6596a84b
+  data.tar.gz: 9127b590ee3001d5ca31cfb97428f9a479b5f598f7376fe004eb4c3613e7c564
 SHA512:
-  metadata.gz: 61030c67b76c62f00b682f135f753d4d82d5fb91e3a5c822e57bf865138befd9b19d574ed9a8c70140fc7ea5aa7a14067a34a69d46ad663165d15849d587d17b
-  data.tar.gz: a2844c3dd6b4a9591cc6bed56edc3a896e8b83b162d7db9f5546573fe8d2e1fcf1bf2338e125a80145835aebfa2f6b7c55bb8f970f5db47e0317f682e242d3fc
+  metadata.gz: fb8304019cfe8d1eed96bf4a7535ea21399b4f6d58e9eb649034a1b8ee6423c12d3b4ea18ee922e8b8351466c1f05f58d96a8999ecf65367a3845bbac2115586
+  data.tar.gz: 3cedfd64107a5fd9639509375eda6b42fda41d04f7339aa2d5a1eae24ec16da62a9f5623e67557e0b8f0b70a979f15c86a01469cdffb9045be3cb2654b36aba3

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -43,6 +43,8 @@
 
 #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
@@ -60,7 +62,7 @@ void Init_tarf_monte_carlo();
 double european_payoff(double, double, int, int, double);
 double get_equivalent_notional(int, double, float);
 double get_equivalent_rebate(int, double, float, int);
-
+bool barrier_variations(int);
 // Prototype for our methods - methods are prefixed by 'method_' here
 VALUE method_box_muller( VALUE );
 VALUE method_run_monte_carlo( VALUE, VALUE );
@@ -99,6 +101,9 @@ double get_equivalent_rebate(int conversion_sign, double rebate, float rate, int
   return (-1 * dir_sign * total);
 }
 
+bool barrier_variations(int KType){
+  return (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE || KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKIN || KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKOUT);
+}
 // main method for running monte carlo simulation from sidekiq worker/outside method
 VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   VALUE MCInputs = rb_ary_shift(args);
@@ -143,6 +148,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") );
@@ -157,6 +163,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) {
@@ -184,10 +191,13 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       LSts_array[leg] = NUM2DBL( rb_ary_entry(LSts, leg) );
       USts_array[leg] = NUM2DBL( rb_ary_entry(USts, leg) );
     }
-  } else if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE) {
+  } else if (barrier_variations(KType)) {
     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,7 +205,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   // first create a 1-D array of pointers, and then, for each array entry, create another 1-D array.
   //
   double **metrics;
-  if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
+  if (barrier_variations(KType)){
     metrics = ( double** ) malloc( BARRIER_DP * sizeof(double*) );
     for( metric = 0; metric < BARRIER_DP; metric++ ) {
       metrics[metric] = ( double* ) malloc( DATAPOINTS * sizeof(double) );
@@ -204,7 +214,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       }
     }
   }
-  else{  
+  else{
     metrics = ( double** ) malloc( DATAPOINTS * sizeof(double*) );
     for( metric = 0; metric < DATAPOINTS; metric++ ) {
       metrics[metric] = ( double* ) malloc( NL * sizeof(double) );
@@ -262,6 +272,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 ));
@@ -385,7 +405,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       //
       // Store spot and spot dash
       //
-      if (!(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE)){
+      if (!(barrier_variations(KType))){
         if( point_pos <= SIM_LIMIT && (sim_count + 2) % INTERVAL == 0 ) {
           // rb_p(rb_str_new2("Leg:Spots"));
           metrics[ point_pos ][ leg ] = Spot;
@@ -403,7 +423,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     // legs loop end
     // start from the Knock value
     double ko_so_far = K, ko_so_far_dash = K;
-    if(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE){
+    if(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKIN){
       if(knockedLeg >= 0){
         double equivalent_notional = get_equivalent_notional(conversion_sign, ( *( Ns_array + knockedLeg ) ), sim[NL-1]);
         sim_pos[NL-1] = european_payoff(( *( Xs_array + (knockedLeg) ) ), sim[NL-1], cp_sign, dir_sign, equivalent_notional);
@@ -419,17 +439,14 @@ 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){
+    } else if(KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE || KType == FX_AMERICAN_BARRIER_WINDOW_KNOCKOUT){
       if(knockedLeg >= 0){
         sim_pos[knockedLeg] = get_equivalent_rebate(rebate_conversion_sign, ( *( Rebate_array + knockedLeg ) ), sim[knockedLeg], dir_sign);
       }
@@ -446,10 +463,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 ) {
@@ -559,7 +578,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       //
       // store and send whichever payoff you want
       //
-      if (!(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE)){
+      if (!(barrier_variations(KType))){
         if( point_pos <= SIM_LIMIT && (sim_count + 2) % INTERVAL == 0 ) {
           // rb_p(rb_str_new2("Leg:Payoffs"));
           metrics[ point_pos + 1 ][ leg ] = (sim_pos[leg] + sim_neg[leg]);
@@ -579,11 +598,11 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     pvs_neg[sim_count]      = sim_neg_sum;
     pvs_neg[sim_count + 1]  = sim_dash_neg_sum;
 
-    
-    if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
+
+    if (barrier_variations(KType)){
       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;
       }
     }
@@ -630,8 +649,8 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   //
   // rb_p(rb_str_new2("Converting metrics"));
   VALUE final_metrics = rb_ary_new();
-  
-  if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
+
+  if (barrier_variations(KType)){
     for(metric = 0; metric < BARRIER_DP; metric++) {
       VALUE leg_metrics = rb_ary_new();
 
@@ -664,7 +683,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) ? BARRIER_DP : DATAPOINTS; 
+  int dp = barrier_variations(KType) ? 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.46"
+  VERSION = "3.52"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.46'
+  version: '3.52'
 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.