The Fast Multipole Method (FMM)
|
|
- Alaina Snow
- 5 years ago
- Views:
Transcription
1 The Fast Multipole Method (FMM)
2 Motivation for FMM Computational Physics Problems involving mutual interactions of N particles Gravitational or Electrostatic forces Collective (but weak) long-range forces cannot be ignored Example Application DNA, proteins, biomollecular processes 200k to 1M atoms History 1977 Pincus, Scheraga Group of distant particles replaced with a Pseudoparticle 1987 Greengard, Rokhlin Rigorous numerical analysis: how many terms are needed Local expansion: Groups interact with groups 1992 Rokhlin: Applied to Helmholz Operator (MoM) FMM Slide 2
3 Need for Fast Algorithms in EM Basic Problem in EM Everything is coupled to everything Problem complexity tends to scale as N 2, where N is the number of unknowns Large EM Problems Radar, Remote Sensing (Large scatterers) Wireless Communications (Environment) Antenna/environment studies (body, vehicles) Complete RF system simulation FMM Slide 3
4 Example of a Fast Algorithm in SP Example: DFT Computation: Complexity: O(N 2 ) Alternative: FFT Partition (odd/even) Recursion (Reuse information) Complexity: O(NlogN) FMM Slide 4
5 Fast Multipole Algorithm Method of Moments (MoM) Get system of the form v = Z i solution i = Z -1 v Coupling matrix Z is a dense NxN matrix Computing inverse O(N 2 ) Fast Multipole Method Exploit structure of Z fast algorithm to compute Z i Combined with conjugate gradient method for inverse Obtain O(N 3/2 ) algorithm Can also do multistage: O(NlogN) algorithm Follow development in Coifman, et. al., The Fast Multipole Method for the Wave Equation: A Pedestrian Prescription, IEEE Antennas and Propagation Magazine FMM Slide 5
6 MoM Problem Scalar wave equation Boundary condition S Application of MoM yields FMM Slide 6
7 Computation of Inverse Small systems (N 100) Larger systems Compute inverse iteratively (conjugate gradient) Repeatedly solve Problem: Matrix vector multiply is O(N 2 ) Can we accelerate this? FMM Slide 7
8 Airports Problem: Connect all airports with flights Option 1: Connect all cities to all cities Problem: Inefficient Bridgeport New York Newark Atlantic City Philadelphia Wilmington Washington Richmond Virgina Beach FMM Slide 8
9 Airports Problem: Connect all airports with flights Option 1: Connect all cities to all cities Problem: Inefficient Bridgeport New York Newark Atlantic City Philadelphia Wilmington Washington Richmond Virgina Beach FMM Slide 9
10 Airports Problem: Connect all airports with flights Option 1: Connect all cities to all cities Problem: Inefficient Option 2: Use hubs All cities served Far fewer connections to deal with! Newark New York Atlantic City Richmond Bridgeport Philadelphia Wilmington Washington Virgina Beach FMM Slide 10
11 Airports Problem: Connect all airports with flights Option 1: Connect all cities to all cities Problem: Inefficient Option 2: Use hubs All cities served Far fewer connections to deal with! Newark New York Atlantic City Richmond Bridgeport Philadelphia Wilmington Washington Virgina Beach FMM Slide 11
12 FFT Analogy Revisited Transformation Matrix Can generate sparse factorization Exploits algebraic properties (exact) More efficient to compute y = W x with FMM Slide 12
13 Sparse Factorization of Z Z can be decomposed as All sparse matrices Exploits analytic properties of Z Approximate, but can control level of approximation (down to precision of computer) Result: O(N 3/2 ) method to compute product of Z with an arbitrary vector Can repeat decomposition (approach O(NlogN)) Not really necessary FMM Slide 13
14 Factoring Z: Identity 1 Consider the Green s function Spherical Bessel Function Legendre Polynomial Spherical Hankel Function Transformed spherical wave from x to x Sum of product of spherical waves Separate long interaction (X) from near interactions (d) FMM Slide 14
15 Factoring Z: Identity 2 Expand spherical factor in d as Substituting, Write j(.) as sum of plane waves (all directions on unit sphere) FMM Slide 15
16 Meaning of Factorization Superimpose all plane waves Phase shift required due to small offset d Complex amplitude of plane wave traveling in k direction from x to x +X Fields at x +X are now a superposition of plane waves Fields at a small distance d from x +X found by shifting the phase of the plane waves Note: L > kd for good accuracy (d should be kept small) FMM Slide 16
17 Back to MoM Result so far Note: Substitute into coupling matrix Interpretation #1 FMM Slide 17
18 Back to MoM Result so far Note: Substitute into coupling matrix Interpretation #1 Source Basis Function (a source with local support) FMM Slide 18
19 Back to MoM Result so far Note: Substitute into coupling matrix Spatial Fourier Transform in x => Transform source to plane waves (k is pw dir) Interpretation #1 FMM Slide 19
20 Back to MoM Result so far Note: Substitute into coupling matrix Interpretation #1 Complex Transfer Coefficient: Attenuation/ phase shift of pw in k dir from x to x +X FMM Slide 20
21 Back to MoM Result so far Note: Substitute into coupling matrix Interpretation #1 Inverse Fourier Transform: pw back to fields in observation space FMM Slide 21
22 Back to MoM Result so far Note: Substitute into coupling matrix Project onto testing functions (Galerkin) Interpretation #1 FMM Slide 22
23 Back to MoM Interpretation #2 Rearrange Let X=X m -X m : X m is coordinate near local support of f n (x) FT of nth weighting function (plane waves received by f n (x) One-to-one coupling of pw s from source to observation FT of n th basis function (plane waves radiated by f n (x) Fields due to current i FMM Slide 23
24 Hub Interpretation Old MoM (far interactions) Apply current to each source basis function Observe field on each weighting function N 2 operations f n (x) V = Z i f n (x ) FMM Slide 24
25 Hub Interpretation Using Factorization Use pw basis functions: Currents to plane waves For single hub, plane waves superimpose 1-to-1 connection of K plane waves from source to dest. At destination, use pw weighting functions Transform from pw to observed fields (integrate over k) V n (k) V n (k) X m X m FMM Slide 25
26 Algorithm Setup 1/4 Divide N basis functions into M groups: M N 1/2 α=1 α=3 α=2 α=4 α=1 α=3 α=2 α=1 α=2 X 1 X 2 X 3 α=4 α=3 α=4 α=5 α=6 α=5 α=6 α=5 α=6 m=1 m=2 m=3 m = Group number α = Local element number n(m, α) = Global element number (like with normal MoM) X m = coordinate at center of mth group FMM Slide 26
27 Algorithm Setup 2/4 Close Interactions Require many terms L in expansion Relatively few close elements for large-scale problem Just use normal dense computation of Z for these For group pairs (m, m ) that are close Close means min{dist[n(m,α),n(m,α )]} = 2π/k Z mαm α is computed using Here we have to handle singularities of Green s function Requires O(N 3/2 ) computations FMM Slide 27
28 Algorithm Setup 3/4 Discretize Directions Form set of K directions Accuracy requires K L 2 > (kd) 2 D = max(d) Compute excitation vectors Fourier transform of basis/weighting functions From partitioning, kd on order of (N/M) 1/2 Requies O(N 2 /M) computations FMM Slide 28
29 Algorithm Setup 4/4 Regions that are not nearby i.e. Z mαm α = 0 Compute the matrix elements FMM Slide 29
30 Fast Matrix-Vector Multiply 1/3 Need to compute Compute KM quantities FT current in mα => plane waves Superimpose => plane waves from mth element Computation: O(N 2 /M) FMM Slide 30
31 Fast Matrix-Vector Multiply 2/3 Compute KM quantities FMM Slide 31
32 Fast Matrix-Vector Multiply 2/3 Compute KM quantities Coupling of plane waves from group m to group m Note: 1-to-1 coupling in k! FMM Slide 32
33 Fast Matrix-Vector Multiply 2/3 Compute KM quantities Observed plane-waves in group m Coupling of plane waves from group m to group m Note: 1-to-1 coupling in k! FMM Slide 33
34 Fast Matrix-Vector Multiply 2/3 Compute KM quantities Observed plane-waves in group m Coupling of plane waves from group m to group m Note: 1-to-1 coupling in k! Note: Just for distant terms Requires O(MN) operations FMM Slide 34
35 Fast Matrix-Vector Multiply 3/3 Put it all together Dense terms Transform pw s in group m back to fields Requires O(N 2 /M) operations FMM Slide 35
36 Conclusion Large-scale Problems in Physics/Engineering Problems involving mutual interactions of N particles Classical methods often lead to O(N 2 ) algorithms Fast Multipole Method Separate near and far interactions (hub analogy) Near interactions captured with exact (dense) operator Far interactions Collect particle interactions into groups Relate groups with each other (not individual particles) Application to MoM in EM Reduced an O(N 2 ) problem to O(N 3/2 ) Significant savings for large scale problems FMM Slide 36
Iterative methods for use with the Fast Multipole Method
Iterative methods for use with the Fast Multipole Method Ramani Duraiswami Perceptual Interfaces and Reality Lab. Computer Science & UMIACS University of Maryland, College Park, MD Joint work with Nail
More informationFast Multipole Accelerated Indirect Boundary Elements for the Helmholtz Equation
Fast Multipole Accelerated Indirect Boundary Elements for the Helmholtz Equation Nail A. Gumerov Ross Adelman Ramani Duraiswami University of Maryland Institute for Advanced Computer Studies and Fantalgo,
More informationAdvanced Surface Based MoM Techniques for Packaging and Interconnect Analysis
Electrical Interconnect and Packaging Advanced Surface Based MoM Techniques for Packaging and Interconnect Analysis Jason Morsey Barry Rubin, Lijun Jiang, Lon Eisenberg, Alina Deutsch Introduction Fast
More informationFast Spherical Filtering in the Broadband FMBEM using a nonequally
Fast Spherical Filtering in the Broadband FMBEM using a nonequally spaced FFT Daniel R. Wilkes (1) and Alec. J. Duncan (1) (1) Centre for Marine Science and Technology, Department of Imaging and Applied
More informationLecture 7: Introduction to HFSS-IE
Lecture 7: Introduction to HFSS-IE 2015.0 Release ANSYS HFSS for Antenna Design 1 2015 ANSYS, Inc. HFSS-IE: Integral Equation Solver Introduction HFSS-IE: Technology An Integral Equation solver technology
More informationArray Decomposition-Fast Multipole Method for finite array analysis
RADIO SCIENCE, VOL 39,, doi:10109/003rs00887, 004 Array Decomposition-Fast Multipole Method for finite array analysis Rick W Kindt 1 and John L Volakis 1 ElectroScience Laboratory, Department of Electrical
More informationEnhanced Characteristic Basis Function Method for Solving the Monostatic Radar Cross Section of Conducting Targets
Progress In Electromagnetics Research M, Vol. 68, 173 180, 2018 Enhanced Characteristic Basis Function Method for Solving the Monostatic Radar Cross Section of Conducting Targets Jinyu Zhu, Yufa Sun *,
More informationA Kernel-independent Adaptive Fast Multipole Method
A Kernel-independent Adaptive Fast Multipole Method Lexing Ying Caltech Joint work with George Biros and Denis Zorin Problem Statement Given G an elliptic PDE kernel, e.g. {x i } points in {φ i } charges
More informationComparison of iteration convergences of SIE and VSIE for solving electromagnetic scattering problems for coated objects
RADIO SCIENCE, VOL. 38, NO. 2, 1028, doi:10.1029/2002rs002610, 2003 Comparison of iteration convergences of SIE and VSIE for solving electromagnetic scattering problems for coated objects Cai-Cheng Lu
More informationEfficient O(N log N) algorithms for scattered data interpolation
Efficient O(N log N) algorithms for scattered data interpolation Nail Gumerov University of Maryland Institute for Advanced Computer Studies Joint work with Ramani Duraiswami February Fourier Talks 2007
More informationMulti-Domain Pattern. I. Problem. II. Driving Forces. III. Solution
Multi-Domain Pattern I. Problem The problem represents computations characterized by an underlying system of mathematical equations, often simulating behaviors of physical objects through discrete time
More informationFast multipole methods for axisymmetric geometries
Fast multipole methods for axisymmetric geometries Victor Churchill New York University May 13, 2016 Abstract Fast multipole methods (FMMs) are one of the main numerical methods used in the solution of
More informationFast Algorithm for Matrix-Vector Multiply of Asymmetric Multilevel Block-Toeplitz Matrices
" Fast Algorithm for Matrix-Vector Multiply of Asymmetric Multilevel Block-Toeplitz Matrices B. E. Barrowes, F. L. Teixeira, and J. A. Kong Research Laboratory of Electronics, MIT, Cambridge, MA 02139-4307
More informationFast Multipole and Related Algorithms
Fast Multipole and Related Algorithms Ramani Duraiswami University of Maryland, College Park http://www.umiacs.umd.edu/~ramani Joint work with Nail A. Gumerov Efficiency by exploiting symmetry and A general
More informationComputational Fluid Dynamics - Incompressible Flows
Computational Fluid Dynamics - Incompressible Flows March 25, 2008 Incompressible Flows Basis Functions Discrete Equations CFD - Incompressible Flows CFD is a Huge field Numerical Techniques for solving
More informationCOSC6365. Introduction to HPC. Lecture 21. Lennart Johnsson Department of Computer Science
Introduction to HPC Lecture 21 Department of Computer Science Most slides from UC Berkeley CS 267 Spring 2011, Lecture 12, Dense Linear Algebra (part 2), Parallel Gaussian Elimination. Jim Demmel Dense
More informationCan Xia 1, 2,WanqingYou 2, and Yufa Sun 2, *
Progress In Electromagnetics Research Letters, Vol. 81, 133 139, 2019 Fast Calculation of Monostatic Radar Cross Section of Conducting Targets Using Hierarchical Characteristic Basis Function Method and
More informationFMM accelerated BEM for 3D Helmholtz equation
FMM accelerated BEM for 3D Helmholtz equation Nail A. Gumerov and Ramani Duraiswami Institute for Advanced Computer Studies University of Maryland, U.S.A. also @ Fantalgo, LLC, U.S.A. www.umiacs.umd.edu/~gumerov
More informationAcceleration of the Multi Level Fast Multipole Method on a GPU
Thesis National Aerospace Laboratory NLR Acceleration of the Multi Level Fast Multipole Method on a GPU Author: Claudia Wagenaar(4090772) Supervisors: H. van der Ven D. van der Heul C. Vuik February 1,
More informationEfficient computation of source magnetic scalar potential
Adv. Radio Sci., 4, 59 63, 2006 Author(s) 2006. This work is licensed under a Creative Commons License. Advances in Radio Science Efficient computation of source magnetic scalar potential W. Hafla, A.
More informationGPU accelerated heterogeneous computing for Particle/FMM Approaches and for Acoustic Imaging
GPU accelerated heterogeneous computing for Particle/FMM Approaches and for Acoustic Imaging Ramani Duraiswami University of Maryland, College Park http://www.umiacs.umd.edu/~ramani With Nail A. Gumerov,
More informationLocal Multilevel Fast Multipole Algorithm for 3D Electromagnetic Scattering
Progress In Electromagnetics Research Symposium 2005, Hangzhou, China, August 22-26 745 Local Multilevel Fast Multipole Algorithm for 3D Electromagnetic Scattering Jun Hu, Zaiping Nie, Lin Lei, and Jun
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationSpherical Microphone Arrays
Spherical Microphone Arrays Acoustic Wave Equation Helmholtz Equation Assuming the solutions of wave equation are time harmonic waves of frequency ω satisfies the homogeneous Helmholtz equation: Boundary
More informationImage Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi
Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi 1. Introduction The choice of a particular transform in a given application depends on the amount of
More informationEM Software & Systems GmbH
EM Software & Systems GmbH Otto-Lilienthal-Straße 36 D-71034 Böblingen GERMANY Telefon +49 7031 714 5200 Telefax +49 7031 714 5249 E-Mail Web info@emss.de http://www.emss.de FEKO Benchmark to handle big
More information3. Image formation, Fourier analysis and CTF theory. Paula da Fonseca
3. Image formation, Fourier analysis and CTF theory Paula da Fonseca EM course 2017 - Agenda - Overview of: Introduction to Fourier analysis o o o o Sine waves Fourier transform (simple examples of 1D
More informationTechnische Universiteit Delft. Acceleration of the Multi Level Fast Multipole Algorithm on a Graphics Processing Unit C. Wagenaar
Technische Universiteit Delft Acceleration of the Multi Level Fast Multipole Algorithm on a Graphics Processing Unit C. Wagenaar Acceleration of the Multi Level Fast Multipole Algorithm on a Graphics
More informationAdvanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude
Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude A. Migukin *, V. atkovnik and J. Astola Department of Signal Processing, Tampere University of Technology,
More informationTHE application of advanced computer architecture and
544 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 45, NO. 3, MARCH 1997 Scalable Solutions to Integral-Equation and Finite-Element Simulations Tom Cwik, Senior Member, IEEE, Daniel S. Katz, Member,
More informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Images from Few Fourier Samples Greg Ongie*, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa SampTA 2015 Washington, D.C. 1. Introduction 2.
More informationSTEPHEN WOLFRAM MATHEMATICADO. Fourth Edition WOLFRAM MEDIA CAMBRIDGE UNIVERSITY PRESS
STEPHEN WOLFRAM MATHEMATICADO OO Fourth Edition WOLFRAM MEDIA CAMBRIDGE UNIVERSITY PRESS Table of Contents XXI a section new for Version 3 a section new for Version 4 a section substantially modified for
More informationAccelerated flow acoustic boundary element solver and the noise generation of fish
Accelerated flow acoustic boundary element solver and the noise generation of fish JUSTIN W. JAWORSKI, NATHAN WAGENHOFFER, KEITH W. MOORED LEHIGH UNIVERSITY, BETHLEHEM, USA FLINOVIA PENN STATE 27 APRIL
More informationOptimum Array Processing
Optimum Array Processing Part IV of Detection, Estimation, and Modulation Theory Harry L. Van Trees WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION Preface xix 1 Introduction 1 1.1 Array Processing
More informationMethod of Moments enhancement technique for the analysis of Sierpinski pre-fractal antennas. Total number of pages: 10
FRACTALCOMS Exploring the limits of Fractal Electrodynamics for the future telecommunication technologies IST-2001-33055 Method of Moments enhancement technique for the analysis of Sierpinski pre-fractal
More informationImplementation of rotation-based operators for Fast Multipole Method in X10 Andrew Haigh
School of Computer Science College of Engineering and Computer Science Implementation of rotation-based operators for Fast Multipole Method in X10 Andrew Haigh Supervisors: Dr Eric McCreath Josh Milthorpe
More informationP052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges
P052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges M. Ayzenberg* (Norwegian University of Science and Technology, A. Aizenberg (Institute of Petroleum Geology and Geophysics,
More informationLong time integrations of a convective PDE on the sphere by RBF collocation
Long time integrations of a convective PDE on the sphere by RBF collocation Bengt Fornberg and Natasha Flyer University of Colorado NCAR Department of Applied Mathematics Institute for Mathematics Applied
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationPrecomputed Radiance Transfer: Theory and Practice
1 Precomputed Radiance Transfer: Peter-Pike Sloan Microsoft Jaakko Lehtinen Helsinki Univ. of Techn. & Remedy Entertainment Jan Kautz MIT 2 Introduction Jan Kautz MIT 3 Introduction We see here an example
More informationExample Lecture 12: The Stiffness Method Prismatic Beams. Consider again the two span beam previously discussed and determine
Example 1.1 Consider again the two span beam previously discussed and determine The shearing force M1 at end B of member B. The bending moment M at end B of member B. The shearing force M3 at end B of
More informationIndependent systems consist of x
5.1 Simultaneous Linear Equations In consistent equations, *Find the solution to each system by graphing. 1. y Independent systems consist of x Three Cases: A. consistent and independent 2. y B. inconsistent
More informationReflectivity modeling for stratified anelastic media
Reflectivity modeling for stratified anelastic media Peng Cheng and Gary F. Margrave Reflectivity modeling ABSTRACT The reflectivity method is widely used for the computation of synthetic seismograms for
More informationHigh-Frequency Algorithmic Advances in EM Tools for Signal Integrity Part 1. electromagnetic. (EM) simulation. tool of the practic-
From January 2011 High Frequency Electronics Copyright 2011 Summit Technical Media, LLC High-Frequency Algorithmic Advances in EM Tools for Signal Integrity Part 1 By John Dunn AWR Corporation Only 30
More informationParallel FFT Program Optimizations on Heterogeneous Computers
Parallel FFT Program Optimizations on Heterogeneous Computers Shuo Chen, Xiaoming Li Department of Electrical and Computer Engineering University of Delaware, Newark, DE 19716 Outline Part I: A Hybrid
More informationA pedestrian introduction to fast multipole methods
SCIENCE CHIN Mathematics. RTICLES. May 2012 Vol. 55 No. 5: 1043 1051 doi: 10.1007/s11425-012-4392-0 pedestrian introduction to fast multipole methods YING Lexing Department of Mathematics and ICES, University
More informationFRESNEL DIFFRACTION AND PARAXIAL WAVE EQUATION. A. Fresnel diffraction
19 IV. FRESNEL DIFFRACTION AND PARAXIAL WAVE EQUATION A. Fresnel diffraction Any physical optical beam is of finite transverse cross section. Beams of finite cross section may be described in terms of
More informationAdvanced modelling of gratings in VirtualLab software. Site Zhang, development engineer Lignt Trans
Advanced modelling of gratings in VirtualLab software Site Zhang, development engineer Lignt Trans 1 2 3 4 Content Grating Order Analyzer Rigorous Simulation of Holographic Generated Volume Grating Coupled
More informationAPPLYING EXTRAPOLATION AND INTERPOLATION METHODS TO MEASURED AND SIMULATED HRTF DATA USING SPHERICAL HARMONIC DECOMPOSITION.
APPLYING EXTRAPOLATION AND INTERPOLATION METHODS TO MEASURED AND SIMULATED HRTF DATA USING SPHERICAL HARMONIC DECOMPOSITION Martin Pollow Institute of Technical Acoustics RWTH Aachen University Neustraße
More informationSynthesis of Pencil Beam Linear Antenna Arrays using Simple FFT/CF/GA Based Technique
International Journal of Engineering & Technology IJET-IJENS Vol:13 No:05 86 Synthesis of Pencil Beam Linear Antenna Arrays using Simple FFT/CF/GA Based Technique B. Eldosouky 1, A. H. Hussein 1, H. H.
More informationUnit II-2. Orthogonal projection. Orthogonal projection. Orthogonal projection. the scalar is called the component of u along v. two vectors u,v are
Orthogonal projection Unit II-2 Orthogonal projection the scalar is called the component of u along v in real ips this may be a positive or negative value in complex ips this may have any complex value
More informationLarge-Scale Full-Wave Simulation
Large-Scale Full-Wave Simulation Sharad Kapur and David Long Integrand Software, Inc. Areas of interest Consistent trends in IC design Increasing operating frequencies Modeling of passive structures (components,
More informationUsing GPUs to compute the multilevel summation of electrostatic forces
Using GPUs to compute the multilevel summation of electrostatic forces David J. Hardy Theoretical and Computational Biophysics Group Beckman Institute for Advanced Science and Technology University of
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationAQA GCSE Maths - Higher Self-Assessment Checklist
AQA GCSE Maths - Higher Self-Assessment Checklist Number 1 Use place value when calculating with decimals. 1 Order positive and negative integers and decimals using the symbols =,, , and. 1 Round to
More informationImproved Multilevel Fast Multipole Method for Higher-Order Discretizations
Improved Multilevel Fast Multipole Method for Higher-Order Discretizations Oscar Borries 1,, Peter Meincke, Erik Jørgensen, Stig Busk Sørensen, Per Christian Hansen 1 1 Technical University of Denmark,
More informationThe meshfree computation of stationary electric current densities in complex shaped conductors using 3D boundary element methods
Boundary Elements and Other Mesh Reduction Methods XXXVII 121 The meshfree computation of stationary electric current densities in complex shaped conductors using 3D boundary element methods A. Buchau
More informationLegal and impossible dependences
Transformations and Dependences 1 operations, column Fourier-Motzkin elimination us use these tools to determine (i) legality of permutation and Let generation of transformed code. (ii) Recall: Polyhedral
More informationIntroduction to Inverse Problems
Introduction to Inverse Problems What is an image? Attributes and Representations Forward vs Inverse Optical Imaging as Inverse Problem Incoherent and Coherent limits Dimensional mismatch: continuous vs
More informationA Broadband Fast Multipole Accelerated Boundary Element Method for the 3D Helmholtz Equation. Abstract
A Broadband Fast Multipole Accelerated Boundary Element Method for the 3D Helmholtz Equation Nail A. Gumerov and Ramani Duraiswami (Dated: 31 December 2007, Revised 22 July 2008, Revised 03 October 2008.)
More informationSimulation Advances for RF, Microwave and Antenna Applications
Simulation Advances for RF, Microwave and Antenna Applications Bill McGinn Application Engineer 1 Overview Advanced Integrated Solver Technologies Finite Arrays with Domain Decomposition Hybrid solving:
More informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationDURING the past decade, several fast algorithms for efficiently. Time Domain Adaptive Integral Method for Surface Integral Equations
2692 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 52, NO. 10, OCTOBER 2004 Time Domain Adaptive Integral Method for Surface Integral Equations Ali E. Yılmaz, Student Member, IEEE, Jian-Ming Jin,
More informationX-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes
X-ray tomography Important application of the Fast Fourier transform: X-ray tomography. Also referred to as CAT scan (Computerized Axial Tomography) X-ray tomography This has revolutionized medical diagnosis.
More informationElectromagnetic Field Transformations for Measurements and Simulations
Progress In Electromagnetics Research, Vol. 151, 127 150, 2015 Electromagnetic Field Transformations for Measurements and Simulations Thomas F. Eibert *, Emre Kılıç, Carlos Lopez, Raimund A. M. Mauermayer,
More informationHFSS Ansys ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary
HFSS 12.0 Ansys 2009 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Comparison of HFSS 11 and HFSS 12 for JSF Antenna Model UHF blade antenna on Joint Strike Fighter Inherent improvements in
More informationGeometric transformations assign a point to a point, so it is a point valued function of points. Geometric transformation may destroy the equation
Geometric transformations assign a point to a point, so it is a point valued function of points. Geometric transformation may destroy the equation and the type of an object. Even simple scaling turns a
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationFAST AND ACCURATE TRANSIENT ULTRASOUND PROPAGATION AND B-MODE IMAGING SIMULATION METHODS. Yi Zhu
FAST AND ACCURATE TRANSIENT ULTRASOUND PROPAGATION AND B-MODE IMAGING SIMULATION METHODS By Yi Zhu A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree
More informationSung-Eui Yoon ( 윤성의 )
CS480: Computer Graphics Curves and Surfaces Sung-Eui Yoon ( 윤성의 ) Course URL: http://jupiter.kaist.ac.kr/~sungeui/cg Today s Topics Surface representations Smooth curves Subdivision 2 Smooth Curves and
More informationGPU-based Distributed Behavior Models with CUDA
GPU-based Distributed Behavior Models with CUDA Courtesy: YouTube, ISIS Lab, Universita degli Studi di Salerno Bradly Alicea Introduction Flocking: Reynolds boids algorithm. * models simple local behaviors
More informationLinear Equation Systems Iterative Methods
Linear Equation Systems Iterative Methods Content Iterative Methods Jacobi Iterative Method Gauss Seidel Iterative Method Iterative Methods Iterative methods are those that produce a sequence of successive
More informationCHAPTER 3 DIFFERENT DOMAINS OF WATERMARKING. domain. In spatial domain the watermark bits directly added to the pixels of the cover
38 CHAPTER 3 DIFFERENT DOMAINS OF WATERMARKING Digital image watermarking can be done in both spatial domain and transform domain. In spatial domain the watermark bits directly added to the pixels of the
More informationAnalysis of Two-dimensional Scattering by a Periodic Array of Conducting Cylinders Using the Method of Auxiliary Sources
PIERS ONLINE, VOL. 4, NO. 5, 8 51 Analysis of Two-dimensional Scattering by a Periodic Array of Conducting Cylinders Using the Method of Auxiliary Sources Naamen Hichem and Taoufik Aguili Ecole Nationale
More informationInterlude: Solving systems of Equations
Interlude: Solving systems of Equations Solving Ax = b What happens to x under Ax? The singular value decomposition Rotation matrices Singular matrices Condition number Null space Solving Ax = 0 under
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More information5. GENERALIZED INVERSE SOLUTIONS
5. GENERALIZED INVERSE SOLUTIONS The Geometry of Generalized Inverse Solutions The generalized inverse solution to the control allocation problem involves constructing a matrix which satisfies the equation
More information29. Diffraction of waves
9. Diffraction of waves Light bends! Diffraction assumptions The Kirchhoff diffraction integral Fresnel Diffraction diffraction from a slit Diffraction Light does not always travel in a straight line.
More informationMEMORY EFFICIENT WDR (WAVELET DIFFERENCE REDUCTION) using INVERSE OF ECHELON FORM by EQUATION SOLVING
Available Online at www.ijcsmc.com International Journal of Computer Science and Mobile Computing A Monthly Journal of Computer Science and Information Technology IJCSMC Vol. 3 Issue. 7 July 2014 pg.512
More informationA Level Set-Based Topology Optimization Method For Three-Dimensional Acoustic Problems Using Fast Multipole Boundary Element Method
9 th World Congress on Structural and Multidisciplinary Optimization June 13-17, 2011, Shizuoka, Japan A Level Set-Based Topology Optimization Method For Three-imensional Acoustic Problems Using Fast Multipole
More informationPanel methods are currently capable of rapidly solving the potential flow equation on rather complex
A Fast, Unstructured Panel Solver John Moore 8.337 Final Project, Fall, 202 A parallel high-order Boundary Element Method accelerated by the Fast Multipole Method is presented in this report. The case
More informationksa 400 Growth Rate Analysis Routines
k-space Associates, Inc., 2182 Bishop Circle East, Dexter, MI 48130 USA ksa 400 Growth Rate Analysis Routines Table of Contents ksa 400 Growth Rate Analysis Routines... 2 1. Introduction... 2 1.1. Scan
More informationEllipsoid Algorithm :Algorithms in the Real World. Ellipsoid Algorithm. Reduction from general case
Ellipsoid Algorithm 15-853:Algorithms in the Real World Linear and Integer Programming II Ellipsoid algorithm Interior point methods First polynomial-time algorithm for linear programming (Khachian 79)
More informationCMSC 714 Lecture 6 MPI vs. OpenMP and OpenACC. Guest Lecturer: Sukhyun Song (original slides by Alan Sussman)
CMSC 714 Lecture 6 MPI vs. OpenMP and OpenACC Guest Lecturer: Sukhyun Song (original slides by Alan Sussman) Parallel Programming with Message Passing and Directives 2 MPI + OpenMP Some applications can
More information1.1 calculator viewing window find roots in your calculator 1.2 functions find domain and range (from a graph) may need to review interval notation
1.1 calculator viewing window find roots in your calculator 1.2 functions find domain and range (from a graph) may need to review interval notation functions vertical line test function notation evaluate
More informationFast Radial Basis Functions for Engineering Applications. Prof. Marco Evangelos Biancolini University of Rome Tor Vergata
Fast Radial Basis Functions for Engineering Applications Prof. Marco Evangelos Biancolini University of Rome Tor Vergata Outline 2 RBF background Fast RBF on HPC Engineering Applications Mesh morphing
More informationGraphics Processing Unit Implementation of Multilevel Plane-Wave Time-Domain Algorithm
IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS, VOL. 13, 2014 1671 Graphics Processing Unit Implementation of Multilevel Plane-Wave Time-Domain Algorithm Yang Liu, Student Member, IEEE, Abdulkadir C. Yücel,
More informationLecture 2: Introduction
Lecture 2: Introduction v2015.0 Release ANSYS HFSS for Antenna Design 1 2015 ANSYS, Inc. Multiple Advanced Techniques Allow HFSS to Excel at a Wide Variety of Applications Platform Integration and RCS
More informationSystem Demonstration of Spiral: Generator for High-Performance Linear Transform Libraries
System Demonstration of Spiral: Generator for High-Performance Linear Transform Libraries Yevgen Voronenko, Franz Franchetti, Frédéric de Mesmay, and Markus Püschel Department of Electrical and Computer
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 14 130307 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Stereo Dense Motion Estimation Translational
More informationPhysics I : Oscillations and Waves Prof. S Bharadwaj Department of Physics & Meteorology Indian Institute of Technology, Kharagpur
Physics I : Oscillations and Waves Prof. S Bharadwaj Department of Physics & Meteorology Indian Institute of Technology, Kharagpur Lecture - 20 Diffraction - I We have been discussing interference, the
More information2D and 3D Far-Field Radiation Patterns Reconstruction Based on Compressive Sensing
Progress In Electromagnetics Research M, Vol. 46, 47 56, 206 2D and 3D Far-Field Radiation Patterns Reconstruction Based on Compressive Sensing Berenice Verdin * and Patrick Debroux Abstract The measurement
More informationMulti-azimuth velocity estimation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 87 Multi-azimuth velocity estimation Robert G. Clapp and Biondo Biondi 1 ABSTRACT It is well known that the inverse problem of estimating interval
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Review of Motion Modelling and Estimation Introduction to Motion Modelling & Estimation Forward Motion Backward Motion Block Motion Estimation Motion
More informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationAim. Structure and matrix sparsity: Part 1 The simplex method: Exploiting sparsity. Structure and matrix sparsity: Overview
Aim Structure and matrix sparsity: Part 1 The simplex method: Exploiting sparsity Julian Hall School of Mathematics University of Edinburgh jajhall@ed.ac.uk What should a 2-hour PhD lecture on structure
More informationELECTROMAGNETIC diffraction by perfectly conducting
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 47, NO. 11, NOVEMBER 1999 1697 Oblique Scattering by a Pair of Conducting Half Planes: TM Case Jong-Won Yu and Noh-Hoon Myung Abstract An exact series
More informationMethod of moments solution of volume integral equations using higher-order hierarchical Legendre basis functions
RADIO SCIENCE, VOL. 39,, doi:10.1029/2004rs003041, 2004 Method of moments solution of volume integral equations using higher-order hierarchical Legendre basis functions O. S. Kim, P. Meincke, and O. Breinbjerg
More informationVirtual EM Inc. Ann Arbor, Michigan, USA
Functional Description of the Architecture of a Special Purpose Processor for Orders of Magnitude Reduction in Run Time in Computational Electromagnetics Tayfun Özdemir Virtual EM Inc. Ann Arbor, Michigan,
More informationNumerical Algorithms
Chapter 10 Slide 464 Numerical Algorithms Slide 465 Numerical Algorithms In textbook do: Matrix multiplication Solving a system of linear equations Slide 466 Matrices A Review An n m matrix Column a 0,0
More information