tarf_monte_carlo 3.7 → 3.15

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: cefc9ef8eb8823816be8c5856f835261f38587509092f4c1f24a680dfe5ed074
-  data.tar.gz: 3f4dc845dfd85eb5a7cd942525ef74fbce147399454acf3ab1d68fee22495e3c
+  metadata.gz: 0f33686b3c3246f31356eb1f1cc2161bb1a23a87f305b4a34e370d6c59584180
+  data.tar.gz: 3576f3fdb01ab53522ee92926d7536994551cc97503a9497c894b904c06582aa
 SHA512:
-  metadata.gz: 259da15728d5b6717a09edf7eba18f28353c27e5c6a2273582a5c1cc07616e2163845b798a4062c0ca44a0dd4540d45a509420e9139059bc21adcdd99671d7d2
-  data.tar.gz: 3b9f7a303b3807c9726591c8c50547065c6d21f1186e0fb70902412845f0a84513b57f3023fedeb724a03092cf87eaf20323a5eaf589b17ffd28d07d8e6f6ed4
+  metadata.gz: 27375bd106d3365323189a6d00ab4398b2d6d94d32fa82ed6ea69dc4d1f050a78d333ce4a694f8fd0a722f5c380b4e76ff946269b6b8cd5f29a02243120cf2d9
+  data.tar.gz: 5465e552c843af156e62bccbd1dbe34779c1f9e06f6596b2f142a817c4465299d635a8a1e8e4be0576e64ceae3f4017613209984aed02a2c0e4a2071425b9356

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,22 @@ 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) {
+      Ps_array[leg] = NUM2DBL( rb_ary_entry(Ps, 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.
@@ -160,16 +187,98 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       Spot = Spot * exp( drift + vSqrdt * eps );
       Spot_dash = Spot_dash * exp( drift + vSqrdt * eps_dash );
 
-      if ( BS == BUY ) {
-        // profit if spot is higher than market price
-        profit_loss = ( Spot - ( *( Xs_array + leg ) ) );
+      // Payoff of Dual strike points is independent of direction
+      if ( KType == DUAL_STRIKE_POINTS ) {
+        if ( Spot < *( LSts_array + leg ) ) {
+          profit_loss = Spot - (*( LSts_array + leg ));
+        } else if(Spot > *( USts_array + leg )){
+          profit_loss = (*( USts_array + leg ) - Spot);
+        } else if(Spot >= *( LSts_array + leg ) && (Spot < *(Ps_array + leg)))  {
+          profit_loss = (*( LSts_array + leg )) - Spot;
+        } else {
+          profit_loss = Spot - (*( USts_array + leg ))
+        }
+        // dash
+        if ( Spot_dash < *( LSts_array + leg ) ) {
+          profit_loss_dash = Spot_dash - (*( LSts_array + leg ));
+        } else if(Spot_dash > *( USts_array + leg )){
+          profit_loss_dash = (*( USts_array + leg ) - Spot_dash);
+        } else if(Spot_dash >= *( LSts_array + leg ) && (Spot_dash < *(Ps_array + leg)))  {
+          profit_loss_dash = (*( LSts_array + leg )) - Spot_dash;
+        } else {
+          profit_loss_dash = Spot_dash - (*( USts_array + leg ))
+        }
+      }
+      else if ( BS == BUY ) {
+        if ( KType == PIVOT_POINTS ) {
+          if ( *( Xs_array + leg ) > Spot && Spot > *( Ps_array + leg ) ) {
+            profit_loss = 0.0;
+          } else {
+            profit_loss = Spot - (*( Xs_array + leg ));
+          }
+          // dash
+          if ( *( Xs_array + leg ) > Spot_dash && Spot_dash > *( Ps_array + leg ) ) {
+            profit_loss_dash = 0.0;
+          } else {
+            profit_loss_dash = Spot_dash - (*( Xs_array + leg ));
+          }
+        } else if ( KType == COLLAR_POINTS ) {
+          if ( *( Xs_array + leg ) > Spot && Spot > *( CFs_array + leg ) ) {
+            profit_loss = 0.0;
+          } else if(Spot <= *( CFs_array + leg )){
+            profit_loss = Spot - (*( CFs_array + leg ));
+          } else {
+            profit_loss = Spot - (*( Xs_array + leg ));
+          }
+          // dash
+          if ( *( Xs_array + leg ) > Spot_dash && Spot_dash > *( CFs_array + leg ) ) {
+            profit_loss_dash = 0.0;
+          } else if(Spot_dash <= *( CFs_array + leg )){
+            profit_loss_dash = Spot_dash - (*( CFs_array + leg ));
+          } else {
+            profit_loss_dash = Spot_dash - (*( Xs_array + leg ));
+          }
+        } else {
         // profit if spot is higher than market price
-        profit_loss_dash = ( Spot_dash - ( *( Xs_array + leg ) ) );
+          profit_loss = ( Spot - ( *( Xs_array + leg ) ) );
+          // 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 if ( KType == COLLAR_POINTS ) {
+          if ( *( Xs_array + leg ) < Spot && Spot < *( CFs_array + leg ) ) {
+            profit_loss = 0.0;
+          } else if(Spot >= *( CFs_array + leg )){
+            profit_loss = ( ( *( CFs_array + leg ) ) - Spot );
+          } else {
+            profit_loss = ( ( *( Xs_array + leg ) ) - Spot );
+          }
+          // dash
+          if ( *( Xs_array + leg ) < Spot_dash && Spot_dash < *( CFs_array + leg ) ) {
+            profit_loss_dash = 0.0;
+          } else if(Spot_dash >= *( CFs_array + leg )){
+            profit_loss_dash = ( ( *( CFs_array + leg ) ) - Spot_dash );
+          } 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 +288,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 +301,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 ) {
@@ -235,7 +343,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 +410,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 +432,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 +452,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 +466,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 +478,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 ));

data/lib/tarf_monte_carlo.rb

@@ -4,11 +4,11 @@ require "tarf_monte_carlo/version"
 # my_malloc/my_malloc will be the shared object built when the gem is
 # installed then copied into lib/my_malloc/.
 
-# begin
-#   require "tarf_monte_carlo/#{RUBY_VERSION[/\d+.\d+/]}/tarf_monte_carlo"
-# rescue LoadError
-#   require "tarf_monte_carlo/tarf_monte_carlo"
-# end
+begin
+  require "tarf_monte_carlo/#{RUBY_VERSION[/\d+.\d+/]}/tarf_monte_carlo"
+rescue LoadError
+  require "tarf_monte_carlo/tarf_monte_carlo"
+end
 
 module TarfMonteCarlo
   # Your code goes here...

data/lib/tarf_monte_carlo/version.rb

@@ -1,7 +1,7 @@
 # gem build tarf_monte_carlo.gemspec
-# gem push -v tarf_monte_carlo-2.3.gem
+# gem push -v tarf_monte_carlo-3.15.gem
 # gem yank tarf_monte_carlo -v 2.3
 
 module TarfMonteCarlo
-  VERSION = "3.7"
+  VERSION = "3.15"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.7'
+  version: '3.15'
 platform: ruby
 authors:
 - Vivek Routh
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2020-04-06 00:00:00.000000000 Z
+date: 2020-09-16 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