Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [numerical-methods]

The tag has no usage guidance.

2
votes
0answers
38 views

BSM PDE and diagonal dominance

I'm using the Thomas Algorithm to solve the Crank-Nicolson FD scheme. However, I've found that the matrix can occasionally fail to be diagonally dominant particularly when flat/local vol is very small....
4
votes
2answers
91 views

Implied volatility in Monte Carlo models

Suppose I want to get the implied volatility for a given option, whose process does not generate a closed-form formula. In that framework, how is the IV calculated, given the fact that bisection ...
2
votes
1answer
69 views

Finite difference methods for (continuously) strike-resettable American options

For simplicity, let us consider an American call/put with a continuously resettable strike price. Current time is $t=0$, maturity is at $t=T$, and the initial strike is $K_0$. We consider a "...
2
votes
0answers
76 views

Architecture of a global pricing library with immutable payoffs

By global pricing library I mean a library handling equity, rate etc, hybrid products having several models (BS, LV, SV, LSV) having several numerical methods (analytic formula, MC, PDE FD/FE) I ...
1
vote
0answers
121 views

What is the best source to get 10 milliseconds time-series data for numerical computation?

I am working with 4th order Runge-Kutta method to compute a second order differential equation. For the best accuracy, I need a 10 milliseconds ohlcv time-series data. I know that I can build it ...
3
votes
0answers
67 views

Is there a more efficient data structure to implement binomial trees than 2d array?

I'm just curious what is the "industry standard" for implementing a binomial tree (if "standards" exist in this case). For simplicity, let's just talk about the simplest trees with recombining nodes. ...
1
vote
1answer
59 views

In search of nice (approx) function forms of the volatility of cumulative simple returns

Let's consider a period $t\in[0,T]$, and let the simple return over year $t$ ($1\le t\le T$) be $r_t$. Assume $r_t$ are iid normal. The cumualative simple return over the whole period $[0,T]$ is $$R_T=...
1
vote
2answers
102 views

Paper recommendation with examples

Could you please advise me some papers/working documents with applications mainly focused on Fixed Income/Financial Engineering (Numerical Methods) as the one in the link below? Preferably on Matlab, ...
1
vote
0answers
39 views

Optimal allocation problem by finite differences

I am attempting to apply implicit finite difference to solve Merton's problem of optimal portfolio allocation for constant parameters. The equation to solve is the Hamilton-Jacobi-Bellman equation: $$...
1
vote
1answer
39 views

How many decimals of accuracy can I expect from FDM and MC (both valuation and risk)

I have implemented some Monte Carlo and FDM code. I can then get greeks by bumping. I am comparing to to exact formulas of price + greeks, and am wondering how many decimals of accuracy I can expect ...
0
votes
0answers
35 views

How to examine the impact of the parameters in the Hull White Model on the yield curve

I want to examine the yield curve resulting from the 2 Factor Hull-White model. Is there any way to examine the influence of the parameters on the yields curve without calibrating the model?
2
votes
1answer
186 views

SABR PDE spot/forward upper boundary condition implementation

When running my Finite Difference code, I observe something odd. Although implementing a classical (non-reverting) SABR model, I initialized the variables such that it should be equal to Black-...
1
vote
1answer
124 views

Finite Difference implicit scheme

I'm trying to solve the following PDE numerically using an implicit FD scheme: \begin{equation} \frac{\sigma_s^2}{2}\frac{\partial^2 V}{\partial S^2} + \rho \sigma_S \sigma_\alpha\frac{\partial^2 V}{\...
2
votes
1answer
127 views

Unable to obtain correct Finite Difference Results

A rather general question regarding a specific problem I am facing with my Matlab implementation of the implicit FD method for this PDE: \begin{equation} \frac{\sigma_s^2}{2}\frac{\partial^2 V}{\...
4
votes
1answer
220 views

Finite Difference method in Matlab for SABR volatility model fails to provide correct option values

Currently, I'm trying to implement a Finite Difference (FD) method in Matlab for my thesis (Quantitative Finance). I implemented the FD method for Black-Scholes already and got correct results. ...
0
votes
0answers
47 views

Computation of an integral containing d ln x (Scale of Market Shocks)

I am trying to implement a Scale of Market Shocks method (SMS) which was presented in a 1999 working document by Olsen & Associates named Introducing a Scale of Market Shocks and later refined in ...
0
votes
1answer
57 views

How to simulate a path through its solution and conditional expectation / variance

Hi I want to simulate in Matlab the following stochastic integral: $ x(t) = x(s) e^{-a(t-s)} + \sigma \int_s^t e^{-a(t-u)} dW_1(u)$ with $E[x(t) \vert F_s] = x(s) e^{-a(t-s)}$ $Var[x(t) \vert ...
0
votes
0answers
26 views

Transform the payoff to be non-zero

Is there any way to transform the basic call option payoff $V(s,0) = \max(s-K,0)$ such that $g(V(s,0))\neq 0$ $\forall s $, where $g()$ is the transform function of the payoff. This is to use in a ...
2
votes
0answers
37 views

Transforming and minimisation of the BS PDE

I'm trying a novel numerical substitution/fitting method to solve the BS PDE, but the issue is that due to the large range of magnitude of prices $V(s,t)\in[10^{-20},10^1]$, when I try to minimise the ...
1
vote
1answer
52 views

4-point Trapezium rule for numerical integration

Background: This is in reference to Mark Joshi's concepts of mathematical finance ch.7 problem 11. Question: We have in the Black-Scholes model: $S_0 = 1, T = 1, \sigma = 0.1, r = 0$. A ...
0
votes
1answer
84 views

Correct way to compute RSI in a moving window

I am trying to compute the RSI using the ticker price. I have troubles at the time to implement the RSI when I want to compute it every 1 minute but using windows to compute the OHLC of 10 minutes. I ...
2
votes
0answers
70 views

Practical precision for Options Pricing

When pricing options, especially in the theoretical literature getting high precision, say up to 8 decimal places is always a competitive goal. Though realistically in a practical setting is such ...
1
vote
0answers
301 views

Brennan-Schwartz algorithm for pricing American options

I'm reading Pricing American Options using LU decomposition by Ikonen and Toivanen (IT). They reference The valuation of American put options by Brennan and Schwartz, and cast it as method that uses ...
1
vote
1answer
113 views

MATLAB exercise on an European call option with time-varying volatility

I have to solve the following exercise: compute and plot the value $V = V(S, t),\ t<T$, ($T=$ maturity) of an European CALL option (with arbitrary $t$, $T$, $K$ (strike price), $r$ (risk-free ...
1
vote
0answers
172 views

Adjoint Algorithmic Differentiation: swap pricing

I have tried to implement an AAD routine to price call options using the Black-Scholes formula, but my greeks are not quite agreeing with the expected ones, so I have decided to start with something a ...
0
votes
1answer
198 views

Euler discretisation error for stochastic volatility model

Given the following model$$dS_t=S_t(\mu dt+\sigma(t,S_t)dW_t)$$ Using Monte Carlo Pricing method, I want to determine the price of the option. However I have been encountered the following problems: ...
0
votes
0answers
51 views

what is the meaning of $U^{n+1/3}$ ADI method

For the ADI in numerical method $$\frac{U^{n+1/3}-U^n}{k/3} = \Delta^2_x U^{n+1/3} + \Delta^2_y U^n + \Delta^2_z U^{n+2/3}$$ $$....$$ $$....$$ don't like $U^{n+1/2} = \dfrac 1 2 (U^{n+1}+U^n),$ I ...
1
vote
1answer
52 views

Extreme cases of normal random numbers and NaN

While trying to implement my version of Euler's method for simulating a SDE in C++, I came up with a problem. It occurs in some cases that the path generated by my method ends up giving values which ...
1
vote
1answer
137 views

How should I develop my coding ability in order to set myself up for a quant role? [closed]

First of all, I apologise if similar questions have already been asked; I've googled around but most similar questions aren't focused on developing specifically quant-friendly programming skills. I'm ...
0
votes
1answer
190 views

How to calibrate an SDE's by finite difference equation?

I would like a general framework for the calibration of the unknown parameters in an arbitrary stochastic differential equation. I have a proposed method that seems reasonable in theory, but is ...
5
votes
2answers
267 views

Smoothing of the payoff function as a terminal condition for numerical option pricing

I am interested in using a 4th order finite difference method in (underlying asset) space to price a European call basket option. I have developed the solver and everything works as expected, except ...
2
votes
1answer
182 views

pricing american put option with fdm

Assume I use some finite difference solver to solve for American type of exercise in BS framework where stock pays dividend discretely. Then at every time iteration, for call option, I firstly adjust ...
6
votes
5answers
311 views

Heston Model Integration Oscillations

Is there a way to reduce oscillations for the numerical integration when evaluating the Heston model. I am pricing a series of 5000 options scattered over the Heston model parameter space and I find ...
5
votes
4answers
208 views

Is it possible to defend a Computer Science master thesis by writing a project about quantitative finance?

What are features and examples of computational finance (financial computing) problems (for thesis project in Master in Computer Science)? Is it possible to defend Computer Science master thesis by ...
1
vote
1answer
260 views

Numerical Methods for Merton Model

The stochastic differential equation for an underlying with jumps in Merton model is: $$d{{S}_{t}}=\mu \,{{S}_{t}}dt+\sigma \,{{S}_{t}}\,d{{W}_{t}}^{P}+(J-1){{S}_{t}}d{{q}_{t}}$$ where $t \quad\,\,\, ...
5
votes
1answer
269 views

Solve Black scholes PDE without using any transformation

I know that one of the methods of solving the black scholes PDE given by : $\frac{\partial V}{\partial t} + \frac{\sigma^2 S^2}{2}\frac{\partial^2V}{\partial S^2} + rS\frac{\partial V}{\partial S} -rV ...
3
votes
1answer
219 views

Numerical Optimizer Matlab Calibration LMM

I am trying to mimimize the following function in order to calibrate the Libor Market Model $$\sum_{i=1}^{n} \left(\sigma_i^{market}-\sigma_i^{Reb}\left(a,b,c,d,\beta\right)/\sqrt{T_i}\right)^2,$$ ...
1
vote
0answers
156 views

School project about Black Scholes with stochastic volatility

In a university project I am looking at Black Scholes model with a stochastic volatility. I’m still not quite sure about my focus (I am in the beginning 'Idea phase'). I want to explain the theory ...
1
vote
0answers
54 views

Jacobian for Newton method for American options by front fixing

In this paper Penalty and front-fixing methods for the numerical solution of American option problems a front fixing method based on Newton is described for an American put option is described. I am ...
1
vote
1answer
109 views

Why can't we use Finite Differences with non-parabolic PDEs?

The title of the question says it all. Why can we only apply the method to parabolic PDEs like the heat equation, and not to ordinary PDEs?
-1
votes
2answers
390 views

Numerical computation of Heston model Integral: Simpsone Rule or Gauss-Legendre Method

I want to price a call option using the Heston model for a given set of parameters. theory from URL: http://elis.sigmath.es.osaka-u.ac.jp/research/Heston-original.pdf The integral equation (18) ...
1
vote
0answers
27 views

Stiffness of numerical methods for SDE

What can I do with stiffness of numerical methods for SDE? I want to use numerical approach for solving SDE in market's scenarios generation. Is there any general approach to handle it?
2
votes
3answers
517 views

Numerical Solution to BS PDE - Digital Option

Here is a relatively simple question about PDE's pricing. Assume that we are within the BS framework and moreover that interest rate is zero. The price $V(t,S_t)$ of the digital is known to be $\Phi(...
1
vote
1answer
237 views

Numerical Solutions for PIDE

I want to solve an exotic options of PIDE by Numerical Methods.I just focus on the integral part of PIDE and want to underestand some tips on numerical solution of how to numerically solve it. Exactly ...
3
votes
0answers
75 views

Time discretisations, FDM vs FEM

I am interested in adaptive mesh methods for numerical solution of PDEs with applications to finance. As part of a school project, I have been pricing vanilla European call and put options using 2D ...
8
votes
1answer
715 views

Speeding up computations: when to use Quasi and standard Monte-Carlo in pricing

I am familiar with the theory of Monte-Carlo techniques in the numerical integration, and recently I have started my experiments with these methods applied to derivatives pricing. I am using ...
5
votes
4answers
8k views

What does “convergence” in Monte Carlo simulation mean?

I have read about convergence in terms of MC simulation for derivative pricing, but I am not clear on what it exactly means. Let us suppose I price an option 100,000 paths twice and both result in the ...
1
vote
2answers
478 views

Numerical delta of Bond Options

I'm trying to calculate the delta for bond Call options. I'm using the vasicek model which gives the following solution for a Zero-coupon bond call option: $Z = N P(t,S) \Phi(d_1) - K P(t,T) \Phi(d_2)...
3
votes
1answer
2k views

Covariance matrix and Cholesky decomposition

I am simulating a spread option with stochastic volatility using Monte Carlo simulation. I have the positive-definite covariance matrix $$ \rho = \left( \begin{array}{cccc} 1 & \rho_{1,2} & \...
0
votes
1answer
521 views

how to make a nonlinear gird where grid points are not equally spaced?

I need to make a grid [0,1] with points that are concentrated close to the edges (close to 0 and 1) while the remaining points in the middle can be equally spaced. The reason for doing this is that I ...