tarf_monte_carlo 3.41 → 3.42

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 613891fdd7b41ca68d9e30b1b3106545dbfdcaab578ea5be585b8e8c33875f7a
-  data.tar.gz: 62de3fcb90c0edfe1298938ed9e351a54a39e42366de3da5a1a01cc230c65712
+  metadata.gz: 401dcffa11be28cb765633d86b3b3351af887311364893a4b4e98bcda42f1f2b
+  data.tar.gz: 8aa5d4b647a4ead25e93d8ec6d5d77f6927d767b368915c367a1557cfcd0d578
 SHA512:
-  metadata.gz: 4da0be66316bb60b507176f283a8d87bc02c83e3b9218fd198812c0847c94410f79ff5c87c4ceae83084546747c711264eca1c98dc104af7511e704a221c11e2
-  data.tar.gz: b2ee2e95b2f23b628e69b6e6b4e334c9453f6325e4374d43dffbaf4cff51217899f56b48f41a7ce2bd3c0abba4c819689233aa4df0fe343d5a644d26d58b5f53
+  metadata.gz: a955c81445337676aa5c291fed82e34fa74b37a999926158b9cd09bb5c6ae14888f483c43dcafc33b04d43f4fc3c016374f81c19ac39aa71933c6746196a0978
+  data.tar.gz: 73b9c4e40243dc65b8f711d515a1ab6196fa4c17991f12d72543cb8918b990e83844f4679e55fb28ebf3bac6857ff0e9cbdd2cf5ef55f4fdeb9dc4fa1b9135f9

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -417,11 +417,11 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       else{
         sim_dash_pos[NL-1] = get_equivalent_rebate(rebate_conversion_sign, ( *( Rebate_array + 0 ) ), sim_dash[NL-1], dir_sign);
       }
-      if ((sim_count + 2) % INTERVAL == 0){
-        metrics[0][sim_count] = sim[NL-1];
-        metrics[1][sim_count] = (knockedLeg >= 0) ? ( *( Xs_array + (knockedLeg) ) ) : 0;
-        metrics[0][sim_count+1] = sim_dash[NL-1];
-        metrics[1][sim_count+1] = (knockedLeg_dash >= 0) ? ( *( Xs_array + (knockedLeg_dash) ) ) : 0;
+      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+1] = sim_dash[NL-1];
+        metrics[1][point_pos+1] = (knockedLeg_dash >= 0) ? ( *( Xs_array + (knockedLeg_dash) ) ) : 0;
       }
     } else if(KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
       if(knockedLeg >= 0){
@@ -439,11 +439,11 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
         sim_dash_pos[NL-1] = european_payoff(( *( Xs_array + 0 ) ), sim_dash[NL-1], cp_sign, dir_sign, equivalent_notional);
       }
 
-      if ((sim_count + 2) % INTERVAL == 0){
-        metrics[0][sim_count] = (knockedLeg >= 0) ? sim[knockedLeg] : sim[NL-1];
-        metrics[1][sim_count] = ( *( Xs_array + 0 ) );
-        metrics[0][sim_count+1] = (knockedLeg_dash >= 0) ? sim[knockedLeg_dash] : sim[NL-1];
-        metrics[1][sim_count+1] = ( *( Xs_array + 0 ) );
+      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+1] = (knockedLeg_dash >= 0) ? sim[knockedLeg_dash] : sim[NL-1];
+        metrics[1][point_pos+1] = ( *( Xs_array + 0 ) );
       }
 
     } else if( KType == ABSOLUTE || KType == PIVOT_ABSOLUTE || KType == COLLAR_ABSOLUTE || KType == DOUBLE_STRIKE_ABSOLUTE ) {
@@ -572,19 +572,22 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     pvs_pos[sim_count + 1]  = sim_dash_pos_sum;
     pvs_neg[sim_count]      = sim_neg_sum;
     pvs_neg[sim_count + 1]  = sim_dash_neg_sum;
-    if ((sim_count + 2) % INTERVAL == 0){
-      metrics[2][sim_count]   = sim_pos_sum;
-      metrics[2][sim_count+1] = sim_dash_pos_sum;
+    if ((point_pos < DATAPOINTS) && (sim_count % INTERVAL) == 0){
+      metrics[2][point_pos]   = sim_pos_sum;
+      metrics[2][point_pos+1] = sim_dash_pos_sum;
+      point_pos += 2;
     }
 
     //
     // increment metric storing point by 4
     // Note: no limit condition here
     //
-    if( (sim_count + 2) % INTERVAL == 0 ) {
-      // rb_p(rb_str_new2("T:"));
-      // rb_p(INT2NUM(point_pos));
-      point_pos += 4;
+    if (!(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE)){
+      if( (sim_count + 2) % INTERVAL == 0 ) {
+        // rb_p(rb_str_new2("T:"));
+        // rb_p(INT2NUM(point_pos));
+        point_pos += 4;
+      }
     }
   }
   // run simulations loop
@@ -618,7 +621,7 @@ 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){
     for(metric = 0; metric < 3; metric++) {
       VALUE leg_metrics = rb_ary_new();
@@ -652,8 +655,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) ? 3 : DATAPOINTS;
-
+  int dp = (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE) ? 3 : 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.41"
+  VERSION = "3.42"
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.41'
+  version: '3.42'
 platform: ruby
 authors:
 - Vivek Routh