DOCSLIB.ORG

65 Numerical Analysis

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text

Load more

e Q (e t o 5 M SectionsSet 1Q (Section 65)MR September 2012 65 NUMERICAL ANALYSIS MR2918625 65-06 FRecent advances in scientific computing and matrix analysis. Proceedings of the International Workshop held at the University of Macau, Macau, December 28{30, 2009. Edited by Xiao-Qing Jin, Hai-Wei Sun and Seak-Weng Vong. International Press, Somerville, MA; Higher Education Press, Beijing, 2011. xii+126 pp. ISBN 978-1-57146-202-2 Contents: Zheng-jian Bai and Xiao-qing Jin [Xiao Qing Jin1], A note on the Ulm-like method for inverse eigenvalue problems (1{7) MR2908437; Che-man Cheng [Che-Man Cheng], Kin-sio Fong [Kin-Sio Fong] and Io-kei Lok [Io-Kei Lok], Another proof for commutators with maximal Frobenius norm (9{14) MR2908438; Wai-ki Ching [Wai-Ki Ching] and Dong-mei Zhu [Dong Mei Zhu1], On high-dimensional Markov chain models for categorical data sequences with applications (15{34) MR2908439; Yan-nei Law [Yan Nei Law], Hwee-kuan Lee [Hwee Kuan Lee], Chao-qiang Liu [Chaoqiang Liu] and Andy M. Yip, An additive variational model for image segmentation (35{48) MR2908440; Hai-yong Liao [Haiyong Liao] and Michael K. Ng, Total variation image restoration with automatic selection of regularization parameters (49{59) MR2908441; Franklin T. Luk and San-zheng Qiao [San Zheng Qiao], Matrices and the LLL algorithm (61{69) MR2908442; Mila Nikolova, Michael K. Ng and Chi-pan Tam [Chi-Pan Tam], A fast nonconvex nonsmooth minimization method for image restoration and reconstruction (71{83) MR2908443; Gang Wu [Gang Wu1], Eigenvalues of certain augmented complex stochastic matrices with applications to PageRank (85{92) MR2908444; Yan Xuan and Fu-rong Lin, Clenshaw-Curtis-rational quadrature rule for Wiener-Hopf equations of the second kind (93{110) MR2908445; Man-chung Yeung [Man-Chung Yeung], On the solution of singular systems by Krylov subspace methods (111{116) MR2908446; Qi- fang Yu, San-zheng Qiao [San Zheng Qiao] and Yi-min Wei, A comparative study of the LLL algorithm (117{126) MR2908447. fThe papers are being reviewed individually.g 65C Probabilistic methods, simulation and stochastic differential equations MR2871782 65C05 11B39 65D30 Bilyk, Dmitriy (1-SC; Columbia, SC); Temlyakov, V. N. [Temlyakov, Vladimir N.] (1-SC; Columbia, SC); Yu, Rui [Yu, Rui1] (1-SC; Columbia, SC) Fibonacci sets and symmetrization in discrepancy theory. (English summary) J. Complexity 28 (2012), no. 1, 18{36. Summary: \We study the Fibonacci sets from the point of view of their quality 1 with respect to discrepancy and numerical integration. Let fbngn=0 be the sequence 2 of Fibonacci numbers. The bn-point Fibonacci set Fn ⊂ [0; 1] is defined as Fn := bn f(µ/bn; fµbn−1=bng)gµ=1, where fxg is the fractional part of a number x 2 R. It is known that cubature formulas based on the Fibonacci set Fn give optimal rate of error of numerical integration for certain classes of functions with mixed smoothness. \We give a Fourier analytic proof of the fact that the symmetrized Fibonacci set 0 Fn = Fn [ f(p1; 1 − p2): (pq; p2) 2 Fng has asymptotically minimal L2 discrepancy. This approach also yields an exact formula for this quantity, which allows us to evaluate the constant in the discrepancy estimates. Numerical computations indicate that these sets 1 e Q (e t o 5 M SectionsSet 1Q (Section 65)MR September 2012 have the smallest currently known L2 discrepancy among two-dimensional point sets. \We also introduce quartered Lp discrepancy, which is a modification of the Lp discrepancy symmetrized with respect to the center of the unit square. We prove that the Fibonacci set Fn has minimal in the sense of order quartered Lp discrepancy for all p 2 (1; 1). This in turn implies that certain two-fold symmetrizations of the Fibonacci set Fn are optimal with respect to the standard Lp discrepancy." MR2871781 65C05 46E22 65D30 Gnewuch, Michael (D-KIEL-II; Kiel) Weighted geometric discrepancies and numerical integration on reproducing kernel Hilbert spaces. (English summary) J. Complexity 28 (2012), no. 1, 2{17. Summary: \We extend the notion of L2-B-discrepancy introduced in [E. Novak and H. Wo´zniakowski, in Analytic number theory, 359{388, Cambridge Univ. Press, Cam- bridge, 2009; MR2508657 (2010d:11088)] to what we shall call weighted geometric L2-discrepancy. This extension enables us to consider weights in order to moderate the importance of different groups of variables, as well as to consider volume measures dif- ferent from the Lebesgue measure and classes of test sets different from measurable subsets of Euclidean spaces. \We relate the weighted geometric L2-discrepancy to numerical integration defined over weighted reproducing kernel Hilbert spaces and settle in this way an open problem posed by Novak and Wo´zniakowski. \Furthermore, we prove an upper bound for the numerical integration error for cubature formulas that use admissible sample points. The set of admissible sample points may actually be a subset of the integration domain of measure zero. We illustrate that particularly in infinite-dimensional numerical integration it is crucial to distinguish between the whole integration domain and the set of those sample points that actually can be used by the algorithms." MR2822407 65C05 60G70 62G20 62G32 90B25 Guyader, Arnaud (F-RENN2-NDM; Rennes); Hengartner, Nicholas [Hengartner, Nicolas W.] (1-LANL-IF; Los Alamos, NM); Matzner-Løber, Eric [Matzner-Løber, Eric]´ (F-RENN2-NDM; Rennes) Simulation and estimation of extreme quantiles and extreme probabilities. (English summary) Appl. Math. Optim. 64 (2011), no. 2, 171{196. The paper presents an efficient algorithm for estimating (i) a tail probability given a quantile, or (ii) a quantile given a tail probability. The algorithm improves the multilevel splitting methods known in the literature. By means of Poisson process tools the exact distribution of the estimated probabilities and quantiles are established. Bernd Heidergott MR2788863 65C05 Lee, Jeong Eun (5-QUT-SM; Brisbane); McVinish, Ross [McVinish, R¨oss] (5-QLD; Brisbane); Mengersen, Kerrie [Mengersen, Kerrie L.] (5-QUT-SM; Brisbane) Population Monte Carlo algorithm in high dimensions. (English summary) Methodol. Comput. Appl. Probab. 13 (2011), no. 2, 369{389. Summary: \The population Monte Carlo algorithm is an iterative importance sampling scheme for solving static problems. We examine the population Monte Carlo algorithm in a simplified setting, a single step of the general algorithm, and study a fundamental 2 e Q (e t o 5 M SectionsSet 1Q (Section 65)MR September 2012 problem that occurs in applying importance sampling to high-dimensional problem. The precision of the computed estimate from the simplified setting is measured by the asymptotic variance of estimate under conditions on the importance function. We demonstrate the exponential growth of the asymptotic variance with the dimension and show that the optimal covariance matrix for the importance function can be estimated in special cases." MR2839113 65C05 93D09 Nurges, Ulo¨ (ES-TALL; Tallinn) Discussion on: \On the generation of random stable polynomials" [MR2839112]. Eur. J. Control 17 (2011), no. 2, 160{161. This paper presents a discussion of [P. S. Shcherbakov and F. Dabbene, Eur. J. Control 17 (2011), no. 2, 145{159; MR2839112] on the generation of stable random polynomials. Two issues regarding the paper are raised. The first concerns the terminology assigned to a set of parameters appearing in a particular method (based on the Schur-Cohn stability test, and called Levinson-Durbin (LD) parameters in [op. cit.]; the authors of the original paper and of this discussion both acknowledge that these parameters go by many names in many different fields). The second relates to the fact that the uniform distribution for LD parameters for high-order polynomials leads to the generation of polynomials whose roots cluster close to the stability boundary. A modified algorithm avoiding this behavior is presented, together with simulation results. fFor further information pertaining to this item see [P. Shcherbakov and F. Dabbene, Eur. J. Control 17 (2011), no. 2, 161; MR2839120].g David C. Saunders MR2839120 65C05 93D05 Shcherbakov, P.; Dabbene, F. [Dabbene, Fabrizio] Final comments by the authors: \On the generation of random stable polynomials" [MR2839112; MR2839113]. Eur. J. Control 17 (2011), no. 2, 161. The authors respond to the two points raised in [U.¨ Nurges, Eur. J. Control 17 (2011), no. 2, 160{161; MR2839113]. They justify their choice of terminology on historical grounds. They raise some issues regarding the algorithm proposed in [op. cit.], and further note that clustering of roots near the stability boundary may not necessarily be a bad thing from the point of view of practical applications, pointing to the motivating examples in [P. S. Shcherbakov and F. Dabbene, Eur. J. Control 17 (2011), no. 2, 145{ 159; MR2839112] for justification. David C. Saunders MR2839112 65C05 93D09 Shcherbakov, Pavel [Shcherbakov, Pavel S.] (RS-AOS-CN; Moscow); Dabbene, Fabrizio (I-TRNP-IEI; Turin) On the generation of random stable polynomials. (English summary) Eur. J. Control 17 (2011), no. 2, 145{159. This paper presents a survey of various techniques for the random generation of stable polynomials. Both the case where the roots must lie in the open unit disk (Schur stable), and the case where they must lie in the left half-plane C− = fs 2 C : Re s < 0g (Hurwitz stable), are considered. The paper begins with an introduction presenting motivating applications. Next, how general random sampling schemes may be adapted to the problem of generating stable random polynomials is discussed. Random generation schemes based on classical tests for stability of polynomials are then presented. Finally, the paper concludes with a list of interesting open problems in the area. Simulation results are presented for many of the schemes discussed in the paper. 3 e Q (e t o 5 M SectionsSet 1Q (Section 65)MR September 2012 fFor further information pertaining to this item see [U.¨ Nurges, Eur.
Recommended publications
  • Numerical Solution of Jump-Diffusion LIBOR Market Models

    Numerical Solution of Jump-Diffusion LIBOR Market Models

    Finance Stochast. 7, 1–27 (2003) c Springer-Verlag 2003 Numerical solution of jump-diffusion LIBOR market models Paul Glasserman1, Nicolas Merener2 1 403 Uris Hall, Graduate School of Business, Columbia University, New York, NY 10027, USA (e-mail: [email protected]) 2 Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027, USA (e-mail: [email protected]) Abstract. This paper develops, analyzes, and tests computational procedures for the numerical solution of LIBOR market models with jumps. We consider, in par- ticular, a class of models in which jumps are driven by marked point processes with intensities that depend on the LIBOR rates themselves. While this formulation of- fers some attractive modeling features, it presents a challenge for computational work. As a first step, we therefore show how to reformulate a term structure model driven by marked point processes with suitably bounded state-dependent intensities into one driven by a Poisson random measure. This facilitates the development of discretization schemes because the Poisson random measure can be simulated with- out discretization error. Jumps in LIBOR rates are then thinned from the Poisson random measure using state-dependent thinning probabilities. Because of discon- tinuities inherent to the thinning process, this procedure falls outside the scope of existing convergence results; we provide some theoretical support for our method through a result establishing first and second order convergence of schemes that accommodates thinning but imposes stronger conditions on other problem data. The bias and computational efficiency of various schemes are compared through numerical experiments. Key words: Interest rate models, Monte Carlo simulation, market models, marked point processes JEL Classification: G13, E43 Mathematics Subject Classification (1991): 60G55, 60J75, 65C05, 90A09 The authors thank Professor Steve Kou for helpful comments and discussions.
  • CURRICULUM VITAE Anatoliy SWISHCHUK Department Of

    CURRICULUM VITAE Anatoliy SWISHCHUK Department Of

    CURRICULUM VITAE Anatoliy SWISHCHUK Department of Mathematics & Statistics, University of Calgary 2500 University Drive NW, Calgary, Alberta, Canada T2N 1N4 Office: MS552 E-mails: [email protected] Tel: +1 (403) 220-3274 (office) home page: http://www.math.ucalgary.ca/~aswish/ Education: • Doctor of Phys. & Math. Sci. (1992, Doctorate, Institute of Mathematics, National Academy of Sciences of Ukraine (NASU), Kiev, Ukraine) • Ph.D. (1984, Institute of Mathematics, NASU, Kiev, Ukraine) • M.Sc., B.Sc. (1974-1979, Kiev State University, Faculty of Mathematics & Mechanics, Probability Theory & Mathematical Statistics Department, Kiev, Ukraine) Work Experience: • Full Professor, Department of Mathematics and Statistics, University of Calgary, Calgary, Canada (April 2012-present) • Co-Director, Mathematical and Computational Finance Laboratory, Department of Math- ematics and Statistics, University of Calgary, Calgary, Canada (October 2004-present) • Associate Professor, Department of Mathematics and Statistics, University of Calgary, Calgary, Canada (July 2006-March 2012) • Assistant Professor, Department of Mathematics and Statistics, University of Calgary, Calgary, Canada (August 2004-June 2006) • Course Director, Department of Mathematics & Statistics, York University, Toronto, ON, Canada (January 2003-June 2004) • Research Associate, Laboratory for Industrial & Applied Mathematics, Department of Mathematics & Statistics, York University, Toronto, ON, Canada (November 1, 2001-July 2004) • Professor, Probability Theory & Mathematical Statistics
  • Research on the Digitial Image Based on Hyper-Chaotic And

    Research on the Digitial Image Based on Hyper-Chaotic And

    Research on digital image watermark encryption based on hyperchaos Item Type Thesis or dissertation Authors Wu, Pianhui Publisher University of Derby Download date 27/09/2021 09:45:19 Link to Item http://hdl.handle.net/10545/305004 UNIVERSITY OF DERBY RESEARCH ON DIGITAL IMAGE WATERMARK ENCRYPTION BASED ON HYPERCHAOS Pianhui Wu Doctor of Philosophy 2013 RESEARCH ON DIGITAL IMAGE WATERMARK ENCRYPTION BASED ON HYPERCHAOS A thesis submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy By Pianhui Wu BSc. MSc. Faculty of Business, Computing and Law University of Derby May 2013 To my parents Acknowledgements I would like to thank sincerely Professor Zhengxu Zhao for his guidance, understanding, patience and most importantly, his friendship during my graduate studies at the University of Derby. His mentorship was paramount in providing a well-round experience consistent with my long-term career goals. I am grateful to many people in Faculty of Business, Computing and Law at the University of Derby for their support and help. I would also like to thank my parents, who have given me huge support and encouragement. Their advice is invaluable. An extra special recognition to my sister whose love and aid have made this thesis possible, and my time in Derby a colorful and wonderful experience. I Glossary AC Alternating Current AES Advanced Encryption Standard CCS Combination Coordinate Space CWT Continue Wavelet Transform BMP Bit Map DC Direct Current DCT Discrete Cosine Transform DWT Discrete Wavelet Transform
  • Introduction to Black's Model for Interest Rate

    Introduction to Black's Model for Interest Rate

    INTRODUCTION TO BLACK'S MODEL FOR INTEREST RATE DERIVATIVES GRAEME WEST AND LYDIA WEST, FINANCIAL MODELLING AGENCY© Contents 1. Introduction 2 2. European Bond Options2 2.1. Different volatility measures3 3. Caplets and Floorlets3 4. Caps and Floors4 4.1. A call/put on rates is a put/call on a bond4 4.2. Greeks 5 5. Stripping Black caps into caplets7 6. Swaptions 10 6.1. Valuation 11 6.2. Greeks 12 7. Why Black is useless for exotics 13 8. Exercises 13 Date: July 11, 2011. 1 2 GRAEME WEST AND LYDIA WEST, FINANCIAL MODELLING AGENCY© Bibliography 15 1. Introduction We consider the Black Model for futures/forwards which is the market standard for quoting prices (via implied volatilities). Black[1976] considered the problem of writing options on commodity futures and this was the first \natural" extension of the Black-Scholes model. This model also is used to price options on interest rates and interest rate sensitive instruments such as bonds. Since the Black-Scholes analysis assumes constant (or deterministic) interest rates, and so forward interest rates are realised, it is difficult initially to see how this model applies to interest rate dependent derivatives. However, if f is a forward interest rate, it can be shown that it is consistent to assume that • The discounting process can be taken to be the existing yield curve. • The forward rates are stochastic and log-normally distributed. The forward rates will be log-normally distributed in what is called the T -forward measure, where T is the pay date of the option.
  • Newton-Krylov-BDDC Solvers for Nonlinear Cardiac Mechanics

    Newton-Krylov-BDDC Solvers for Nonlinear Cardiac Mechanics

    Newton-Krylov-BDDC solvers for nonlinear cardiac mechanics Item Type Article Authors Pavarino, L.F.; Scacchi, S.; Zampini, Stefano Citation Newton-Krylov-BDDC solvers for nonlinear cardiac mechanics 2015 Computer Methods in Applied Mechanics and Engineering Eprint version Post-print DOI 10.1016/j.cma.2015.07.009 Publisher Elsevier BV Journal Computer Methods in Applied Mechanics and Engineering Rights NOTICE: this is the author’s version of a work that was accepted for publication in Computer Methods in Applied Mechanics and Engineering. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Computer Methods in Applied Mechanics and Engineering, 18 July 2015. DOI:10.1016/ j.cma.2015.07.009 Download date 07/10/2021 05:34:35 Link to Item http://hdl.handle.net/10754/561071 Accepted Manuscript Newton-Krylov-BDDC solvers for nonlinear cardiac mechanics L.F. Pavarino, S. Scacchi, S. Zampini PII: S0045-7825(15)00221-2 DOI: http://dx.doi.org/10.1016/j.cma.2015.07.009 Reference: CMA 10661 To appear in: Comput. Methods Appl. Mech. Engrg. Received date: 13 December 2014 Revised date: 3 June 2015 Accepted date: 8 July 2015 Please cite this article as: L.F. Pavarino, S. Scacchi, S. Zampini, Newton-Krylov-BDDC solvers for nonlinear cardiac mechanics, Comput. Methods Appl. Mech. Engrg. (2015), http://dx.doi.org/10.1016/j.cma.2015.07.009 This is a PDF file of an unedited manuscript that has been accepted for publication.
  • A High-Order Boris Integrator

    A High-Order Boris Integrator

    A high-order Boris integrator Mathias Winkela,∗, Robert Speckb,a, Daniel Ruprechta aInstitute of Computational Science, University of Lugano, Switzerland. bJ¨ulichSupercomputing Centre, Forschungszentrum J¨ulich,Germany. Abstract This work introduces the high-order Boris-SDC method for integrating the equations of motion for electrically charged particles in an electric and magnetic field. Boris-SDC relies on a combination of the Boris-integrator with spectral deferred corrections (SDC). SDC can be considered as preconditioned Picard iteration to com- pute the stages of a collocation method. In this interpretation, inverting the preconditioner corresponds to a sweep with a low-order method. In Boris-SDC, the Boris method, a second-order Lorentz force integrator based on velocity-Verlet, is used as a sweeper/preconditioner. The presented method provides a generic way to extend the classical Boris integrator, which is widely used in essentially all particle-based plasma physics simulations involving magnetic fields, to a high-order method. Stability, convergence order and conserva- tion properties of the method are demonstrated for different simulation setups. Boris-SDC reproduces the expected high order of convergence for a single particle and for the center-of-mass of a particle cloud in a Penning trap and shows good long-term energy stability. Keywords: Boris integrator, time integration, magnetic field, high-order, spectral deferred corrections (SDC), collocation method 1. Introduction Often when modeling phenomena in plasma physics, for example particle dynamics in fusion vessels or particle accelerators, an externally applied magnetic field is vital to confine the particles in the physical device [1, 2]. In many cases, such as instabilities [3] and high-intensity laser plasma interaction [4], the magnetic field even governs the microscopic evolution and drives the phenomena to be studied.
  • Strong Form Meshfree Collocation Method for Higher Order and Nonlinear Pdes in Engineering Applications

    Strong Form Meshfree Collocation Method for Higher Order and Nonlinear Pdes in Engineering Applications

    University of Colorado, Boulder CU Scholar Civil Engineering Graduate Theses & Dissertations Civil, Environmental, and Architectural Engineering Spring 1-1-2019 Strong Form Meshfree Collocation Method for Higher Order and Nonlinear Pdes in Engineering Applications Andrew Christian Beel University of Colorado at Boulder, [email protected] Follow this and additional works at: https://scholar.colorado.edu/cven_gradetds Part of the Civil Engineering Commons, and the Partial Differential Equations Commons Recommended Citation Beel, Andrew Christian, "Strong Form Meshfree Collocation Method for Higher Order and Nonlinear Pdes in Engineering Applications" (2019). Civil Engineering Graduate Theses & Dissertations. 471. https://scholar.colorado.edu/cven_gradetds/471 This Thesis is brought to you for free and open access by Civil, Environmental, and Architectural Engineering at CU Scholar. It has been accepted for inclusion in Civil Engineering Graduate Theses & Dissertations by an authorized administrator of CU Scholar. For more information, please contact [email protected]. Strong Form Meshfree Collocation Method for Higher Order and Nonlinear PDEs in Engineering Applications Andrew Christian Beel B.S., University of Colorado Boulder, 2019 A thesis submitted to the Faculty of the Graduate School of the University of Colorado Boulder in partial fulfillment of the requirement for the degree of Master of Science Department of Civil, Environmental and Architectural Engineering 2019 This thesis entitled: Strong Form Meshfree Collocation Method for Higher Order and Nonlinear PDEs in Engineering Applications written by Andrew Christian Beel has been approved for the Department of Civil, Environmental and Architectural Engineering Dr. Jeong-Hoon Song (Committee Chair) Dr. Ronald Pak Dr. Victor Saouma Dr. Richard Regueiro Date The final copy of this thesis has been examined by the signatories, and we find that both the content and the form meet acceptable presentation standards of scholarly work in the above mentioned discipline.
  • Numerical Methods for Wave Equations: Part I: Smooth Solutions

    Numerical Methods for Wave Equations: Part I: Smooth Solutions

    WPPII Computational Fluid Dynamics I Numerical Methods for Wave Equations: Part I: Smooth Solutions Instructor: Hong G. Im University of Michigan Fall 2001 Computational Fluid Dynamics I WPPII Outline Solution Methods for Wave Equation Part I • Method of Characteristics • Finite Volume Approach and Conservative Forms • Methods for Continuous Solutions - Central and Upwind Difference - Stability, CFL Condition - Various Stable Methods Part II • Methods for Discontinuous Solutions - Burgers Equation and Shock Formation - Entropy Condition - Various Numerical Schemes WPPII Computational Fluid Dynamics I Method of Characteristics WPPII Computational Fluid Dynamics I 1st Order Wave Equation ∂f ∂f + c = 0 ∂t ∂x The characteristics for this equation are: dx df = c; = 0; f dt dt t f x WPPII Computational Fluid Dynamics I 1-D Wave Equation (2nd Order Hyperbolic PDE) ∂ 2 f ∂ 2 f − c2 = 0 ∂t 2 ∂x2 Define ∂f ∂f v = ; w = ; ∂t ∂x which leads to ∂v ∂w − c2 = 0 ∂t ∂x ∂w ∂v − = 0 ∂t ∂x WPPII Computational Fluid Dynamics I In matrix form v 0 − c2 v t + x = 0 or u + Au = 0 − t x wt 1 0 wx Can it be transformed into the form v λ 0v t + x = 0 or u + λu = 0 ? λ t x wt 0 wx Find the eigenvalue, eigenvector ()AT − λI q = 0 − λ −1 l 1 = 0 − 2 − λ c l2 WPPII Computational Fluid Dynamics I Eigenvalue dx AT − λI = 0; λ2 − c2 = 0; λ = = ±c dt Eigenvector 1 For λ = +c, − cl − l = 0; q = 1 2 1 − c 1 For λ = −c, cl − l = 0; q = 1 2 2 c The solution (v,w) is governed by ODE’s along the characteristic lines dx / dt = ±c WPPII
  • Pricing Bermudan Swaptions on the LIBOR Market Model Using the Stochastic Grid Bundling Method

    Pricing Bermudan Swaptions on the LIBOR Market Model Using the Stochastic Grid Bundling Method

    Pricing Bermudan Swaptions on the LIBOR Market Model using the Stochastic Grid Bundling Method Stef Maree,∗ Jacques du Toity Abstract We examine using the Stochastic Grid Bundling Method (SGBM) to price a Bermu- dan swaption driven by a one-factor LIBOR Market Model (LMM). Using a well- known approximation formula from the finance literature, we implement SGBM with one basis function and show that it is around six times faster than the equivalent Longstaff–Schwartz method. The two methods agree in price to one basis point, and the SGBM path estimator gives better (higher) prices than the Longstaff–Schwartz prices. A closer examination shows that inaccuracies in the approximation formula introduce a small bias into the SGBM direct estimator. 1 Introduction The LIBOR Market Model (LMM) is an interest rate market model. It owes much of its popularity to the fact that it is consistent with Black’s pricing formulas for simple interest rate derivatives such as caps and swaptions, see, for example, [3]. Under LMM all forward rates are modelled as individual processes, which typically results in a high dimensional model. Consequently when pricing more complex instruments under LMM, Monte Carlo simulation is usually used. When pricing products with early-exercise features, such as Bermudan swaptions, a method is required to determine the early-exercise policy. The most popular approach to this is the Longstaff–Schwartz Method (LSM) [8]. When time is discrete (as is usually the case in numerical schemes) Bellman’s backward induction principle says that the price of the option at any exercise time is the maximum of the spot payoff and the so-called contin- uation value.
  • Workshop High-Order Numerical Approximation for Partial Differential Equations

    Workshop High-Order Numerical Approximation for Partial Differential Equations

    Workshop High-Order Numerical Approximation for Partial Differential Equations Hausdorff Center for Mathematics University of Bonn February 6-10, 2012 Organizers Alexey Chernov (University of Bonn) Christoph Schwab (ETH Zurich) Program overview Time Monday Tuesday Wednesday Thursday Friday 8:00 Registration 8:40 Opening 9:00 Stephan Sloan Ainsworth Demkowicz Buffa 10:00 Canuto Nobile Sherwin Houston Beir~ao da Veiga 10:30 Reali 11:00 Coffee break 11:30 Quarteroni Gittelson Sch¨oberl Sch¨otzau Sangalli 12:00 Bespalov 12:30 Lunch break 14:30 Dauge Xiu Wihler D¨uster Excursions: 15:00 Yosibash Arithmeum 15:30 Melenk Litvinenko Hesthaven or Coffee break 16:00 Maischak Rank Botanic 16:30 Coffee break Coffee break Gardens 17:00 Maday Hackbusch Hiptmair {18:00 Free Time 20:00 Dinner 2 Detailed program Monday, 9:00{9:55 Ernst P. Stephan hp-adaptive DG-FEM for Parabolic Obstacle Problems Monday, 10:00{10:55 Claudio Canuto Adaptive Fourier-Galerkin Methods Monday, 11:30{12:25 Alfio Quarteroni Discontinuous approximation of elastodynamics equations Monday, 14:30{15:25 Monique Dauge Weighted analytic regularity in polyhedra Monday, 15:30{16:25 Jens Markus Melenk Stability and convergence of Galerkin discretizations of the Helmholtz equation Monday, 17:00{17:55 Yvon Maday A generalized empirical interpolation method based on moments and application Tuesday, 9:00{9:55 Ian H. Sloan PDE with random coefficients as a problem in high dimensional integration Tuesday, 10:00{10:55 Fabio Nobile Stochastic Polynomial approximation of PDEs with random coefficients Tuesday,
  • Preconditioning the Coarse Problem of BDDC Methods—Three-Level, Algebraic Multigrid, and Vertex-Based Preconditioners

    Preconditioning the Coarse Problem of BDDC Methods—Three-Level, Algebraic Multigrid, and Vertex-Based Preconditioners

    Electronic Transactions on Numerical Analysis. Volume 51, pp. 432–450, 2019. ETNA Kent State University and Copyright c 2019, Kent State University. Johann Radon Institute (RICAM) ISSN 1068–9613. DOI: 10.1553/etna_vol51s432 PRECONDITIONING THE COARSE PROBLEM OF BDDC METHODS— THREE-LEVEL, ALGEBRAIC MULTIGRID, AND VERTEX-BASED PRECONDITIONERS∗ AXEL KLAWONNyz, MARTIN LANSERyz, OLIVER RHEINBACHx, AND JANINE WEBERy Abstract. A comparison of three Balancing Domain Decomposition by Constraints (BDDC) methods with an approximate coarse space solver using the same software building blocks is attempted for the first time. The comparison is made for a BDDC method with an algebraic multigrid preconditioner for the coarse problem, a three-level BDDC method, and a BDDC method with a vertex-based coarse preconditioner. It is new that all methods are presented and discussed in a common framework. Condition number bounds are provided for all approaches. All methods are implemented in a common highly parallel scalable BDDC software package based on PETSc to allow for a simple and meaningful comparison. Numerical results showing the parallel scalability are presented for the equations of linear elasticity. For the first time, this includes parallel scalability tests for a vertex-based approximate BDDC method. Key words. approximate BDDC, three-level BDDC, multilevel BDDC, vertex-based BDDC AMS subject classifications. 68W10, 65N22, 65N55, 65F08, 65F10, 65Y05 1. Introduction. During the last decade, approximate variants of the BDDC (Balanc- ing Domain Decomposition by Constraints) and FETI-DP (Finite Element Tearing and Interconnecting-Dual-Primal) methods have become popular for the solution of various linear and nonlinear partial differential equations [1,8,9, 12, 14, 15, 17, 19, 21, 24, 25].
  • Projective Synchronization of Chaotic Systems Via Backstepping Design

    Projective Synchronization of Chaotic Systems Via Backstepping Design

    International Journal of Applied Mathematical Research, 1 (4) (2012) 531-540 ©Science Publishing Corporation www.sciencepubco.com/index.php/IJBAS Projective Synchronization of Chaotic Systems Via Backstepping Design 1 Anindita Tarai (Poria), 2 Mohammad Ali Khan 1 Department of Mathematics, Aligunj R.R.B. High School Midnapore (West), West Bengal, India E-mail: [email protected] 2 Department of Mathematics, Garhbeta Ramsundar Vidyabhavan, Garhbeta, Midnapore (West), West Bengal, India E-mail: [email protected] Abstract Chaos synchronization of discrete dynamical systems are investigated. An algorithm is proposed for projective synchronization of chaotic 2D Duffing map and chaotic Tinkerbell map. The control law was derived from the Lyapunov stability theory. The numerical simulation results are presented to verify the effectiveness of the proposed algorithm Keywords: Lyapunov function, Projective synchronization, Backstepping Design, Duffing map and Tinkerbell map. 1 Introduction Adjacent chaotic trajectories governed by the same nonlinear systems may evolve into a state utterly uncorrelated, but in 1990 Pecora and Carrol [1] shown that it could be synchronized through a proper coupling. Since their seminal paper in 1990, chaos synchronization is an interesting research topic of great attention. Hayes et.al. (1993)[2] have studied to some potential applications in secure communication. Blasius et.al. [3] have observed complex dynamics and phase synchronization in spatially extended ecological systems in 1999 and system identification was investigated by Kocarev et.al. in 1996 [4]. Different forms of synchronization phenomena have been observed in a variety of chaotic systems, such as identical synchronization [1]. In 1996 Rosenblum et.al. [5] have studied phase synchronization of chaotic oscillators.