escript basics Lutz Gross School of Earth Sciences The University of Queensland 1/46 13/03/13
|
|
- Byron Johns
- 5 years ago
- Views:
Transcription
1 escript basics Lutz Gross School of Earth Sciences The University of Queensland 1/46 13/03/13
2 escript Sponsors AuScope National Collaborative Research Infrastructure Strategy Australian Geophysical Observing System (AGOS) (EIF) School of Earth Sciences at the University of Queensland. 2/46 13/03/13
3 Escript? Software for solving partial differential equations Coupled, time-dependent, non-linear Easy to use uses PDE terminology programming in python based on the finite element method parallelized for multi-cores and clusters 3/46 13/03/13
4 Porous Media Flow: Perth Basin 4/46 13/03/13
5 Stress Localization 5/46 13/03/13
6 Mantel Convection 6/46 13/03/13
7 How to get escript Distributed via launchpad: visit FAQ current version some binary versions available Documentation Cookbook (for beginner) User's guide (the bible) Inversion cookbook Reference guide 7/46 13/03/13
8 Python? Programming language Interactive or scripts Object oriented Open source Standard in Linux Tons of useful modules numpy: linear algebra SciPy, obspy, matplotlib, sympy Tutorial: 8/46 13/03/13
9 Python example $ python Python (#62, Sep , 19:20:46) [MSC v bit (Intel)] on win32 Type "help", "copyright", "credits" or "license" for more information. >>> print `Hello World` Hello World >> x=6 >> print x+9 15 >> s=`hello`+` `+`World` # concatenate strings >> print s Hello World >> q=[ 1, `H`, x, s ] # a list of objects >> print q[2] 6 >> print q[2:] [6, 'Hello World'] >> CRTL^D # get out of here $ 9/46 13/03/13
10 escript Data Handling: esys.escript PDEs: esys.escript.linearpdes Finite Element solver: Unstructured: esys.finley Grids: esys.ripley Inversion: esys.downunder Data port: esys.weipa Geometries: esys.pycad 10/46 13/03/13
11 Problem: Temperature Distribution T=0 (l 0,l 1 ) heat flux=0 heat flux=0 x 1 (0,0) T=T bot heat source x 0 11/46 13/03/13
12 1D Temperature Diffusion x 0 ( κ T x 0 ) =Q T temperature κ thermal conductivity (may depend on x) Q heat source 12/46 13/03/13
13 Escript's Notation for derivatives With Z, 0 = Z x 0 (κ T,0 ), 0 =Q 13/46 13/03/13
14 2D Temperature Diffusion ( κ T )=Q or x 0 ( κ T x 0 ) x 1 ( κ T x 1 ) =Q 14/46 13/03/13
15 2D Temperature Diffusion (cont.) x 0 ( κ T x 0 ) x 1 ( κ T x 1 ) =Q or (κ T,0 ), 0 (κ T, 1 ),1 =Q 15/46 13/03/13
16 Insulation Boundary Condition heat flux=0 heat flux=0 boundary heat flux=κ T n =κ n T =κ n T i, i (summation over i) 16/46 13/03/13
17 Flux Boundary Conditions (l 0,l 1 ) outer normal field: n= (n i ) (0,0) Ω n Apply at face x 0 =0 and x 0 =l 0 heat flux=κ n i T, i =0 Called Neuman or natural boundary condition 17/46 13/03/13
18 Temperature Boundary Condition T=0 T=T bot (known value) Also called Dirichlet boundary condition or constraint 18/46 13/03/13
19 How to define a Constraint T=T D =0 appropriate T D over entire domain (l 0,l 1 ) (0,0) Ω T=T D =T bot Constraint is given as: T=T D at x 1 =0 and x 1 =l 1 For instance : T D (x 0, x 1 )= T bot l 1 (l 1 x 1 ) 19/46 13/03/13
20 PDE for Temperature Diffusion (κ T,i ), i =Q κ n i T,i =0 at faces x 0 =0 and x 0 =l 0 T =T D at faces x 1 =0 and x 1 =l 1 (summation over i) 20/46 13/03/13
21 Domain in escript - region where to solve the PDE - a discretization method (0,l 1 ) (l 0,l 1 ) (0,0) (l 0,0) 21/46 13/03/13
22 Finite Element Mesh in finley Node Element Quadrature Point Temperature by values on Nodes Gradients by value on quadrature points. 22/46 13/03/13
23 Create a finley domain # get the tools we want to use from esys.escript import * from esys.finley import Rectangle # generate n0 x n1 elements over # [0,l0] x [0,l1] L0=10. L1=5. mydomain=rectangle(l0=l0,l1=l1,n0=80,n1=20) print dimension =,mydomain.getdim() x=mydomain.getx() print x # coordinates of FEM nodes as list or summary 23/46 13/03/13
24 How do I run this? 1. create a file myprog.py, eg. >> gedit myprog.py & 2. enter statements into file and save 3. run python >> run-escript myprog.py 24/46 13/03/13
25 Write to file from esys.escript import * from esys.finley import Rectangle from esys.weipa import * L0=10. L1=5. T_bot=100. mydomain=rectangle(l0=l0,l1=l1,n0=80,n1=40) x=mydomain.getx() T_D=T_bot/L1*(L1-x[1]) # save T_D for visualization with a VTK tool savevtk('u.vtu',t=t_d) 25/46 13/03/13
26 Run & Visualization >> run-escript myprog.py >> ls myprog.py u.vtu >> visit & # fire up visualization Or any tool supporting VTK files: VisIt, paraview, mayavi2 26/46 13/03/13
27 Visualization with VisIt Click Select data file Click 27/46 13/03/13
28 VisIt (cont.) 1. Click 3. Click to draw 2. Select data set 'T' 28/46 13/03/13
29 Visualization with VisIt (cont.) 29/46 13/03/13
30 Tell me more about T_D! escript data objects: T_D, x Values at data points >> print T_D.getShape(), x.getshape() () (2,) >>print T_D.getRank(), x.getrank() 0 1 >> print x.getfunctionspace() Finley_Nodes >> print T_D.getFunctionSpace() # no surprise! Finley_Nodes >> print grad(t_d).getfunctionspace() Finley_Elements >> print grad(t_d).getshape(), grad(x).getshape() (2,), (2,2) 30/46 13/03/13
31 PDE for Temperature Diffusion (κ T,i ), i =Q κ n i T,i =0 at x 0 =0 and x 0 =l 0 T =T D at x 1 =0 and x 1 =l 1 31/46 13/03/13
32 Steady PDEs in Escript from esys.escript.linearpdes import LinearPDE from esys.finley import Rectangle mydomain=rectangle(l0=l0,l1=l1,n0=80,n1=40) # create a PDE: mypde=linearpde(mydomain) # set the coefficients of the PDE mypde.setvalue(?????) # solve the PDE T = mypde.getsolution() 32/46 13/03/13
33 LinearPDE Interface (simplified) Interface to a single linear PDE for u ( A ij u, j ),i (B i u),i +C j u, j +D u= X j, j +Y natural BC: n i A ij u, j +n i B i u=n j X j + y contraints:u=r where q>0 Summation over i and j 33/46 13/03/13
34 LinearPDE for Diffusion Identify terms: ( A ij u,j ),i ( B i u ),i +C j u,j +Du= ( X,j ),j +Y (κt,i ),i =Q 34/46 13/03/13
35 LinearPDE for Diffusion (cont.) Natural boundary condition: n i A ij u,j +n i B i u=n j X,j +y n i κt,i =0 35/46 13/03/13
36 LinearPDE for Diffusion u=r where q> 0 T=T D where x 1 =0 or x 1 =l 1 36/46 13/03/13
37 LinearPDE for Diffusion (cont.) from esys.escript.linearpdes import LinearPDE from esys.finley import Rectangle mydomain=rectangle(l0=l0,l1=l1,n0=80,n1=40) # create a PDE: mypde=linearpde(mydomain) # set the coefficients of the PDE mypde.setvalue(a=?, Y=?, q=?, r=?, y=0) # solve the PDE T=mypde.getSolution() Not required! 37/46 13/03/13
38 Diffusion Term A (κt, i ),i =(κ δ ij T, j ),i =( A ij T, j ), i Kronecker symbol δ ij =[ ] = { A ij =κ δ ij =[ 1 if i= j } 0 else κ 0 ] 0 κ sum over i and j Check natural boundary conditions!!! 38/46 13/03/13
39 LinearPDE for Diffusion (cont.) from esys.escript.linearpdes import LinearPDE from esys.finley import Rectangle mydomain=rectangle(l0=l0,l1=l1,n0=80,n1=40) # create a PDE: mypde=linearpde(mydomain) # set the coefficients of the PDE mypde.setvalue(a=k*kronecker(mydomain)) Anisotropy is easy. 39/46 13/03/13
40 Source Q Ω Qc Q ( x )={ Q c for x x c <r } 0 otherwise 40/46 13/03/13
41 Source Q (continued) y = y y y 2 2 =esys.escript.length ( y ) from esys.escript import * xc=[7.,2.] r=0.8 Qc=100. x=mydomain.getx() Q=Qc*whereNegative(length(x-xc)-r) mypde.setvalue(y=q) 41/46 13/03/13
42 How to set q Make q one at top and bottom and zero elsewhere. # get the coordinates of points in the domain x=mydomain.getx() # f1=1 where x[1]==0 f1=wherezero(x[1]) f2==1 # f2=1 where x[1]==l1 x[1]==l1 f2=wherezero(x[1]-l1) q=f2+f1 mypde.setvalue(r=t_d,\ Ω q=q) f1==1 x[1]==0 42/46 13/03/13
43 Put it all together from esys.escript import * from esys.escript.linearpdes import LinearPDE from esys.finley import Rectangle from esys.weipa import savevtk L0=10.;L1=5. T_bot=100.; k=1. xc=[7.,2.]; r=0.8; Qc=100. mydomain=rectangle(l0=l0,l1=l1,n0=80,n1=40) x=mydomain.getx() T_D=T_bot/L1*(L1-x[1]) Q=Qc*whereNegative(length(x-xc)-r) mypde=linearpde(mydomain) mypde.setvalue(a=k*kronecker(mydomain),y=q, r=t_d, \ q=wherezero(x[1])+wherezero(x[1]-l1)) T=mypde.getSolution() savevtk('u.vtu', T=T, flux=-k*grad(t)) 43/46 13/03/13
44 Result 44/46 13/03/13
45 Some Remarks This PDE problem is symmetric: use mypde.setsymmetryon() Speeds up calculation Use mypde.settolerance() to control accuracy of linear solver for piecewise constant PDE coefficients use tagging ( mesh generation esys.pycad) 45/46 13/03/13
46 The End (for now) Questions? 46/46 13/03/13
Finite Element Method
In Finite Difference Methods: the solution domain is divided into a grid of discrete points or nodes the PDE is then written for each node and its derivatives replaced by finite-divided differences In
More informationPDE Project Course 4. An Introduction to DOLFIN and Puffin
PDE Project Course 4. An Introduction to DOLFIN and Puffin Anders Logg logg@math.chalmers.se Department of Computational Mathematics PDE Project Course 03/04 p. 1 Lecture plan DOLFIN Overview Input / output
More informationApplication of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem
Application of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem Malgorzata A. Jankowska 1, Andrzej Marciniak 2 and Tomasz Hoffmann 2 1 Poznan University
More informationd n U i dx n dx n δ n U i
Last time Taylor s series on equally spaced nodes Forward difference d n U i d n n U i h n + 0 h Backward difference d n U i d n n U i h n + 0 h Centered difference d n U i d n δ n U i or 2 h n + 0 h2
More informationFinancial Computing with Python
Introduction to Financial Computing with Python Matthieu Mariapragassam Why coding seems so easy? But is actually not Sprezzatura : «It s an art that doesn t seem to be an art» - The Book of the Courtier
More informationAbout Weak Form Modeling
Weak Form Modeling About Weak Form Modeling Do not be misled by the term weak; the weak form is very powerful and flexible. The term weak form is borrowed from mathematics. The distinguishing characteristics
More informationAs an example, we consider the following PDE with one variable; Finite difference method is one of numerical method for the PDE.
7. Introduction to the numerical integration of PDE. As an example, we consider the following PDE with one variable; Finite difference method is one of numerical method for the PDE. Accuracy requirements
More informationUsing Optimal Time Step Selection to Boost the Accuracy of FD Schemes for Variable-Coefficient PDEs
Using Optimal Time Step Selection to Boost the Accuracy of FD Schemes for Variable-Coefficient PDEs Kevin T. Chu and James V. Lambers Abstract We extend the technique of optimal time step OTS selection
More informationCS476/676 Mar 6, Today s Topics. American Option: early exercise curve. PDE overview. Discretizations. Finite difference approximations
CS476/676 Mar 6, 2019 1 Today s Topics American Option: early exercise curve PDE overview Discretizations Finite difference approximations CS476/676 Mar 6, 2019 2 American Option American Option: PDE Complementarity
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More informationChapter 20: An Introduction to ADI and Splitting Schemes
Chapter 20: An Introduction to ADI and Splitting Schemes 20.1INTRODUCTION AND OBJECTIVES In this chapter we discuss how to apply finite difference schemes to approximate the solution of multidimensional
More information(Refer Slide Time: 01:17)
Computational Electromagnetics and Applications Professor Krish Sankaran Indian Institute of Technology Bombay Lecture 06/Exercise 03 Finite Difference Methods 1 The Example which we are going to look
More informationBoundary conditions for options
Boundary conditions for options Boundary conditions for options can refer to the non-arbitrage conditions that option prices has to satisfy. If these conditions are broken, arbitrage can exist. to the
More informationReduced models for sparse grid discretizations of the multi-asset Black-Scholes equation
Reduced models for sparse grid discretizations of the multi-asset Black-Scholes equation The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationBoise State University GEOS 397 Fall 2016
Homework #3: Hill-slope evolution and plotting Due: 5:00 PM 09/16/16 Please read the following questions carefully and make sure to answer the problems completely. In your MATLAB script(s), please include
More information1 The Hull-White Interest Rate Model
Abstract Numerical Implementation of Hull-White Interest Rate Model: Hull-White Tree vs Finite Differences Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 30 April 2002 We implement the
More information32.4. Parabolic PDEs. Introduction. Prerequisites. Learning Outcomes
Parabolic PDEs 32.4 Introduction Second-order partial differential equations (PDEs) may be classified as parabolic, hyperbolic or elliptic. Parabolic and hyperbolic PDEs often model time dependent processes
More informationResearch Article Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation
Applied Mathematics Volume 1, Article ID 796814, 1 pages doi:11155/1/796814 Research Article Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation Zhongdi
More informationMonte Carlo Simulations in the Teaching Process
Monte Carlo Simulations in the Teaching Process Blanka Šedivá Department of Mathematics, Faculty of Applied Sciences University of West Bohemia, Plzeň, Czech Republic CADGME 2018 Conference on Digital
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More informationPDE Project Course 1. Adaptive finite element methods
PDE Project Course 1. Adaptive finite element methods Anders Logg logg@math.chalmers.se Department of Computational Mathematics PDE Project Course 03/04 p. 1 Lecture plan Introduction to FEM FEM for Poisson
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 2, Mikhail Zaslavsky 2 University of Michigan, Ann
More informationA local RBF method based on a finite collocation approach
Boundary Elements and Other Mesh Reduction Methods XXXVIII 73 A local RBF method based on a finite collocation approach D. Stevens & H. Power Department of Mechanical Materials and Manufacturing Engineering,
More informationIntroduction to Numerical PDEs
Introduction to Numerical PDEs Varun Shankar February 16, 2016 1 Introduction In this chapter, we will introduce a general classification scheme for linear second-order PDEs, and discuss when they have
More informationHOW TO GUIDE. The FINANCE module
HOW TO GUIDE The FINANCE module Copyright and publisher: EMD International A/S Niels Jernes vej 10 9220 Aalborg Ø Denmark Phone: +45 9635 44444 e-mail: emd@emd.dk web: www.emd.dk About energypro energypro
More informationOn the use of time step prediction
On the use of time step prediction CODE_BRIGHT TEAM Sebastià Olivella Contents 1 Introduction... 3 Convergence failure or large variations of unknowns... 3 Other aspects... 3 Model to use as test case...
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 3, Mikhail Zaslavsky 3 University of Michigan, Ann
More informationPICOF, Palaiseau, April 2-4, 2012
The Sobolev gradient regularization strategy for optical tomography coupled with a finite element formulation of the radiative transfer equation Fabien Dubot, Olivier Balima, Yann Favennec, Daniel Rousse
More informationExact shape-reconstruction by one-step linearization in EIT
Exact shape-reconstruction by one-step linearization in EIT Bastian von Harrach harrach@ma.tum.de Department of Mathematics - M1, Technische Universität München, Germany Joint work with Jin Keun Seo, Yonsei
More informationOutline. 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 informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
More informationChapter DIFFERENTIAL EQUATIONS: PHASE SPACE, NUMERICAL SOLUTIONS
Chapter 10 10. DIFFERENTIAL EQUATIONS: PHASE SPACE, NUMERICAL SOLUTIONS Abstract Solving differential equations analytically is not always the easiest strategy or even possible. In these cases one may
More informationConnected Mathematics 2, 6 th and 7th Grade Units 2009 Correlated to: Washington Mathematics Standards (Grade 6)
Grade 6 6.1. Core Content: Multiplication and division of fractions and decimals (Numbers, Operations, Algebra) 6.1.A Compare and order non-negative fractions, decimals, and integers using the number line,
More informationfor 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 informationProbabilistic Meshless Methods for Bayesian Inverse Problems. Jon Cockayne July 8, 2016
Probabilistic Meshless Methods for Bayesian Inverse Problems Jon Cockayne July 8, 2016 1 Co-Authors Chris Oates Tim Sullivan Mark Girolami 2 What is PN? Many problems in mathematics have no analytical
More informationMulti-Asset Options. A Numerical Study VILHELM NIKLASSON FRIDA TIVEDAL. Master s thesis in Engineering Mathematics and Computational Science
Multi-Asset Options A Numerical Study Master s thesis in Engineering Mathematics and Computational Science VILHELM NIKLASSON FRIDA TIVEDAL Department of Mathematical Sciences Chalmers University of Technology
More informationF.2 Factoring Trinomials
1 F.2 Factoring Trinomials In this section, we discuss factoring trinomials. We start with factoring quadratic trinomials of the form 2 + bbbb + cc, then quadratic trinomials of the form aa 2 + bbbb +
More informationSeminar: Efficient Monte Carlo Methods for Uncertainty Quantification
Seminar: Efficient Monte Carlo Methods for Uncertainty Quantification Elisabeth Ullmann Lehrstuhl für Numerische Mathematik (M2) TU München Elisabeth Ullmann (TU München) Efficient Monte Carlo for UQ 1
More informationAmerican Equity Option Valuation Practical Guide
Valuation Practical Guide John Smith FinPricing Summary American Equity Option Introduction The Use of American Equity Options Valuation Practical Guide A Real World Example American Option Introduction
More informationHIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR 1D PARABOLIC EQUATIONS. Ahmet İzmirlioğlu. BS, University of Pittsburgh, 2004
HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR D PARABOLIC EQUATIONS by Ahmet İzmirlioğlu BS, University of Pittsburgh, 24 Submitted to the Graduate Faculty of Art and Sciences in partial fulfillment of
More informationA GAUSS Implementation of Non Uniform Grids for PDE The PDE library
A GAUSS Implementation of Non Uniform Grids for PDE The PDE library Jérome Bodeau, Gaël Riboulet and Thierry Roncalli Groupe de Recherche Opérationnelle Bercy-Expo Immeuble Bercy Sud 4è étage 90quai de
More informationTopic #1: Evaluating and Simplifying Algebraic Expressions
John Jay College of Criminal Justice The City University of New York Department of Mathematics and Computer Science MAT 105 - College Algebra Departmental Final Examination Review Topic #1: Evaluating
More informationHeinz W. Engl. Industrial Mathematics Institute Johannes Kepler Universität Linz, Austria
Some Identification Problems in Finance Heinz W. Engl Industrial Mathematics Institute Johannes Kepler Universität Linz, Austria www.indmath.uni-linz.ac.at Johann Radon Institute for Computational and
More informationMAC Learning Objectives. Learning Objectives (Cont.)
MAC 1140 Module 12 Introduction to Sequences, Counting, The Binomial Theorem, and Mathematical Induction Learning Objectives Upon completing this module, you should be able to 1. represent sequences. 2.
More informationUsing radial basis functions for option pricing
Using radial basis functions for option pricing Elisabeth Larsson Division of Scientific Computing Department of Information Technology Uppsala University Actuarial Mathematics Workshop, March 19, 2013,
More informationManual Npv Formula Excel 2010 Example
Manual Npv Formula Excel 2010 Example As the example spreadsheet embedded below shows, the NPV is by its nature an Excel 2010 and 2013 Basics Demystifying The Excel Pro-Forma: What It Is And file below
More informationElastic demand solution methods
solution methods CE 392C October 6, 2016 REVIEW We ve added another consistency relationship: Route choices No vehicle left behind Link performance functions Shortest path OD matrix Travel times Demand
More informationNumerical software & tools for the actuarial community
Numerical software & tools for the actuarial community John Holden john.holden@nag.co.uk 20 th March 203 The Actuarial Profession Staple Inn Hall Experts in numerical algorithms and HPC services Agenda
More informationReview of Beginning Algebra MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
Review of Beginning Algebra MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Classify as an expression or an equation. 1) 2x + 9 1) A) Expression B)
More informationLattice-Theoretic Framework for Data-Flow Analysis. Defining Available Expressions Analysis. Reality Check! Reaching Constants
Lattice-Theoretic Framework for Data-Flow Analysis Defining Available Expressions Analysis Last time Generalizing data-flow analysis Today Finish generalizing data-flow analysis Reaching Constants introduction
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #1 March 3, 2014 Name: SOLUTION ID#: PS.:
More informationRemarks 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 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 informationHow to Model Multi-Party Service in SoaML? Written Date : September 4, 2013
Written Date : September 4, 2013 Multi-party service contracts refer to services that involve the participation of more than two participants, and with interaction between one another. An example would
More informationFinancial Modelling Using Discrete Stochastic Calculus
Preprint typeset in JHEP style - HYPER VERSION Financial Modelling Using Discrete Stochastic Calculus Eric A. Forgy, Ph.D. E-mail: eforgy@yahoo.com Abstract: In the present report, a review of discrete
More informationProject 1: Double Pendulum
Final Projects Introduction to Numerical Analysis II http://www.math.ucsb.edu/ atzberg/winter2009numericalanalysis/index.html Professor: Paul J. Atzberger Due: Friday, March 20th Turn in to TA s Mailbox:
More informationModels of the ocean: which ocean?
Models of the ocean: which ocean? Anne Marie Treguier, CNRS, Laboratoire de Physique des Océans, Brest, France Part 1: Some general statements - ocean models - convergence of solution - diffusion equation
More informationJAYARAM COLLEGE OF ENGINEERING AND TECHNOLOGY DEPARTMENT OF INFORMATION TECHNOLOGY
JAYARAM COLLEGE OF ENGINEERING AND TECHNOLOGY DEPARTMENT OF INFORMATION TECHNOLOGY Two Mark Question for Student s Reference 1. Define software project management. Software Project Management has key ideas
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationVasicek Model in Interest Rates
Vasicek Model in Interest Rates In case of Stocks we can model them as a diffusion process where they can end up anywhere in the universe of Prices, but Interest Rates are not so. Interest rates can not
More informationAEM Computational Fluid Dynamics Instructor: Dr. M. A. R. Sharif
AEM 620 - Computational Fluid Dynamics Instructor: Dr. M. A. R. Sharif Numerical Solution Techniques for 1-D Parabolic Partial Differential Equations: Transient Flow Problem by Parshant Dhand September
More informationOptimizing 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 informationWhat can we do with numerical optimization?
Optimization motivation and background Eddie Wadbro Introduction to PDE Constrained Optimization, 2016 February 15 16, 2016 Eddie Wadbro, Introduction to PDE Constrained Optimization, February 15 16, 2016
More informationPDE Methods for the Maximum Drawdown
PDE Methods for the Maximum Drawdown Libor Pospisil, Jan Vecer Columbia University, Department of Statistics, New York, NY 127, USA April 1, 28 Abstract Maximum drawdown is a risk measure that plays an
More informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationON AN IMPLEMENTATION OF BLACK SCHOLES MODEL FOR ESTIMATION OF CALL- AND PUT-OPTION VIA PROGRAMMING ENVIRONMENT MATHEMATICA
Доклади на Българската академия на науките Comptes rendus de l Académie bulgare des Sciences Tome 66, No 5, 2013 MATHEMATIQUES Mathématiques appliquées ON AN IMPLEMENTATION OF BLACK SCHOLES MODEL FOR ESTIMATION
More informationComputational Finance Finite Difference Methods
Explicit finite difference method Computational Finance Finite Difference Methods School of Mathematics 2018 Today s Lecture We now introduce the final numerical scheme which is related to the PDE solution.
More informationUsing condition numbers to assess numerical quality in HPC applications
Using condition numbers to assess numerical quality in HPC applications Marc Baboulin Inria Saclay / Université Paris-Sud, France INRIA - Illinois Petascale Computing Joint Laboratory 9th workshop, June
More informationScenario 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 information3 Ways to Write Ratios
RATIO & PROPORTION Sec 1. Defining Ratio & Proportion A RATIO is a comparison between two quantities. We use ratios everyday; one Pepsi costs 50 cents describes a ratio. On a map, the legend might tell
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationFitob R Package: Manual & Tutorial
Fitob R Package: Manual & Tutorial Janos Benk, benkjanos@gmail.com 2012.11.14 Contents 1 Introduction 1 1.1 Mathematical and Numerical Background............. 2 2 User Interface 4 2.1 Inputing a Financial
More informationMagnet Resonance Electrical Impedance Tomography (MREIT): convergence of the Harmonic B z Algorithm
Magnet Resonance Electrical Impedance Tomography (MREIT): convergence of the Harmonic B z Algorithm Dominik Garmatter garmatter@math.uni-frankfurt.de Group for Numerics of PDEs, Goethe University Frankfurt,
More informationInfinitely Many Solutions to the Black-Scholes PDE; Physics Point of View
CBS 2018-05-23 1 Infinitely Many Solutions to the Black-Scholes PDE; Physics Point of View 서울대학교물리학과 2018. 05. 23. 16:00 (56 동 106 호 ) 최병선 ( 경제학부 ) 최무영 ( 물리천문학부 ) CBS 2018-05-23 2 Featuring: 최병선 Pictures
More informationFINITE DIFFERENCE METHODS
FINITE DIFFERENCE METHODS School of Mathematics 2013 OUTLINE Review 1 REVIEW Last time Today s Lecture OUTLINE Review 1 REVIEW Last time Today s Lecture 2 DISCRETISING THE PROBLEM Finite-difference approximations
More informationExact shape-reconstruction by one-step linearization in EIT
Exact shape-reconstruction by one-step linearization in EIT Bastian von Harrach harrach@math.uni-mainz.de Zentrum Mathematik, M1, Technische Universität München, Germany Joint work with Jin Keun Seo, Yonsei
More informationEnergy Systems under Uncertainty: Modeling and Computations
Energy Systems under Uncertainty: Modeling and Computations W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Systems Analysis 2015, November 11 13, IIASA (Laxenburg,
More informationThe Intermediate Value Theorem states that if a function g is continuous, then for any number M satisfying. g(x 1 ) M g(x 2 )
APPM/MATH 450 Problem Set 5 s This assignment is due by 4pm on Friday, October 25th. You may either turn it in to me in class or in the box outside my office door (ECOT 235). Minimal credit will be given
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationFINITE MATH LECTURE NOTES. c Janice Epstein 1998, 1999, 2000 All rights reserved.
FINITE MATH LECTURE NOTES c Janice Epstein 1998, 1999, 2000 All rights reserved. August 27, 2001 Chapter 1 Straight Lines and Linear Functions In this chapter we will learn about lines - how to draw them
More informationLecture 4: Divide and Conquer
Lecture 4: Divide and Conquer Divide and Conquer Merge sort is an example of a divide-and-conquer algorithm Recall the three steps (at each level to solve a divideand-conquer problem recursively Divide
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Program Advisor: Howard Elman Department of Computer Science May 5, 2016 Tengfei
More informationFinite Difference Approximation of Hedging Quantities in the Heston model
Finite Difference Approximation of Hedging Quantities in the Heston model Karel in t Hout Department of Mathematics and Computer cience, University of Antwerp, Middelheimlaan, 22 Antwerp, Belgium Abstract.
More informationSection 7C Finding the Equation of a Line
Section 7C Finding the Equation of a Line When we discover a linear relationship between two variables, we often try to discover a formula that relates the two variables and allows us to use one variable
More informationA NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS. Burhaneddin İZGİ
A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS Burhaneddin İZGİ Department of Mathematics, Istanbul Technical University, Istanbul, Turkey
More informationNet lift and return maximization. Victor D. Zurkowski Analytics Consultant Metrics and Analytics CIBC National Collection
Net lift and return maximization Victor D. Zurkowski Analytics Consultant Metrics and Analytics CIBC National Collection Page 2 Page 3 Could I have been wrong all along? Page 4 There has been recent mentions
More informationLab12_sol. November 21, 2017
Lab12_sol November 21, 2017 1 Sample solutions of exercises of Lab 12 Suppose we want to find the current prices of an European option so that the error at the current stock price S(0) = S 0 was less that
More informationNUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE
Trends in Mathematics - New Series Information Center for Mathematical Sciences Volume 13, Number 1, 011, pages 1 5 NUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE YONGHOON
More informationImplementing hybrid PDE solvers
Implementing hybrid PDE solvers George Sarailidis and Manolis Vavalis ECE Department, University of Thessaly November 14, 2015 Outline Objectives State-of-the-art Basic Implementations Experimentation
More information2. ANALYTICAL TOOLS. E(X) = P i X i = X (2.1) i=1
2. ANALYTICAL TOOLS Goals: After reading this chapter, you will 1. Know the basic concepts of statistics: expected value, standard deviation, variance, covariance, and coefficient of correlation. 2. Use
More informationPartial Differential Equations of Fluid Dynamics
Partial Differential Equations of Fluid Dynamics Ville Vuorinen,D.Sc.(Tech.) 1 1 Department of Energy Technology, Internal Combustion Engine Research Group Department of Energy Technology Outline Introduction
More informationFactoring is the process of changing a polynomial expression that is essentially a sum into an expression that is essentially a product.
Ch. 8 Polynomial Factoring Sec. 1 Factoring is the process of changing a polynomial expression that is essentially a sum into an expression that is essentially a product. Factoring polynomials is not much
More informationFinite Potential Well
Finite Potential Well These notes are provided as a supplement to the text and a replacement for the lecture on 11/16/17. Make sure you fill in the steps outlined in red. The finite potential well problem
More informationLecture 2: The Neoclassical Growth Model
Lecture 2: The Neoclassical Growth Model Florian Scheuer 1 Plan Introduce production technology, storage multiple goods 2 The Neoclassical Model Three goods: Final output Capital Labor One household, with
More informationCommon Core Georgia Performance Standards
A Correlation of Pearson Connected Mathematics 2 2012 to the Common Core Georgia Performance s Grade 6 FORMAT FOR CORRELATION TO THE COMMON CORE GEORGIA PERFORMANCE STANDARDS (CCGPS) Subject Area: K-12
More informationThe Software Engineering Discipline. Computer Aided Software Engineering (CASE) tools. Chapter 7: Software Engineering
Chapter 7: Software Engineering Computer Science: An Overview Tenth Edition by J. Glenn Brookshear Chapter 7: Software Engineering 7.1 The Software Engineering Discipline 7.2 The Software Life Cycle 7.3
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 13. The Black-Scholes PDE Steve Yang Stevens Institute of Technology 04/25/2013 Outline 1 The Black-Scholes PDE 2 PDEs in Asset Pricing 3 Exotic
More informationLecture 20 Heat Equation and Parabolic Problems
.539 Lecture 2 Heat Equation and Parabolic Problems Prof. Dean Wang. Classification of Differential Equations In most texts the classification is given for a linear second-order differential equation in
More informationMaking the Link between Actuaries and Data Science
Making the Link between Actuaries and Data Science Simon Lee, Cecilia Chow, Thibault Imbert AXA Asia 2 nd ASHK General Insurance & Data Analytics Seminar Friday 7 October 2016 1 Agenda Data Driving Insurers
More informationNo-arbitrage theorem for multi-factor uncertain stock model with floating interest rate
Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer
More information