# Tag Info

3

I am not sure if I understood your question correctly but I will try to answer it anyway. If you have a standard normal random vector $z \sim N(\mathbb{0},I_n)$ (where $z,0 \in \mathbb{R}^{n\times1}$ and $I_n \in \mathbb{R}^{n\times n}$ is the identity matrix) and you want to transform it into a multivariate normal $x \sim N(\mu,\Sigma)$ you do it the ...

3

The Papageorgiou paper is presumably referring specifically to quasi-random sequences used in path generation. Researchers had noticed that, in high dimensions, QR sequences tend to have good space coverage for the first couple of dimensions: but terrible coverage for the latter dimensions: (Plots here are points 101-200 from a 32-dimensional QR ...

3

The method described in Hallerbach (2004) always worked well for me. We derive an estimator for Black-Scholes-Merton implied volatility that, when compared to the familiar Corrado & Miller [JBaF, 1996] estimator, has substantially higher approximation accuracy and extends over a wider region of moneyness.

2

You can find a brief but useful explanation of Brownian bridge techniques in Andersen and Piterbarg (page 125), which includes references for further reading. It's probably the best place to start. They discuss valuing barrier options specifically, and discuss the performance issues mentioned here. Later (pg 647), they use Brownian bridges in constructing ...

2

There is a qualitative shift in the shape of the density. When V is small it is monotone decaying. When V is large it looks more like a Gaussian. Another reason he uses two schemes is that he wants match two moments of the density. When V is small, the moment matching equations for the quadratic Gaussian are unsolvable. When V is large they are unsolvable ...

2

Peter Jaeckel wrote a paper just on how to solve this problem: By Implication (July 2006; Wilmott, pages 60-66, November 2006). Probably the most complicated trivial issue in financial mathematics: how to compute Black's implied volatility robustly, simply, efficiently, and fast downloadable from jaeckel.org In my experience the most important thing is to ...

1

Presumably you are trying to use a finite difference method to solve a differential equation. The non-uniformity of the grid has an impact on accuracy. Hence, it is useful to include a parameter in the grid-generation algorithm that controls the rate at which the spacing increases away from the boundary. There are many approaches for generating non-uniform ...

1

Below is the root search algorithm code I wrote in college. This is written in octave. It's simple to understand and re-write in C++. Develop numerical methods algos as a separate module and integrate with your pricing and other code I want to WARN you to re-check for bugs. It always converges for my objective functions First function is Dekker method ...

1

Bracketing methods such as Bisection and Regula Falsi are always known to converge but they are very slow. Newton Raphson and secant methods are fast (quadratic convergence) but has convergence problems. Google for Newton Raphson convergence pitfalls. Classical ones such as"Trapped in local minima", "Diverge instead of converge" etc Algorithms such as ...

1

Use your total wealth allocated to the trades as denominator. Total wealth allocated would include all collateral. In this way you (or your broker) make sure that the denominator is always positive. Presumably this would also reflect what you really want to track. The only problem that remains is what amount of your wealth needs to be allocated. But this is ...

1

1) Brownian Bridge is used in Quasi Monte Carlo pricing of asian options to reexpress paths in a basis where few selected components/subspaces bring the most contribution, so as to align these to the best distributed dimensions/subspaces of a low discrepancy sequence. This allows for better coverage and thus faster convergence and paths amount reduction. ...

1

The only special function needed for computing Black-Scholes option prices is the cumulative normal function ("N" or "Phi") or equivalently the error function ("erf"). These are very widely available with good standard library implementations. The erf function in single and double precision is part of the c99 and c++11 math standard libraries. For your ...

Only top voted, non community-wiki answers of a minimum length are eligible