tarf_monte_carlo 3.13 → 3.15.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: bf901445c70d80f370053ffa90ba2cc6fb3e12b9ad7c372777cae63b92ceecdb
-  data.tar.gz: 04331203b779e146db0e48a7eb9d1a0452c31aee17311c5f3c863d2104bc080b
+  metadata.gz: 9f59add36257db0337d6c1a72238af5d328892cb59959117797130e5b90b0f31
+  data.tar.gz: d78b75d035a82c487bdc15797e3589ff4c15f33f9180eea8f8d211b847c22fc6
 SHA512:
-  metadata.gz: b912836d2d1a64d62c9b1d394be9ca89c5dbeb245312ce18d9490751bf61fd1fc9ade361a60446d1a0b1df7928882e5e4d0691b0e8c7528733f9bdfc49ccf6b2
-  data.tar.gz: '08248d29944ff3fb8981d7bda9506bd056869575e54cb1ec6aa3be403b945cbe1b11c6bcb1cc42193f5617330d1be6633f4f1c2ab331cd2cc70d4cf455ef7df1'
+  metadata.gz: 675daf6c71ee357cc6eca37080acc93fa5211dc9737104c92ef0adc1c9d8bd8b61275d681fa64b1fd3da8832867b1f8e879eea8198d4b29dd06d1188eb645ba7
+  data.tar.gz: bbc2a4bfe53e589db3225c86bec1449adee93e2f9022f6b62775eca692b6a62a832988e627be1a0edf35478aa4a876a260735d1d0288eb3fac03023eaf3ae2c7

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -28,9 +28,19 @@
 #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 DOUBLE_STRIKE_POINTS 6
+
+#define PIVOT_ABSOLUTE 7
+#define COLLAR_ABSOLUTE 8
+#define DOUBLE_STRIKE_ABSOLUTE 9
+
+#define PIVOT_LEGS 10
+#define COLLAR_LEGS 11
+#define DOUBLE_STRIKE_LEGS 12
+
 #define DATAPOINTS 200 // data for plotting
 #define INTERVAL 50
 #define SIM_LIMIT 196 // 196 + 4 = 200 simulations nedded
@@ -126,16 +136,17 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     DFs_array[leg] = NUM2DBL( rb_ary_entry(DFs, leg) );
   }
   // extra tarf structures
-  if ( KType == PIVOT_POINTS )  {
+  if ( KType == PIVOT_POINTS || KType == PIVOT_ABSOLUTE || KType == PIVOT_LEGS )  {
     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 || KType == COLLAR_ABSOLUTE || KType == COLLAR_LEGS ) {
     for (leg = 0; leg < NL; ++leg) {
       CFs_array[leg] = NUM2DBL( rb_ary_entry(CFs, leg) );
     }
-  } else if ( KType == DUAL_STRIKE_POINTS ) {
+  } else if ( KType == DOUBLE_STRIKE_POINTS || KType == DOUBLE_STRIKE_ABSOLUTE || KType == DOUBLE_STRIKE_LEGS) {
     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) );
     }
@@ -186,13 +197,64 @@ 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 == DOUBLE_STRIKE_POINTS || KType == DOUBLE_STRIKE_ABSOLUTE || KType == DOUBLE_STRIKE_LEGS ) {
+        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 || KType == PIVOT_ABSOLUTE || KType == PIVOT_LEGS ) {
+          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 || KType == COLLAR_ABSOLUTE || KType == COLLAR_LEGS ) {
+          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 ( KType == PIVOT_POINTS || KType == PIVOT_ABSOLUTE || KType == PIVOT_LEGS ) {
           if ( *( Xs_array + leg ) < Spot && Spot < *( Ps_array + leg ) ) {
             profit_loss = 0.0;
           } else {
@@ -204,6 +266,22 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
           } else {
             profit_loss_dash = ( ( *( Xs_array + leg ) ) - Spot_dash );
           }
+        } else if ( KType == COLLAR_POINTS || KType == COLLAR_ABSOLUTE || KType == COLLAR_LEGS ) {
+          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 );
@@ -234,7 +312,7 @@ 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 ) {
+    if( KType == ABSOLUTE || KType == PIVOT_ABSOLUTE || KType == COLLAR_ABSOLUTE || KType == DOUBLE_STRIKE_ABSOLUTE ) {
       for( leg = 0; leg < NL; ++leg ) {
         // simulation normal
         if( ko_so_far > 0.0 ) {
@@ -274,7 +352,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
           }
         }
       }
-    } else if( KType == POINTS || KType == PIVOT_POINTS || KType == COLLAR_POINTS || KType == DUAL_STRIKE_POINTS ) {
+    } else if( KType == POINTS || KType == PIVOT_POINTS || KType == COLLAR_POINTS || KType == DOUBLE_STRIKE_POINTS ) {
       for( leg = 0; leg < NL; ++leg ) {
         // simulation normal
         if ( ko_so_far > 0.0 ) {
@@ -308,7 +386,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
           }
         }
       }
-    } else if( KType == LEGS ) {
+    } else if( KType == LEGS || KType == PIVOT_LEGS || KType == COLLAR_LEGS || KType == DOUBLE_STRIKE_LEGS ) {
       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.3.gem
 # gem yank tarf_monte_carlo -v 2.3
 
 module TarfMonteCarlo
-  VERSION = "3.13"
+  VERSION = "3.15.3"
 end

metadata

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