Near Real-Time Risk Simulation of Complex Portfolios on Heterogeneous Computing Systems with OpenCL

Size: px
Start display at page:

Download "Near Real-Time Risk Simulation of Complex Portfolios on Heterogeneous Computing Systems with OpenCL"

Transcription

1 Near Real-Time Risk Simulation of Complex Portfolios on Heterogeneous Computing Systems with OpenCL Javier Alejandro Varela, Norbert Wehn Microelectronic Systems Design Research Group University of Kaiserslautern, Germany Javier Alejandro Varela, Norbert Wehn Near Real-Time Risk Simulation of Complex Portfolios on Heterogeneous Computing Systems with OpenCL. In Proceedings of IWOCL 17, Toronto, Canada, May 16-18, 2017, 10 pages. DOI:

2 Summary I. Overview: Portfolio Risk II. OpenCL + Workload Allocation III. Algorithmic Optimizations / Numerical Scheme IV. Minimizing Device Global Memory V. Conclusion 2

3 Summary I. Overview: Portfolio Risk II. OpenCL + Workload Allocation III. Algorithmic Optimizations / Numerical Scheme IV. Minimizing Device Global Memory V. Conclusion 3

4 Portfolio Composition: Overview: Portfolio Risk Stocks (S 1, S 2 ) Corporate Bond Foreign Currency Options: European Asian (on S 1 ) What is the RISK at which this portfolio is exposed? European Barrier (on S 2 ) American Vanilla (on S 1 ) American 2D Max (on S 1,S 2 ) Cash (no simulation required) 4

5 Overview: Portfolio Risk What is Value-at-Risk (VaR)? I am α percent certain there will not be a loss of more than VaR USD in the next N days (Hull, Options Futures and Other Derivatives). VaR of level α Portfolio loss function L α - quantile Var α L t+n inf l R P L t+n > l Loss empirical distribution function F L l P L t+n > l where α is as high as 95%, 99%, 99.5%. Industry standard (Basel II and Basel III) What is Expected Shortfall (ES or cvar)? ES α L t+n 1 1 α α (1 α) = inf l R F L l α 1 Var γ L t+n dγ 5

6 Overview: Portfolio Risk How to compute Value-at-Risk (VaR)? 1. Historical simulation: describes future changes based on empirical distribution of observed past data. 2. Variance-Covariance method: first-order approximation of the loss function L. assumes normally distributed returns. 3. Monte Carlo (MC) method: simulates the loss function L. Does not rely on historical data. No need for approximation. No assumption of a normal distribution. Computationally intensive problem Heterogeneous problem (many different algorithms involved) Exploit OpenCL 6

7 Overview: Portfolio Risk Portfolio Simulation: External Internal Stocks (S 1, S 2 ) Corporate Bond Foreign Currency Options: European Asian (on S 1 ) European Barrier (on S 2 ) American Vanilla (on S 1 ) American 2D Max (on S 1,S 2 ) Cash (no simulation required) Two Models: Black-Scholes (BS) Heston t t + N 7

8 Overview: Portfolio Risk 252 steps (1008 steps for european Barrier) External Internal 32k paths 32k paths 1 step (BS model) 4 steps (Heston model) t t + N Single-precision floating point operations 8

9 Sequence of steps: Overview: Portfolio Risk Asset 1 Asset N P&L portf Sorted P&L portf Simulated Scenarios (.h 1 ) + + (.h N ) Portf. value - today = Sorting (ascending order) VaR ES (mean) where: h i = holding of Asset i (for i=1,,n) 9

10 Summary I. Overview: Portfolio Risk II. OpenCL + Workload Allocation III. Algorithmic Optimizations / Numerical Scheme IV. Minimizing Device Global Memory V. Conclusion 10

11 OpenCL + Workload Allocation Asset 1 Asset N P&L portf Sorted P&L portf Simulated Scenarios (.h 1 )+ +( Portf..h N ) - value = today Sorting (ascending order) VaR ES (mean) where: h i = holding of Asset i (for i=1,,n) SuperMicro Superserver 7048GR-TR Intel Xeon E v3 CPU (Host) OpenCL Devices: CPU Xeon PHI (dual-socket motherboard) Intel Xeon E v3 Intel Xeon Phi 7120P Individual OpenCL Command Queues GPU 0 GPU 1 (shared PCIe link) Nvidia K80 11

12 OpenCL + Workload Allocation ILP formulation (Matlab s intlinprog) min x f T. x subject to x i are integers A. x b A eq. x b eq lb x lu Set of Kernels: K = k Stocks2Corr, k Bond, k F.Currency, Set of Devices: D = d CPU, d PHI, d GPU0, Binary decision variables: x k,d, where: 0 x k,d 1 Total runtime: x T, where: 0 x T 12

13 OpenCL + Workload Allocation ILP formulation (continued) Constraints: Single assignment of each kernel among all devices: k K, d D x k,d = 1 Device runtime: d D, k K t k,d. x k,d x T 0 Device global memory: d D, m k,d. x k,d M d k K Cost function: f T. x = x T 13

14 OpenCL + Workload Allocation Workload Allocation with ILP 1) Profile each kernel on every device 2) ILP CPU (Host) OpenCL Devices: CPU Xeon PHI GPU 0 GPU 1 14

15 OpenCL + Workload Allocation Kernels Runtime CPU (Host) OpenCL Devices: CPU Xeon PHI GPU 0 GPU 1 15

16 Summary I. Overview: Portfolio Risk II. OpenCL + Workload Allocation III. Algorithmic Optimizations / Numerical Scheme IV. Minimizing Device Global Memory V. Conclusion 16

17 Algorithmic Optimizations Random Numbers Reuse Approach Nested MC-Based Risk Measurement of Complex Portfolios: Acceleration and Energy Efficiency S. Desmettre, R. Korn, J. Varela, N. Wehn. Risks Vol. 4, no. 4, pages 36, October, (All parameters in the internal simulation remain constant) Incorrect Different option prices Correct 17

18 Algorithmic Optimizations Paths-Reuse Approach (+Normalization) (Note: All parameters in the internal simulation remain constant) 18

19 Algorithmic Optimizations Paths-Reuse: European-style Options 0Path Values [USD] Normalized S0=1 Kernel P1 Internal Simulation Maturity Maturity Max Value / Path Min Value / Path Simulation Steps Extracted Stochastic Information Option Price [USD] Kernel P2 Information-Reuse Approach 19

20 Algorithmic Optimizations Paths-Reuse: American-style Options Kernel P1 Internal Simulation 0Path Values [USD] Normalized S0=1 Simulation Steps Maturity Option price [USD] Kernel P2 Backward Step-By-Step Decision Process (Regression) 20

21 Algorithmic Optimizations Runtime: Paths-/Info-Reuse (under BS model) 4.5x 6.8x 1. RNs-reuse 2. Info-reuse 3. Paths-reuse 1000x 1000x European Options American Options 21

22 Algorithmic Optimizations Runtime: Paths-/Info-Reuse (under BS model) 9.2x 12.0x 1. RNs-reuse 2. Info-reuse 3. Paths-reuse 1500x 2500x European Options American Options 22

23 Numerical Scheme Interpolation: Kernel: AmericanVanillaBT AmericanVanillaLS P1 P2 Fast, accurate and inelegant valuation of American options. Adriaan Joubert and L.C.G. Rogers In Proceedings of the Numerical Methods Workshop at the Isaac Newton Institute, April 1995, L.C.G. Rogers and D. Talay (Eds.). CambridgeS. Kernel: Interpolate 1D Option Price [USD] Extension: American2DMax P1 P2 Interpolate 2D Option Price [USD] 23

24 Numerical Scheme Runtime: Interpolation (American-style Options) 10000x 1. RNs-reuse 2. Interpolation+Paths-reuse 3. Interpolation+Surface Load&Transfer 1000x

25 Summary I. Overview: Portfolio Risk II. OpenCL + Workload Allocation III. Algorithmic Optimizations / Numerical Scheme IV. Minimizing Device Global Memory V. Conclusion 25

26 Minimizing Device Global Memory Device Global Memory Cost (Kernels) 26

27 Minimizing Device Global Memory Autoaggregation Asset 1 Asset N P&L portf Sorted P&L portf Simulated Scenarios (.h 1 )+ +( Portf..h N ) - value = today Sorting (ascending order) VaR ES (mean) CPU (Host) where: h i = holding of Asset i (for i=1,,n) OpenCL Devices: CPU Xeon PHI Device global memory Accum. Minimizes Host workload Minimizes data transfers via PCIe GPU 0 GPU 1 Single output vector per device 27

28 Summary I. Overview: Portfolio Risk II. OpenCL + Workload Allocation III. Algorithmic Optimizations / Numerical Scheme IV. Minimizing Device Global Memory V. Conclusion 28

29 I. Overview: Portfolio Risk II. Conclusion Simulated Scenarios Asset 1 OpenCL + Workload Allocation (.h 1 )+ +( Asset N P&L portf Portf..h N ) - value = today Sorting (ascending order) Sorted P&L portf VaR ES (mean) CPU (Host) OpenCL Devices: CPU Xeon PHI GPU 0 GPU 1 III. Algorithmic Optimizations / Numerical Scheme 4.5x 6.8x 10000x 1000x 1000x 1000x IV. Minimizing Device Global Memory Portability (Numerical Results) 29

30 Questions? Thanks for your attention 30

Financial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA

Financial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA Financial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA Rajesh Bordawekar and Daniel Beece IBM T. J. Watson Research Center 3/17/2015 2014 IBM Corporation

More information

Contents Critique 26. portfolio optimization 32

Contents Critique 26. portfolio optimization 32 Contents Preface vii 1 Financial problems and numerical methods 3 1.1 MATLAB environment 4 1.1.1 Why MATLAB? 5 1.2 Fixed-income securities: analysis and portfolio immunization 6 1.2.1 Basic valuation of

More information

Economic Scenario Generator: Applications in Enterprise Risk Management. Ping Sun Executive Director, Financial Engineering Numerix LLC

Economic Scenario Generator: Applications in Enterprise Risk Management. Ping Sun Executive Director, Financial Engineering Numerix LLC Economic Scenario Generator: Applications in Enterprise Risk Management Ping Sun Executive Director, Financial Engineering Numerix LLC Numerix makes no representation or warranties in relation to information

More information

Accelerating Financial Computation

Accelerating Financial Computation Accelerating Financial Computation Wayne Luk Department of Computing Imperial College London HPC Finance Conference and Training Event Computational Methods and Technologies for Finance 13 May 2013 1 Accelerated

More information

Fast Convergence of Regress-later Series Estimators

Fast Convergence of Regress-later Series Estimators Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser

More information

F1 Acceleration for Montecarlo: financial algorithms on FPGA

F1 Acceleration for Montecarlo: financial algorithms on FPGA F1 Acceleration for Montecarlo: financial algorithms on FPGA Presented By Liang Ma, Luciano Lavagno Dec 10 th 2018 Contents Financial problems and mathematical models High level synthesis Optimization

More information

Barrier Option. 2 of 33 3/13/2014

Barrier Option. 2 of 33 3/13/2014 FPGA-based Reconfigurable Computing for Pricing Multi-Asset Barrier Options RAHUL SRIDHARAN, GEORGE COOKE, KENNETH HILL, HERMAN LAM, ALAN GEORGE, SAAHPC '12, PROCEEDINGS OF THE 2012 SYMPOSIUM ON APPLICATION

More information

The Evaluation of American Compound Option Prices under Stochastic Volatility. Carl Chiarella and Boda Kang

The Evaluation of American Compound Option Prices under Stochastic Volatility. Carl Chiarella and Boda Kang The Evaluation of American Compound Option Prices under Stochastic Volatility Carl Chiarella and Boda Kang School of Finance and Economics University of Technology, Sydney CNR-IMATI Finance Day Wednesday,

More information

Asian Option Pricing: Monte Carlo Control Variate. A discrete arithmetic Asian call option has the payoff. S T i N N + 1

Asian Option Pricing: Monte Carlo Control Variate. A discrete arithmetic Asian call option has the payoff. S T i N N + 1 Asian Option Pricing: Monte Carlo Control Variate A discrete arithmetic Asian call option has the payoff ( 1 N N + 1 i=0 S T i N K ) + A discrete geometric Asian call option has the payoff [ N i=0 S T

More information

Accelerated Option Pricing Multiple Scenarios

Accelerated Option Pricing Multiple Scenarios Accelerated Option Pricing in Multiple Scenarios 04.07.2008 Stefan Dirnstorfer (stefan@thetaris.com) Andreas J. Grau (grau@thetaris.com) 1 Abstract This paper covers a massive acceleration of Monte-Carlo

More information

Domokos Vermes. Min Zhao

Domokos Vermes. Min Zhao Domokos Vermes and Min Zhao WPI Financial Mathematics Laboratory BSM Assumptions Gaussian returns Constant volatility Market Reality Non-zero skew Positive and negative surprises not equally likely Excess

More information

GRAPHICAL ASIAN OPTIONS

GRAPHICAL ASIAN OPTIONS GRAPHICAL ASIAN OPTIONS MARK S. JOSHI Abstract. We discuss the problem of pricing Asian options in Black Scholes model using CUDA on a graphics processing unit. We survey some of the issues with GPU programming

More information

Implementing Models in Quantitative Finance: Methods and Cases

Implementing Models in Quantitative Finance: Methods and Cases Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1

More information

Market Risk Analysis Volume IV. Value-at-Risk Models

Market Risk Analysis Volume IV. Value-at-Risk Models Market Risk Analysis Volume IV Value-at-Risk Models Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume IV xiii xvi xxi xxv xxix IV.l Value

More information

Portfolio selection with multiple risk measures

Portfolio selection with multiple risk measures Portfolio selection with multiple risk measures Garud Iyengar Columbia University Industrial Engineering and Operations Research Joint work with Carlos Abad Outline Portfolio selection and risk measures

More information

GPU-Accelerated Quant Finance: The Way Forward

GPU-Accelerated Quant Finance: The Way Forward GPU-Accelerated Quant Finance: The Way Forward NVIDIA GTC Express Webinar Gerald A. Hanweck, Jr., PhD CEO, Hanweck Associates, LLC Hanweck Associates, LLC 30 Broad St., 42nd Floor New York, NY 10004 www.hanweckassoc.com

More information

Scenario-Based Value-at-Risk Optimization

Scenario-Based Value-at-Risk Optimization Scenario-Based Value-at-Risk Optimization Oleksandr Romanko Quantitative Research Group, Algorithmics Incorporated, an IBM Company Joint work with Helmut Mausser Fields Industrial Optimization Seminar

More information

Computer Exercise 2 Simulation

Computer Exercise 2 Simulation Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing

More information

High Performance and Low Power Monte Carlo Methods to Option Pricing Models via High Level Design and Synthesis

High Performance and Low Power Monte Carlo Methods to Option Pricing Models via High Level Design and Synthesis High Performance and Low Power Monte Carlo Methods to Option Pricing Models via High Level Design and Synthesis Liang Ma, Fahad Bin Muslim, Luciano Lavagno Department of Electronics and Telecommunication

More information

Analytics in 10 Micro-Seconds Using FPGAs. David B. Thomas Imperial College London

Analytics in 10 Micro-Seconds Using FPGAs. David B. Thomas Imperial College London Analytics in 10 Micro-Seconds Using FPGAs David B. Thomas dt10@imperial.ac.uk Imperial College London Overview 1. The case for low-latency computation 2. Quasi-Random Monte-Carlo in 10us 3. Binomial Trees

More information

Algorithmic Differentiation of a GPU Accelerated Application

Algorithmic Differentiation of a GPU Accelerated Application of a GPU Accelerated Application Numerical Algorithms Group 1/31 Disclaimer This is not a speedup talk There won t be any speed or hardware comparisons here This is about what is possible and how to do

More information

Financial Mathematics and Supercomputing

Financial Mathematics and Supercomputing GPU acceleration in early-exercise option valuation Álvaro Leitao and Cornelis W. Oosterlee Financial Mathematics and Supercomputing A Coruña - September 26, 2018 Á. Leitao & Kees Oosterlee SGBM on GPU

More information

IEOR E4602: Quantitative Risk Management

IEOR E4602: Quantitative Risk Management IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Market risk measurement in practice

Market risk measurement in practice Lecture notes on risk management, public policy, and the financial system Allan M. Malz Columbia University 2018 Allan M. Malz Last updated: October 23, 2018 2/32 Outline Nonlinearity in market risk Market

More information

Stochastic Grid Bundling Method

Stochastic Grid Bundling Method Stochastic Grid Bundling Method GPU Acceleration Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee London - December 17, 2015 A. Leitao &

More information

Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing

Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014

More information

Assessing Solvency by Brute Force is Computationally Tractable

Assessing Solvency by Brute Force is Computationally Tractable O T Y H E H U N I V E R S I T F G Assessing Solvency by Brute Force is Computationally Tractable (Applying High Performance Computing to Actuarial Calculations) E D I N B U R M.Tucker@epcc.ed.ac.uk Assessing

More information

Running Financial Risk Management Applications on FPGA in the Amazon Cloud

Running Financial Risk Management Applications on FPGA in the Amazon Cloud Running Financial Risk Management Applications on FPGA in the Amazon Cloud Javier Alejandro Varela, Norbert Wehn Microelectronic Systems Design Research Group, University of Kaiserslautern, 67663 Kaiserslautern,

More information

Chapter 5 Portfolio. O. Afonso, P. B. Vasconcelos. Computational Economics: a concise introduction

Chapter 5 Portfolio. O. Afonso, P. B. Vasconcelos. Computational Economics: a concise introduction Chapter 5 Portfolio O. Afonso, P. B. Vasconcelos Computational Economics: a concise introduction O. Afonso, P. B. Vasconcelos Computational Economics 1 / 22 Overview 1 Introduction 2 Economic model 3 Numerical

More information

2.1 Mathematical Basis: Risk-Neutral Pricing

2.1 Mathematical Basis: Risk-Neutral Pricing Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t

More information

Outline. GPU for Finance SciFinance SciFinance CUDA Risk Applications Testing. Conclusions. Monte Carlo PDE

Outline. GPU for Finance SciFinance SciFinance CUDA Risk Applications Testing. Conclusions. Monte Carlo PDE Outline GPU for Finance SciFinance SciFinance CUDA Risk Applications Testing Monte Carlo PDE Conclusions 2 Why GPU for Finance? Need for effective portfolio/risk management solutions Accurately measuring,

More information

Handbook of Financial Risk Management

Handbook of Financial Risk Management Handbook of Financial Risk Management Simulations and Case Studies N.H. Chan H.Y. Wong The Chinese University of Hong Kong WILEY Contents Preface xi 1 An Introduction to Excel VBA 1 1.1 How to Start Excel

More information

Computational Finance in CUDA. Options Pricing with Black-Scholes and Monte Carlo

Computational Finance in CUDA. Options Pricing with Black-Scholes and Monte Carlo Computational Finance in CUDA Options Pricing with Black-Scholes and Monte Carlo Overview CUDA is ideal for finance computations Massive data parallelism in finance Highly independent computations High

More information

RISKMETRICS. Dr Philip Symes

RISKMETRICS. Dr Philip Symes 1 RISKMETRICS Dr Philip Symes 1. Introduction 2 RiskMetrics is JP Morgan's risk management methodology. It was released in 1994 This was to standardise risk analysis in the industry. Scenarios are generated

More information

Numerix Pricing with CUDA. Ghali BOUKFAOUI Numerix LLC

Numerix Pricing with CUDA. Ghali BOUKFAOUI Numerix LLC Numerix Pricing with CUDA Ghali BOUKFAOUI Numerix LLC What is Numerix? Started in 1996 Roots in pricing exotic derivatives Sophisticated models CrossAsset product Excel and SDK for pricing Expanded into

More information

Pricing Early-exercise options

Pricing Early-exercise options Pricing Early-exercise options GPU Acceleration of SGBM method Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee Lausanne - December 4, 2016

More information

The Fundamental Review of the Trading Book: from VaR to ES

The Fundamental Review of the Trading Book: from VaR to ES The Fundamental Review of the Trading Book: from VaR to ES Chiara Benazzoli Simon Rabanser Francesco Cordoni Marcus Cordi Gennaro Cibelli University of Verona Ph. D. Modelling Week Finance Group (UniVr)

More information

HIGH PERFORMANCE COMPUTING IN THE LEAST SQUARES MONTE CARLO APPROACH. GILLES DESVILLES Consultant, Rationnel Maître de Conférences, CNAM

HIGH PERFORMANCE COMPUTING IN THE LEAST SQUARES MONTE CARLO APPROACH. GILLES DESVILLES Consultant, Rationnel Maître de Conférences, CNAM HIGH PERFORMANCE COMPUTING IN THE LEAST SQUARES MONTE CARLO APPROACH GILLES DESVILLES Consultant, Rationnel Maître de Conférences, CNAM Introduction Valuation of American options on several assets requires

More information

EC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods

EC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions

More information

Machine Learning for Quantitative Finance

Machine Learning for Quantitative Finance Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing

More information

An Energy Efficient FPGA Accelerator for Monte Carlo Option Pricing with the Heston Model

An Energy Efficient FPGA Accelerator for Monte Carlo Option Pricing with the Heston Model 2011 International Conference on Reconfigurable Computing and FPGAs An Energy Efficient FPGA Accelerator for Monte Carlo Option Pricing with the Heston Model Christian de Schryver, Ivan Shcherbakov, Frank

More information

Hedging Strategy Simulation and Backtesting with DSLs, GPUs and the Cloud

Hedging Strategy Simulation and Backtesting with DSLs, GPUs and the Cloud Hedging Strategy Simulation and Backtesting with DSLs, GPUs and the Cloud GPU Technology Conference 2013 Aon Benfield Securities, Inc. Annuity Solutions Group (ASG) This document is the confidential property

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Simulation Efficiency and an Introduction to Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Assicurazioni Generali: An Option Pricing Case with NAGARCH

Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: Business Snapshot Find our latest analyses and trade ideas on bsic.it Assicurazioni Generali SpA is an Italy-based insurance

More information

Approximations of Stochastic Programs. Scenario Tree Reduction and Construction

Approximations of Stochastic Programs. Scenario Tree Reduction and Construction Approximations of Stochastic Programs. Scenario Tree Reduction and Construction W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany www.mathematik.hu-berlin.de/~romisch

More information

Fast and accurate pricing of discretely monitored barrier options by numerical path integration

Fast and accurate pricing of discretely monitored barrier options by numerical path integration Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:

More information

Multistage risk-averse asset allocation with transaction costs

Multistage risk-averse asset allocation with transaction costs Multistage risk-averse asset allocation with transaction costs 1 Introduction Václav Kozmík 1 Abstract. This paper deals with asset allocation problems formulated as multistage stochastic programming models.

More information

PRICING AMERICAN OPTIONS WITH LEAST SQUARES MONTE CARLO ON GPUS. Massimiliano Fatica, NVIDIA Corporation

PRICING AMERICAN OPTIONS WITH LEAST SQUARES MONTE CARLO ON GPUS. Massimiliano Fatica, NVIDIA Corporation PRICING AMERICAN OPTIONS WITH LEAST SQUARES MONTE CARLO ON GPUS Massimiliano Fatica, NVIDIA Corporation OUTLINE! Overview! Least Squares Monte Carlo! GPU implementation! Results! Conclusions OVERVIEW!

More information

Calculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the

Calculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the VaR Pro and Contra Pro: Easy to calculate and to understand. It is a common language of communication within the organizations as well as outside (e.g. regulators, auditors, shareholders). It is not really

More information

Executive Summary: A CVaR Scenario-based Framework For Minimizing Downside Risk In Multi-Asset Class Portfolios

Executive Summary: A CVaR Scenario-based Framework For Minimizing Downside Risk In Multi-Asset Class Portfolios Executive Summary: A CVaR Scenario-based Framework For Minimizing Downside Risk In Multi-Asset Class Portfolios Axioma, Inc. by Kartik Sivaramakrishnan, PhD, and Robert Stamicar, PhD August 2016 In this

More information

Stochastic Programming in Gas Storage and Gas Portfolio Management. ÖGOR-Workshop, September 23rd, 2010 Dr. Georg Ostermaier

Stochastic Programming in Gas Storage and Gas Portfolio Management. ÖGOR-Workshop, September 23rd, 2010 Dr. Georg Ostermaier Stochastic Programming in Gas Storage and Gas Portfolio Management ÖGOR-Workshop, September 23rd, 2010 Dr. Georg Ostermaier Agenda Optimization tasks in gas storage and gas portfolio management Scenario

More information

Computational Finance. Computational Finance p. 1

Computational Finance. Computational Finance p. 1 Computational Finance Computational Finance p. 1 Outline Binomial model: option pricing and optimal investment Monte Carlo techniques for pricing of options pricing of non-standard options improving accuracy

More information

Exotic Derivatives & Structured Products. Zénó Farkas (MSCI)

Exotic Derivatives & Structured Products. Zénó Farkas (MSCI) Exotic Derivatives & Structured Products Zénó Farkas (MSCI) Part 1: Exotic Derivatives Over the counter products Generally more profitable (and more risky) than vanilla derivatives Why do they exist? Possible

More information

New Developments in MATLAB for Computational Finance Kevin Shea, CFA Principal Software Developer MathWorks

New Developments in MATLAB for Computational Finance Kevin Shea, CFA Principal Software Developer MathWorks New Developments in MATLAB for Computational Finance Kevin Shea, CFA Principal Software Developer MathWorks 2014 The MathWorks, Inc. 1 Who uses MATLAB in Financial Services? The top 15 assetmanagement

More information

New GPU Pricing Library

New GPU Pricing Library New GPU Pricing Library! Client project for Bank Sarasin! Highly regarded sustainable Swiss private bank! Founded 1841! Core business! Asset management! Investment advisory! Investment funds! Structured

More information

- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t

- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t - 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label

More information

HPC IN THE POST 2008 CRISIS WORLD

HPC IN THE POST 2008 CRISIS WORLD GTC 2016 HPC IN THE POST 2008 CRISIS WORLD Pierre SPATZ MUREX 2016 STANFORD CENTER FOR FINANCIAL AND RISK ANALYTICS HPC IN THE POST 2008 CRISIS WORLD Pierre SPATZ MUREX 2016 BACK TO 2008 FINANCIAL MARKETS

More information

MFE/3F Questions Answer Key

MFE/3F Questions Answer Key MFE/3F Questions Download free full solutions from www.actuarialbrew.com, or purchase a hard copy from www.actexmadriver.com, or www.actuarialbookstore.com. Chapter 1 Put-Call Parity and Replication 1.01

More information

Master s in Financial Engineering Foundations of Buy-Side Finance: Quantitative Risk and Portfolio Management. > Teaching > Courses

Master s in Financial Engineering Foundations of Buy-Side Finance: Quantitative Risk and Portfolio Management.  > Teaching > Courses Master s in Financial Engineering Foundations of Buy-Side Finance: Quantitative Risk and Portfolio Management www.symmys.com > Teaching > Courses Spring 2008, Monday 7:10 pm 9:30 pm, Room 303 Attilio Meucci

More information

INTEREST RATES AND FX MODELS

INTEREST RATES AND FX MODELS INTEREST RATES AND FX MODELS 7. Risk Management Andrew Lesniewski Courant Institute of Mathematical Sciences New York University New York March 8, 2012 2 Interest Rates & FX Models Contents 1 Introduction

More information

Risk Reward Optimisation for Long-Run Investors: an Empirical Analysis

Risk Reward Optimisation for Long-Run Investors: an Empirical Analysis GoBack Risk Reward Optimisation for Long-Run Investors: an Empirical Analysis M. Gilli University of Geneva and Swiss Finance Institute E. Schumann University of Geneva AFIR / LIFE Colloquium 2009 München,

More information

Risk Management. Exercises

Risk Management. Exercises Risk Management Exercises Exercise Value at Risk calculations Problem Consider a stock S valued at $1 today, which after one period can be worth S T : $2 or $0.50. Consider also a convertible bond B, which

More information

Math 623 (IOE 623), Winter 2008: Final exam

Math 623 (IOE 623), Winter 2008: Final exam Math 623 (IOE 623), Winter 2008: Final exam Name: Student ID: This is a closed book exam. You may bring up to ten one sided A4 pages of notes to the exam. You may also use a calculator but not its memory

More information

Equity correlations implied by index options: estimation and model uncertainty analysis

Equity correlations implied by index options: estimation and model uncertainty analysis 1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to

More information

SPEED UP OF NUMERIC CALCULATIONS USING A GRAPHICS PROCESSING UNIT (GPU)

SPEED UP OF NUMERIC CALCULATIONS USING A GRAPHICS PROCESSING UNIT (GPU) SPEED UP OF NUMERIC CALCULATIONS USING A GRAPHICS PROCESSING UNIT (GPU) NIKOLA VASILEV, DR. ANATOLIY ANTONOV Eurorisk Systems Ltd. 31, General Kiselov str. BG-9002 Varna, Bulgaria Phone +359 52 612 367

More information

Callable Libor exotic products. Ismail Laachir. March 1, 2012

Callable Libor exotic products. Ismail Laachir. March 1, 2012 5 pages 1 Callable Libor exotic products Ismail Laachir March 1, 2012 Contents 1 Callable Libor exotics 1 1.1 Bermudan swaption.............................. 2 1.2 Callable capped floater............................

More information

Remarks on stochastic automatic adjoint differentiation and financial models calibration

Remarks on stochastic automatic adjoint differentiation and financial models calibration arxiv:1901.04200v1 [q-fin.cp] 14 Jan 2019 Remarks on stochastic automatic adjoint differentiation and financial models calibration Dmitri Goloubentcev, Evgeny Lakshtanov Abstract In this work, we discuss

More information

Accelerating Quantitative Financial Computing with CUDA and GPUs

Accelerating Quantitative Financial Computing with CUDA and GPUs Accelerating Quantitative Financial Computing with CUDA and GPUs NVIDIA GPU Technology Conference San Jose, California Gerald A. Hanweck, Jr., PhD CEO, Hanweck Associates, LLC Hanweck Associates, LLC 30

More information

MATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, Student Name (print):

MATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, Student Name (print): MATH4143 Page 1 of 17 Winter 2007 MATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, 2007 Student Name (print): Student Signature: Student ID: Question

More information

Asset Allocation Model with Tail Risk Parity

Asset Allocation Model with Tail Risk Parity Proceedings of the Asia Pacific Industrial Engineering & Management Systems Conference 2017 Asset Allocation Model with Tail Risk Parity Hirotaka Kato Graduate School of Science and Technology Keio University,

More information

Computer Exercise 2 Simulation

Computer Exercise 2 Simulation Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Spring 2010 Computer Exercise 2 Simulation This lab deals with

More information

Statistical Models and Methods for Financial Markets

Statistical Models and Methods for Financial Markets Tze Leung Lai/ Haipeng Xing Statistical Models and Methods for Financial Markets B 374756 4Q Springer Preface \ vii Part I Basic Statistical Methods and Financial Applications 1 Linear Regression Models

More information

Optimizing Modular Expansions in an Industrial Setting Using Real Options

Optimizing Modular Expansions in an Industrial Setting Using Real Options Optimizing Modular Expansions in an Industrial Setting Using Real Options Abstract Matt Davison Yuri Lawryshyn Biyun Zhang The optimization of a modular expansion strategy, while extremely relevant in

More information

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society

More information

for Finance Python Yves Hilpisch Koln Sebastopol Tokyo O'REILLY Farnham Cambridge Beijing

for Finance Python Yves Hilpisch Koln Sebastopol Tokyo O'REILLY Farnham Cambridge Beijing Python for Finance Yves Hilpisch Beijing Cambridge Farnham Koln Sebastopol Tokyo O'REILLY Table of Contents Preface xi Part I. Python and Finance 1. Why Python for Finance? 3 What Is Python? 3 Brief History

More information

Efficient Lifetime Portfolio Sensitivities: AAD Versus Longstaff-Schwartz Compression Chris Kenyon

Efficient Lifetime Portfolio Sensitivities: AAD Versus Longstaff-Schwartz Compression Chris Kenyon Efficient Lifetime Portfolio Sensitivities: AAD Versus Longstaff-Schwartz Compression Chris Kenyon 26.03.2014 Contact: Chris.Kenyon@lloydsbanking.com Acknowledgments & Disclaimers Joint work with Andrew

More information

Efficient Reconfigurable Design for Pricing Asian Options

Efficient Reconfigurable Design for Pricing Asian Options Efficient Reconfigurable Design for Pricing Asian Options Anson H.T. Tse, David B. Thomas, K.H. Tsoi, Wayne Luk Department of Computing Imperial College London, UK {htt08,dt10,khtsoi,wl}@doc.ic.ac.uk ABSTRACT

More information

Stochastic Dual Dynamic Programming

Stochastic Dual Dynamic Programming 1 / 43 Stochastic Dual Dynamic Programming Operations Research Anthony Papavasiliou 2 / 43 Contents [ 10.4 of BL], [Pereira, 1991] 1 Recalling the Nested L-Shaped Decomposition 2 Drawbacks of Nested Decomposition

More information

Statistical Methods in Financial Risk Management

Statistical Methods in Financial Risk Management Statistical Methods in Financial Risk Management Lecture 1: Mapping Risks to Risk Factors Alexander J. McNeil Maxwell Institute of Mathematical Sciences Heriot-Watt University Edinburgh 2nd Workshop on

More information

Reinforcement Learning and Simulation-Based Search

Reinforcement Learning and Simulation-Based Search Reinforcement Learning and Simulation-Based Search David Silver Outline 1 Reinforcement Learning 2 3 Planning Under Uncertainty Reinforcement Learning Markov Decision Process Definition A Markov Decision

More information

Reverse Sensitivity Testing: What does it take to break the model? Silvana Pesenti

Reverse Sensitivity Testing: What does it take to break the model? Silvana Pesenti Reverse Sensitivity Testing: What does it take to break the model? Silvana Pesenti Silvana.Pesenti@cass.city.ac.uk joint work with Pietro Millossovich and Andreas Tsanakas Insurance Data Science Conference,

More information

MFE/3F Questions Answer Key

MFE/3F Questions Answer Key MFE/3F Questions Download free full solutions from www.actuarialbrew.com, or purchase a hard copy from www.actexmadriver.com, or www.actuarialbookstore.com. Chapter 1 Put-Call Parity and Replication 1.01

More information

A Study on Optimal Limit Order Strategy using Multi-Period Stochastic Programming considering Nonexecution Risk

A Study on Optimal Limit Order Strategy using Multi-Period Stochastic Programming considering Nonexecution Risk Proceedings of the Asia Pacific Industrial Engineering & Management Systems Conference 2018 A Study on Optimal Limit Order Strategy using Multi-Period Stochastic Programming considering Nonexecution Ris

More information

MONTE-CARLO SIMULATION CALCULATION OF VAR (VALUE-AT-RISK) & CVAR (CONDITIONAL VALUE-AT-RISK)

MONTE-CARLO SIMULATION CALCULATION OF VAR (VALUE-AT-RISK) & CVAR (CONDITIONAL VALUE-AT-RISK) MONTE-CARLO SIMULATION CALCULATION OF VAR (VALUE-AT-RISK) & CVAR (CONDITIONAL VALUE-AT-RISK) PRESENTER: SANJOY ROY 15-APR-2018 TERMINOLOGY V-a-R (Value-At-Risk) How much can one expect to lose Parameters

More information

Worst-case-expectation approach to optimization under uncertainty

Worst-case-expectation approach to optimization under uncertainty Worst-case-expectation approach to optimization under uncertainty Wajdi Tekaya Joint research with Alexander Shapiro, Murilo Pereira Soares and Joari Paulo da Costa : Cambridge Systems Associates; : Georgia

More information

Distributed Computing in Finance: Case Model Calibration

Distributed Computing in Finance: Case Model Calibration Distributed Computing in Finance: Case Model Calibration Global Derivatives Trading & Risk Management 19 May 2010 Techila Technologies, Tampere University of Technology juho.kanniainen@techila.fi juho.kanniainen@tut.fi

More information

"Pricing Exotic Options using Strong Convergence Properties

Pricing Exotic Options using Strong Convergence Properties Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike

More information

Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)

Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 2 Random number generation January 18, 2018

More information

RunnING Risk on GPUs. Answering The Computational Challenges of a New Environment. Tim Wood Market Risk Management Trading - ING Bank

RunnING Risk on GPUs. Answering The Computational Challenges of a New Environment. Tim Wood Market Risk Management Trading - ING Bank RunnING Risk on GPUs Answering The Computational Challenges of a New Environment Tim Wood Market Risk Management Trading - ING Bank Nvidia GTC Express September 19 th 2012 www.ing.com ING Bank Part of

More information

Stochastic Receding Horizon Control for Dynamic Option Hedging

Stochastic Receding Horizon Control for Dynamic Option Hedging Stochastic Receding Horizon Control for Dynamic Option Hedging Alberto Bemporad http://www.dii.unisi.it/~bemporad University of Siena Department of Information Engineering joint work with L. Bellucci and

More information

Financial Times Series. Lecture 6

Financial Times Series. Lecture 6 Financial Times Series Lecture 6 Extensions of the GARCH There are numerous extensions of the GARCH Among the more well known are EGARCH (Nelson 1991) and GJR (Glosten et al 1993) Both models allow for

More information

NUMERICAL AND SIMULATION TECHNIQUES IN FINANCE

NUMERICAL AND SIMULATION TECHNIQUES IN FINANCE NUMERICAL AND SIMULATION TECHNIQUES IN FINANCE Edward D. Weinberger, Ph.D., F.R.M Adjunct Assoc. Professor Dept. of Finance and Risk Engineering edw2026@nyu.edu Office Hours by appointment This half-semester

More information

About the Risk Quantification of Technical Systems

About the Risk Quantification of Technical Systems About the Risk Quantification of Technical Systems Magda Schiegl ASTIN Colloquium 2013, The Hague Outline Introduction / Overview Fault Tree Analysis (FTA) Method of quantitative risk analysis Results

More information

Portfolio Optimization. Prof. Daniel P. Palomar

Portfolio Optimization. Prof. Daniel P. Palomar Portfolio Optimization Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics Fall 2018-19, HKUST, Hong

More information

Many-core Accelerated LIBOR Swaption Portfolio Pricing

Many-core Accelerated LIBOR Swaption Portfolio Pricing 2012 SC Companion: High Performance Computing, Networking Storage and Analysis Many-core Accelerated LIBOR Swaption Portfolio Pricing Jörg Lotze, Paul D. Sutton, Hicham Lahlou Xcelerit Dunlop House, Fenian

More information

A new breed of Monte Carlo to meet FRTB computational challenges

A new breed of Monte Carlo to meet FRTB computational challenges A new breed of Monte Carlo to meet FRTB computational challenges 10/01/2017 Adil REGHAI Acknowledgement & Disclaimer Thanks to Abdelkrim Lajmi, Antoine Kremer, Luc Mathieu, Carole Camozzi, José Luu, Rida

More information

The risk/return trade-off has been a

The risk/return trade-off has been a Efficient Risk/Return Frontiers for Credit Risk HELMUT MAUSSER AND DAN ROSEN HELMUT MAUSSER is a mathematician at Algorithmics Inc. in Toronto, Canada. DAN ROSEN is the director of research at Algorithmics

More information

Scenario reduction and scenario tree construction for power management problems

Scenario reduction and scenario tree construction for power management problems Scenario reduction and scenario tree construction for power management problems N. Gröwe-Kuska, H. Heitsch and W. Römisch Humboldt-University Berlin Institute of Mathematics Page 1 of 20 IEEE Bologna POWER

More information

Calibrating to Market Data Getting the Model into Shape

Calibrating to Market Data Getting the Model into Shape Calibrating to Market Data Getting the Model into Shape Tutorial on Reconfigurable Architectures in Finance Tilman Sayer Department of Financial Mathematics, Fraunhofer Institute for Industrial Mathematics

More information