HILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT
|
|
- Emery Davidson
- 5 years ago
- Views:
Transcription
1 Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is a Lovey Boundary Eement Research Too Matab ibrary for h-adaptive Gaerkin BEM currenty, owest-order eements for 2D Lapacian P 0 for fuxes S 1 for traces research code for FWF project P21732 overview on current state of the art starting point for further investigations free for academic use Vienna University of Technoogy Institute for Anaysis and Scientific Computing Features of HILBERT HILBERT? C impementation of integra operators via MEX interface V (P 0 Γ P0 Γ ) K(S 1 Γ P0 Γ ) W(S 1 Γ S1 Γ ) N(P 1 Ω P0 Γ ) remaining codes in Matab fuy vectorized different error estimators different marking strategies oca mesh-refinement demo fies and adaptive agorithms for weaky-singuar integra equation hypersinguar integra equation symmetric integra formuation of mixed BVP with/without voume force
2 Current Restrictions of HILBERT Γ = Ω piecewise affine boundary ony owest-order eements ony canonica bases characteristic functions for P 0 hat functions for S 1 Weaky-Singuar Integra Equation no matrix compression, i.e., dense matrices direct soution of inear systems by Matab backsash operator Outine Mode Probem 1 HILBERT? 2 Weaky-Singuar Integra Equation 3 Introduction to HILBERT 4 Concusion & Outook Dirichet Probem u = 0 in Ω u = g on Γ Weaky-Singuar Integra Equation V φ = (K + 1/2)g H := H 1/2 (Γ) energy space equivaent energy norm on H induced by V
3 Variationa Formuation Continuous Formuation Seek φ H := H 1/2 (Γ) s.t. V φ,ψ = (K + 1/2)g,ψ for a ψ H T partition of Γ Gaerkin Formuation Seek Φ X := P 0 (T ) s.t. V Φ,Ψ = (K + 1/2)g,Ψ for a Ψ X Perturbed Gaerkin Formuation 2/2 Matrix Formuation of Perturbed Gaerkin Formuation Compute coefficient vector x of Φ by soving Vx = where g is noda vector of g (K M ) g mass-type matrix M anayticay computabe and sparse integra operator matrices V, K anayticay computabe drawback: RHS hard to compute in genera Perturbed Gaerkin Formuation 1/2 Contro Data Approximation Reminder: Continuous Probem assume additiona reguarity g H 1 (Γ) G S 1 (T ) noda interpoant of g Perturbed Gaerkin Formuation Seek Φ X s.t. V φ = (K + 1/2)g Reminder: Perturbed Continuous Probem V φ = (K + 1/2)G V Φ,Ψ = (K + 1/2)G,Ψ for a Ψ X Lemma (Aurada, Godenits, P. 09) φ φ g G H 1/2 (Γ) h1/2 (g G ) L 2 (Γ) =: osc
4 Contro Gaerkin Error Loca Refinement Indicators Φ Gaerkin soution w.r.t. T Φ Gaerkin soution w.r.t. T L 2 -projection Π onto X = P 0 (T ) Lemma (Ferraz-Leite, P. 08) There hods η η µ µ η for error estimators η = Φ Φ η = Φ Π Φ µ = h 1/2 ( Φ Φ ) L 2 (Γ) µ = h 1/2 ( Φ Π Φ ) L 2 (Γ) ρ (T) 2 := h 1/2 ( Φ Π Φ ) 2 L 2 (T) + h1/2 (g G ) 2 L 2 (T) φ Φ 2 (η + osc ) 2 µ 2 + osc2 = T T ρ (T) 2 everything works aso with µ instead of µ Contro Gaerkin and Data Approximation Error Saturation Assumption for Non-Perturbed Probem Φ Gaerkin soution w.r.t. T and non-perturbed RHS Φ Gaerkin soution w.r.t. T and non-perturbed RHS φ Φ q φ Φ Coroary (Aurada, Godenits, P. 09) η φ Φ + osc saturation assumption = φ Φ η + osc Adaptive Mesh-Refinement initia mesh T 0 adaptivity parameter 0 < θ < 1 Adaptive Agorithm 1 compute discrete soution Φ X = P 0 ( T ) 2 for a T T, compute ρ (T) 3 find (minima) set M T s.t. θ T T ρ (T) 2 T M ρ (T) 2 4 refine (at east) marked eements T M to obtain T +1 ensure κ(t +1 ) 2 κ(t 0 ) 5 increase counter + 1 and iterate
5 Convergence of Adaptive Scheme Data Structure Theorem (Aurada, Godenits, P. 09) Use refinement procedure of [AGP 09] Then, concept of estimator reduction appies, whence im ( µ2 + osc2 ) = 0 Saturation assumption for non-perturbed probem impies im φ Φ = 0 = im φ Φ T = {T 1,...,T N } partition N = {z 1,...,z N } nodes coordinates (N 2) doube-array z j = (x,y) R 2 corresponds to coordinates(j,:) = [x, y ] eements (N 2) int-array T i = conv{z j,z k } corresponds to eements(i,:) = [j, k ] Buid Gaerkin Data Introduction to HILBERT buidv assembes V = V (P 0 P 0 ) use of symmetry outer integration over eement with ength(t j ) ength(t k ) outer integration by Gaussian quadrature for far-fied entries buidk assembes K = K(S 1 P 0 ) buidsymmrhs assembes b for perturbed RHS (K + 1/2)G
6 Loca Mesh-Refinement (h h/2)-error Estimators 1/2 markeements returns index vector for marked eements M T from θ T T ρ (T) 2 T M ρ (T) 2 even mutipe trianguations, e.g., FEM-BEM, mixed BVP refineboundarymesh agorithm from [Aurada, Godenits, P. 09] compexity O(N og N) for oca mesh-refinement compexity O(N) for uniform mesh-refinement computeestspeta returns η 2 = Φ Φ 2. computeestspetatide returns η 2 = Φ Π Φ 2. Data Osciations (h h/2)-error Estimators 2/2 computeoscdirichet returns v R N with osc 2 õsc2 = N j=1 v j, where v j ength(t j ) (g G ) 2 L 2 (T j ) computeestspmu returns v R N with µ 2 = N j=1 v j, where v j = ength(t j ) Φ Φ 2 L 2 (T j ) quadrature ocated at nodes and midpoints exact for g S 2 osc = O(h 3/2 ) and osc õsc = O(h 5/2 ) for smooth g computeestspmutide returns ṽ R N with µ 2 = N j=1 ṽj, where ṽ j = ength(t j ) Φ Π Φ 2 L 2 (T j )
7 Adaptive Mesh-Refining Agorithm 1 whie s i z e ( eements,1) < nemax 2 3 % bui d uniformy r ef i ned mesh 4 [ co ordin ates f i ne, eements fine, father2son ]... 5 = refineboundarymesh ( coordinates, eements ) ; 6 7 % compute fine mesh sou t i on 8 V fine = buidv ( c oordinates f ine, e ements f i ne ) ; 9 b f i ne = buidsymmrhs( coo r dina t es f i ne, eements ) ; 10 x f i n e = V fine \ b f i ne ; % compute (h h/2) e r r o r estimator t i de mu 13 mu tide = computeestspmutide( coordinates, eements, father2son, x f i n e ) ; 15 % compute data o s c i a t i o n s 16 o sc f i ne = computeoscdirichet ( coo rd ina t es f in e, eements ) ; 17 osc = o s c f i n e ( father2son ( :, 1) ) + o s c f i n e ( father2son ( :, 2 ) ) ; % mark eements f or refinement 20 marked = markeements ( theta, mu tide + osc ) ; % generate new mesh 23 [ coordinates, eements ] = refineboundarymesh ( coordinates, eements, marked ) ; end estimators Numerica Experiment 2/2 O(N 2/3 ) err + osc (unif) O(N 11/10 ) eµ (unif) osc (unif) err + osc (adap) O(N 3/2 ) eµ (adap) osc (adap) number of eements Numerica Experiment 1/2 Adaptivity Pays 1/ adaptive pays w.r.t. number of eements observe optima convergence behaviour does adaptivity pay w.r.t. computationa time? := unif () (T 0 ) +assembe +sove t (adap) 0 := assembe +sove +estimate +mark +refine t (unif) t (adap) +1 := t (adap) +... known exact soution φ generic singuarity of φ at reentrant corner (green) singuarity of g at uppermost corner (red)
8 Adaptivity Pays 2/2 HILBERT Reease One estimators err + osc (unif) eµ (unif) osc (unif) err + osc (adap) eµ (adap) osc (adap) computationa time [s] reease date: Friday 13. November 2009 avaiabe at: Reease One contains owest-order BEM for 2D Lapacian Dirichet probem Neumann probem mixed Dirichet-Neumann probem with/without voume forces Reease One incudes (h h/2)-error estimators for a probems oca mesh-refinement for boundary/voume meshes adaptive agorithms & demo fies severa visuaization toos fu documentation of a Matab codes (about 60 pages) Short-Term Extensions Concusion & Outook further error estimators, e.g., weighted-residua error estimators two-eve error estimators Faermann residua error estimator Steinbach-Schuz error estimator adaptive FEM-BEM couping Lamé equation
9 Thank You for Listening Vienna University of Technoogy Institute for Anaysis and Scientific Computing
HILBERT-A MATLAB Implementation of Adaptive 2D-BEM
ASC Report No. 24/2011 HILBERT-A MATLAB Impementation of Adaptive 2D-BEM Markus Aurada, Michae Ebner, Michae Feisch, Samue Ferraz-Leite, Petra Godenits, Michae Karkuik, Markus Mayr, Dirk Praetorius Institute
More informationTechnische Universität Graz
Technische Universität Graz Convergence of adaptive 3D BEM for weaky singuar integra equations based on isotropic mesh-refinement M. Karkuik, G. Of, D. Praetorius Berichte aus dem Institut für Numerische
More informationSmoothers for ecient multigrid methods in IGA
Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third
More informationConvergence of adaptive 3D BEM for weakly singular integral equations based on isotropic mesh-refinement
ASC Report No. 20/2012 Convergence of adaptive 3D BEM for weaky singuar integra equations based on isotropic mesh-refinement M. Karkuik, G. Of, and D. Praetorius Institute for Anaysis and Scientific Computing
More informationM. Aurada 1,M.Feischl 1, J. Kemetmüller 1,M.Page 1 and D. Praetorius 1
ESAIM: M2AN 47 (2013) 1207 1235 DOI: 10.1051/m2an/2013069 ESAIM: Mathematica Modeing and Numerica Anaysis www.esaim-m2an.org EACH H 1/2 STABLE PROJECTION YIELDS CONVERGENCE AND QUASI OPTIMALITY OF ADAPTIVE
More informationConvergence of Adaptive FEM for some Elliptic Obstacle Problem with Inhomogeneous Dirichlet Data
ASC Report No. 33/2010 Convergence of Adaptive FEM for some Eiptic Obstace Probem with Inhomogeneous Dirichet Data Michae Feisch, Marcus Page, Dirk Praetorius Institute for Anaysis and Scientific Computing
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationNumerical methods for PDEs FEM - abstract formulation, the Galerkin method
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona
More informationIntroduction to PDEs and Numerical Methods Tutorial 10. Finite Element Analysis
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Introduction to PDEs and Numerica Methods Tutoria. Finite Eement Anaysis Dr. Noemi Friedman, 3..2 FROM STRONG FORM TO WEAK FORM inhomogeneous
More informationMultigrid Method for Elliptic Control Problems
J OHANNES KEPLER UNIVERSITÄT LINZ Netzwerk f ür Forschung, L ehre und Praxis Mutigrid Method for Eiptic Contro Probems MASTERARBEIT zur Erangung des akademischen Grades MASTER OF SCIENCE in der Studienrichtung
More informationConvergence of adaptive BEM and adaptive FEM-BEM coupling for estimators without h-weighting factor
ASC Report No. 16/2014 Convergence of adaptive BEM and adaptive FEM-BEM couping for estimators without h-weighting factor M. Feisch, T. Führer, G. Mitscha-Eib, D. Praetorius, and E.P. Stephan Institute
More informationOn 2D newest vertex bisection: Optimality of mesh-closure and H 1 -stability of L 2 -projection
ASC Report No. 10/2012 On 2D newest vertex bisection: Optimaity of mesh-cosure and H 1 -stabiity of L 2 -projection Michae Karkuik, David Pavicek, and Dirk Praetorius Institute for Anaysis and Scientific
More informationAdaptive Boundary Element Methods Part 2: ABEM. Dirk Praetorius
CENTRAL School on Analysis and Numerics for PDEs November 09-12, 2015 Adaptive Boundary Element Methods Part 2: ABEM Dirk Praetorius TU Wien Institute for Analysis and Scientific Computing Outline 1 Boundary
More informationarxiv: v1 [math.na] 24 Jan 2019
ADAPTIVE BOUNDARY ELEMENT METHODS FOR THE COMPUTATION OF THE ELECTROSTATIC CAPACITY ON COMPLEX POLYHEDRA TIMO BETCKE, ALEXANDER HABERL, AND DIRK PRAETORIUS arxiv:1901.08393v1 [math.na] 24 Jan 2019 Abstract.
More informationAnalysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case
Technica Report PC-04-00 Anaysis of Emerson s Mutipe Mode Interpoation Estimation Agorithms: The MIMO Case João P. Hespanha Dae E. Seborg University of Caifornia, Santa Barbara February 0, 004 Anaysis
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,
More informationAn extension of the MAC scheme to locally refined meshes : convergence analysis for the full tensor time dependent Navier Stokes equations
An extension of the MAC scheme to ocay refined meshes : convergence anaysis for the fu tensor time dependent Navier Stokes equations Eric Chénier, Robert Eymard, Thierry Gaouët and Raphaèe Herbin Apri
More informationA two-level Schwarz preconditioner for heterogeneous problems
A two-eve Schwarz preconditioner for heterogeneous probems V. Doean 1, F. Nataf 2, R. Scheich 3 and N. Spiane 2 1 Introduction Coarse space correction is essentia to achieve agorithmic scaabiity in domain
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationNumerische Mathematik
Numer. Math. (1999) 84: 97 119 Digita Object Identifier (DOI) 10.1007/s002119900096 Numerische Mathematik c Springer-Verag 1999 Mutigrid methods for a parameter dependent probem in prima variabes Joachim
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationFFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection
FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary
More informationSolving Maxwell s Equations Using the Ultra Weak Variational Formulation
Soving Maxwe s Equations Using the Utra Weak Variationa Formuation T. Huttunen, M. Mainen and P. Monk Department of Appied Physics, University of Kuopio, P.O.Box 1627, 7211 Kuopio, Finand Department of
More informationON THE STRONG CONVERGENCE OF GRADIENTS IN STABILIZED DEGENERATE CONVEX MINIMIZATION PROBLEMS
SIAM J NUMER ANAL Vo 47, No 6, pp 4569 4580 c 010 Society for Industria and Appied Mathematics ON THE STRONG CONVERGENCE O GRADIENTS IN STABILIZED DEGENERATE CONVEX MINIMIZATION PROBLEMS WOLGANG BOIGER
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationPreprint BUW-SC 11/1. Bergische Universität Wuppertal. Fachbereich C Mathematik und Naturwissenschaften. Mathematik
Preprint BUW-SC 11/1 Bergische Universität Wupperta Fachbereich C Mathematik und Naturwissenschaften Mathematik A. Brandt, J. Brannick, K. Kah, and I. Livshits Bootstrap AMG Januar 2011 http://www-ai.math.uni-wupperta.de/scicomp/
More informationNumerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet
Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport
More informationA nonlinear multigrid for imaging electrical conductivity and permittivity at low frequency
INSTITUTE OF PHYSICS PUBLISHING INVERSE PROBLEMS Inverse Probems 17 (001) 39 359 www.iop.org/journas/ip PII: S066-5611(01)14169-9 A noninear mutigrid for imaging eectrica conductivity and permittivity
More information18-660: Numerical Methods for Engineering Design and Optimization
8-660: Numerica Methods for Engineering esign and Optimization in i epartment of ECE Carnegie Meon University Pittsburgh, PA 523 Side Overview Conjugate Gradient Method (Part 4) Pre-conditioning Noninear
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationConvergence and quasi-optimality of adaptive finite element methods for harmonic forms
Noname manuscript No. (wi be inserted by the editor) Convergence and quasi-optimaity of adaptive finite eement methods for harmonic forms Aan Demow 1 the date of receipt and acceptance shoud be inserted
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationLocal defect correction for time-dependent problems
Loca defect correction for time-dependent probems Minero, R. DOI: 10.6100/IR608579 Pubished: 01/01/2006 Document Version Pubisher s PDF, aso known as Version of Record (incudes fina page, issue and voume
More informationFurther analysis of multilevel Monte Carlo methods for elliptic PDEs with random coefficients
Further anaysis of mutieve Monte Caro methods for eiptic PDEs with random coefficients A. L. Teckentrup, R. Scheich, M. B. Gies, and E. Umann Abstract We consider the appication of mutieve Monte Caro methods
More informationA UNIVERSAL METRIC FOR THE CANONICAL BUNDLE OF A HOLOMORPHIC FAMILY OF PROJECTIVE ALGEBRAIC MANIFOLDS
A UNIERSAL METRIC FOR THE CANONICAL BUNDLE OF A HOLOMORPHIC FAMILY OF PROJECTIE ALGEBRAIC MANIFOLDS DROR AROLIN Dedicated to M Saah Baouendi on the occasion of his 60th birthday 1 Introduction In his ceebrated
More informationASummaryofGaussianProcesses Coryn A.L. Bailer-Jones
ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe
More informationc 2011 Society for Industrial and Applied Mathematics BOOTSTRAP AMG
SIAM J. SCI. COMPUT. Vo. 0, No. 0, pp. 000 000 c 2011 Society for Industria and Appied Mathematics BOOTSTRAP AMG A. BRANDT, J. BRANNICK, K. KAHL, AND I. LIVSHITS Abstract. We deveop an agebraic mutigrid
More informationPrimal and dual active-set methods for convex quadratic programming
Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:
More informationA MULTI-ORDER DISCONTINUOUS GALERKIN MONTE CARLO METHOD FOR HYPERBOLIC PROBLEMS WITH STOCHASTIC PARAMETERS
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 A MULTI-ORDER DISCONTINUOUS GALERKIN MONTE CARLO METHOD FOR HYPERBOLIC PROBLEMS WITH
More informationAcoustic-Structure Simulation of Exhaust Systems by Coupled FEM and Fast BEM
Acoustic-Structure Simuation of Exhaust Systems by Couped FEM and Fast BEM Lothar Gau, Michae Junge Institute of Appied and Experimenta Mechanics, University of Stuttgart Pfaffenwadring 9, 755 Stuttgart,
More informationBCCS TECHNICAL REPORT SERIES
BCCS TCHNICAL RPORT SRIS The Crouzeix-Raviart F on Nonmatching Grids with an Approximate Mortar Condition Taa Rahman, Petter Bjørstad and Xuejun Xu RPORT No. 17 March 006 UNIFOB the University of Bergen
More information3. LS-DYNA Anwenderforum, Bamberg 2004 Netzfreie Verfahren. Method (EFG) in LS-DYNA Implementation and Applications
3. Anwenderforum, Bamberg 004 Netzfreie Verfahren VERMORE SOFTWARE TECHNOOGY CORPORATON GM RESEARCH & DEVEOPMENT CENTER Eement-Free Free-Gaerkin Method (EFG) in mpementation and Appications C.T. Wu and
More informationAerodynamic Design of Turbine Blades Using an Adjoint Equation Method
43rd AIAA Aerospace Sciences Meeting and Exhibit, January 10-13, 2005, Reno, Nevada Aerodynamic Design of Turbine ades Using an Adjoint Equation Method Hsiao-Yuan Wu, and Feng Liu Department of Mechanica
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vo. 36, No. 2, pp. 58 606 c 999 Society for Industria and Appied Mathematics OVERLAPPING NONMATCHING GRID MORTAR ELEMENT METHODS FOR ELLIPTIC PROBLEMS XIAO-CHUAN CAI, MAKSYMILIAN DRYJA,
More informationConvergence of adaptive stochastic Galerkin FEM
ASC Report No. 31/2018 Convergence of adaptive stochastic Gaerkin FEM A. Bespaov, D. Praetorius, L. Rocchi, and M. Ruggeri Institute for Anaysis and Scientific Computing Vienna University of Technoogy
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More informationSubstructuring Preconditioners for the Bidomain Extracellular Potential Problem
Substructuring Preconditioners for the Bidomain Extraceuar Potentia Probem Mico Pennacchio 1 and Vaeria Simoncini 2,1 1 IMATI - CNR, via Ferrata, 1, 27100 Pavia, Itay mico@imaticnrit 2 Dipartimento di
More informationA Novel Learning Method for Elman Neural Network Using Local Search
Neura Information Processing Letters and Reviews Vo. 11, No. 8, August 2007 LETTER A Nove Learning Method for Eman Neura Networ Using Loca Search Facuty of Engineering, Toyama University, Gofuu 3190 Toyama
More informationPreconditioned Locally Harmonic Residual Method for Computing Interior Eigenpairs of Certain Classes of Hermitian Matrices
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Preconditioned Locay Harmonic Residua Method for Computing Interior Eigenpairs of Certain Casses of Hermitian Matrices Vecharynski, E.; Knyazev,
More informationTitleCryptanalysis of the Quaternion Rai. IEICE Transactions on Fundamentals.
TiteCryptanaysis of the Quaternion Rai Author(s Hashimoto, Yasufumi Citation IEICE Transactions on Fundamentas Communications and Computer Science Issue Date 205-0-0 URL http://hd.hande.net/20.500.2000/
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationIntroduction to LMTO method
1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More informationDistributed average consensus: Beyond the realm of linearity
Distributed average consensus: Beyond the ream of inearity Usman A. Khan, Soummya Kar, and José M. F. Moura Department of Eectrica and Computer Engineering Carnegie Meon University 5 Forbes Ave, Pittsburgh,
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationA Robust Multigrid Method for Isogeometric Analysis using Boundary Correction. C. Hofreither, S. Takacs, W. Zulehner. G+S Report No.
A Robust Mutigrid Method for Isogeometric Anaysis using Boundary Correction C. Hofreither, S. Takacs, W. Zuehner G+S Report No. 33 Juy 2015 A Robust Mutigrid Method for Isogeometric Anaysis using Boundary
More informationOptimal Control of Assembly Systems with Multiple Stages and Multiple Demand Classes 1
Optima Contro of Assemby Systems with Mutipe Stages and Mutipe Demand Casses Saif Benjaafar Mohsen EHafsi 2 Chung-Yee Lee 3 Weihua Zhou 3 Industria & Systems Engineering, Department of Mechanica Engineering,
More informationAbsolute Value Preconditioning for Symmetric Indefinite Linear Systems
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Absoute Vaue Preconditioning for Symmetric Indefinite Linear Systems Vecharynski, E.; Knyazev, A.V. TR2013-016 March 2013 Abstract We introduce
More informationA. Distribution of the test statistic
A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch
More informationSUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS
ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia
More informationFEMAG: A High Performance Parallel Finite Element Toolbox for Electromagnetic Computations
Internationa Journa of Energy and Power Engineering 2016; 5(1-1): 57-64 Pubished onine November 28, 2015 (http://www.sciencepubishinggroup.com/j/ijepe) doi: 10.11648/j.ijepe.s.2016050101.19 ISSN: 2326-957X
More informationModel Solutions (week 4)
CIV-E16 (17) Engineering Computation and Simuation 1 Home Exercise 6.3 Mode Soutions (week 4) Construct the inear Lagrange basis functions (noda vaues as degrees of freedom) of the ine segment reference
More informationHomotopy Perturbation Method for Solving Partial Differential Equations of Fractional Order
Int Journa of Math Anaysis, Vo 6, 2012, no 49, 2431-2448 Homotopy Perturbation Method for Soving Partia Differentia Equations of Fractiona Order A A Hemeda Department of Mathematics, Facuty of Science
More informationHigh Spectral Resolution Infrared Radiance Modeling Using Optimal Spectral Sampling (OSS) Method
High Spectra Resoution Infrared Radiance Modeing Using Optima Spectra Samping (OSS) Method J.-L. Moncet and G. Uymin Background Optima Spectra Samping (OSS) method is a fast and accurate monochromatic
More informationThe Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation
The Symmetric Antipersymmetric Soutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B 2 + + A X B C Its Optima Approximation Ying Zhang Member IAENG Abstract A matrix A (a ij) R n n is said to be symmetric
More informationDISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE
DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationOn norm equivalence between the displacement and velocity vectors for free linear dynamical systems
Kohaupt Cogent Mathematics (5 : 95699 http://dxdoiorg/8/33835595699 COMPUTATIONAL SCIENCE RESEARCH ARTICLE On norm equivaence between the dispacement and veocity vectors for free inear dynamica systems
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationNonlinear Gaussian Filtering via Radial Basis Function Approximation
51st IEEE Conference on Decision and Contro December 10-13 01 Maui Hawaii USA Noninear Gaussian Fitering via Radia Basis Function Approximation Huazhen Fang Jia Wang and Raymond A de Caafon Abstract This
More informationNEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION
NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION Hsiao-Chang Chen Dept. of Systems Engineering University of Pennsyvania Phiadephia, PA 904-635, U.S.A. Chun-Hung Chen
More informationCASCADIC MULTILEVEL METHODS FOR FAST NONSYMMETRIC BLUR- AND NOISE-REMOVAL. Dedicated to Richard S. Varga on the occasion of his 80th birthday.
CASCADIC MULTILEVEL METHODS FOR FAST NONSYMMETRIC BLUR- AND NOISE-REMOVAL S. MORIGI, L. REICHEL, AND F. SGALLARI Dedicated to Richard S. Varga on the occasion of his 80th birthday. Abstract. Image deburring
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationCoupling of LWR and phase transition models at boundary
Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of
More informationWavelet Galerkin Solution for Boundary Value Problems
Internationa Journa of Engineering Research and Deveopment e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Voume, Issue 5 (May 24), PP.2-3 Waveet Gaerkin Soution for Boundary Vaue Probems D. Pate, M.K.
More informationc 2015 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vo. 53, No. 1, pp. 43 62 c 2015 Society for Industria and Appied Mathematics CONVERGENCE AND OPTIMALITY OF ADAPTIVE LEAST SQUARES FINITE ELEMENT METHODS CARSTEN CARSTENSEN AND EUN-JAE
More informationOn the Number of Limit Cycles for Discontinuous Generalized Liénard Polynomial Differential Systems
Internationa Journa of Bifurcation and Chaos Vo. 25 No. 10 2015 1550131 10 pages c Word Scientific Pubishing Company DOI: 10.112/S02181271550131X On the Number of Limit Cyces for Discontinuous Generaized
More informationA Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem
Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,
More informationIntroduction to Riemann Solvers
CO 5 BOLD WORKSHOP 2, 2 4 June Introduction to Riemann Sovers Oskar Steiner . Three major advancements in the numerica treatment of the hydrodynamic equations Three major progresses in computationa fuid
More informationHigh-order approximations to the Mie series for electromagnetic scattering in three dimensions
Proceedings of the 9th WSEAS Internationa Conference on Appied Mathematics Istanbu Turkey May 27-29 2006 (pp199-204) High-order approximations to the Mie series for eectromagnetic scattering in three dimensions
More informationarxiv: v4 [math.na] 25 Aug 2014
BIT manuscript No. (wi be inserted by the editor) A muti-eve spectra deferred correction method Robert Speck Danie Ruprecht Matthew Emmett Michae Minion Matthias Boten Rof Krause arxiv:1307.1312v4 [math.na]
More informationHigher dimensional PDEs and multidimensional eigenvalue problems
Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u
More informationA Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)
A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,
More informationMultiscale Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Problems UNCORRECTED PROOF
AQ Mutiscae Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Probems Yachin Efendiev, Juan Gavis, Raytcho Lazarov, Svetozar Margenov 2, 4 and Jun Ren 5 Department of Mathematics, TAMU,
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Heilbronner Str. 10, D Berlin - Wilmersdorf
Konrad-Zuse-Zentrum für Informationstechnik Berin Heibronner Str. 10, D-10711 Berin - Wimersdorf Raf Kornhuber Monotone Mutigrid Methods for Eiptic Variationa Inequaities II Preprint SC 93{19 (Marz 1994)
More informationScalable Spectrum Allocation for Large Networks Based on Sparse Optimization
Scaabe Spectrum ocation for Large Networks ased on Sparse Optimization innan Zhuang Modem R&D Lab Samsung Semiconductor, Inc. San Diego, C Dongning Guo, Ermin Wei, and Michae L. Honig Department of Eectrica
More informationP-MULTIGRID PRECONDITIONERS APPLIED TO HIGH-ORDER DG AND HDG DISCRETIZATIONS
6th European Conference on Computationa Mechanics (ECCM 6) 7th European Conference on Computationa Fuid Dynamics (ECFD 7) 11-15 June 2018, Gasgow, UK P-MULTIGRID PRECONDITIONERS APPLIED TO HIGH-ORDER DG
More informationA proposed nonparametric mixture density estimation using B-spline functions
A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),
More informationTwo Kinds of Parabolic Equation algorithms in the Computational Electromagnetics
Avaiabe onine at www.sciencedirect.com Procedia Engineering 9 (0) 45 49 0 Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More informationEfficient Visual-Inertial Navigation using a Rolling-Shutter Camera with Inaccurate Timestamps
Efficient Visua-Inertia Navigation using a Roing-Shutter Camera with Inaccurate Timestamps Chao X. Guo, Dimitrios G. Kottas, Ryan C. DuToit Ahmed Ahmed, Ruipeng Li and Stergios I. Roumeiotis Mutipe Autonomous
More informationA Solution to the 4-bit Parity Problem with a Single Quaternary Neuron
Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria
More informationAN A POSTERIORI ANALYSIS OF C 0 INTERIOR PENALTY METHODS FOR THE OBSTACLE PROBLEM OF CLAMPED KIRCHHOFF PLATES
AN A POSERIORI ANALYSIS OF C 0 INERIOR PENALY MEHODS FOR HE OBSACLE PROBLEM OF CLAMPED KIRCHHOFF PLAES SUSANNE C. BRENNER, JOSCHA GEDICKE, LI-YENG SUNG, AND YI ZHANG Abstract. We deveop an a posteriori
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationProceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds.
Proceedings of the 2012 Winter Simuation Conference C. aroque, J. Himmespach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds. TIGHT BOUNDS FOR AMERICAN OPTIONS VIA MUTIEVE MONTE CARO Denis Beomestny Duisburg-Essen
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationThe Binary Space Partitioning-Tree Process Supplementary Material
The inary Space Partitioning-Tree Process Suppementary Materia Xuhui Fan in Li Scott. Sisson Schoo of omputer Science Fudan University ibin@fudan.edu.cn Schoo of Mathematics and Statistics University of
More informationEXACT CLOSED FORM FORMULA FOR SELF INDUC- TANCE OF CONDUCTOR OF RECTANGULAR CROSS SECTION
Progress In Eectromagnetics Research M, Vo. 26, 225 236, 22 EXACT COSED FORM FORMUA FOR SEF INDUC- TANCE OF CONDUCTOR OF RECTANGUAR CROSS SECTION Z. Piatek, * and B. Baron 2 Czestochowa University of Technoogy,
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More information