tarf_monte_carlo 3.6 → 3.14

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: bb1130c50f083ee54950ca2a4d8325d462762da404fd52c8cd0a67534f1ea410
-  data.tar.gz: 5cfb4c9642286a14869655c2ba082add5226d076be8998c78c5df2eccbbd702f
+  metadata.gz: d8d2906e29ad90b4dc0ff480eec64fed2a8d7aaca0e5933de6fa1d672185a38b
+  data.tar.gz: d8afc4b56ea885bcd67299273ab887b2c97b6e73e6d218729a9e39b54433a240
 SHA512:
-  metadata.gz: a1bc8757d3a7945059c5784f6180e596576a828a8283bc03ffa85243425f2300136d10cf5a69e4480c79ede4be88010c2f73171f8750182bbf6da7d2b0b3716b
-  data.tar.gz: d7805cc57138165266e6e63b1ec062fad176ec598610c5317d21e9d756e097c060fe8b9101403df5ca91240619b316f46db2f47f433cb1aa80e4658985e0d8c7
+  metadata.gz: ae61b13786e7deb50903881364d793c96db0d453f34f54dd6dd4040a4eb3f86cfaeb6217be6433cbc8f9827c7ba12caabcfc2953edf3a18088a17eceda03fafc
+  data.tar.gz: 546c462c9ccb5038e9ebb385aec5ccc1f66745d011d72f55c68afa9ec0a83847b0d8beea21bc531610453bd473e82e3a6792c0dc9138125247c32afd2228539d

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 ) {
@@ -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 ));

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

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

metadata

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