ALGORITHMS FOR THE EVOLUTION OF SURFACES
|
|
- Howard Todd Davidson
- 5 years ago
- Views:
Transcription
1 ALGORITHMS FOR THE EVOLUTION OF SURFACES MATTHEW STONE Abstract. We develop algorithms for simulating various curvature dependent motions of interfaces. These motions are described by a class of partial differential equations for which a notion of weak solution the viscosity solution has been developed that is unique and persists beyond topological changes that may occur during the evolution. A comparison principle that holds for these evolutions is the essential tool in the theory. It is also known that to establish the convergence to the viscosity solution of a consistent approximation scheme for such an evolution it is sufficient to show that it respects the comparison principle in which case it is called monotone. A familiar example that comes up in many applications is motion by mean curvature. We try to extend recent consistent and monotone numerical schemes for motion by mean curvature to evolution laws that have a more general dependence on the principal curvatures.. Introduction In [] Esedoglu Ruuth and Tsai ERT developed an algorithm for approximating the motion by certain functions of mean curvature of an interface. They adapt a similar method used by Merriman Bence and Osher MBO in [2] which alternated two steps: Convolution and simple thresholding. The MBO algorithm proceeded as follows: Let Σ R N be a region whose border denoted Σ is to be evolved via motion by mean curvature. For any time step δt > 0 the algorithm gives approximations { Σ n } to motion by mean curvature where Σ n is an approximation of the curve at time t = n δt. These approximations are generated inductively by the following with Σ 0 = Σ: Convolution: Form u : R N R as ux = G t Σn x where G t x is the N-dimensional Gaussian kernel x G t x = 2 e 4t 4πt N/2 2 Thresholding: Compute the next approximation with { Σ n+ = x : ux } 2 This algorithm while computationally inexpensive has its flaws. Most proently unless grid size is refined along with the time step the approximate motion generated by the algorithm will get stuck [][2].
2 2 MATTHEW STONE Thus Esedoglu Ruuth and Tsai replace the thresholding step with another efficient procedure: constructing the signed distance function to the interface and using formulas for approximations of the mean curvature in terms of convolutions with the Gaussian. Due to the Lipschitz continuity of signed distance functions this eliated the inaccuracies found in the MBO algorithm. This paper explores similar methods to the ERT algorithm investigating functions not only of mean curvature but of its components the principal curvatures. Due to certain problems encountered with maintaining the monotonicity of the algorithm our final goal is to approximate motion by fκ + κ 2 where f is any Lipschitz function that is increasing in each coordinate κ > κ 2 are the principal curvatures. The function x + is defined as 0 for x 0 and x otherwise. Similarly x is defined as 0 for x 0 and x otherwise... Things that need to be done. Prove theorems. 2. Expansions for the distance function In this section we first write down the Taylor expansion of the signed distance function d Σ x along a unital direction v in the neighborhood of a point p Σ on the smooth boundary Σ of a set Σ. We work in R 3 where we write x = x y z. This expansion allows us to obtain a Taylor expansion for the convolution of d Σ x y z with the Gaussian kernel G t v x y z which is defined as: G t v x = x v 2 e 4t. 4πt N/2 This will allow us to express the expansion coefficients in terms of the second fundamental form and thus using the convolution for every direction v will allow us to express the principle curvatures of Σ. 2.. Expansion for a smooth interface. First let us recall a few well known properties of the signed distance function. The following will be taken directly from []. For the following d s will denote dσ s d xs will denote 2 d Σ x s and so on. The taylor expansion at s = 0 is given by: 2. dx y z s = dx y z + s v = dx y z + d s x y z 0s + 2 d ssx y z 0s 2 + d xs x y z 0xs + d ys x y z 0ys + d zs x y z 0zs + 6 d sssx y z 0s d ssxx y z 0s 2 x + 2 d ssyx y z 0s 2 y + 2 d sszx y z 0s 2 z + d sxy x y z 0sxy + d sxz x y z 0sxz + d syz x y z 0syz + higher order terms We can substitute the expansion into the convolution integral G t v x y z σdx y z s σdσ R
3 ALGORITHMS FOR THE EVOLUTION OF SURFACES 3 to get a Taylor expansion for the convolution G t v dx y z s at x y z s = The terms we need are: c G t v x y z 0 = c s G t v x y z 0 = 0 s 2 G t v x y z 0 = 2t s 3 G t v x y z 0 = 0 Using these we arrive at the following expansion: Proposition 2.. Convolution of the signed distance function d with the Gaussian kernel G t v has the following expansion 2.2 d G t v x y z 0 = dx y z + 2 d ssx y z 0 2t + 2 d ssxx y z 0 2tx + 2 d ssyx y z 0 2ty + 2 d sszx y z 0 2tz + higher order terms Provided that each of x y z are Ot we get 2.3 d G t v x y z 0 = dx y z + d ss x y z 0t + Ot 2 Thus we find By its definition Finally giving G t v d d = d ss x y z 0t + Ot 2 d ss x y z 0 = D 2 d v v. G t v d d = D 2 d v v t + Ot 2 For any signed distance function d the eigenvalues of D 2 d are given by the principle curvatures with eigenvectors along the principle directions plus an eigenvalue of 0 with corresponding eigenvector of the normal to the surface. Thus v = D 2 d v v that is the greatest eigenvalue is either κ or 0 depending which is larger. That is v = D 2 d v v. Similarly v = D 2 d v v = κ 2. This leads to the conclusions that: 2.4 v = G t v d d = κ + t + Ot v = G t v d d = κ 2 t + Ot 2 Using these two equations for any Lipschitz function f 2.6 f t G t v d d v = t G t v d d = fκ + κ 2 + Ot v = This leads to the following algorithm:
4 4 MATTHEW STONE Algorithm 2.2. Given the initial set Σ 0 through its signed distance function d 0 x and a time step δt > 0 generate the sets Σ j via their signed distance functions d j x at subsequent discrete times t = jδt by alternating the following steps: Form the function 2.7 Ax := d j +δtf Mδt G Mδt v d d v = 2 Construct distance function d j+ by 2.8 d j+ x = RedistA. Mδt G Mδt v d d v = At the jth step of the algorithm the set Σ j can be recovered through the relation Σ j = {x : d j x > 0} So using Equation 2.6 gives the 0-level set of which has moved with the desired speed proving the consistency of the algorithm. Finally we must show the monotonicity of our algorithm assug f is Lipschitz and increasing: Proposition 2.3. If M L f then algorithm 2.7 and 2.8 is monotone for any choice of time step size δt > 0. Proof. Let Σ Σ 2 be two sets satisfying Σ Σ 2 and d x d 2 x be signed distance functions to Σ Σ 2 respectively. Firstly we have x d x d 2 x Using the same notation as in the algorithm let A x = d x+δtf A 2 x = d 2 x+δtf Then we calculate: 2.9 A 2 A = d 2 d + δt f Mδt G Mδt v d d v = Mδt G Mδt v d d v = Mδt G Mδt v d 2 d 2 v = Mδt G Mδt v d 2 d 2 v = [ f Mδt G Mδt v d 2 d 2 v = Mδt G Mδt v d d v = Mδt G Mδt v d 2 d 2 v = Mδt G Mδt v d d v = ]
5 ALGORITHMS FOR THE EVOLUTION OF SURFACES d 2 d v G Mδt v d 2 G Mδt v d G Mδt v d 2 d 2 G Mδt v d d and v = v = G Mδt v d 2 d 2 G Mδt v d d v = v = Substituting 2.0 into 2.9 using the fact that f is increasing in both variables: 2. A 2 A d 2 d + δt f [ f Mδt G Mδt v d 2 d 2 v = Mδt G Mδt v d 2 d v = Mδt G Mδt v d 2 d 2 v = Mδt G Mδt v d 2 d v = ] Using the fact that f is Lipschitz with constant L f and M L f : 2.2 A 2 A d 2 d δt 2Lf d 2 d 0 Mδt Thus we have that our algorithm is both consistent and monotone sufficient conditions to show the convergence of the algorithm to the desired motion as δt 0 References [] S. Esedoglu S. Ruuth and R. Tsai Diffusion generated motion using signed distance functions Journal of Computational Physics [2] B. Merriman J.K. Bence and S. Osher Motion of multiple junctions: a level set approach Journal of Computational Physics Department of Mathematics Yale University New Haven CT address: matthew.i.stone@yale.edu
A generalized MBO diffusion generated motion for constrained harmonic maps
A generalized MBO diffusion generated motion for constrained harmonic maps Dong Wang Department of Mathematics, University of Utah Joint work with Braxton Osting (U. Utah) Workshop on Modeling and Simulation
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationAlgorithms for Multiphase Partitioning
Algorithms for Multiphase Partitioning by Matthew Jacobs A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Mathematics) in the University of Michigan
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationCorner detection: the basic idea
Corner detection: the basic idea At a corner, shifting a window in any direction should give a large change in intensity flat region: no change in all directions edge : no change along the edge direction
More informationINTEREST POINTS AT DIFFERENT SCALES
INTEREST POINTS AT DIFFERENT SCALES Thank you for the slides. They come mostly from the following sources. Dan Huttenlocher Cornell U David Lowe U. of British Columbia Martial Hebert CMU Intuitively, junctions
More informationIntroduction to gradient descent
6-1: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction to gradient descent Derivation and intuitions Hessian 6-2: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction Our
More informationAn improved threshold dynamics method for wetting dynamics
An improved threshold dynamics method for wetting dynamics Dong Wang a, Xiao-Ping Wang b,, Xianmin Xu c,d a Department of Mathematics, University of Utah, Salt Lake City, UT, 8, USA. b Department of Mathematics,
More informationClass Meeting # 2: The Diffusion (aka Heat) Equation
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 2: The Diffusion (aka Heat) Equation The heat equation for a function u(, x (.0.). Introduction
More informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
More informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More informationMax-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig Threshold Dynamics for Networks with Arbitrary Surface Tensions (revised version: January 203) by Selim Esedoglu and Felix Otto Preprint
More informationSTAT 200C: High-dimensional Statistics
STAT 200C: High-dimensional Statistics Arash A. Amini May 30, 2018 1 / 57 Table of Contents 1 Sparse linear models Basis Pursuit and restricted null space property Sufficient conditions for RNS 2 / 57
More informationFIXED POINT ITERATIONS
FIXED POINT ITERATIONS MARKUS GRASMAIR 1. Fixed Point Iteration for Non-linear Equations Our goal is the solution of an equation (1) F (x) = 0, where F : R n R n is a continuous vector valued mapping in
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms
ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:
More informationEcon Lecture 14. Outline
Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations
More informationLinear Algebra Practice Problems
Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,
More informationECOLE POLYTECHNIQUE. Approximation of the anisotropic mean curvature ow CENTRE DE MATHÉMATIQUES APPLIQUÉES UMR CNRS A. Chambolle, M.
ECOLE POLYTECHNIQUE CENTRE DE MATHÉMATIQUES APPLIQUÉES UMR CNRS 7641 91128 PALAISEAU CEDEX (FRANCE). Tél: 01 69 33 41 50. Fax: 01 69 33 30 11 http://www.cmap.polytechnique.fr/ Approximation of the anisotropic
More informationGood Problems. Math 641
Math 641 Good Problems Questions get two ratings: A number which is relevance to the course material, a measure of how much I expect you to be prepared to do such a problem on the exam. 3 means of course
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationIntroduction to Vector Functions
Introduction to Vector Functions Limits and Continuity Philippe B. Laval KSU Today Philippe B. Laval (KSU) Vector Functions Today 1 / 14 Introduction Until now, the functions we studied took a real number
More informationHeat Source Identification Based on L 1 Optimization
Heat Source Identification Based on L 1 Optimization Yingying Li, Stanley Osher and Richard Tsai August 27, 2009 Yingying Li, Stanley Osher and Richard Tsai () L 1 Heat Source Identification Aug 27, 2009
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationarxiv: v1 [math.oc] 22 Dec 2018
AN EFFICIENT THRESHOLD DYNAMICS METHOD FOR TOPOLOGY OPTIMIZATION FOR FLUIDS HUANGXIN CHEN, HAITAO LENG, DONG WANG, AND XIAO-PING WANG arxiv:1812.09437v1 [math.oc] 22 Dec 2018 Abstract. We propose an efficient
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationSingular Diffusion Equations With Nonuniform Driving Force. Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009
Singular Diffusion Equations With Nonuniform Driving Force Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009 1 Contents 0. Introduction 1. Typical Problems 2. Variational Characterization
More informationEdge Detection. CS 650: Computer Vision
CS 650: Computer Vision Edges and Gradients Edge: local indication of an object transition Edge detection: local operators that find edges (usually involves convolution) Local intensity transitions are
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationEvolution of semiclassical Wigner function (the higher dimensio
Evolution of semiclassical Wigner function (the higher dimensional case) Workshop on Fast Computations in Phase Space, WPI-Vienna, November 2008 Dept. Appl. Math., Univ. Crete & IACM-FORTH 1 2 3 4 5 6
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationSeptember Math Course: First Order Derivative
September Math Course: First Order Derivative Arina Nikandrova Functions Function y = f (x), where x is either be a scalar or a vector of several variables (x,..., x n ), can be thought of as a rule which
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationLecture 10: October 27, 2016
Mathematical Toolkit Autumn 206 Lecturer: Madhur Tulsiani Lecture 0: October 27, 206 The conjugate gradient method In the last lecture we saw the steepest descent or gradient descent method for finding
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationEcon 204 Differential Equations. 1 Existence and Uniqueness of Solutions
Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationNonlinear Diffusion. 1 Introduction: Motivation for non-standard diffusion
Nonlinear Diffusion These notes summarize the way I present this material, for my benefit. But everything in here is said in more detail, and better, in Weickert s paper. 1 Introduction: Motivation for
More informationWhat is a Space Curve?
What is a Space Curve? A space curve is a smooth map γ : I R R 3. In our analysis of defining the curvature for space curves we will be able to take the inclusion (γ, 0) and have that the curvature of
More informationNonparametric Density Estimation
Nonparametric Density Estimation Econ 690 Purdue University Justin L. Tobias (Purdue) Nonparametric Density Estimation 1 / 29 Density Estimation Suppose that you had some data, say on wages, and you wanted
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationNonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University
Nonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University this presentation derived from that presented at the Pan-American Advanced
More informationPostprint.
http://www.diva-portal.org Postprint This is the accepted version of a chapter published in Domain Decomposition Methods in Science and Engineering XXI. Citation for the original published chapter: Gander,
More informationCAAM 336 DIFFERENTIAL EQUATIONS IN SCI AND ENG Examination 1
CAAM 6 DIFFERENTIAL EQUATIONS IN SCI AND ENG Examination Instructions: Time limit: uninterrupted hours There are four questions worth a total of 5 points Please do not look at the questions until you begin
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationTHE problem of simulating the motion of interfaces with
Simulation of Triple Junction Motion ith Arbitrary Surface Tensions Nur Shofianah, Rhudaina Z. Mohammad, Karel Svadlenka Abstract We simulate triple junction motion that is given by the gradient flo of
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationTHE INVERSE FUNCTION THEOREM FOR LIPSCHITZ MAPS
THE INVERSE FUNCTION THEOREM FOR LIPSCHITZ MAPS RALPH HOWARD DEPARTMENT OF MATHEMATICS UNIVERSITY OF SOUTH CAROLINA COLUMBIA, S.C. 29208, USA HOWARD@MATH.SC.EDU Abstract. This is an edited version of a
More informationSurface Evolution under Curvature Flows
Journal of Visual Communication and Image Representation 13, 65 81 (2002) doi:10.1006/jvci.2001.0476, available online at http://www.idealibrary.com on Surface Evolution under Curvature Flows Conglin Lu,
More informationMath 5588 Final Exam Solutions
Math 5588 Final Exam Solutions Prof. Jeff Calder May 9, 2017 1. Find the function u : [0, 1] R that minimizes I(u) = subject to u(0) = 0 and u(1) = 1. 1 0 e u(x) u (x) + u (x) 2 dx, Solution. Since the
More informationLinear ODE s with periodic coefficients
Linear ODE s with periodic coefficients 1 Examples y = sin(t)y, solutions Ce cos t. Periodic, go to 0 as t +. y = 2 sin 2 (t)y, solutions Ce t sin(2t)/2. Not periodic, go to to 0 as t +. y = (1 + sin(t))y,
More informationE[X n ]= dn dt n M X(t). ). What is the mgf? Solution. Found this the other day in the Kernel matching exercise: 1 M X (t) =
Chapter 7 Generating functions Definition 7.. Let X be a random variable. The moment generating function is given by M X (t) =E[e tx ], provided that the expectation exists for t in some neighborhood of
More informationLECTURE 7. k=1 (, v k)u k. Moreover r
LECTURE 7 Finite rank operators Definition. T is said to be of rank r (r < ) if dim T(H) = r. The class of operators of rank r is denoted by K r and K := r K r. Theorem 1. T K r iff T K r. Proof. Let T
More informationA finite element level set method for anisotropic mean curvature flow with space dependent weight
A finite element level set method for anisotropic mean curvature flow with space dependent weight Klaus Deckelnick and Gerhard Dziuk Centre for Mathematical Analysis and Its Applications, School of Mathematical
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationTopics in Curve and Surface Implicitization
Topics in Curve and Surface Implicitization David A. Cox Amherst College PASI 2009 p.1/41 Outline Curves: Moving Lines & µ-bases Moving Curve Ideal & the Rees Algebra Adjoint Curves PASI 2009 p.2/41 Outline
More informationChapter 5. The Second Fundamental Form
Chapter 5. The Second Fundamental Form Directional Derivatives in IR 3. Let f : U IR 3 IR be a smooth function defined on an open subset of IR 3. Fix p U and X T p IR 3. The directional derivative of f
More informationWorksheet 9. Topics: Taylor series; using Taylor polynomials for approximate computations. Polar coordinates.
ATH 57H Spring 0 Worksheet 9 Topics: Taylor series; using Taylor polynomials for approximate computations. Polar coordinates.. Let f(x) = +x. Find f (00) (0) - the 00th derivative of f at point x = 0.
More informationElements of Linear Algebra, Topology, and Calculus
Appendix A Elements of Linear Algebra, Topology, and Calculus A.1 LINEAR ALGEBRA We follow the usual conventions of matrix computations. R n p is the set of all n p real matrices (m rows and p columns).
More information3 Compact Operators, Generalized Inverse, Best- Approximate Solution
3 Compact Operators, Generalized Inverse, Best- Approximate Solution As we have already heard in the lecture a mathematical problem is well - posed in the sense of Hadamard if the following properties
More informationStructurally Stable Singularities for a Nonlinear Wave Equation
Structurally Stable Singularities for a Nonlinear Wave Equation Alberto Bressan, Tao Huang, and Fang Yu Department of Mathematics, Penn State University University Park, Pa. 1682, U.S.A. e-mails: bressan@math.psu.edu,
More informationLearning Eigenfunctions: Links with Spectral Clustering and Kernel PCA
Learning Eigenfunctions: Links with Spectral Clustering and Kernel PCA Yoshua Bengio Pascal Vincent Jean-François Paiement University of Montreal April 2, Snowbird Learning 2003 Learning Modal Structures
More informationExamples of numerics in commutative algebra and algebraic geo
Examples of numerics in commutative algebra and algebraic geometry MCAAG - JuanFest Colorado State University May 16, 2016 Portions of this talk include joint work with: Sandra Di Rocco David Eklund Michael
More informationRemarks on Extremization Problems Related To Young s Inequality
Remarks on Extremization Problems Related To Young s Inequality Michael Christ University of California, Berkeley University of Wisconsin May 18, 2016 Part 1: Introduction Young s convolution inequality
More information6. Brownian Motion. Q(A) = P [ ω : x(, ω) A )
6. Brownian Motion. stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ, P) and a real valued stochastic process can be defined
More informationDifferential Geometry and Lie Groups with Applications to Medical Imaging, Computer Vision and Geometric Modeling CIS610, Spring 2008
Differential Geometry and Lie Groups with Applications to Medical Imaging, Computer Vision and Geometric Modeling CIS610, Spring 2008 Jean Gallier Department of Computer and Information Science University
More informationUnit Speed Curves. Recall that a curve Α is said to be a unit speed curve if
Unit Speed Curves Recall that a curve Α is said to be a unit speed curve if The reason that we like unit speed curves that the parameter t is equal to arc length; i.e. the value of t tells us how far along
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationThe Fourier Transform Method
The Fourier Transform Method R. C. Trinity University Partial Differential Equations April 22, 2014 Recall The Fourier transform The Fourier transform of a piecewise smooth f L 1 (R) is ˆf(ω) = F(f)(ω)
More informationLecture 9 Metric spaces. The contraction fixed point theorem. The implicit function theorem. The existence of solutions to differenti. equations.
Lecture 9 Metric spaces. The contraction fixed point theorem. The implicit function theorem. The existence of solutions to differential equations. 1 Metric spaces 2 Completeness and completion. 3 The contraction
More informationIntegration over curves and surfaces defined by the closest point mapping
Integration over curves and surfaces defined by the closest point mapping Catherine Kublik and Richard Tsai Abstract We propose a new formulation using the closest point mapping for integrating over smooth
More informationLecture 1: Systems of linear equations and their solutions
Lecture 1: Systems of linear equations and their solutions Course overview Topics to be covered this semester: Systems of linear equations and Gaussian elimination: Solving linear equations and applications
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationSpectral Regularization
Spectral Regularization Lorenzo Rosasco 9.520 Class 07 February 27, 2008 About this class Goal To discuss how a class of regularization methods originally designed for solving ill-posed inverse problems,
More informationEconomics 204 Summer/Fall 2010 Lecture 10 Friday August 6, 2010
Economics 204 Summer/Fall 2010 Lecture 10 Friday August 6, 2010 Diagonalization of Symmetric Real Matrices (from Handout Definition 1 Let δ ij = { 1 if i = j 0 if i j A basis V = {v 1,..., v n } of R n
More informationAnalysis II: The Implicit and Inverse Function Theorems
Analysis II: The Implicit and Inverse Function Theorems Jesse Ratzkin November 17, 2009 Let f : R n R m be C 1. When is the zero set Z = {x R n : f(x) = 0} the graph of another function? When is Z nicely
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationExamination paper for TMA4145 Linear Methods
Department of Mathematical Sciences Examination paper for TMA4145 Linear Methods Academic contact during examination: Franz Luef Phone: 40614405 Examination date: 5.1.016 Examination time (from to): 09:00-13:00
More informationStatistical Pattern Recognition
Statistical Pattern Recognition Feature Extraction Hamid R. Rabiee Jafar Muhammadi, Alireza Ghasemi, Payam Siyari Spring 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Agenda Dimensionality Reduction
More informationMathematica Bohemica
Mathematica Bohemica Miroslav Kolář; Michal Beneš; Daniel Ševčovič Computational studies of conserved mean-curvature flow Mathematica Bohemica, Vol. 39 (4), No. 4, 677--684 Persistent URL: http://dml.cz/dmlcz/4444
More informationStatistics 3657 : Moment Generating Functions
Statistics 3657 : Moment Generating Functions A useful tool for studying sums of independent random variables is generating functions. course we consider moment generating functions. In this Definition
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More information10. Smooth Varieties. 82 Andreas Gathmann
82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It
More informationAuction Dynamics: A Volume Constrained MBO Scheme
Auction Dynamics: A Volume Constrained MBO Scheme Matt Jacobs Ekaterina Merkurjev Selim Esedoḡlu October 22, 2017 Abstract We show how auction algorithms, originally developed for the assigment problem,
More information1 The Differential Geometry of Surfaces
1 The Differential Geometry of Surfaces Three-dimensional objects are bounded by surfaces. This section reviews some of the basic definitions and concepts relating to the geometry of smooth surfaces. 1.1
More information1 Matrices and Systems of Linear Equations. a 1n a 2n
March 31, 2013 16-1 16. Systems of Linear Equations 1 Matrices and Systems of Linear Equations An m n matrix is an array A = (a ij ) of the form a 11 a 21 a m1 a 1n a 2n... a mn where each a ij is a real
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 7 I. A SIMPLE EXAMPLE OF ANGULAR MOMENTUM ADDITION Given two spin-/ angular momenta, S and S, we define S S S The problem is to find the eigenstates of the
More informationNon-Convex Optimization. CS6787 Lecture 7 Fall 2017
Non-Convex Optimization CS6787 Lecture 7 Fall 2017 First some words about grading I sent out a bunch of grades on the course management system Everyone should have all their grades in Not including paper
More informationAHAHA: Preliminary results on p-adic groups and their representations.
AHAHA: Preliminary results on p-adic groups and their representations. Nate Harman September 16, 2014 1 Introduction and motivation Let k be a locally compact non-discrete field with non-archimedean valuation
More informationLECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationRegularization via Spectral Filtering
Regularization via Spectral Filtering Lorenzo Rosasco MIT, 9.520 Class 7 About this class Goal To discuss how a class of regularization methods originally designed for solving ill-posed inverse problems,
More informationEstimates for the density of functionals of SDE s with irregular drift
Estimates for the density of functionals of SDE s with irregular drift Arturo KOHATSU-HIGA a, Azmi MAKHLOUF a, a Ritsumeikan University and Japan Science and Technology Agency, Japan Abstract We obtain
More informationMarkov Chain Monte Carlo Inference. Siamak Ravanbakhsh Winter 2018
Graphical Models Markov Chain Monte Carlo Inference Siamak Ravanbakhsh Winter 2018 Learning objectives Markov chains the idea behind Markov Chain Monte Carlo (MCMC) two important examples: Gibbs sampling
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More information