P016 Toward Gauss-Newton and Exact Newton Optimization for Full Waveform Inversion
|
|
- Audra Carroll
- 5 years ago
- Views:
Transcription
1 P016 Toward Gauss-Newton and Exact Newton Optiization for Full Wavefor Inversion L. Métivier* ISTerre, R. Brossier ISTerre, J. Virieux ISTerre & S. Operto Géoazur SUMMARY Full Wavefor Inversion FWI applications classically rely on efficient first-order optiization schees, as the steepest descent or the nonlinear conjugate gradient optiization. However, second-order inforation provided by the Hessian atrix is proven to give a useful help in the scaling of the FWI proble and in the speed-up of the optiization. In this study, we propose an efficient atrix-free Hessian-vector foralis, that should allow to tackle Gauss-Newton GN and Exact-Newton EN optiization for large and realistic FWI targets. Our ethod relies on general second order adjoint forulas, based on a Lagrangian foralis. These forulas yield the possibility of coputing Hessianvector products at the cost of 2 forward siulations per shot. In this context, the coputational cost per shot of one GN or one EN nonlinear iteration aounts to the resolution of 2 forward siulations for the coputation of the gradient plus 2 forward siulations per inner linear conjugate gradient iteration. A nuerical test is provided, ephasizing the possible iproveent of the resolution when accounting for the exact Hessian in the inversion algorith.
2 Introduction Full Wavefor Inversion FWI is becoing an efficient tool to derive high resolution quantitative odels of the subsurface paraeters. The ethod relies on the iniization, through an iterative procedure, of the residual between recorded data and synthetic data coputed by solving the two-way wave equation in a subsurface odel. The growth of available coputational resources and recent developents of the ethod akes now possible applications to 2D and 3D data in the acoustic approxiation see for exaple Prieux et al., 2011; Plessix et al., 2010 and even in the elastic approxiation Brossier et al., Most of the FWI applications rely on fast optiization schees as preconditioned steepest descent or preconditioned conjugate-gradient ethods PCG. Second-order inforation provided by the Hessian is often neglected in FWI, due to the high coputational cost to build this atrix and solve the noral equation syste. However, a significant iproveent of the results can be obtained using this inforation: Pratt et al have shown the iproved resolution of Gauss-Newton ethod copared to the steepest descent one in a canonical application. Hu et al have shown results iproveent provided by a non-diagonal truncated Hessian in PCG. Brossier et al have shown the estiated Hessian s ipact of a quasi-newton l-bfgs Nocedal, 1980, on iage resolution and convergence speed copared to PCG. Epanoeritakis et al. 2008; Fichtner and Trapert 2011 have also discussed the interest of Hessian for inversion and uncertainty estiation, and also the prohibitive cost of Hessian coputation and storage. In this study, we develop the atheatical fraework to propose an efficient atrix-free Hessian-vector product algorith for FWI, and give an illustration of the interest of accounting for the exact Hessian in the inversion process. The final ai is to tackle Gauss-Newton and Full-Newton ethod for large scale applications. Our developent relies on a general second-order adjoint-state forula, valid either in the tie or in the frequency doain. Proble stateent We consider the forward proble equation Spu = ϕ, 1 where p denotes the subsurface paraeters odel space, Sp denotes the linear forward proble operator corresponding to the two-way wave equation discretization 1, ϕ is the source vector, and u is the wavefield vector. In the following, we consider that p and u have N coponents 2. These notations are general and can be applied either in the tie doain or in the frequency doain. The FWI proble is expressed as the least-square proble in f p = 1 p 2 N s R s u s p d s 2, 2 where d s and u s p are respectively the recorded dataset and the solution of the forward proble associated with the source ϕ s, N s is the total nuber of sources in the dataset, and R s is a restriction operator that aps the wavefield u s to the receivers locations. The isfit gradient is f p = R J s R s R s u s p d s 3 1 Note that the operator Sp depends non-linearly on p 2 This siplification is satisfied whenever u is discretized on the sae grid than p. However, the forulas can be straightforwardly extended to the situation where p and u do not have the sae nuber of coponents.
3 where J s p denotes the Jacobian atrix p u s p. The coplex conjugate transpose operator is denoted by the sybol and the real part application by R. The Hessian operator is Hp = R J s R s R s J s + [R s R s u s p d s ] j H s j, H s j = 2 ppu s j p, 4 and its Gauss-Newton approxiation is Bp = R J s R s R s J s 5 In our study, the proble 2 is solved using a Newton or a Gauss-Newton algorith, a local iterative optiization ethod that coputes a sequence p k fro an initial guess p 0 using the update forula where d k is the solution of p k+1 = p k + α k d k, 6 Hp k d k = f p k, or Bp k d k = f p k in the Gauss-Newton approxiation, 7 and α k is coputed through a globalization procedure linesearch or trust region. Explicit coputation and storage of atrices Jp, Hp and Bp is prohibitive for large scale probles and realistic FWI applications. Hence, the equation 7 should be solved using atrix-free linear iterative solvers, as proposed in Epanoeritakis et al. 2008; Fichtner and Trapert 2011, such as the conjugate gradient CG ethod. This requires to copute efficiently both the gradient f p and Hessian-vector products Hpv or Bpv for arbitrary vectors v. While the coputation of f p can be perfored efficiently through the classical adjoint-state forula Plessix, 2006 where λp the adjoint state is the solution of f p i = R pi Spup,λp, i = 1,..., 8 Sp λ = R d Ru, 9 we propose general second-order adjoint-state forula for the coputation of Hpv and Bpv related to the Lagrangian foralis and the nonlinear constrained optiization theory. Coputing Hessian vector products For the sake of clarity, we consider in the following that N s = 1 and we drop index s. Forulas for N s > 1 are directly obtained by suation. First, we define the function g v p such that g v p = f p,v = R J R Rup d,v, 10 where up is the solution of 1. We have g v p = Hpv. We define the Lagrangian function associated with the functional g v p L v p,u,α,λ,µ = R R Ru d,α + R Spu ϕ, µ + R Spα + v j p j Spu,λ 11 The Lagrangian L v is coposed of three ters: the first one accounts for the function g v, the second one accounts for the constraints on the wavefield u, solution of the forward proble, the third one accounts for the constraints on the first-order derivatives of the wavefield u with respect to p. For ũ and α such that Spũ = ϕ, Sp α = v j p j Spũ, 12
4 we have g v p = p L v p,ũ, α,λ,µ + u L v p,ũ, α,λ,µ p ũp + α L v p,ũ, α,λ,µ p αp. 13 We define λ and µ such that We have and Hpv i = R u L v p,ũ, α, λ, µ Sp µ = R R α pi Spũ, µ + = 0, α L v p,ũ, α, λ, µ = v j p j Sp λ, Sp λ = R Rũ d 15 p pi Sp α, λ + v j j pi Sp ũ, λ,i = 1,...,. 16 Note that λ corresponds to the adjoint state defined for the coputation of f. In addition, it can be proved that the coputation of Bpv aounts to setting λ to 0 in equations 15 and 16. Bpv i = R pi Spũ, µ,i = 1,..., with Sp µ = R R α, 17 The coputation of one atrix vector product Hpv or Bpv thus requires to solve one additional forward proble for α and one additional adjoint proble for µ, as reported in Epanoeritakis et al. 2008; Fichtner and Trapert For N s > 1, the overall coputation cost is ultiplied by N s. Nuerical results We consider the estiation of the pressure wave velocity v p using a 2D acoustic FWI algorith in the frequency doain. The exact odel v p is coposed of a hoogeneous background v 0 p = 1500.s 1 defined on a 2000 length square, and two inclusions where v p = 4500.s 1. PML are all around. The two inclusions are distant only fro 40 see fig.1. The high velocity contrast generates high aplitude ulti-scattered waves. We use a full acquisition with four lines of 29 sources placed each 50 on each side of the doain. Each source is associated with four lines of 29 receivers placed each 50 all around the doain. The odel is discretized over a grid with a spatial step of 20. The initial guess is the background v 0 p. We use one dataset corresponding to the frequency of 5 Hz. The associated average wavelength is λ 300. At this frequency, the distance between the two inclusions is saller than the expected resolution of classical FWI algoriths based on l-bfgs or Gauss-Newton approxiation. We copare the results obtained perforing 50 iterations of the l-bfgs algorith, 20 iterations of Gauss-Newton inversion and 20 iterations of Exact-Newton inversion 3. The corresponding estiated odels are plot in Figure 1. Figure 1 Acoustic FWI for pressure wave velocity. Fro left to right, exact odel, l-bfgs result, Gauss-Newton result, Exact-Newton result. 3 This corresponds approxiately to the sae coputation cost for 2.5 conjugate gradient inner iterations per Newton nonlinear iteration
5 These results ephasize the role of the second-order part of the Hessian, which is neglected in the Gauss-Newton approxiation, and hardly estiated in the l-bfgs approxiation. As entioned by different authors, this part of the Hessian allows to account for double scattering during the inversion Pratt et al., 1998; Virieux and Operto, An enhanceent of the resolution is obtained: only when using the exact Hessian, the two inclusions are identified. In this case, the noralized isfit is In the case of the Gauss-Newton inversion, the final noralized isfit is around 0.25: neglecting the double scattered waves prevent the optiizer to converge. In the case of the l-bfgs inversion, the final noralized isfit is around 10 3 : the l-bfgs Hessian approxiation allows to account for a part of the scattered wavefield. Conclusions and perspectives The preliinary results we have obtained illustrate how accounting for the Hessian effect can iprove the FWI results. The second-order adjoint forula is an efficient tool to ipleent Gaus-Newton or Exact-Newton algoriths in a atrix-free fashion. The next step consists now in developing globalization ethods either based on linesearch or trust region that render the overall ethod copetitive with the l-bfgs ethod. Two ain probles have to be tackled: the first one is the definition of a proper criterion to stop the inner conjugate gradient iterations, in order to iniize the coputation cost. The second one is the developent of efficient preconditioners to speed-up the convergence of the inner conjugate gradient loop. Aong several possibilities, we are interested in the Newton Steihaug algorith Steihaug, 1983, based on the trust-region ethod, that provides naturally a stopping criterion for the inner iterations. In addition, for each resolution of the inner linear systes, a l-bfgs preconditioner can be coputed siultaneously, which is applied to the inner linear syste raised by the next outer nonlinear iteration. Acknowledgeents This research was funded by the SEISCOPE consortiu sponsored by BP, CGG-VERITAS, ENI, EXXON MOBIL, PETROBRAS, SAUDI ARAMCO, SHELL, STATOIL and TOTAL. The linear systes were solved with the MUMPS package. This work was perfored by accessing to the high-perforance coputing facilities of CIMENT Université de Grenoble, France and to the HPC resources of GENCI- CINES under Grant References Brossier, R., Operto, S. and Virieux, J. [2009] Seisic iaging of coplex onshore structures by 2D elastic frequency-doain full-wavefor inversion. Geophysics, 746, WCC63 WCC76, doi: / Epanoeritakis, I., Akçelik, V., Ghattas, O. and Bielak, J. [2008] A Newton-CG ethod for large-scale threediensional elastic full wavefor seisic inversion. Inverse Probles, 24, Fichtner, A. and Trapert, J. [2011] Hessian kernels of seisic data functionals based upon adjoint techniques. Geophysical Journal International, 1852, , ISSN X, doi: /j x x. Hu, W., Abubakar, A., Habashy, T.M. and Liu, J. [2011] Preconditioned non-linear conjugate gradient ethod for frequency doain full-wavefor seisic inversion. Geophysical Prospecting, 593, , ISSN , doi: /j x. Nocedal, J. [1980] Updating Quasi-Newton Matrices With Liited Storage. Matheatics of Coputation, 35151, Plessix, R.E. [2006] A review of the adjoint-state ethod for coputing the gradient of a functional with geophysical applications. Geophysical Journal International, 1672, Plessix, R.E., Baeten, G., de Maag, J.W., Klaassen, M., Rujie, Z. and Zhifei, T. [2010] Application of acoustic full wavefor inversion to a low-frequency large-offset land data set. SEG Technical Progra Expanded Abstracts, 291, , doi: / Pratt, R.G., Shin, C. and Hicks, G.J. [1998] Gauss-Newton and full Newton ethods in frequency-space seisic wavefor inversion. Geophysical Journal International, 133, Prieux, V. et al. [2011] On the footprint of anisotropy on isotropic full wavefor inversion: the Valhall case study. Geophysical Journal International, 187, , doi:doi: /j X x, Steihaug, T. [1983] The conjugate gradient ethod and trust regions in large scale optiization. SIAM Journal on Nuerical Analysis, 20, Virieux, J. and Operto, S. [2009] An overview of full wavefor inversion in exploration geophysics. Geophysics, 746, WCC127 WCC152.
The truncated Newton method for Full Waveform Inversion
The truncated Newton method for Full Waveform Inversion Ludovic Métivier, Romain Brossier, Jean Virieux, Stéphane Operto To cite this version: Ludovic Métivier, Romain Brossier, Jean Virieux, Stéphane
More informationOptimization schemes for Full Waveform Inversion: the preconditioned truncated Newton method
Optimization schemes for Full Waveform Inversion: the preconditioned truncated Newton method Ludovic Métivier, Romain Brossier, Jean Virieux, Stéphane Operto To cite this version: Ludovic Métivier, Romain
More informationSuppress Parameter Cross-talk for Elastic Full-waveform Inversion: Parameterization and Acquisition Geometry
Suress Paraeter Cross-talk for Elastic Full-wavefor Inversion: Paraeterization and Acquisition Geoetry Wenyong Pan and Kris Innanen CREWES Project, Deartent of Geoscience, University of Calgary Suary Full-wavefor
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationA Nonlinear Sparsity Promoting Formulation and Algorithm for Full Waveform Inversion
A Nonlinear Sparsity Prooting Forulation and Algorith for Full Wavefor Inversion Aleksandr Aravkin, Tristan van Leeuwen, Jaes V. Burke 2 and Felix Herrann Dept. of Earth and Ocean sciences University of
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More informationSUMMARY REVIEW OF THE FREQUENCY DOMAIN L2 FWI-HESSIAN
Efficient stochastic Hessian estimation for full waveform inversion Lucas A. Willemsen, Alison E. Malcolm and Russell J. Hewett, Massachusetts Institute of Technology SUMMARY In this abstract we present
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More information3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University
More informationVariations on Backpropagation
2 Variations on Backpropagation 2 Variations Heuristic Modifications Moentu Variable Learning Rate Standard Nuerical Optiization Conjugate Gradient Newton s Method (Levenberg-Marquardt) 2 2 Perforance
More informationW011 Full Waveform Inversion for Detailed Velocity Model Building
W011 Full Waveform Inversion for Detailed Velocity Model Building S. Kapoor* (WesternGeco, LLC), D. Vigh (WesternGeco), H. Li (WesternGeco) & D. Derharoutian (WesternGeco) SUMMARY An accurate earth model
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationP032 3D Seismic Diffraction Modeling in Multilayered Media in Terms of Surface Integrals
P032 3D Seisic Diffraction Modeling in Multilayered Media in Ters of Surface Integrals A.M. Aizenberg (Institute of Geophysics SB RAS, M. Ayzenberg* (Norwegian University of Science & Technology, H.B.
More informationINTRODUCTION. Residual migration has proved to be a useful tool in imaging and in velocity analysis.
Stanford Exploration Project, Report, June 3, 999, pages 5 59 Short Note On Stolt prestack residual igration Paul Sava keywords: Stolt, residual igration INTRODUCTION Residual igration has proved to be
More informationŞtefan ŞTEFĂNESCU * is the minimum global value for the function h (x)
7Applying Nelder Mead s Optiization Algorith APPLYING NELDER MEAD S OPTIMIZATION ALGORITHM FOR MULTIPLE GLOBAL MINIMA Abstract Ştefan ŞTEFĂNESCU * The iterative deterinistic optiization ethod could not
More informationSharp Time Data Tradeoffs for Linear Inverse Problems
Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used
More informationA projected Hessian for full waveform inversion
CWP-679 A projected Hessian for full waveform inversion Yong Ma & Dave Hale Center for Wave Phenomena, Colorado School of Mines, Golden, CO 80401, USA (c) Figure 1. Update directions for one iteration
More informationGeneralized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,
Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationKernel Methods and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic
More informationRecovering Data from Underdetermined Quadratic Measurements (CS 229a Project: Final Writeup)
Recovering Data fro Underdeterined Quadratic Measureents (CS 229a Project: Final Writeup) Mahdi Soltanolkotabi Deceber 16, 2011 1 Introduction Data that arises fro engineering applications often contains
More informationOn Constant Power Water-filling
On Constant Power Water-filling Wei Yu and John M. Cioffi Electrical Engineering Departent Stanford University, Stanford, CA94305, U.S.A. eails: {weiyu,cioffi}@stanford.edu Abstract This paper derives
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationANALYSIS OF REFLECTOR AND HORN ANTENNAS USING MULTILEVEL FAST MULTIPOLE ALGORITHM
European Congress on Coputational Methods in Applied Sciences and Engineering ECCOMAS 2 Barcelona, 11-14 Septeber 2 ECCOMAS ANALYSIS OF REFLECTOR AND HORN ANTENNAS USING MULTILEVEL FAST MULTIPOLE ALGORITHM
More informationThe linear sampling method and the MUSIC algorithm
INSTITUTE OF PHYSICS PUBLISHING INVERSE PROBLEMS Inverse Probles 17 (2001) 591 595 www.iop.org/journals/ip PII: S0266-5611(01)16989-3 The linear sapling ethod and the MUSIC algorith Margaret Cheney Departent
More informationPattern Recognition and Machine Learning. Artificial Neural networks
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lessons 7 20 Dec 2017 Outline Artificial Neural networks Notation...2 Introduction...3 Key Equations... 3 Artificial
More informationQ ESTIMATION WITHIN A FORMATION PROGRAM q_estimation
Foration Attributes Progra q_estiation Q ESTIMATION WITHIN A FOMATION POGAM q_estiation Estiating Q between stratal slices Progra q_estiation estiate seisic attenuation (1/Q) on coplex stratal slices using
More informationA model reduction approach to numerical inversion for a parabolic partial differential equation
Inverse Probles Inverse Probles 30 (204) 250 (33pp) doi:0.088/0266-56/30/2/250 A odel reduction approach to nuerical inversion for a parabolic partial differential equation Liliana Borcea, Vladiir Drusin
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationDepartment of Electronic and Optical Engineering, Ordnance Engineering College, Shijiazhuang, , China
6th International Conference on Machinery, Materials, Environent, Biotechnology and Coputer (MMEBC 06) Solving Multi-Sensor Multi-Target Assignent Proble Based on Copositive Cobat Efficiency and QPSO Algorith
More informationSymbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm
Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat
More informationSimultaneous estimation of wavefields & medium parameters
Simultaneous estimation of wavefields & medium parameters reduced-space versus full-space waveform inversion Bas Peters, Felix J. Herrmann Workshop W- 2, The limit of FWI in subsurface parameter recovery.
More informationCHAPTER 8 CONSTRAINED OPTIMIZATION 2: SEQUENTIAL QUADRATIC PROGRAMMING, INTERIOR POINT AND GENERALIZED REDUCED GRADIENT METHODS
CHAPER 8 CONSRAINED OPIMIZAION : SEQUENIAL QUADRAIC PROGRAMMING, INERIOR POIN AND GENERALIZED REDUCED GRADIEN MEHODS 8. Introduction In the previous chapter we eained the necessary and sufficient conditions
More informationIntelligent Systems: Reasoning and Recognition. Artificial Neural Networks
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley MOSIG M1 Winter Seester 2018 Lesson 7 1 March 2018 Outline Artificial Neural Networks Notation...2 Introduction...3 Key Equations... 3 Artificial
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationAnalytical solution for the electric potential in arbitrary anisotropic layered media applying the set of Hankel transforms of integer order
Geophysical Prospecting, 2006, 54, 651 661 Analytical solution for the electric potential in arbitrary anisotropic layered edia applying the set of Hankel transfors of integer order E. Pervago 1,A.Mousatov
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationComparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations
International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah
More informationStability Analysis of the Matrix-Free Linearly Implicit 2 Euler Method 3 UNCORRECTED PROOF
1 Stability Analysis of the Matrix-Free Linearly Iplicit 2 Euler Method 3 Adrian Sandu 1 andaikst-cyr 2 4 1 Coputational Science Laboratory, Departent of Coputer Science, Virginia 5 Polytechnic Institute,
More informationUBC-GIF: Capabilities for EM Modelling and Inversion of LSBB data
The University of British Colubia Geophysical Inversion Facility slide 1 UBC-GIF: Capabilities for EM Modelling and Inversion of LSBB data Douglas W. Oldenburg Departent of Earth and Ocean Sciences June
More informationLow-complexity, Low-memory EMS algorithm for non-binary LDPC codes
Low-coplexity, Low-eory EMS algorith for non-binary LDPC codes Adrian Voicila,David Declercq, François Verdier ETIS ENSEA/CP/CNRS MR-85 954 Cergy-Pontoise, (France) Marc Fossorier Dept. Electrical Engineering
More informationIterative Linear Solvers and Jacobian-free Newton-Krylov Methods
Eric de Sturler Iterative Linear Solvers and Jacobian-free Newton-Krylov Methods Eric de Sturler Departent of Matheatics, Virginia Tech www.ath.vt.edu/people/sturler/index.htl sturler@vt.edu Efficient
More informationPattern Recognition and Machine Learning. Artificial Neural networks
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2016 Lessons 7 14 Dec 2016 Outline Artificial Neural networks Notation...2 1. Introduction...3... 3 The Artificial
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informationRESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS
BIT Nuerical Matheatics 43: 459 466, 2003. 2003 Kluwer Acadeic Publishers. Printed in The Netherlands 459 RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS V. SIMONCINI Dipartiento di
More informationA method to determine relative stroke detection efficiencies from multiplicity distributions
A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,
More informationAn improved self-adaptive harmony search algorithm for joint replenishment problems
An iproved self-adaptive harony search algorith for joint replenishent probles Lin Wang School of Manageent, Huazhong University of Science & Technology zhoulearner@gail.co Xiaojian Zhou School of Manageent,
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationFull-Waveform Inversion with Gauss- Newton-Krylov Method
Full-Waveform Inversion with Gauss- Newton-Krylov Method Yogi A. Erlangga and Felix J. Herrmann {yerlangga,fherrmann}@eos.ubc.ca Seismic Laboratory for Imaging and Modeling The University of British Columbia
More informationInteractive Markov Models of Evolutionary Algorithms
Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary
More informationUncertainty quantification for Wavefield Reconstruction Inversion
Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),
More informationA MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION
A eshsize boosting algorith in kernel density estiation A MESHSIZE BOOSTING ALGORITHM IN KERNEL DENSITY ESTIMATION C.C. Ishiekwene, S.M. Ogbonwan and J.E. Osewenkhae Departent of Matheatics, University
More informationREDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION
ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationGLOBALLY CONVERGENT LEVENBERG-MARQUARDT METHOD FOR PHASE RETRIEVAL
GLOBALLY CONVERGENT LEVENBERG-MARQUARDT METHOD FOR PHASE RETRIEVAL CHAO MA, XIN LIU, AND ZAIWEN WEN Abstract. In this paper, we consider a nonlinear least squares odel for the phase retrieval proble. Since
More informationUNIVERSITY OF TRENTO ON THE USE OF SVM FOR ELECTROMAGNETIC SUBSURFACE SENSING. A. Boni, M. Conci, A. Massa, and S. Piffer.
UIVRSITY OF TRTO DIPARTITO DI IGGRIA SCIZA DLL IFORAZIO 3823 Povo Trento (Italy) Via Soarive 4 http://www.disi.unitn.it O TH US OF SV FOR LCTROAGTIC SUBSURFAC SSIG A. Boni. Conci A. assa and S. Piffer
More informationApproximate- vs. full-hessian in FWI: 1D analytical and numerical experiments
Approximate- vs. full-hessian in FWI: 1D analytical and numerical experiments Raul Cova and Kris Innanen ABSTRACT Feasibility of using Full Waveform Inversion (FWI) to build velocity models has been increasing
More informationMSEC MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL SOLUTION FOR MAINTENANCE AND PERFORMANCE
Proceeding of the ASME 9 International Manufacturing Science and Engineering Conference MSEC9 October 4-7, 9, West Lafayette, Indiana, USA MSEC9-8466 MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL
More informationEnvelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Damping Characteristics
Copyright c 2007 Tech Science Press CMES, vol.22, no.2, pp.129-149, 2007 Envelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Daping Characteristics D. Moens 1, M.
More information2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationAn Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles
An Inverse Interpolation Method Utilizing In-Flight Strain Measureents for Deterining Loads and Structural Response of Aerospace Vehicles S. Shkarayev and R. Krashantisa University of Arizona, Tucson,
More informationInspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information
Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub
More informationMulti-Scale/Multi-Resolution: Wavelet Transform
Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the
More informationA Generalized Permanent Estimator and its Application in Computing Multi- Homogeneous Bézout Number
Research Journal of Applied Sciences, Engineering and Technology 4(23): 5206-52, 202 ISSN: 2040-7467 Maxwell Scientific Organization, 202 Subitted: April 25, 202 Accepted: May 3, 202 Published: Deceber
More informationFull-waveform inversion application in different geological settings Denes Vigh*, Jerry Kapoor and Hongyan Li, WesternGeco
Full-waveform inversion application in different geological settings Denes Vigh*, Jerry Kapoor and Hongyan Li, WesternGeco Summary After the synthetic data inversion examples, real 3D data sets have been
More informationStochastic Subgradient Methods
Stochastic Subgradient Methods Lingjie Weng Yutian Chen Bren School of Inforation and Coputer Science University of California, Irvine {wengl, yutianc}@ics.uci.edu Abstract Stochastic subgradient ethods
More informationTopic 5a Introduction to Curve Fitting & Linear Regression
/7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline
More informationIAENG International Journal of Computer Science, 42:2, IJCS_42_2_06. Approximation Capabilities of Interpretable Fuzzy Inference Systems
IAENG International Journal of Coputer Science, 4:, IJCS_4 6 Approxiation Capabilities of Interpretable Fuzzy Inference Systes Hirofui Miyajia, Noritaka Shigei, and Hiroi Miyajia 3 Abstract Many studies
More informationpaper prepared for the 1996 PTRC Conference, September 2-6, Brunel University, UK ON THE CALIBRATION OF THE GRAVITY MODEL
paper prepared for the 1996 PTRC Conference, Septeber 2-6, Brunel University, UK ON THE CALIBRATION OF THE GRAVITY MODEL Nanne J. van der Zijpp 1 Transportation and Traffic Engineering Section Delft University
More informationThis model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.
CS 493: Algoriths for Massive Data Sets Feb 2, 2002 Local Models, Bloo Filter Scribe: Qin Lv Local Models In global odels, every inverted file entry is copressed with the sae odel. This work wells when
More informationEstimation of ADC Nonlinearities from the Measurement in Input Voltage Intervals
Estiation of ADC Nonlinearities fro the Measureent in Input Voltage Intervals M. Godla, L. Michaeli, 3 J. Šaliga, 4 R. Palenčár,,3 Deptartent of Electronics and Multiedia Counications, FEI TU of Košice,
More informationSupervised Baysian SAR image Classification Using The Full Polarimetric Data
Supervised Baysian SAR iage Classification Using The Full Polarietric Data (1) () Ziad BELHADJ (1) SUPCOM, Route de Raoued 3.5 083 El Ghazala - TUNSA () ENT, BP. 37, 100 Tunis Belvedere, TUNSA Abstract
More informationBernoulli Wavelet Based Numerical Method for Solving Fredholm Integral Equations of the Second Kind
ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol., No., 6, pp.-9 Bernoulli Wavelet Based Nuerical Method for Solving Fredhol Integral Equations of the Second Kind S. C. Shiralashetti*,
More informationIN modern society that various systems have become more
Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto
More information2D Laplace-Domain Waveform Inversion of Field Data Using a Power Objective Function
Pure Appl. Geophys. Ó 213 Springer Basel DOI 1.17/s24-13-651-4 Pure and Applied Geophysics 2D Laplace-Domain Waveform Inversion of Field Data Using a Power Objective Function EUNJIN PARK, 1 WANSOO HA,
More informationOn the Analysis of the Quantum-inspired Evolutionary Algorithm with a Single Individual
6 IEEE Congress on Evolutionary Coputation Sheraton Vancouver Wall Centre Hotel, Vancouver, BC, Canada July 16-1, 6 On the Analysis of the Quantu-inspired Evolutionary Algorith with a Single Individual
More informationFast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials
Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter
More informationAnisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys. 5 00 007 04 PII: S00-770867-6 Anisotropic reference edia and the possible linearized approxiations for phase
More informationNBN Algorithm Introduction Computational Fundamentals. Bogdan M. Wilamoswki Auburn University. Hao Yu Auburn University
NBN Algorith Bogdan M. Wilaoswki Auburn University Hao Yu Auburn University Nicholas Cotton Auburn University. Introduction. -. Coputational Fundaentals - Definition of Basic Concepts in Neural Network
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationA general forulation of the cross-nested logit odel Michel Bierlaire, Dpt of Matheatics, EPFL, Lausanne Phone: Fax:
A general forulation of the cross-nested logit odel Michel Bierlaire, EPFL Conference paper STRC 2001 Session: Choices A general forulation of the cross-nested logit odel Michel Bierlaire, Dpt of Matheatics,
More informationCS Lecture 13. More Maximum Likelihood
CS 6347 Lecture 13 More Maxiu Likelihood Recap Last tie: Introduction to axiu likelihood estiation MLE for Bayesian networks Optial CPTs correspond to epirical counts Today: MLE for CRFs 2 Maxiu Likelihood
More informationNumerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk
More informationESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics
ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents
More information146 THE LEADING EDGE February 2016 Special Section: Imaging/inversion: Estimating the earth model
A downscaling strategy from FWI to microscale reservoir properties from high-resolution images B. Dupuy 1, A. Asnaashari 2, R. Brossier 3, S. Garambois 3, L. Métivier 3, A. Ribodetti 4, and J. Virieux
More informationThe Methods of Solution for Constrained Nonlinear Programming
Research Inventy: International Journal Of Engineering And Science Vol.4, Issue 3(March 2014), PP 01-06 Issn (e): 2278-4721, Issn (p):2319-6483, www.researchinventy.co The Methods of Solution for Constrained
More informationWavefield Reconstruction Inversion (WRI) a new take on wave-equation based inversion Felix J. Herrmann
Wavefield Reconstruction Inversion (WRI) a new take on wave-equation based inversion Felix J. Herrmann SLIM University of British Columbia van Leeuwen, T and Herrmann, F J (2013). Mitigating local minima
More informationGREY FORECASTING AND NEURAL NETWORK MODEL OF SPORT PERFORMANCE
Journal of heoretical and Applied Inforation echnology 3 st March 03 Vol 49 No3 005-03 JAI & LLS All rights reserved ISSN: 99-8645 wwwatitorg E-ISSN: 87-395 GREY FORECASING AND NEURAL NEWORK MODEL OF SPOR
More informationOn the approximation of Feynman-Kac path integrals
On the approxiation of Feynan-Kac path integrals Stephen D. Bond, Brian B. Laird, and Benedict J. Leikuhler University of California, San Diego, Departents of Matheatics and Cheistry, La Jolla, CA 993,
More informationEE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng
EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel
More informationEnsemble Based on Data Envelopment Analysis
Enseble Based on Data Envelopent Analysis So Young Sohn & Hong Choi Departent of Coputer Science & Industrial Systes Engineering, Yonsei University, Seoul, Korea Tel) 82-2-223-404, Fax) 82-2- 364-7807
More information+ -d-t-' )=1. = vpi. Aportaciones Matematicas Comunicaciones 17 (1996) 5-10.
Aportaciones Mateaticas Counicaciones 17 (1996) 5-10. 1. A suary of the proble Much of the processing that is used in the petroleu industry requires the consideration of a large nuber of cheical reactions.
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationTracking using CONDENSATION: Conditional Density Propagation
Tracking using CONDENSATION: Conditional Density Propagation Goal Model-based visual tracking in dense clutter at near video frae rates M. Isard and A. Blake, CONDENSATION Conditional density propagation
More informationGeneralized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials
Copyright 7 Tech Science Press CMC vol.53 no.4 pp.37-34 7 Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials S.D. Akbarov and M. Negin 3 Abstract: Dispersion of
More informationOptimal Resource Allocation in Multicast Device-to-Device Communications Underlaying LTE Networks
1 Optial Resource Allocation in Multicast Device-to-Device Counications Underlaying LTE Networks Hadi Meshgi 1, Dongei Zhao 1 and Rong Zheng 2 1 Departent of Electrical and Coputer Engineering, McMaster
More information