tarf_monte_carlo 3.5 → 3.13

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 67b1ad29c02a6eaa79c27b2bcd092d8d58736de856228479553c01275c74aa23
-  data.tar.gz: a0f8a59cc28cec075cdaf3181aa41df49068efa7ef53613fb840d7582f91194d
+  metadata.gz: bf901445c70d80f370053ffa90ba2cc6fb3e12b9ad7c372777cae63b92ceecdb
+  data.tar.gz: 04331203b779e146db0e48a7eb9d1a0452c31aee17311c5f3c863d2104bc080b
 SHA512:
-  metadata.gz: 71691b20d4de8c2c218280998dbd67a2c89d829001b46181186560ed15f8e644b1cc6e9c83b737c4f59e77adfe22e65bc34f2184d9e8d9f1f4d59efff8a2193e
-  data.tar.gz: c8b637dfab5f1153e02bd61ff4e48742cd15507fd3dba45af07bad6a63224ce962ed8afbc88011d241cb82a919f94901fd76fbaa6a8e6f192beb3518c2fa8888
+  metadata.gz: b912836d2d1a64d62c9b1d394be9ca89c5dbeb245312ce18d9490751bf61fd1fc9ade361a60446d1a0b1df7928882e5e4d0691b0e8c7528733f9bdfc49ccf6b2
+  data.tar.gz: '08248d29944ff3fb8981d7bda9506bd056869575e54cb1ec6aa3be403b945cbe1b11c6bcb1cc42193f5617330d1be6633f4f1c2ab331cd2cc70d4cf455ef7df1'

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -28,6 +28,9 @@
 #define ABSOLUTE 1 // knockout by absolute
 #define POINTS 2 // knockout by points
 #define LEGS 3 // knockout by legs
+#define PIVOT_POINTS 4
+#define COLLAR_POINTS 5
+#define DUAL_STRIKE_POINTS 6
 #define DATAPOINTS 200 // data for plotting
 #define INTERVAL 50
 #define SIM_LIMIT 196 // 196 + 4 = 200 simulations nedded
@@ -67,8 +70,8 @@ VALUE method_box_muller( VALUE self ) {
 // 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);
-  rb_p(rb_str_new2("MC Inputs:"));
-  rb_p(MCInputs);
+  // rb_p(rb_str_new2("MC Inputs:"));
+  // rb_p(MCInputs);
 
   // seed value for rand() function
   srand( time(0) );
@@ -87,15 +90,19 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   int Ko_compare_mult   = NUM2INT( rb_hash_aref(MCInputs, rb_str_new2("multiplier")) );
 
   // assign leg specific attributes
-  double *pvs_pos   = ( double* ) malloc( SCount * sizeof(double) );
-  double *pvs_neg   = ( double* ) malloc( SCount * sizeof(double) );
-  double *Ls_array  = ( double* ) malloc( NL * sizeof(double) );
-  double *Ts_array  = ( double* ) malloc( NL * sizeof(double) );
-  double *Ns_array  = ( double* ) malloc( NL * sizeof(double) );
-  double *Xs_array  = ( double* ) malloc( NL * sizeof(double) );
-  double *vs_array  = ( double* ) malloc( NL * sizeof(double) );
-  double *bs_array  = ( double* ) malloc( NL * sizeof(double) );
-  double *DFs_array = ( double* ) malloc( NL * sizeof(double) );
+  double *pvs_pos     = ( double* ) malloc( SCount * sizeof(double) );
+  double *pvs_neg     = ( double* ) malloc( SCount * sizeof(double) );
+  double *Ls_array    = ( double* ) malloc( NL * sizeof(double) );
+  double *Ts_array    = ( double* ) malloc( NL * sizeof(double) );
+  double *Ns_array    = ( double* ) malloc( NL * sizeof(double) );
+  double *Xs_array    = ( double* ) malloc( NL * sizeof(double) );
+  double *vs_array    = ( double* ) malloc( NL * sizeof(double) );
+  double *bs_array    = ( double* ) malloc( NL * sizeof(double) );
+  double *DFs_array   = ( double* ) malloc( NL * sizeof(double) );
+  double *Ps_array    = ( double* ) malloc( NL * sizeof(double) );
+  double *CFs_array   = ( double* ) malloc( NL * sizeof(double) );
+  double *LSts_array  = ( double* ) malloc( NL * sizeof(double) );
+  double *USts_array  = ( double* ) malloc( NL * sizeof(double) );
 
   VALUE Ls = rb_hash_aref(MCInputs, rb_str_new2("leverage_ratios") );
   VALUE Ts = rb_hash_aref(MCInputs, rb_str_new2("expiration_times") );
@@ -104,6 +111,10 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   VALUE vs = rb_hash_aref(MCInputs, rb_str_new2("volatilities") );
   VALUE bs = rb_hash_aref(MCInputs, rb_str_new2("carry_rates") );
   VALUE DFs = rb_hash_aref(MCInputs, rb_str_new2("discount_factors") );
+  VALUE Ps = rb_hash_aref(MCInputs, rb_str_new2("pivots") );
+  VALUE CFs = rb_hash_aref(MCInputs, rb_str_new2("caps") );
+  VALUE LSts = rb_hash_aref(MCInputs, rb_str_new2("lower_strikes") );
+  VALUE USts = rb_hash_aref(MCInputs, rb_str_new2("upper_strikes") );
 
   for (leg = 0; leg < NL; ++leg) {
     Ls_array[leg] = NUM2DBL( rb_ary_entry(Ls, leg) );
@@ -114,6 +125,21 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     bs_array[leg] = NUM2DBL( rb_ary_entry(bs, leg) );
     DFs_array[leg] = NUM2DBL( rb_ary_entry(DFs, leg) );
   }
+  // extra tarf structures
+  if ( KType == PIVOT_POINTS )  {
+    for (leg = 0; leg < NL; ++leg) {
+      Ps_array[leg] = NUM2DBL( rb_ary_entry(Ps, leg) );
+    }
+  } else if ( KType == COLLAR_POINTS ) {
+    for (leg = 0; leg < NL; ++leg) {
+      CFs_array[leg] = NUM2DBL( rb_ary_entry(CFs, leg) );
+    }
+  } else if ( KType == DUAL_STRIKE_POINTS ) {
+    for (leg = 0; leg < NL; ++leg) {
+      LSts_array[leg] = NUM2DBL( rb_ary_entry(LSts, leg) );
+      USts_array[leg] = NUM2DBL( rb_ary_entry(USts, leg) );
+    }
+  }
 
   //
   // first create a 1-D array of pointers, and then, for each array entry, create another 1-D array.
@@ -166,10 +192,24 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
         // profit if spot is higher than market price
         profit_loss_dash = ( Spot_dash - ( *( Xs_array + leg ) ) );
       } else if ( BS == SELL ) {
-        // profit if spot is lower than market price
-        profit_loss = ( ( *( Xs_array + leg ) ) - Spot );
-        // profit if spot is lower than market price
-        profit_loss_dash = ( ( *( Xs_array + leg ) ) - Spot_dash );
+        if ( KType == PIVOT_POINTS ) {
+          if ( *( Xs_array + leg ) < Spot && Spot < *( Ps_array + leg ) ) {
+            profit_loss = 0.0;
+          } else {
+            profit_loss = ( ( *( Xs_array + leg ) ) - Spot );
+          }
+          // dash
+          if ( *( Xs_array + leg ) < Spot_dash && Spot_dash < *( Ps_array + leg ) ) {
+            profit_loss_dash = 0.0;
+          } else {
+            profit_loss_dash = ( ( *( Xs_array + leg ) ) - Spot_dash );
+          }
+        } else {
+          // profit if spot is lower than market price
+          profit_loss = ( ( *( Xs_array + leg ) ) - Spot );
+          // profit if spot is lower than market price
+          profit_loss_dash = ( ( *( Xs_array + leg ) ) - Spot_dash );
+        }
       }
 
       sim[leg] = profit_loss;
@@ -179,7 +219,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       // Store spot and spot dash
       //
       if( point_pos <= SIM_LIMIT && (sim_count + 2) % INTERVAL == 0 ) {
-        rb_p(rb_str_new2("Leg:Spots"));
+        // rb_p(rb_str_new2("Leg:Spots"));
         metrics[ point_pos ][ leg ] = Spot;
         metrics[ point_pos + 2 ][ leg ] = Spot_dash;
       }
@@ -192,7 +232,6 @@ 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 == ABSOLUTE ) {
@@ -205,7 +244,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
             double temp_payoff_ko_ccy = temp_payoff * pow( *( Xs_array + leg ), Ko_compare_mult );
 
             if( temp_payoff_ko_ccy > ko_so_far ) {
-              sim_pos[leg] = ko_so_far * pow( *( Xs_array + leg ), Ko_compare_mult );
+              sim_pos[leg] = ko_so_far * pow( *( Xs_array + leg ), (-1 * Ko_compare_mult) );
               ko_so_far = 0.0;
             } else if( temp_payoff_ko_ccy <= ko_so_far ) {
               sim_pos[leg] = temp_payoff; // take the payoff in +ve's
@@ -224,7 +263,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
             double temp_payoff_ko_ccy = temp_payoff_dash * pow( *( Xs_array + leg ), Ko_compare_mult );
 
             if( temp_payoff_ko_ccy > ko_so_far_dash ) {
-              sim_dash_pos[leg] = ko_so_far_dash * pow( *( Xs_array + leg ), Ko_compare_mult );
+              sim_dash_pos[leg] = ko_so_far_dash * pow( *( Xs_array + leg ), (-1 * Ko_compare_mult) );
               ko_so_far_dash = 0.0;
             } else if( temp_payoff_ko_ccy <= ko_so_far_dash ) {
               sim_dash_pos[leg] = temp_payoff_dash; // take the payoff in +ve's
@@ -235,7 +274,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
           }
         }
       }
-    } else if( KType == POINTS ) {
+    } else if( KType == POINTS || KType == PIVOT_POINTS || KType == COLLAR_POINTS || KType == DUAL_STRIKE_POINTS ) {
       for( leg = 0; leg < NL; ++leg ) {
         // simulation normal
         if ( ko_so_far > 0.0 ) {
@@ -302,7 +341,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       // store and send whichever payoff you want
       //
       if( point_pos <= SIM_LIMIT && (sim_count + 2) % INTERVAL == 0 ) {
-        rb_p(rb_str_new2("Leg:Payoffs"));
+        // rb_p(rb_str_new2("Leg:Payoffs"));
         metrics[ point_pos + 1 ][ leg ] = (sim_pos[leg] + sim_neg[leg]);
         metrics[ point_pos + 3 ][ leg ] = (sim_dash_pos[leg] + sim_dash_neg[leg]);
       }
@@ -324,18 +363,18 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     // Note: no limit condition here
     //
     if( (sim_count + 2) % INTERVAL == 0 ) {
-      rb_p(rb_str_new2("T:"));
-      rb_p(INT2NUM(point_pos));
+      // rb_p(rb_str_new2("T:"));
+      // rb_p(INT2NUM(point_pos));
       point_pos += 4;
     }
   }
   // run simulations loop
-  rb_p(rb_str_new2("Simulation loop done"));
+  // rb_p(rb_str_new2("Simulation loop done"));
 
   //
   // sum all final +ve and -ve payoffs
   //
-  rb_p(rb_str_new2("Calculating final PVs"));
+  // rb_p(rb_str_new2("Calculating final PVs"));
   for(sim_count = 0; sim_count < SCount; ++sim_count) {
     pvs_pos_sum += *(pvs_pos + sim_count);
     pvs_neg_sum += *(pvs_neg + sim_count);
@@ -344,7 +383,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   //
   // free dynamically alloted heap memory
   //
-  rb_p(rb_str_new2("Free allocated arrays"));
+  // rb_p(rb_str_new2("Free allocated arrays"));
   free(pvs_pos);
   free(pvs_neg);
   free(Ls_array);
@@ -358,7 +397,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   //
   // return both payoffs
   //
-  rb_p(rb_str_new2("Converting metrics"));
+  // rb_p(rb_str_new2("Converting metrics"));
   VALUE final_metrics = rb_ary_new();
   for(metric = 0; metric < DATAPOINTS; metric++) {
     VALUE leg_metrics = rb_ary_new();
@@ -370,7 +409,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     rb_ary_push(final_metrics, leg_metrics);
   }
 
-  rb_p(rb_str_new2("Generating output hash"));
+  // rb_p(rb_str_new2("Generating output hash"));
   VALUE final_pvs = rb_hash_new();
   rb_hash_aset(final_pvs, rb_str_new2("positive_pv"), DBL2NUM( pvs_pos_sum / SCount ));
   rb_hash_aset(final_pvs, rb_str_new2("negative_pv"), DBL2NUM( pvs_neg_sum / SCount ));
@@ -402,4 +441,3 @@ void Init_tarf_monte_carlo() {
   rb_define_singleton_method(cTarfMonteCarlo, "box_muller", method_box_muller, 0);
   rb_define_singleton_method(cTarfMonteCarlo, "run_monte_carlo", method_run_monte_carlo, -2);
 }
-

data/lib/tarf_monte_carlo/version.rb

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

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.5'
+  version: '3.13'
 platform: ruby
 authors:
 - Vivek Routh
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2020-04-05 00:00:00.000000000 Z
+date: 2020-09-15 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -116,6 +116,7 @@ files:
 - bin/setup
 - ext/tarf_monte_carlo/extconf.rb
 - ext/tarf_monte_carlo/tarf_monte_carlo.c
+- lib/.DS_Store
 - lib/tarf_monte_carlo.rb
 - lib/tarf_monte_carlo/version.rb
 - tarf_monte_carlo.gemspec