A Heuristic Method for Statistical Digital Circuit Sizing

Size: px
Start display at page:

Download "A Heuristic Method for Statistical Digital Circuit Sizing"

Transcription

1 A Heuristic Method for Statistical Digital Circuit Sizing Stephen Boyd Seung-Jean Kim Dinesh Patil Mark Horowitz Microlithography 06 2/23/06

2 Statistical variation in digital circuits growing in importance as devices shrink modeling still open many sources: environmental, process parameter variation, lithography intrachip, interchip variation distributions, correlations not well known, change as process matures Microlithography 06 2/23/06 1

3 Statistical digital circuit sizing standard design approaches: margining, guardbanding, design over corners statistical design explicitly takes statistical variation into account (combines circuit design with design for manufacturing, yield optimization, design centering,... ) statistical design is very hard problem (even for small circuits) this talk: a (relatively) simple heuristic method for statistical design that appears to work well Microlithography 06 2/23/06 2

4 Outline A quick example Digital circuit sizing: models and optimization New method for statistical digital circuit sizing Digital circuit sizing example Conclusions and future work Microlithography 06 2/23/06 3

5 A Quick Example

6 Example: Ladner-Fisher 32-bit adder 64 inputs, 33 outputs, 451 gates, 3214 paths, max depth 8 simplified RC delay model design variables: 451 scale factors for gates cycle time T cycle is max path delay minimize cycle time subject to limits on area, min/max scale factor Microlithography 06 2/23/06 4

7 Optimization results (no statistical variation) path delays with optimized & uniform scale factors (same total area) PSfrag replacements 4000 # of paths 2000 optimal uniform path delay Microlithography 06 2/23/06 5

8 Statistical variation in gate delay PSfrag replacements simple Pelgrom model; larger gates have less (relative) variation in delay min sized gate has 10% variation x = 5 probability x = 2 x = delay Microlithography 06 2/23/06 6

9 Effects of statistical variation on nominal optimal design T cycle PDF estimated via Monte Carlo PSfrag replacements nominal T cycle Q.95 (T cycle ) cycle time PDF cycle time Microlithography 06 2/23/06 7

10 Why isn t T cycle PDF centered around nominal value? T cycle is max of 3214 random path delays max of RVs behaves differently from sum of RVs in sum, negative and positive deviations tend to cancel out; PDF is centered, has smaller relative variation in max, large deviation of any leads to large value; PDF is shifted, skewed to right, has large relative deviation Microlithography 06 2/23/06 8

11 PSfrag replacements PDF of sum of random variables Z = M i=1 X i, X i N (1, 0.1) independent PDF M = PDF M = PDF M = Microlithography 06 2/23/06 9

12 PSfrag replacements PDF of max of random variables Z = max{x 1,..., X M }, X i N (1, 0.1) independent PDF M = 1 PDF M = 10 PDF M = Microlithography 06 2/23/06 10

13 Simple worst-case design use slow model for all gates, e.g., 1.2D i gives same design can we do better? Microlithography 06 2/23/06 11

14 Statistically robust design via new method same circuit, uncertainty model, and constraints PSfrag replacements distribution of cycle time robust design nominal optimal design cycle time Microlithography 06 2/23/06 12

15 Statistically robust design via new method nominal delay E D σ D Q.95 (D) nominal optimal robust same circuit, uncertainty model, and constraints compared to nominal optimal design, some gates are upsized, others are downsized Microlithography 06 2/23/06 13

16 Nominal vs. statistical robust designs 4000 PSfrag replacements # of paths path delay Microlithography 06 2/23/06 14

17 Path delay mean/std. dev. scatter plots path delay std. dev. path delay std. dev. PSfrag replacements nominal optimal design mean path delay robust design mean path delay Microlithography 06 2/23/06 15

18 Area/delay trade-off analysis PSfrag replacements A max 7000 robust design nominal optimal design nominal cycle time Microlithography 06 2/23/06 16

19 Area/delay trade-off analysis PSfrag replacements A max 7000 robust design nominal optimal design % cycle time Microlithography 06 2/23/06 17

20 Digital Circuit Sizing: Models

21 R A R B R C input flip flops 1 Gate scaling combinational logic block output flip flops 4 6 in 2 out C X C Y clock combinational logic; circuit topology & gate types given gate sizes (scale factors x i 1) to be determined scale factors affect total circuit area, power and delay Microlithography 06 2/23/06 18

22 PSfrag replacements R i RC gate delay model V dd C in i C in i R i C int i C L i input & intrinsic capacitances, driving resistance, load capacitance C in i = C in i x i, C int i = C int i x i, R i = R i /x i, C L i = C in j j FO(i) RC gate delay: D i = 0.69R i (C L i + C int i ) Microlithography 06 2/23/06 19

23 Path and circuit delay 1 PSfrag replacements delay of a path: sum of delays of gates on path circuit delay (cycle time): maximum delay over all paths Microlithography 06 2/23/06 20

24 Area & power total circuit area: A = x 1 Ā x n Ā n total power is P = P dyn + P stat n dynamic power P dyn = f i (Ci L i=1 f i is gate switching frequency static (leakage) power P stat = I leak i n i=1 + C int i )V 2 dd Ii leak V dd is leakage current (average over input states) Microlithography 06 2/23/06 21

25 Parameters used in example model parameters: gate type Cin Cint R Ā INV NAND NOR AOI OAI time unit is τ, delay of min-size inverter ( = 1) area (total width) unit is width of NMOS in min-size inverter Microlithography 06 2/23/06 22

26 Statistical variation in threshold voltage we focus on statistical variation in threshold voltage V th (can also model variations in other parameters, e.g., t ox, L eff,... ) Pelgrom model: σ Vth = σ Vth x 1/2 where σ 2 V th is V th variance for unit scaled gate larger gates have less V th variation Microlithography 06 2/23/06 23

27 alpha-power law model: Statistical gate delay model D V dd (V dd V th ) α (α 1.3) for small variation in V th, σ D D σ V th = α(v dd V th ) 1 σ Vth x 0.5 D V th gate scaling affects mean delay and relative variation differently relative variation decreases as gate scale factor increases: σ D /D x 0.5 Microlithography 06 2/23/06 24

28 Statistical variation in gate delay PSfrag replacements 10% relative variation for min sized gate (σ D /D = 0.1) inverter driving C L = 4 x = 5 probability x = 2 x = delay Microlithography 06 2/23/06 25

29 inverter driving C L = 4 PSfrag replacements Statistical variation in gate delay 10τ µ + 3σ delay 5τ µ µ 3σ scale factor Microlithography 06 2/23/06 26

30 Statistical leakage power model leakage current (V ) I leak xe V th/v 0 linearization does not give accurate prediction of E I leak, σ I leak exact values for V th Gaussian: E I leak = I leak,nom e σ2 V th /(2V 2 0 x), σ I leak = (e σ2 V th /(V 2 0 x) 1) 1/2 E I leak I leak,nom is leakage current when statistical variation is ignored Microlithography 06 2/23/06 27

31 Effects of statistical variation on leakage power V th N ( V th, 0.15 V th ), V th = 0.25, V 0 = 0.04 x = 1 PSfrag replacements PDF x = 2 x = 5 I leak xe V th /V 0 Microlithography 06 2/23/06 28

32 Statistical variation in leakage power 1.7 PSfrag replacements E I leak /I leak,nom scale factor Microlithography 06 2/23/06 29

33 Digital Circuit Sizing: Optimization

34 Basic gate scaling problem (no statistical variation) minimize D subject to P P max, A A max 1 x i, i = 1,..., n a geometric program (GP); can be solved efficiently extensions/variations: minimize area, power, or some combination maximize clock frequency subject to area, power limits add other constraints optimal trade-off of area, power, delay Microlithography 06 2/23/06 30

35 Statistical parameter variation now model gate delay & power as random variables circuit performance measures P, D become random variables P, D distributions of P, D depend on gate scalings x i for fixed design, can estimate PDFs of P, D via Monte Carlo PSfrag replacements frequency cycle time D Microlithography 06 2/23/06 31

36 Statistical design measure random performance measures by 95% quantile (say) minimize Q.95 (D) subject to Q.95 (P) P max, A A max 1 x i, i = 1,..., n extremely difficult stochastic optimization problem; almost no analytic/exact results but, simple heuristic method works well Microlithography 06 2/23/06 32

37 The New Method

38 Statistical power constraint total power is sum of gate powers E P = n E P i i=1 if n is large and P 1,..., P m are independent (enough), P n E P i i=1 can use E P P max as reasonable approximation of Q.95 (P) P max Microlithography 06 2/23/06 33

39 define surrogate gate delays Surrogate gate delay PSfrag replacements D i (x) = D i (x) + κ i σ i (x) κ i σ i (x) is margin on gate delay (κ i is typically 2) µ + κσ gate delay µ gives more margin to smaller gates scale factor Microlithography 06 2/23/06 34

40 Interpretation of gate delay margins margins κ i σ i (x) take statistical gate delay variation into account κ i related to Prob (D i µ i + κ i σ i ) Chebyshev inequality: Prob (D i µ i + κ i σ i ) κ2 i 1 + κ 2 i if D i is Gaussian Prob (D i µ i + κ i σ i ) = 1 2π κ i e t 2 /2 dt Microlithography 06 2/23/06 35

41 Heuristic for statistical design use modified (leakage) power model taking into account statistical variation use surrogate gate delays D i (x) = D i (x) + κ i σ i (x) now solve resulting (deterministic) gate scaling problem verify statistical performance via Monte Carlo analysis (can update κ i s and repeat) Microlithography 06 2/23/06 36

42 Digital Circuit Sizing Example

43 Statistically robust design via new method same circuit, uncertainty model, and constraints PSfrag replacements distribution of cycle time robust design nominal optimal design cycle time Microlithography 06 2/23/06 37

44 Path delay mean/std. dev. scatter plots path delay std. dev. path delay std. dev. PSfrag replacements nominal optimal design mean path delay robust design mean path delay Microlithography 06 2/23/06 38

45 PSfrag replacements Comparison of nominal optimal and robust designs # of gates 400 nominal optimal design # of gates robust design scale factor Microlithography 06 2/23/06 39

46 Comparison of nominal optimal and robust designs PSfrag replacements 32 scale factor (robust design) scale factor (nominal optimal design) Microlithography 06 2/23/06 40

47 Effect of margin coefficients 52 PSfrag replacements 95% cycle time κ Microlithography 06 2/23/06 41

48 Sensitivity to model assumptions question: how sensitive is robust design to our model of process variation? distribution shape correlation between gates Pelgrom model of variance vs. scale factor answer: not very Microlithography 06 2/23/06 42

49 PSfrag replacements Simulation with uniform gate delay distributions distribution of cycle time robust design nominal optimal design cycle time compared with Gaussian gate delays: nominal optimal design not quite as bad; robust design still quite good Microlithography 06 2/23/06 43

50 PSfrag replacements Simulation with correlated gate delays connected gates have delays that are 30% correlated robust design distribution of cycle time nominal optimal design cycle time nominal optimal not as bad; but robust design still quite good Microlithography 06 2/23/06 44

51 Conclusions and Future Work

52 Conclusions statistically robust design is subtle; cannot be done by hand exact or direct methods will not work well computationally intractable depend on details of statistical models heuristic method is relatively simple, scales well, gives good designs reduces problem to a deterministic one Microlithography 06 2/23/06 45

53 References Boyd, Kim, and Mohan, DATE Tutorial 2005 Geometric programming and its applications to EDA problems Boyd, Kim, Patil, and Horowitz, SPIE ML 2006 A heuristic method for statistical digital circuit sizing Kim, Boyd, Patil, and Horowitz, Optimization and Engineering, 2006 A heuristic for optimizing stochastic activity networks with applications to statistical digital circuit sizing Boyd, Kim, Patil, and Horowitz, Operations Research, 2005 Digital circuit optimization via geometric programming Patil, Yun, Kim, Cheung, Horowitz, and Boyd, ISQED 2005 A new method for design of robust digital circuits all available from boyd/research.html Microlithography 06 2/23/06 46

54 References (continued) Mani, Devgan, Orshansky, DAC 2005 An efficient algorithm for statistical minimization of total power under timing yield constraints Satish, Ravindran, Moskewicz, Chinnery, and Keutzer, UCB tech. report, 2005 Evaluating the effectiveness of statistical gate sizing for power optimization Bhardwaj and Vrudhula, DAC 2005 Leakage minimization of nano-scale circuits in the presence of systematic and random variations Microlithography 06 2/23/06 47

Statistical Static Timing Analysis: How simple can we get?

Statistical Static Timing Analysis: How simple can we get? Statistical Static Timing Analysis: How simple can we get? Chirayu Amin, Noel Menezes *, Kip Killpack *, Florentin Dartu *, Umakanta Choudhury *, Nagib Hakim *, Yehea Ismail ECE Department Northwestern

More information

The good, the bad and the statistical

The good, the bad and the statistical The good, the bad and the statistical Noel Menezes Strategic CAD Labs Design and Technology Solutions Intel Corp. Acknowledgements Keith Bowman Yossi Abulafia Steve Burns Mahesh Ketkar Vivek De Jim Tschanz

More information

EE115C Spring 2013 Digital Electronic Circuits. Lecture 19: Timing Analysis

EE115C Spring 2013 Digital Electronic Circuits. Lecture 19: Timing Analysis EE115C Spring 2013 Digital Electronic Circuits Lecture 19: Timing Analysis Outline Timing parameters Clock nonidealities (skew and jitter) Impact of Clk skew on timing Impact of Clk jitter on timing Flip-flop-

More information

Investment strategies and risk management for participating life insurance contracts

Investment strategies and risk management for participating life insurance contracts 1/20 Investment strategies and risk for participating life insurance contracts and Steven Haberman Cass Business School AFIR Colloquium Munich, September 2009 2/20 & Motivation Motivation New supervisory

More information

CPSC 540: Machine Learning

CPSC 540: Machine Learning CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}

More information

CPSC 540: Machine Learning

CPSC 540: Machine Learning CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x

More information

Robust Optimization Applied to a Currency Portfolio

Robust Optimization Applied to a Currency Portfolio Robust Optimization Applied to a Currency Portfolio R. Fonseca, S. Zymler, W. Wiesemann, B. Rustem Workshop on Numerical Methods and Optimization in Finance June, 2009 OUTLINE Introduction Motivation &

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

Portfolio Management and Optimal Execution via Convex Optimization

Portfolio Management and Optimal Execution via Convex Optimization Portfolio Management and Optimal Execution via Convex Optimization Enzo Busseti Stanford University April 9th, 2018 Problems portfolio management choose trades with optimization minimize risk, maximize

More information

Dynamic Replication of Non-Maturing Assets and Liabilities

Dynamic Replication of Non-Maturing Assets and Liabilities Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland

More information

Lecture 8: Skew Tolerant Design (including Dynamic Circuit Issues)

Lecture 8: Skew Tolerant Design (including Dynamic Circuit Issues) Lecture 8: Skew Tolerant Design (including Dynamic Circuit Issues) Computer Systems Laboratory Stanford University horowitz@stanford.edu Copyright 2007 by Mark Horowitz w/ material from David Harris 1

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

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

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

Optimal Portfolio Selection Under the Estimation Risk in Mean Return

Optimal Portfolio Selection Under the Estimation Risk in Mean Return Optimal Portfolio Selection Under the Estimation Risk in Mean Return by Lei Zhu A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master of Mathematics

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

Loss Simulation Model Testing and Enhancement

Loss Simulation Model Testing and Enhancement Loss Simulation Model Testing and Enhancement Casualty Loss Reserve Seminar By Kailan Shang Sept. 2011 Agenda Research Overview Model Testing Real Data Model Enhancement Further Development Enterprise

More information

EE266 Homework 5 Solutions

EE266 Homework 5 Solutions EE, Spring 15-1 Professor S. Lall EE Homework 5 Solutions 1. A refined inventory model. In this problem we consider an inventory model that is more refined than the one you ve seen in the lectures. The

More information

Hedging with Life and General Insurance Products

Hedging with Life and General Insurance Products Hedging with Life and General Insurance Products June 2016 2 Hedging with Life and General Insurance Products Jungmin Choi Department of Mathematics East Carolina University Abstract In this study, a hybrid

More information

Conover Test of Variances (Simulation)

Conover Test of Variances (Simulation) Chapter 561 Conover Test of Variances (Simulation) Introduction This procedure analyzes the power and significance level of the Conover homogeneity test. This test is used to test whether two or more population

More information

Stochastic Optimization Methods in Scheduling. Rolf H. Möhring Technische Universität Berlin Combinatorial Optimization and Graph Algorithms

Stochastic Optimization Methods in Scheduling. Rolf H. Möhring Technische Universität Berlin Combinatorial Optimization and Graph Algorithms Stochastic Optimization Methods in Scheduling Rolf H. Möhring Technische Universität Berlin Combinatorial Optimization and Graph Algorithms More expensive and longer... Eurotunnel Unexpected loss of 400,000,000

More information

Modern Methods of Data Analysis - SS 2009

Modern Methods of Data Analysis - SS 2009 Modern Methods of Data Analysis Lecture II (7.04.09) Contents: Characterize data samples Characterize distributions Correlations, covariance Reminder: Average of a Sample arithmetic mean of data set: weighted

More information

Linda Allen, Jacob Boudoukh and Anthony Saunders, Understanding Market, Credit and Operational Risk: The Value at Risk Approach

Linda Allen, Jacob Boudoukh and Anthony Saunders, Understanding Market, Credit and Operational Risk: The Value at Risk Approach P1.T4. Valuation & Risk Models Linda Allen, Jacob Boudoukh and Anthony Saunders, Understanding Market, Credit and Operational Risk: The Value at Risk Approach Bionic Turtle FRM Study Notes Reading 26 By

More information

Diploma in Business Administration Part 2. Quantitative Methods. Examiner s Suggested Answers

Diploma in Business Administration Part 2. Quantitative Methods. Examiner s Suggested Answers Cumulative frequency Diploma in Business Administration Part Quantitative Methods Examiner s Suggested Answers Question 1 Cumulative Frequency Curve 1 9 8 7 6 5 4 3 1 5 1 15 5 3 35 4 45 Weeks 1 (b) x f

More information

Financial Risk Management

Financial Risk Management Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given

More information

KARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI

KARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI 88 P a g e B S ( B B A ) S y l l a b u s KARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI Course Title : STATISTICS Course Number : BA(BS) 532 Credit Hours : 03 Course 1. Statistical

More information

ECO 317 Economics of Uncertainty Fall Term 2009 Tuesday October 6 Portfolio Allocation Mean-Variance Approach

ECO 317 Economics of Uncertainty Fall Term 2009 Tuesday October 6 Portfolio Allocation Mean-Variance Approach ECO 317 Economics of Uncertainty Fall Term 2009 Tuesday October 6 ortfolio Allocation Mean-Variance Approach Validity of the Mean-Variance Approach Constant absolute risk aversion (CARA): u(w ) = exp(

More information

INSTITUTE AND FACULTY OF ACTUARIES SUMMARY

INSTITUTE AND FACULTY OF ACTUARIES SUMMARY INSTITUTE AND FACULTY OF ACTUARIES SUMMARY Specimen 2019 CP2: Actuarial Modelling Paper 2 Institute and Faculty of Actuaries TQIC Reinsurance Renewal Objective The objective of this project is to use random

More information

Eco504 Spring 2010 C. Sims FINAL EXAM. β t 1 2 φτ2 t subject to (1)

Eco504 Spring 2010 C. Sims FINAL EXAM. β t 1 2 φτ2 t subject to (1) Eco54 Spring 21 C. Sims FINAL EXAM There are three questions that will be equally weighted in grading. Since you may find some questions take longer to answer than others, and partial credit will be given

More information

Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods

Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods ISOPE 2010 Conference Beijing, China 24 June 2010 Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods Xi Ying Zhang, Zhi Ping Cheng, Jer-Fang Wu and Chee Chow Kei ABS 1 Main Contents

More information

ASC Topic 718 Accounting Valuation Report. Company ABC, Inc.

ASC Topic 718 Accounting Valuation Report. Company ABC, Inc. ASC Topic 718 Accounting Valuation Report Company ABC, Inc. Monte-Carlo Simulation Valuation of Several Proposed Relative Total Shareholder Return TSR Component Rank Grants And Index Outperform Grants

More information

Economics 883: The Basic Diffusive Model, Jumps, Variance Measures. George Tauchen. Economics 883FS Spring 2015

Economics 883: The Basic Diffusive Model, Jumps, Variance Measures. George Tauchen. Economics 883FS Spring 2015 Economics 883: The Basic Diffusive Model, Jumps, Variance Measures George Tauchen Economics 883FS Spring 2015 Main Points 1. The Continuous Time Model, Theory and Simulation 2. Observed Data, Plotting

More information

CS 774 Project: Fall 2009 Version: November 27, 2009

CS 774 Project: Fall 2009 Version: November 27, 2009 CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS

More information

Overnight Index Rate: Model, calibration and simulation

Overnight Index Rate: Model, calibration and simulation Research Article Overnight Index Rate: Model, calibration and simulation Olga Yashkir and Yuri Yashkir Cogent Economics & Finance (2014), 2: 936955 Page 1 of 11 Research Article Overnight Index Rate: Model,

More information

Gamma. The finite-difference formula for gamma is

Gamma. The finite-difference formula for gamma is Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S) + P (S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gammas

More information

ECONOMIA DEGLI INTERMEDIARI FINANZIARI AVANZATA MODULO ASSET MANAGEMENT LECTURE 6

ECONOMIA DEGLI INTERMEDIARI FINANZIARI AVANZATA MODULO ASSET MANAGEMENT LECTURE 6 ECONOMIA DEGLI INTERMEDIARI FINANZIARI AVANZATA MODULO ASSET MANAGEMENT LECTURE 6 MVO IN TWO STAGES Calculate the forecasts Calculate forecasts for returns, standard deviations and correlations for the

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

ROBUST OPTIMIZATION OF MULTI-PERIOD PRODUCTION PLANNING UNDER DEMAND UNCERTAINTY. A. Ben-Tal, B. Golany and M. Rozenblit

ROBUST OPTIMIZATION OF MULTI-PERIOD PRODUCTION PLANNING UNDER DEMAND UNCERTAINTY. A. Ben-Tal, B. Golany and M. Rozenblit ROBUST OPTIMIZATION OF MULTI-PERIOD PRODUCTION PLANNING UNDER DEMAND UNCERTAINTY A. Ben-Tal, B. Golany and M. Rozenblit Faculty of Industrial Engineering and Management, Technion, Haifa 32000, Israel ABSTRACT

More information

What is Cyclical in Credit Cycles?

What is Cyclical in Credit Cycles? What is Cyclical in Credit Cycles? Rui Cui May 31, 2014 Introduction Credit cycles are growth cycles Cyclicality in the amount of new credit Explanations: collateral constraints, equity constraints, leverage

More information

Reverse Engineering and Prevention Techniques for Physical Unclonable Functions Using Side Channels

Reverse Engineering and Prevention Techniques for Physical Unclonable Functions Using Side Channels Reverse Engineering and Prevention Techniques for Physical Unclonable Functions Using Side Channels Sheng Wei * James B. Wendt * Ani Nahapetian * Miodrag Potkonjak * * University of California, Los Angeles

More information

MODELLING OPTIMAL HEDGE RATIO IN THE PRESENCE OF FUNDING RISK

MODELLING OPTIMAL HEDGE RATIO IN THE PRESENCE OF FUNDING RISK MODELLING OPTIMAL HEDGE RATIO IN THE PRESENCE O UNDING RISK Barbara Dömötör Department of inance Corvinus University of Budapest 193, Budapest, Hungary E-mail: barbara.domotor@uni-corvinus.hu KEYWORDS

More information

Price Impact and Optimal Execution Strategy

Price Impact and Optimal Execution Strategy OXFORD MAN INSTITUE, UNIVERSITY OF OXFORD SUMMER RESEARCH PROJECT Price Impact and Optimal Execution Strategy Bingqing Liu Supervised by Stephen Roberts and Dieter Hendricks Abstract Price impact refers

More information

Multi-Period Trading via Convex Optimization

Multi-Period Trading via Convex Optimization Multi-Period Trading via Convex Optimization Stephen Boyd Enzo Busseti Steven Diamond Ronald Kahn Kwangmoo Koh Peter Nystrup Jan Speth Stanford University & Blackrock City University of Hong Kong September

More information

Puttable Bond and Vaulation

Puttable Bond and Vaulation and Vaulation Dmitry Popov FinPricing http://www.finpricing.com Summary Puttable Bond Definition The Advantages of Puttable Bonds Puttable Bond Payoffs Valuation Model Selection Criteria LGM Model LGM

More information

The Optimization Process: An example of portfolio optimization

The Optimization Process: An example of portfolio optimization ISyE 6669: Deterministic Optimization The Optimization Process: An example of portfolio optimization Shabbir Ahmed Fall 2002 1 Introduction Optimization can be roughly defined as a quantitative approach

More information

Lecture 8: Skew Tolerant Domino Clocking

Lecture 8: Skew Tolerant Domino Clocking Lecture 8: Skew Tolerant Domino Clocking Computer Systems Laboratory Stanford University horowitz@stanford.edu Copyright 2001 by Mark Horowitz (Original Slides from David Harris) 1 Introduction Domino

More information

How to Consider Risk Demystifying Monte Carlo Risk Analysis

How to Consider Risk Demystifying Monte Carlo Risk Analysis How to Consider Risk Demystifying Monte Carlo Risk Analysis James W. Richardson Regents Professor Senior Faculty Fellow Co-Director, Agricultural and Food Policy Center Department of Agricultural Economics

More information

Support Vector Machines: Training with Stochastic Gradient Descent

Support Vector Machines: Training with Stochastic Gradient Descent Support Vector Machines: Training with Stochastic Gradient Descent Machine Learning Spring 2018 The slides are mainly from Vivek Srikumar 1 Support vector machines Training by maximizing margin The SVM

More information

The Role of ERM in Reinsurance Decisions

The Role of ERM in Reinsurance Decisions The Role of ERM in Reinsurance Decisions Abbe S. Bensimon, FCAS, MAAA ERM Symposium Chicago, March 29, 2007 1 Agenda A Different Framework for Reinsurance Decision-Making An ERM Approach for Reinsurance

More information

Final Exam Suggested Solutions

Final Exam Suggested Solutions University of Washington Fall 003 Department of Economics Eric Zivot Economics 483 Final Exam Suggested Solutions This is a closed book and closed note exam. However, you are allowed one page of handwritten

More information

Bivariate Birnbaum-Saunders Distribution

Bivariate Birnbaum-Saunders Distribution Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators

More information

Statistics 431 Spring 2007 P. Shaman. Preliminaries

Statistics 431 Spring 2007 P. Shaman. Preliminaries Statistics 4 Spring 007 P. Shaman The Binomial Distribution Preliminaries A binomial experiment is defined by the following conditions: A sequence of n trials is conducted, with each trial having two possible

More information

Portfolio Risk Management and Linear Factor Models

Portfolio Risk Management and Linear Factor Models Chapter 9 Portfolio Risk Management and Linear Factor Models 9.1 Portfolio Risk Measures There are many quantities introduced over the years to measure the level of risk that a portfolio carries, and each

More information

Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Solutions to Final Exam

Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Solutions to Final Exam Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (30 pts) Answer briefly the following questions. 1. Suppose that

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

1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and

1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and CHAPTER 13 Solutions Exercise 1 1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and (13.82) (13.86). Also, remember that BDT model will yield a recombining binomial

More information

An Empirical Examination of the Electric Utilities Industry. December 19, Regulatory Induced Risk Aversion in. Contracting Behavior

An Empirical Examination of the Electric Utilities Industry. December 19, Regulatory Induced Risk Aversion in. Contracting Behavior An Empirical Examination of the Electric Utilities Industry December 19, 2011 The Puzzle Why do price-regulated firms purchase input coal through both contract Figure and 1(a): spot Contract transactions,

More information

Optimal Trade Execution: Mean Variance or Mean Quadratic Variation?

Optimal Trade Execution: Mean Variance or Mean Quadratic Variation? Optimal Trade Execution: Mean Variance or Mean Quadratic Variation? Peter Forsyth 1 S. Tse 2 H. Windcliff 2 S. Kennedy 2 1 Cheriton School of Computer Science University of Waterloo 2 Morgan Stanley New

More information

Ho Ho Quantitative Portfolio Manager, CalPERS

Ho Ho Quantitative Portfolio Manager, CalPERS Portfolio Construction and Risk Management under Non-Normality Fiduciary Investors Symposium, Beijing - China October 23 rd 26 th, 2011 Ho Ho Quantitative Portfolio Manager, CalPERS The views expressed

More information

Decision Support Models 2012/2013

Decision Support Models 2012/2013 Risk Analysis Decision Support Models 2012/2013 Bibliography: Goodwin, P. and Wright, G. (2003) Decision Analysis for Management Judgment, John Wiley and Sons (chapter 7) Clemen, R.T. and Reilly, T. (2003).

More information

Moving PUFs out of the lab

Moving PUFs out of the lab Moving PUFs out of the lab Patrick Schaumont 2/3/2012 Research results by Abhranil Maiti, Jeff Casarona, Luke McHale, Logan McDougall, Vikash Gunreddy, Michael Cantrell What is a Physical Unclonable Function?

More information

SAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax:

SAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax: ProSINTAP - A Probabilistic Program for Safety Evaluation Peter Dillström SAQ / SINTAP / 09 SAQ KONTROLL AB Box 49306, 100 29 STOCKHOLM, Sweden Tel: +46 8 617 40 00; Fax: +46 8 651 70 43 June 1999 Page

More information

Lecture 17: More on Markov Decision Processes. Reinforcement learning

Lecture 17: More on Markov Decision Processes. Reinforcement learning Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture

More information

Hedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo

Hedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor

More information

Dynamic Portfolio Choice II

Dynamic Portfolio Choice II Dynamic Portfolio Choice II Dynamic Programming Leonid Kogan MIT, Sloan 15.450, Fall 2010 c Leonid Kogan ( MIT, Sloan ) Dynamic Portfolio Choice II 15.450, Fall 2010 1 / 35 Outline 1 Introduction to Dynamic

More information

Simulation of delta hedging of an option with volume uncertainty. Marc LE DU, Clémence ALASSEUR EDF R&D - OSIRIS

Simulation of delta hedging of an option with volume uncertainty. Marc LE DU, Clémence ALASSEUR EDF R&D - OSIRIS Simulation of delta hedging of an option with volume uncertainty Marc LE DU, Clémence ALASSEUR EDF R&D - OSIRIS Agenda 1. Introduction : volume uncertainty 2. Test description: a simple option 3. Results

More information

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation

More information

MAFS Computational Methods for Pricing Structured Products

MAFS Computational Methods for Pricing Structured Products MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 )

More information

A Stochastic Reserving Today (Beyond Bootstrap)

A Stochastic Reserving Today (Beyond Bootstrap) A Stochastic Reserving Today (Beyond Bootstrap) Presented by Roger M. Hayne, PhD., FCAS, MAAA Casualty Loss Reserve Seminar 6-7 September 2012 Denver, CO CAS Antitrust Notice The Casualty Actuarial Society

More information

A Physical Unclonable Function based on Capacitor Mismatch in a Charge-Redistribution SAR-ADC

A Physical Unclonable Function based on Capacitor Mismatch in a Charge-Redistribution SAR-ADC A Physical Unclonable Function based on Capacitor Mismatch in a Charge-Redistribution SAR-ADC Qianying Tang, Won Ho Choi, Luke Everson, Keshab K. Parhi and Chris H. Kim University of Minnesota Department

More information

Optimal Stochastic Recovery for Base Correlation

Optimal Stochastic Recovery for Base Correlation Optimal Stochastic Recovery for Base Correlation Salah AMRAOUI - Sebastien HITIER BNP PARIBAS June-2008 Abstract On the back of monoline protection unwind and positive gamma hunting, spreads of the senior

More information

Simulation. Decision Models

Simulation. Decision Models Lecture 9 Decision Models Decision Models: Lecture 9 2 Simulation What is Monte Carlo simulation? A model that mimics the behavior of a (stochastic) system Mathematically described the system using a set

More information

36106 Managerial Decision Modeling Monte Carlo Simulation in Excel: Part IV

36106 Managerial Decision Modeling Monte Carlo Simulation in Excel: Part IV 36106 Managerial Decision Modeling Monte Carlo Simulation in Excel: Part IV Kipp Martin University of Chicago Booth School of Business November 29, 2017 Reading and Excel Files 2 Reading: Handout: Optimal

More information

Version A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.

Version A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise. Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x

More information

Fat tails and 4th Moments: Practical Problems of Variance Estimation

Fat tails and 4th Moments: Practical Problems of Variance Estimation Fat tails and 4th Moments: Practical Problems of Variance Estimation Blake LeBaron International Business School Brandeis University www.brandeis.edu/~blebaron QWAFAFEW May 2006 Asset Returns and Fat Tails

More information

Results for option pricing

Results for option pricing Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform

More information

Measures of Center. Mean. 1. Mean 2. Median 3. Mode 4. Midrange (rarely used) Measure of Center. Notation. Mean

Measures of Center. Mean. 1. Mean 2. Median 3. Mode 4. Midrange (rarely used) Measure of Center. Notation. Mean Measure of Center Measures of Center The value at the center or middle of a data set 1. Mean 2. Median 3. Mode 4. Midrange (rarely used) 1 2 Mean Notation The measure of center obtained by adding the values

More information

Portfolio Management Under Epistemic Uncertainty Using Stochastic Dominance and Information-Gap Theory

Portfolio Management Under Epistemic Uncertainty Using Stochastic Dominance and Information-Gap Theory Portfolio Management Under Epistemic Uncertainty Using Stochastic Dominance and Information-Gap Theory D. Berleant, L. Andrieu, J.-P. Argaud, F. Barjon, M.-P. Cheong, M. Dancre, G. Sheble, and C.-C. Teoh

More information

Integrating Reserve Variability and ERM:

Integrating Reserve Variability and ERM: Integrating Reserve Variability and ERM: Mark R. Shapland, FCAS, FSA, MAAA Jeffrey A. Courchene, FCAS, MAAA International Congress of Actuaries 30 March 4 April 2014 Washington, DC What are the Issues?

More information

Review: Population, sample, and sampling distributions

Review: Population, sample, and sampling distributions Review: Population, sample, and sampling distributions A population with mean µ and standard deviation σ For instance, µ = 0, σ = 1 0 1 Sample 1, N=30 Sample 2, N=30 Sample 100000000000 InterquartileRange

More information

Chapter 7. Confidence Intervals and Sample Sizes. Definition. Definition. Definition. Definition. Confidence Interval : CI. Point Estimate.

Chapter 7. Confidence Intervals and Sample Sizes. Definition. Definition. Definition. Definition. Confidence Interval : CI. Point Estimate. Chapter 7 Confidence Intervals and Sample Sizes 7. Estimating a Proportion p 7.3 Estimating a Mean µ (σ known) 7.4 Estimating a Mean µ (σ unknown) 7.5 Estimating a Standard Deviation σ In a recent poll,

More information

CHAPTER 8. Confidence Interval Estimation Point and Interval Estimates

CHAPTER 8. Confidence Interval Estimation Point and Interval Estimates CHAPTER 8. Confidence Interval Estimation Point and Interval Estimates A point estimate is a single number, a confidence interval provides additional information about the variability of the estimate Lower

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

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

The histogram should resemble the uniform density, the mean should be close to 0.5, and the standard deviation should be close to 1/ 12 =

The histogram should resemble the uniform density, the mean should be close to 0.5, and the standard deviation should be close to 1/ 12 = Chapter 19 Monte Carlo Valuation Question 19.1 The histogram should resemble the uniform density, the mean should be close to.5, and the standard deviation should be close to 1/ 1 =.887. Question 19. The

More information

Resource Planning with Uncertainty for NorthWestern Energy

Resource Planning with Uncertainty for NorthWestern Energy Resource Planning with Uncertainty for NorthWestern Energy Selection of Optimal Resource Plan for 213 Resource Procurement Plan August 28, 213 Gary Dorris, Ph.D. Ascend Analytics, LLC gdorris@ascendanalytics.com

More information

Budget Management In GSP (2018)

Budget Management In GSP (2018) Budget Management In GSP (2018) Yahoo! March 18, 2018 Miguel March 18, 2018 1 / 26 Today s Presentation: Budget Management Strategies in Repeated auctions, Balseiro, Kim, and Mahdian, WWW2017 Learning

More information

Practical example of an Economic Scenario Generator

Practical example of an Economic Scenario Generator Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application

More information

Non-Deterministic Search

Non-Deterministic Search Non-Deterministic Search MDP s 1 Non-Deterministic Search How do you plan (search) when your actions might fail? In general case, how do you plan, when the actions have multiple possible outcomes? 2 Example:

More information

Foreign Competition and Banking Industry Dynamics: An Application to Mexico

Foreign Competition and Banking Industry Dynamics: An Application to Mexico Foreign Competition and Banking Industry Dynamics: An Application to Mexico Dean Corbae Pablo D Erasmo 1 Univ. of Wisconsin FRB Philadelphia June 12, 2014 1 The views expressed here do not necessarily

More information

Supply Contracts with Financial Hedging

Supply Contracts with Financial Hedging Supply Contracts with Financial Hedging René Caldentey Martin Haugh Stern School of Business NYU Integrated Risk Management in Operations and Global Supply Chain Management: Risk, Contracts and Insurance

More information

Value at Risk Ch.12. PAK Study Manual

Value at Risk Ch.12. PAK Study Manual Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and

More information

Efficient Rebalancing of Taxable Portfolios

Efficient Rebalancing of Taxable Portfolios Efficient Rebalancing of Taxable Portfolios Sanjiv R. Das & Daniel Ostrov 1 Santa Clara University @JOIM La Jolla, CA April 2015 1 Joint work with Dennis Yi Ding and Vincent Newell. Das and Ostrov (Santa

More information

Context Power analyses for logistic regression models fit to clustered data

Context Power analyses for logistic regression models fit to clustered data . Power Analysis for Logistic Regression Models Fit to Clustered Data: Choosing the Right Rho. CAPS Methods Core Seminar Steve Gregorich May 16, 2014 CAPS Methods Core 1 SGregorich Abstract Context Power

More information

Callable Bond and Vaulation

Callable Bond and Vaulation and Vaulation Dmitry Popov FinPricing http://www.finpricing.com Summary Callable Bond Definition The Advantages of Callable Bonds Callable Bond Payoffs Valuation Model Selection Criteria LGM Model LGM

More information

Log-Robust Portfolio Management

Log-Robust Portfolio Management Log-Robust Portfolio Management Dr. Aurélie Thiele Lehigh University Joint work with Elcin Cetinkaya and Ban Kawas Research partially supported by the National Science Foundation Grant CMMI-0757983 Dr.

More information

Simple Robust Hedging with Nearby Contracts

Simple Robust Hedging with Nearby Contracts Simple Robust Hedging with Nearby Contracts Liuren Wu and Jingyi Zhu Baruch College and University of Utah October 22, 2 at Worcester Polytechnic Institute Wu & Zhu (Baruch & Utah) Robust Hedging with

More information

Simulation Wrap-up, Statistics COS 323

Simulation Wrap-up, Statistics COS 323 Simulation Wrap-up, Statistics COS 323 Today Simulation Re-cap Statistics Variance and confidence intervals for simulations Simulation wrap-up FYI: No class or office hours Thursday Simulation wrap-up

More information

Estimation of dynamic term structure models

Estimation of dynamic term structure models Estimation of dynamic term structure models Greg Duffee Haas School of Business, UC-Berkeley Joint with Richard Stanton, Haas School Presentation at IMA Workshop, May 2004 (full paper at http://faculty.haas.berkeley.edu/duffee)

More information

Efficient Valuation of Large Variable Annuity Portfolios

Efficient Valuation of Large Variable Annuity Portfolios Efficient Valuation of Large Variable Annuity Portfolios Emiliano A. Valdez joint work with Guojun Gan University of Connecticut Seminar Talk at Wisconsin School of Business University of Wisconsin Madison,

More information