tarf_monte_carlo 3.23 → 3.28

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: e31ee8faa0b4ae02a93252953f455286d901887af6983874cb32c388842d9038
-  data.tar.gz: 7ea19ee7948374d024720267a5cb4aa0d532d1acfb0ed7e56fa4401b006af2d7
+  metadata.gz: 61612dbfd673a328672091d8e64ec0ebd5251f0e1d0904938e4fcfa1e083ff84
+  data.tar.gz: 22c0461dc088058378702be27ad055f600f97eac6c5392008d9c5783ff15cea0
 SHA512:
-  metadata.gz: 60943ac5835b69779f2c3fc89536113b3c99a85413e4e51c6bbbc60d2d3c8266cf1d2a27798c89ae6f8e759ad8eace8825233945173f6c3c57ba06a53fdcbca8
-  data.tar.gz: c9ff7608e59e0d94b2c2d4d585e5b2a0e28d9064d4e28c547634d2e72da3d7196b912cf89820c594f7fffcf6c3b2f0e910788f1bd41fda1228092ae1c70f332f
+  metadata.gz: 0f3224cfb1535727083cb0223bf06eca402e43d50f4bd6eb3adccc7cf44cd47f5ddf16cf8882b57ca7ae564014f805cb8a3c5e3137edd967862e5b5653d9e0dd
+  data.tar.gz: 77e3ed4c219a3a797f1f03b8fa663eea5bb954672824a5154610fe1012fe25c1c0a9981881b272e678ae1579ea32ced99566539bb5f016255d1b06f8815d1a2e

data/ext/tarf_monte_carlo/Makefile

@@ -0,0 +1,264 @@
+
+SHELL = /bin/sh
+
+# V=0 quiet, V=1 verbose.  other values don't work.
+V = 0
+Q1 = $(V:1=)
+Q = $(Q1:0=@)
+ECHO1 = $(V:1=@ :)
+ECHO = $(ECHO1:0=@ echo)
+NULLCMD = :
+
+#### Start of system configuration section. ####
+
+srcdir = .
+topdir = /Users/vivektripathi/.rvm/rubies/ruby-2.5.1/include/ruby-2.5.0
+hdrdir = $(topdir)
+arch_hdrdir = /Users/vivektripathi/.rvm/rubies/ruby-2.5.1/include/ruby-2.5.0/x86_64-darwin18
+PATH_SEPARATOR = :
+VPATH = $(srcdir):$(arch_hdrdir)/ruby:$(hdrdir)/ruby
+prefix = $(DESTDIR)/Users/vivektripathi/.rvm/rubies/ruby-2.5.1
+rubysitearchprefix = $(rubylibprefix)/$(sitearch)
+rubyarchprefix = $(rubylibprefix)/$(arch)
+rubylibprefix = $(libdir)/$(RUBY_BASE_NAME)
+exec_prefix = $(prefix)
+vendorarchhdrdir = $(vendorhdrdir)/$(sitearch)
+sitearchhdrdir = $(sitehdrdir)/$(sitearch)
+rubyarchhdrdir = $(rubyhdrdir)/$(arch)
+vendorhdrdir = $(rubyhdrdir)/vendor_ruby
+sitehdrdir = $(rubyhdrdir)/site_ruby
+rubyhdrdir = $(includedir)/$(RUBY_VERSION_NAME)
+vendorarchdir = $(vendorlibdir)/$(sitearch)
+vendorlibdir = $(vendordir)/$(ruby_version)
+vendordir = $(rubylibprefix)/vendor_ruby
+sitearchdir = $(sitelibdir)/$(sitearch)
+sitelibdir = $(sitedir)/$(ruby_version)
+sitedir = $(rubylibprefix)/site_ruby
+rubyarchdir = $(rubylibdir)/$(arch)
+rubylibdir = $(rubylibprefix)/$(ruby_version)
+sitearchincludedir = $(includedir)/$(sitearch)
+archincludedir = $(includedir)/$(arch)
+sitearchlibdir = $(libdir)/$(sitearch)
+archlibdir = $(libdir)/$(arch)
+ridir = $(datarootdir)/$(RI_BASE_NAME)
+mandir = $(datarootdir)/man
+localedir = $(datarootdir)/locale
+libdir = $(exec_prefix)/lib
+psdir = $(docdir)
+pdfdir = $(docdir)
+dvidir = $(docdir)
+htmldir = $(docdir)
+infodir = $(datarootdir)/info
+docdir = $(datarootdir)/doc/$(PACKAGE)
+oldincludedir = $(SDKROOT)/usr/include
+includedir = $(prefix)/include
+localstatedir = $(prefix)/var
+sharedstatedir = $(prefix)/com
+sysconfdir = $(prefix)/etc
+datadir = $(datarootdir)
+datarootdir = $(prefix)/share
+libexecdir = $(exec_prefix)/libexec
+sbindir = $(exec_prefix)/sbin
+bindir = $(exec_prefix)/bin
+archdir = $(rubyarchdir)
+
+
+CC = gcc
+CXX = g++
+LIBRUBY = $(LIBRUBY_SO)
+LIBRUBY_A = lib$(RUBY_SO_NAME)-static.a
+LIBRUBYARG_SHARED = -l$(RUBY_SO_NAME)
+LIBRUBYARG_STATIC = -l$(RUBY_SO_NAME)-static -framework Foundation
+empty =
+OUTFLAG = -o $(empty)
+COUTFLAG = -o $(empty)
+CSRCFLAG = $(empty)
+
+RUBY_EXTCONF_H = 
+cflags   = $(optflags) $(debugflags) $(warnflags)
+cxxflags = $(optflags) $(debugflags) $(warnflags)
+optflags = -O3
+debugflags = -ggdb3
+warnflags = -Wall -Wextra -Wno-unused-parameter -Wno-parentheses -Wno-long-long -Wno-missing-field-initializers -Wno-tautological-compare -Wno-parentheses-equality -Wno-constant-logical-operand -Wno-self-assign -Wunused-variable -Wimplicit-int -Wpointer-arith -Wwrite-strings -Wdeclaration-after-statement -Wshorten-64-to-32 -Wimplicit-function-declaration -Wdivision-by-zero -Wdeprecated-declarations -Wextra-tokens
+CCDLFLAGS = -fno-common
+CFLAGS   = $(CCDLFLAGS) $(cflags)  -fno-common -pipe $(ARCH_FLAG)
+INCFLAGS = -I. -I$(arch_hdrdir) -I$(hdrdir)/ruby/backward -I$(hdrdir) -I$(srcdir)
+DEFS     = 
+CPPFLAGS =   -I/usr/local/opt/libyaml/include -I/usr/local/opt/libksba/include -I/usr/local/opt/readline/include -I/usr/local/opt/zlib/include -I/usr/local/opt/openssl@1.1/include -D_XOPEN_SOURCE -D_DARWIN_C_SOURCE -D_DARWIN_UNLIMITED_SELECT -D_REENTRANT $(DEFS) $(cppflags)
+CXXFLAGS = $(CCDLFLAGS) $(cxxflags) $(ARCH_FLAG)
+ldflags  = -L. -fstack-protector -L/usr/local/lib -L/usr/local/opt/libyaml/lib -L/usr/local/opt/libksba/lib -L/usr/local/opt/readline/lib -L/usr/local/opt/zlib/lib -L/usr/local/opt/openssl@1.1/lib
+dldflags = -Wl,-undefined,dynamic_lookup -Wl,-multiply_defined,suppress -L/usr/local/opt/libyaml/lib -L/usr/local/opt/libksba/lib -L/usr/local/opt/readline/lib -L/usr/local/opt/zlib/lib -L/usr/local/opt/openssl@1.1/lib 
+ARCH_FLAG = 
+DLDFLAGS = $(ldflags) $(dldflags) $(ARCH_FLAG)
+LDSHARED = $(CC) -dynamic -bundle
+LDSHAREDXX = $(CXX) -dynamic -bundle
+AR = libtool -static
+EXEEXT = 
+
+RUBY_INSTALL_NAME = $(RUBY_BASE_NAME)
+RUBY_SO_NAME = ruby.2.5.1
+RUBYW_INSTALL_NAME = 
+RUBY_VERSION_NAME = $(RUBY_BASE_NAME)-$(ruby_version)
+RUBYW_BASE_NAME = rubyw
+RUBY_BASE_NAME = ruby
+
+arch = x86_64-darwin18
+sitearch = $(arch)
+ruby_version = 2.5.0
+ruby = $(bindir)/$(RUBY_BASE_NAME)
+RUBY = $(ruby)
+ruby_headers = $(hdrdir)/ruby.h $(hdrdir)/ruby/backward.h $(hdrdir)/ruby/ruby.h $(hdrdir)/ruby/defines.h $(hdrdir)/ruby/missing.h $(hdrdir)/ruby/intern.h $(hdrdir)/ruby/st.h $(hdrdir)/ruby/subst.h $(arch_hdrdir)/ruby/config.h
+
+RM = rm -f
+RM_RF = $(RUBY) -run -e rm -- -rf
+RMDIRS = rmdir -p
+MAKEDIRS = /usr/local/opt/coreutils/bin/gmkdir -p
+INSTALL = /usr/local/opt/coreutils/bin/ginstall -c
+INSTALL_PROG = $(INSTALL) -m 0755
+INSTALL_DATA = $(INSTALL) -m 644
+COPY = cp
+TOUCH = exit >
+
+#### End of system configuration section. ####
+
+preload = 
+libpath = . $(libdir) /usr/local/opt/libyaml/lib /usr/local/opt/libksba/lib /usr/local/opt/readline/lib /usr/local/opt/zlib/lib /usr/local/opt/openssl@1.1/lib
+LIBPATH =  -L. -L$(libdir) -L/usr/local/opt/libyaml/lib -L/usr/local/opt/libksba/lib -L/usr/local/opt/readline/lib -L/usr/local/opt/zlib/lib -L/usr/local/opt/openssl@1.1/lib
+DEFFILE = 
+
+CLEANFILES = mkmf.log
+DISTCLEANFILES = 
+DISTCLEANDIRS = 
+
+extout = 
+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
+SRCS = $(ORIG_SRCS) 
+OBJS = tarf_monte_carlo (1).o tarf_monte_carlo.o
+HDRS = 
+LOCAL_HDRS = 
+TARGET = tarf_monte_carlo
+TARGET_NAME = tarf_monte_carlo
+TARGET_ENTRY = Init_$(TARGET_NAME)
+DLLIB = $(TARGET).bundle
+EXTSTATIC = 
+STATIC_LIB = 
+
+TIMESTAMP_DIR = .
+BINDIR        = $(bindir)
+RUBYCOMMONDIR = $(sitedir)$(target_prefix)
+RUBYLIBDIR    = $(sitelibdir)$(target_prefix)
+RUBYARCHDIR   = $(sitearchdir)$(target_prefix)
+HDRDIR        = $(rubyhdrdir)/ruby$(target_prefix)
+ARCHHDRDIR    = $(rubyhdrdir)/$(arch)/ruby$(target_prefix)
+TARGET_SO_DIR =
+TARGET_SO     = $(TARGET_SO_DIR)$(DLLIB)
+CLEANLIBS     = $(TARGET_SO) 
+CLEANOBJS     = *.o  *.bak
+
+all:    $(DLLIB)
+static: $(STATIC_LIB)
+.PHONY: all install static install-so install-rb
+.PHONY: clean clean-so clean-static clean-rb
+
+clean-static::
+clean-rb-default::
+clean-rb::
+clean-so::
+clean: clean-so clean-static clean-rb-default clean-rb
+		-$(Q)$(RM) $(CLEANLIBS) $(CLEANOBJS) $(CLEANFILES) .*.time
+
+distclean-rb-default::
+distclean-rb::
+distclean-so::
+distclean-static::
+distclean: clean distclean-so distclean-static distclean-rb-default distclean-rb
+		-$(Q)$(RM) Makefile $(RUBY_EXTCONF_H) conftest.* mkmf.log
+		-$(Q)$(RM) core ruby$(EXEEXT) *~ $(DISTCLEANFILES)
+		-$(Q)$(RMDIRS) $(DISTCLEANDIRS) 2> /dev/null || true
+
+realclean: distclean
+install: install-so install-rb
+
+install-so: $(DLLIB) $(TIMESTAMP_DIR)/.sitearchdir.-.tarf_monte_carlo.time
+	$(INSTALL_PROG) $(DLLIB) $(RUBYARCHDIR)
+clean-static::
+	-$(Q)$(RM) $(STATIC_LIB)
+install-rb: pre-install-rb do-install-rb install-rb-default
+install-rb-default: pre-install-rb-default do-install-rb-default
+pre-install-rb: Makefile
+pre-install-rb-default: Makefile
+do-install-rb:
+do-install-rb-default:
+pre-install-rb-default:
+	@$(NULLCMD)
+$(TIMESTAMP_DIR)/.sitearchdir.-.tarf_monte_carlo.time:
+	$(Q) $(MAKEDIRS) $(@D) $(RUBYARCHDIR)
+	$(Q) $(TOUCH) $@
+
+site-install: site-install-so site-install-rb
+site-install-so: install-so
+site-install-rb: install-rb
+
+.SUFFIXES: .c .m .cc .mm .cxx .cpp .o .S
+
+.cc.o:
+	$(ECHO) compiling $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -c $(CSRCFLAG)$<
+
+.cc.S:
+	$(ECHO) translating $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -S $(CSRCFLAG)$<
+
+.mm.o:
+	$(ECHO) compiling $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -c $(CSRCFLAG)$<
+
+.mm.S:
+	$(ECHO) translating $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -S $(CSRCFLAG)$<
+
+.cxx.o:
+	$(ECHO) compiling $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -c $(CSRCFLAG)$<
+
+.cxx.S:
+	$(ECHO) translating $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -S $(CSRCFLAG)$<
+
+.cpp.o:
+	$(ECHO) compiling $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -c $(CSRCFLAG)$<
+
+.cpp.S:
+	$(ECHO) translating $(<)
+	$(Q) $(CXX) $(INCFLAGS) $(CPPFLAGS) $(CXXFLAGS) $(COUTFLAG)$@ -S $(CSRCFLAG)$<
+
+.c.o:
+	$(ECHO) compiling $(<)
+	$(Q) $(CC) $(INCFLAGS) $(CPPFLAGS) $(CFLAGS) $(COUTFLAG)$@ -c $(CSRCFLAG)$<
+
+.c.S:
+	$(ECHO) translating $(<)
+	$(Q) $(CC) $(INCFLAGS) $(CPPFLAGS) $(CFLAGS) $(COUTFLAG)$@ -S $(CSRCFLAG)$<
+
+.m.o:
+	$(ECHO) compiling $(<)
+	$(Q) $(CC) $(INCFLAGS) $(CPPFLAGS) $(CFLAGS) $(COUTFLAG)$@ -c $(CSRCFLAG)$<
+
+.m.S:
+	$(ECHO) translating $(<)
+	$(Q) $(CC) $(INCFLAGS) $(CPPFLAGS) $(CFLAGS) $(COUTFLAG)$@ -S $(CSRCFLAG)$<
+
+$(TARGET_SO): $(OBJS) Makefile
+	$(ECHO) linking shared-object tarf_monte_carlo/$(DLLIB)
+	-$(Q)$(RM) $(@)
+	$(Q) $(LDSHARED) -o $@ $(OBJS) $(LIBPATH) $(DLDFLAGS) $(LOCAL_LIBS) $(LIBS)
+	$(Q) $(POSTLINK)
+
+
+
+$(OBJS): $(HDRS) $(ruby_headers)

data/ext/tarf_monte_carlo/tarf_monte_carlo.c

@@ -41,6 +41,9 @@
 #define COLLAR_LEGS 11
 #define DOUBLE_STRIKE_LEGS 12
 
+#define FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE 13
+#define FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE 14
+
 #define DATAPOINTS 200 // data for plotting
 #define INTERVAL 50
 #define SIM_LIMIT 196 // 196 + 4 = 200 simulations nedded
@@ -53,6 +56,9 @@ VALUE cTarfMonteCarlo = Qnil;
 
 // Prototype for the initialization method - Ruby calls this, not you
 void Init_tarf_monte_carlo();
+double european_payoff(double, double, int, int, double);
+double get_equivalent_notional(int, double, float);
+double get_equivalent_rebate(int, double, float);
 
 // Prototype for our methods - methods are prefixed by 'method_' here
 VALUE method_box_muller( VALUE );
@@ -77,11 +83,25 @@ VALUE method_box_muller( VALUE self ) {
   return DBL2NUM(eps);
 }
 
+double european_payoff(double strike, double spot, int cp_sign, int dir_sign, double notional){
+  double diff = (spot - strike) * cp_sign;
+  diff = diff > 0 ? diff : 0;
+  return (diff * notional * dir_sign);
+}
+
+double get_equivalent_notional(int conversion_sign, double notional, float rate){
+  return (conversion_sign == 1) ? (notional / rate) : notional
+}
+
+double get_equivalent_rebate(int conversion_sign, double rebate, float rate){
+  return (conversion_sign == 1) ? (rebate * rate) : rebate
+}
+
 // 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(MCInputs);
 
   // seed value for rand() function
   srand( time(0) );
@@ -99,6 +119,9 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   double S              = NUM2DBL( rb_hash_aref(MCInputs, rb_str_new2("spot_rate")) );
   int Ko_compare_mult   = NUM2INT( rb_hash_aref(MCInputs, rb_str_new2("multiplier")) );
   int ConvertNotional   = NUM2INT( rb_hash_aref(MCInputs, rb_str_new2("convert_notional")) );
+  int dir_sign          = NUM2INT( rb_hash_aref(MCInputs, rb_str_new2("direction_sign")) );
+  int dir2_sign         = NUM2INT( rb_hash_aref(MCInputs, rb_str_new2("direction2_sign")) );
+  int cp_sign           = NUM2INT( rb_hash_aref(MCInputs, rb_str_new2("callput_sign")) );
 
   // assign leg specific attributes
   double *pvs_pos     = ( double* ) malloc( SCount * sizeof(double) );
@@ -114,6 +137,9 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   double *CFs_array   = ( double* ) malloc( NL * sizeof(double) );
   double *LSts_array  = ( double* ) malloc( NL * sizeof(double) );
   double *USts_array  = ( double* ) malloc( NL * sizeof(double) );
+  double *TempNs_array = ( double* ) malloc( NL * sizeof(double) );
+  double *Barrier_array = ( double* ) malloc( NL * sizeof(double) );
+  double *Rebate_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") );
@@ -126,6 +152,8 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
   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") );
+  VALUE Brs = rb_hash_aref(MCInputs, rb_str_new2("barriers") );
+  VALUE Rbts = rb_hash_aref(MCInputs, rb_str_new2("rebates") );
 
   for (leg = 0; leg < NL; ++leg) {
     Ls_array[leg] = NUM2DBL( rb_ary_entry(Ls, leg) );
@@ -135,6 +163,7 @@ 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) );
+    TempNs_array[leg] = *(Ns_array + leg);
   }
   // extra tarf structures
   if ( KType == PIVOT_POINTS || KType == PIVOT_ABSOLUTE || KType == PIVOT_LEGS )  {
@@ -151,6 +180,11 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       LSts_array[leg] = NUM2DBL( rb_ary_entry(LSts, leg) );
       USts_array[leg] = NUM2DBL( rb_ary_entry(USts, leg) );
     }
+  } else if (KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE) {
+    for (leg = 0; leg < NL; ++leg) {
+      Barrier_array[leg]  = NUM2DBL( rb_ary_entry(Brs, leg) );
+      Rebate_array[leg]   = NUM2DBL( rb_ary_entry(Rbts, leg) );
+    }
   }
 
   //
@@ -182,6 +216,8 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       sim_dash_neg[leg] = 0.0;
     }
 
+    int knockedLeg = -1;
+    int knockedLeg_dash = -1;
     // legs loop start
     for( leg = 0; leg < NL; ++leg ) {
       double eps, eps_dash, drift, vSqrdt, profit_loss, profit_loss_dash;
@@ -198,8 +234,20 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
       Spot = Spot * exp( drift + vSqrdt * eps );
       Spot_dash = Spot_dash * exp( drift + vSqrdt * eps_dash );
 
-      // Payoff of Dual strike points is independent of direction
-      if ( KType == DOUBLE_STRIKE_POINTS || KType == DOUBLE_STRIKE_ABSOLUTE || KType == DOUBLE_STRIKE_LEGS ) {
+      if(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE || KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
+        profit_loss = Spot;
+        profit_loss_dash = Spot_dash;
+
+        double knockin = (Spot - *(Barrier_array + leg)) * dir2_sign;
+        double knockin_dash = (Spot_dash - *(Barrier_array + leg)) * dir2_sign;
+
+        if(knockedLeg == -1 && knockin > 0){
+          knockedLeg = leg;
+        }
+        if(knockedLeg_dash == -1 && knockin_dash > 0){
+          knockedLeg_dash = leg;
+        }
+      } else 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 )){
@@ -211,10 +259,10 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
         }
 
         if (ConvertNotional == 1){
-          if(Spot < *( Ps_array + leg )){
-            Ns_array[leg] = *(Ns_array + leg) / *( LSts_array + leg );
+          if(Spot < *( Ps_array + leg ) || Spot_dash < *( Ps_array + leg )){
+            Ns_array[leg] = *(TempNs_array + leg) / *( LSts_array + leg );
           } else {
-            Ns_array[leg] = *(Ns_array + leg) / *( USts_array + leg );
+            Ns_array[leg] = *(TempNs_array + leg) / *( USts_array + leg );
           }
         }
 
@@ -246,6 +294,10 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
             profit_loss = 0.0;
           } else if(Spot <= *( CFs_array + leg )){
             profit_loss = Spot - (*( CFs_array + leg ));
+
+            if (ConvertNotional == 1){
+              Ns_array[leg] = *(TempNs_array + leg) / *( CFs_array + leg );
+            }
           } else {
             profit_loss = Spot - (*( Xs_array + leg ));
           }
@@ -254,6 +306,10 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
             profit_loss_dash = 0.0;
           } else if(Spot_dash <= *( CFs_array + leg )){
             profit_loss_dash = Spot_dash - (*( CFs_array + leg ));
+
+            if (ConvertNotional == 1){
+              Ns_array[leg] = *(TempNs_array + leg) / *( CFs_array + leg );
+            }
           } else {
             profit_loss_dash = Spot_dash - (*( Xs_array + leg ));
           }
@@ -281,6 +337,10 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
             profit_loss = 0.0;
           } else if(Spot >= *( CFs_array + leg )){
             profit_loss = ( ( *( CFs_array + leg ) ) - Spot );
+
+            if (ConvertNotional == 1){
+              Ns_array[leg] = *(TempNs_array + leg) / *( CFs_array + leg );
+            }
           } else {
             profit_loss = ( ( *( Xs_array + leg ) ) - Spot );
           }
@@ -289,6 +349,10 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
             profit_loss_dash = 0.0;
           } else if(Spot_dash >= *( CFs_array + leg )){
             profit_loss_dash = ( ( *( CFs_array + leg ) ) - Spot_dash );
+
+            if (ConvertNotional == 1){
+              Ns_array[leg] = *(TempNs_array + leg) / *( CFs_array + leg );
+            }
           } else {
             profit_loss_dash = ( ( *( Xs_array + leg ) ) - Spot_dash );
           }
@@ -322,7 +386,37 @@ 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 || KType == PIVOT_ABSOLUTE || KType == COLLAR_ABSOLUTE || KType == DOUBLE_STRIKE_ABSOLUTE ) {
+    if(KType == FX_AMERICAN_BARRIER_KNOCKIN_DISCRETE){
+      if(knockedLeg >= 0){
+        double equivalent_notional = get_equivalent_notional(conversion_sign, ( *( Ns_array + knockedLeg ) ), sim[NL-1]);
+        sim_pos[NL-1] = european_payoff(( *( Xs_array + (knockedLeg) ) ), sim[NL-1], cp_sign, dir_sign, equivalent_notional);
+      }
+      else{
+        sim_pos[NL-1] = get_equivalent_rebate(conversion_sign, ( *( Rebate_array + 0 ) ), sim[NL-1]);
+      }
+      if(knockedLeg_dash >= 0 ){
+        double equivalent_notional = get_equivalent_notional(conversion_sign, ( *( Ns_array + knockedLeg_dash ) ), sim_dash[NL-1]);
+        sim_dash_pos[NL-1] = european_payoff(( *( Xs_array + (knockedLeg_dash) ) ), sim_dash[NL-1], cp_sign, dir_sign, equivalent_notional);
+      }
+      else{
+        sim_dash_pos[NL-1] = get_equivalent_rebate(conversion_sign, ( *( Rebate_array + 0 ) ), sim_dash[NL-1]);
+      }
+    } else if(KType == FX_AMERICAN_BARRIER_KNOCKOUT_DISCRETE){
+      if(knockedLeg >= 0){
+        sim_pos[knockedLeg] = get_equivalent_rebate(conversion_sign, ( *( Rebate_array + knockedLeg ) ), sim[knockedLeg]);
+      }
+      else{
+        double equivalent_notional = get_equivalent_notional(conversion_sign, ( *( Ns_array + 0 ) ), sim[NL-1]);
+        sim_pos[NL-1] = european_payoff(( *( Xs_array + 0 ) ), sim[NL-1], cp_sign, dir_sign, equivalent_notional);
+      }
+      if(knockedLeg_dash >= 0 ){
+        sim_dash_pos[knockedLeg_dash] = get_equivalent_rebate(conversion_sign, ( *( Rebate_array + knockedLeg_dash ) ), sim[knockedLeg_dash]);
+      }
+      else{
+        double equivalent_notional = get_equivalent_notional(conversion_sign, ( *( Ns_array + 0 ) ), sim_dash[NL-1]);
+        sim_dash_pos[NL-1] = european_payoff(( *( Xs_array + 0 ) ), sim_dash[NL-1], cp_sign, dir_sign, equivalent_notional);
+      }
+    } else 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 ) {
@@ -365,7 +459,7 @@ VALUE method_run_monte_carlo( VALUE self, VALUE args ) {
     } else if( KType == POINTS || KType == PIVOT_POINTS || KType == COLLAR_POINTS || KType == DOUBLE_STRIKE_POINTS ) {
       for( leg = 0; leg < NL; ++leg ) {
         // simulation normal
-        rb_p(DBL2NUM(*(Ns_array + leg )));
+        // rb_p(DBL2NUM(*(Ns_array + leg )));
         if ( ko_so_far > 0.0 ) {
           if( sim[leg] >= 0.0 ) {
             // knock out condition
@@ -529,4 +623,4 @@ void Init_tarf_monte_carlo() {
   // void rb_define_singleton_method(VALUE klass, const char *name, VALUE (*func)(), int argc)
   rb_define_singleton_method(cTarfMonteCarlo, "box_muller", method_box_muller, 0);
   rb_define_singleton_method(cTarfMonteCarlo, "run_monte_carlo", method_run_monte_carlo, -2);
-}
+}

data/lib/tarf_monte_carlo/version.rb

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

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: tarf_monte_carlo
 version: !ruby/object:Gem::Version
-  version: '3.23'
+  version: '3.28'
 platform: ruby
 authors:
 - Vivek Routh
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2020-09-24 00:00:00.000000000 Z
+date: 2020-10-03 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -114,6 +114,7 @@ files:
 - Rakefile
 - bin/console
 - bin/setup
+- ext/tarf_monte_carlo/Makefile
 - ext/tarf_monte_carlo/extconf.rb
 - ext/tarf_monte_carlo/tarf_monte_carlo.c
 - lib/.DS_Store