tarf_monte_carlo 3.38 → 3.43

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: f6cb2b1cffe291ac211ef3517133affc4d62c2368f815a08a502a56650aee788
-  data.tar.gz: b17cfb2cce1123ab4c363d08a2eb214d14530d11249ed757c52e938bebacb38f
+  metadata.gz: 88cc80746dba47969e5cdc166b6510be889fd4c090ec1415cfdc066f47e32ebc
+  data.tar.gz: 4c7ab741260ab16d84a21659176211a466585d61be1c5adbdda62503173afa95
 SHA512:
-  metadata.gz: fb8b79a26f30802f2016981c4a96d68c193c794a76bd5b7b2f4f24d8805d0afc72150440b17b182898508a2cbc25675c7bba3bbc67288edec0aa8b8270c78bbc
-  data.tar.gz: bbe9daf65b7c7197fe76a4f2542ea422340399e2d8c46e4169d303f1da8475f2eef33e19b2e2f706fd3c3ef79de588fc3ebce5228c95be0afcc8bd6610522282
+  metadata.gz: 421f79a3b144d3066b04576f79b2041b24d3643322c76c4a06cca0be951069cff5e8d999cc3eaa6e2dcf62b81a321c7c8c5018aeefb1b3269d86d751d39251c7
+  data.tar.gz: 35b4d0df1af4fa5ca3ab88617c9516543a1c6127bf49287f21f19d00c1afa8bc13489fdb67bf808bcd8ac8d1348544a75fe3132f90f5f1d7979cdb7beae47a00

data/ext/tarf_monte_carlo/Makefile

@@ -135,9 +135,9 @@ extout_prefix =
 target_prefix = /tarf_monte_carlo
 LOCAL_LIBS = 
 LIBS = $(LIBRUBYARG_SHARED)  -lpthread -ldl -lobjc 
-ORIG_SRCS = tarf_monte_carlo (1).c tarf_monte_carlo.c
+ORIG_SRCS = tarf_monte_carlo.c
 SRCS = $(ORIG_SRCS) 
-OBJS = tarf_monte_carlo (1).o tarf_monte_carlo.o
+OBJS = tarf_monte_carlo.o
 HDRS = 
 LOCAL_HDRS = 
 TARGET = tarf_monte_carlo

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -194,13 +194,23 @@ 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;
-  double data[SCount][5];
-  metrics = ( double** ) malloc( DATAPOINTS * sizeof(double*) );
-  for( metric = 0; metric < DATAPOINTS; metric++ ) {
-    metrics[metric] = ( double* ) malloc( NL * sizeof(double) );
+  if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
+    metrics = ( double** ) malloc( 3 * sizeof(double*) );
+    for( metric = 0; metric < 3; metric++ ) {
+      metrics[metric] = ( double* ) malloc( DATAPOINTS * sizeof(double) );
+      for(leg = 0; leg < DATAPOINTS; leg++) {
+        metrics[metric][leg] = 0.0;
+      }
+    }
+  }
+  else{  
+    metrics = ( double** ) malloc( DATAPOINTS * sizeof(double*) );
+    for( metric = 0; metric < DATAPOINTS; metric++ ) {
+      metrics[metric] = ( double* ) malloc( NL * sizeof(double) );
 
-    for(leg = 0; leg < NL; leg++) {
-      metrics[metric][leg] = 0.0;
+      for(leg = 0; leg < NL; leg++) {
+        metrics[metric][leg] = 0.0;
+      }
     }
   }
 
@@ -210,18 +220,6 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     double Spot = S, Spot_dash = S;
     double sim[NL], sim_pos[NL], sim_neg[NL], sim_dash[NL], sim_dash_pos[NL], sim_dash_neg[NL];
 
-    data[sim_count][0] = 0;
-    data[sim_count][1] = 0;
-    data[sim_count][2] = 0;
-    data[sim_count][3] = 0;
-    data[sim_count][4] = 0;
-
-    data[sim_count + 1][0] = 0;
-    data[sim_count + 1][1] = 0;
-    data[sim_count + 1][2] = 0;
-    data[sim_count + 1][3] = 0;
-    data[sim_count + 1][4] = 0;
-
     // reset all simulation leg specific results
     for ( leg = 0; leg < NL; ++leg ) {
       sim[leg] = 0.0;
@@ -386,10 +384,12 @@ 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"));
-        metrics[ point_pos ][ leg ] = Spot;
-        metrics[ point_pos + 2 ][ leg ] = Spot_dash;
+      if (!(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE)){
+        if( point_pos <= SIM_LIMIT && (sim_count + 2) % INTERVAL == 0 ) {
+          // rb_p(rb_str_new2("Leg:Spots"));
+          metrics[ point_pos ][ leg ] = Spot;
+          metrics[ point_pos + 2 ][ leg ] = Spot_dash;
+        }
       }
 
       // if excuted UNCHAINED method
@@ -417,33 +417,35 @@ 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 ((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){
         sim_pos[knockedLeg] = get_equivalent_rebate(rebate_conversion_sign, ( *( Rebate_array + knockedLeg ) ), sim[knockedLeg], dir_sign);
-        data[sim_count][1] = sim_pos[knockedLeg];
-        data[sim_count][3] = sim_pos[knockedLeg];
-        data[sim_count][5] = sim[knockedLeg];
       }
       else{
         double equivalent_notional = get_equivalent_notional(conversion_sign, ( *( Ns_array + 0 ) ), sim[NL-1]);
-        data[sim_count][0] = equivalent_notional;
         sim_pos[NL-1] = european_payoff(( *( Xs_array + 0 ) ), sim[NL-1], cp_sign, dir_sign, equivalent_notional);
-        data[sim_count][2] = sim_pos[NL-1];
-        data[sim_count][5] = sim[NL-1];
       }
       if(knockedLeg_dash >= 0 ){
         sim_dash_pos[knockedLeg_dash] = get_equivalent_rebate(rebate_conversion_sign, ( *( Rebate_array + knockedLeg_dash ) ), sim[knockedLeg_dash], dir_sign);
-        data[sim_count + 1][1] = sim_dash_pos[knockedLeg_dash];
-        data[sim_count + 1][3] = sim_dash_pos[knockedLeg_dash];
-        data[sim_count + 1][5] = sim[knockedLeg_dash];
       }
       else{
         double equivalent_notional = get_equivalent_notional(conversion_sign, ( *( Ns_array + 0 ) ), sim_dash[NL-1]);
-        data[sim_count + 1][0] = equivalent_notional;
         sim_dash_pos[NL-1] = european_payoff(( *( Xs_array + 0 ) ), sim_dash[NL-1], cp_sign, dir_sign, equivalent_notional);
-        data[sim_count + 1][2] = sim_dash_pos[NL-1];
-        data[sim_count + 1][5] = sim[NL-1];
       }
+
+      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 ) {
       for( leg = 0; leg < NL; ++leg ) {
         // simulation normal
@@ -551,10 +553,12 @@ 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"));
-        metrics[ point_pos + 1 ][ leg ] = (sim_pos[leg] + sim_neg[leg]);
-        metrics[ point_pos + 3 ][ leg ] = (sim_dash_pos[leg] + sim_dash_neg[leg]);
+      if (!(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE)){
+        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]);
+          metrics[ point_pos + 3 ][ leg ] = (sim_dash_pos[leg] + sim_dash_neg[leg]);
+        }
       }
 
       sim_pos_sum += sim_pos[leg] * ( *( DFs_array + leg ) );
@@ -568,15 +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 ((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
@@ -610,14 +621,28 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   //
   // rb_p(rb_str_new2("Converting metrics"));
   VALUE final_metrics = rb_ary_new();
-  for(metric = 0; metric < SCount; metric++) {
-    VALUE leg_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();
+
+      for(leg = 0; leg < DATAPOINTS; leg++) {
+        rb_ary_push( leg_metrics, DBL2NUM( metrics[metric][leg] ) );
+      }
 
-    for(leg = 0; leg < 5; leg++) {
-      rb_ary_push( leg_metrics, DBL2NUM( data[metric][leg] ) );
+      rb_ary_push(final_metrics, leg_metrics);
     }
+  }
+  else{
+    for(metric = 0; metric < DATAPOINTS; metric++) {
+      VALUE leg_metrics = rb_ary_new();
 
-    rb_ary_push(final_metrics, leg_metrics);
+      for(leg = 0; leg < NL; leg++) {
+        rb_ary_push( leg_metrics, DBL2NUM( metrics[metric][leg] ) );
+      }
+
+      rb_ary_push(final_metrics, leg_metrics);
+    }
   }
 
   // rb_p(rb_str_new2("Generating output hash"));
@@ -630,7 +655,8 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   //
   // free your arrays from memory once you're done using them
   //
-  for(metric = 0; metric < DATAPOINTS; metric++) {
+  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] );
   }
   free(metrics);

data/lib/tarf_monte_carlo/version.rb

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

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.38'
+  version: '3.43'
 platform: ruby
 authors:
 - Vivek Routh
@@ -140,7 +140,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.0.6
+rubygems_version: 3.0.8
 signing_key: 
 specification_version: 4
 summary: Monte Carlo Simulation.