A Heuristic Method for Statistical Digital Circuit Sizing
|
|
- Edith Lawrence
- 6 years ago
- Views:
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? Chirayu Amin, Noel Menezes *, Kip Killpack *, Florentin Dartu *, Umakanta Choudhury *, Nagib Hakim *, Yehea Ismail ECE Department Northwestern
More informationThe 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 informationEE115C 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 informationInvestment 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 informationCPSC 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 informationCPSC 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 informationRobust 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 informationMFE/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 informationPortfolio 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 informationDynamic 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 informationLecture 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 informationEquity 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 informationOptimal 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 informationAsian 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 informationOptimal 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 informationMFE/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 informationLoss 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 informationEE266 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 informationHedging 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 informationConover 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 informationStochastic 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 informationModern 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 informationLinda 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 informationDiploma 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 informationFinancial 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 informationKARACHI 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 informationECO 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 informationINSTITUTE 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 informationEco504 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 informationDynamic 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 informationASC 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 informationEconomics 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 informationCS 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 informationOvernight 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 informationGamma. 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 informationECONOMIA 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 informationIEOR 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 informationROBUST 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 informationWhat 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 informationReverse 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 informationMODELLING 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 informationPrice 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 informationMulti-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 informationPuttable 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 informationThe 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 informationLecture 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 informationHow 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 informationSupport 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 informationThe 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 informationFinal 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 informationBivariate 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 informationStatistics 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 informationPortfolio 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 informationGraduate 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 informationA 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 information1. 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 informationAn 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 informationOptimal 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 informationHo 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 informationDecision 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 informationMoving 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 informationSAQ 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 informationLecture 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 informationHedging 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 informationDynamic 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 informationSimulation 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 informationLimit 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 informationMAFS 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 informationA 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 informationA 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 informationOptimal 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 informationSimulation. 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 information36106 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 informationVersion 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 informationFat 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 informationResults 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 informationMeasures 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 informationPortfolio 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 informationIntegrating 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 informationReview: 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 informationChapter 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 informationCHAPTER 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 informationMarket 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 informationStochastic 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 informationThe 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 informationResource 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 informationBudget 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 informationPractical 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 informationNon-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 informationForeign 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 informationSupply 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 informationValue 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 informationEfficient 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 informationContext 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 informationCallable 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 informationLog-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 informationSimple 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 informationSimulation 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 informationEstimation 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 informationEfficient 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