tarf_monte_carlo 3.12 → 3.15.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 92b03e70d933ae7e53c9b8816bae025f315a455413684e542d03fb66b391fdb5
-  data.tar.gz: d18d1ac7aa941de3ebf8e5557635e8301a11002b468055a94c303d569edf9fdf
+  metadata.gz: bc23d4cbeb2b027663b1e8f6a2a3bcb559f4c63874bb2751358ce2e11e76412c
+  data.tar.gz: 102bee43e62b5bff5e30e478b5bf030d5b5530c1bea05f49376bc2511d6c5716
 SHA512:
-  metadata.gz: d54912e68221632f24689f1c69c88b8ca44daed63589025adb69910c4952b3819cf26d1ba9bf18e48a3db33bb017c39f7733494144a336633fd8a1da169557d1
-  data.tar.gz: 2a729f8b99b24bd1d7b7e9dea1d940d08bf0867efc22cb1cc18629f232bec75dd252d8360f82f11f5443af1026845293a13426193e7535af3f2c41df4de991d6
+  metadata.gz: c635e944069ef99ea9f7b3cb1d5d18c07b0427288c91bd7d94a221dc1ec47eaed356aada77b8e51a73a3ef18ee20b9573717269702426f73f6efdb58e8e85fee
+  data.tar.gz: ad0c9836a15abc620d10070c08f5cdce970543a58aad47b311151915d39d3f540e7968b02ad63c85eeace9d5bcce17af8b6a1bf7d48d81dd1d927b1ba8476eeb

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -124,11 +124,19 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     vs_array[leg] = NUM2DBL( rb_ary_entry(vs, leg) );
     bs_array[leg] = NUM2DBL( rb_ary_entry(bs, leg) );
     DFs_array[leg] = NUM2DBL( rb_ary_entry(DFs, leg) );
-    if (KType == PIVOT_POINTS) {
+  }
+  // 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){
+    }
+  } 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){
+    }
+  } 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) );
     }
@@ -179,27 +187,90 @@ 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 = Spot - (*( LSts_array + leg ));
+        } else {
+          profit_loss = (*( USts_array + leg )) - Spot;
+        }
+        // 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 = Spot_dash - (*( LSts_array + leg ));
+        } else {
+          profit_loss_dash = (*( USts_array + leg )) - Spot_dash;
+        }
+      } 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 ) {
         if ( KType == PIVOT_POINTS ) {
-          if ( Spot <= *( Xs_array + leg ) && Spot < *( Ps_array + leg ) ) {
-            profit_loss = *( Xs_array + leg ) - Spot;
-          } else if ( Spot > *( Xs_array + leg ) && Spot >= *( Ps_array + leg ) ) {
-            profit_loss = *( Xs_array + leg ) - Spot;
-          } else {
+          if ( *( Xs_array + leg ) < Spot && Spot < *( Ps_array + leg ) ) {
             profit_loss = 0.0;
+          } else {
+            profit_loss = ( ( *( Xs_array + leg ) ) - Spot );
           }
           // dash
-          if ( Spot_dash <= *( Xs_array + leg ) && Spot_dash < *( Ps_array + leg ) ) {
-            profit_loss_dash = *( Xs_array + leg ) - Spot_dash;
-          } else if ( Spot_dash > *( Xs_array + leg ) && Spot_dash >= *( Ps_array + leg ) ) {
-            profit_loss_dash = *( 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
@@ -229,7 +300,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 ) {
@@ -272,7 +342,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 ) {

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.2.gem
 # gem yank tarf_monte_carlo -v 2.3
 
 module TarfMonteCarlo
-  VERSION = "3.12"
+  VERSION = "3.15.2"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.12'
+  version: 3.15.2
 platform: ruby
 authors:
 - Vivek Routh
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2020-09-15 00:00:00.000000000 Z
+date: 2020-09-16 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler