Stability Analysis of the Matrix-Free Linearly Implicit 2 Euler Method 3 UNCORRECTED PROOF
|
|
- Randolph Turner
- 5 years ago
- Views:
Transcription
1 1 Stability Analysis of the Matrix-Free Linearly Iplicit 2 Euler Method 3 Adrian Sandu 1 andaikst-cyr Coputational Science Laboratory, Departent of Coputer Science, Virginia 5 Polytechnic Institute, Blacksburg, VA, 24060, USA. E-ail sandu@cs.vt.edu 6 2 National Center for Atospheric Research, 1850 Table Mesa Drive, Boulder CO, E-ail aik@ucar.edu. 8 Suary. Iplicit tie stepping ethods are useful for the siulation of large scale PDE 9 systes because they avoid the tie step liitations iposed by explicit stability conditions. 10 To alleviate the challenges posed by coputational and eory constraints, any applica- 11 tions solve the resulting linear systes by iterative ethods where the Jacobian-vector prod- 12 ucts are approxiated by finite differences. This paper explains the relation between a linearly 13 iplicit Euler ethod, solved using a Jacobian-free Krylov ethod, and explicit Runge-Kutta 14 ethods. The case with preconditioning is equivalent to a Rosenbrock-W ethod where the 15 approxiate Jacobian, inverted at each stage, corresponds directly to the preconditioner. The 16 accuracy of the resulting Runge-Kutta ethods can be controlled by constraining the Krylov 17 solution. Nuerical experients confir the theoretical findings Introduction 19 Large systes of tie dependent partial differential equations PDEs), arising in 20 ulti-physics siulations, are often discretized using the ethod of lines approach. 21 The independent tie and space nuerical schees allow the coupling of ultiple 22 physics odules, and provide axiu flexibility in choosing appropriate algo- 23 riths. After the sei-discretization in space the syste of PDEs is reduced to a 24 syste of ordinary differential equations ODEs) 25 y = f y), t 0 t t f, yt 0 )=y 0. 1) Here yt) R d is the solution vector and y 0 the initial condition. We denote the 26 Jacobian of the ODE function by Jy) = f y y) R d d, and the identity atrix by 27 I R d d. 28 Stability requireents e.g., the CFL condition for discretized hyperbolic PDEs) 29 liit the tie steps allowable by explicit tie discretizations of 1). When the fastest 30 tie scales in the syste 1) are short, e.g., in the presence of fast waves, the stability 31 condition iposes tie steps uch saller then those required to achieve the target 32 R. Bank et al. eds.), Doain Decoposition Methods in Science and Engineering XX, Lecture Notes in Coputational Science and Engineering 91, DOI / , Springer-Verlag Berlin Heidelberg 2013 Page 427
2 Adrian Sandu and Aik St-Cyr accuracy. The step size liitation by linear stability conditions is referred to as stiff- 33 ness. In order to overcoe this coputational inefficiency, it is desirable to use i- 34 plicit, unconditionally stable discretizations which allow arbitrarily large tie steps 35 [2]. Iplicit ethods have a high cost per step due to the need to solve a non)linear 36 syste of equations. 37 To reduce the coputational and eory costs of direct linear syste solvers, 38 and to aid parallelization, iterative Krylov space ethods are eployed. Further- 39 ore, atrix-free ipleentations approxiate Jacobian vector products by finite 40 differences [4]. This approach avoids additional coding for the Jacobian, preserves 41 the parallel scalability of the explicit odel, and has becoe popular in any appli- 42 cations, e.g., [1, 5, 6]. The hope is that the properties of the iplicit tie discretiza- 43 tion reain unaltered, provided that the iterative solutions are carried out to sufficient 44 accuracy. We show here that the atrix-free approach does alter the properties of the 45 underlying iplicit tie stepping ethod. 46 This study treats a linearly iplicit ethod, together with the Krylov subspace 47 iterations for solving the linear syste, as a single nuerical schee. The analysis 48 reveals that atrix-free ipleentations of linearly iplicit ethods are equivalent 49 to explicit Runge Kutta ethods. Consequently, the unconditional stability property 50 of the base ethod is lost. When preconditioning is used, the atrix-free iplicit 51 ethods are equivalent to Rosenbrock-W ROS-W) ethods where the approxiate 52 Jacobians correspond directly to the preconditioners The Matrix-Free Linearly Iplicit Euler Method 54 Consider the linearly iplicit Euler LIE) ethod applied to 1) 55 I ΔtJyn ) ) w = f y n ), y n+1 = y n + Δt w. 2) When the linear syste is solved exactly odulo roundoff errors) by LU factoriza- 56 tion the ethod 2) is unconditionally stable, and thus suitable for the solution of 57 stiff systes. For any PDEs sei-discretized in the ethod of lines fraework, 58 however, the diension of the linear syste 2) is very large, and the coputational 59 and eory costs associated with a direct solution are prohibitive. Moreover, the 60 construction of the explicit Jacobian atrix J is difficult when the space discretiza- 61 tion is based on a doain decoposition approach. To alleviate these probles, a 62 popular approach is to solve 2) by atrix-free iterative ethods. We seek to analyze 63 the ipact that this approxiate solutions have on the stability and accuracy of the 64 iplicit tie stepping schee. Our approach is to treat the original discretization 65 2) together with the iterations as a single nuerical ethod applied to solve the 66 ODE 1). 67 To be specific, we solve the linear syste in 2) by a Krylov space ethod. The 68 initial guess is y n+1 = y n, i.e., w = 0. After iterations the following -diensional 69 Krylov space is built: 70 { K = span f y n ),..., I ΔtJy n ) ) } 1 f yn ). Page 428
3 Stability Analysis of the Matrix-Free Linearly Iplicit Euler Method In the atrix-free approach, the basis is constructed recursively and the Jacobian- 71 vector products are approxiated by finite differences 72 l i = l i 1 Δt ε 1 f y n + ε l i 1 )+Δt ε 1 l 1, i = 2,...,. 3) We assue that the sae scaling factor ε is used to copute the finite differences in 73 all iterations. The analysis can be easily extended to the case where a different ε is 74 used in each iteration.) Denote 75 k 1 = f y n ); k i = f y n + εl i 1 ), i = 2,,. 4) The recurrence 3) can be expressed in ters of k i as: 76 k i = f y n + Δt Δt 1 ε +i 2) ) i 1 k 1 Δt k j ), i = 2,...,. 5) The solution w = α i l i K can be expressed in ters of k i s: 77 ) ) w = α i + Δt ε 1 i 1)α i k 1 Δt ε 1 α j k i. 6) i=2 j=2 i=2 j=i Equations 5)and6), together with the relation y n+1 = y n +Δt w, are copared with 78 the -stage explicit Runge Kutta ERK) ethod [3] 79 i 1 k i = f y n + a ij k j ), i = 1,...,; y n+1 = y n + Δt b i k i. j=1 The coparison reveals the following. 80 Theore 1. The atrix-free LIE 2) ethod is equivalent to an explicit Runge Kutta 81 ethod. The nuber of Krylov iterations defines the nuber of Runge Kutta stages. 82 Equations 5)and6) define the coefficients of the ERK ethod: 83 a i,1 = Δt 1 ε +i 2); a i, j = 1, for i = 2,,, j = 2,,i 1; b 1 = j + Δt ε j=1α 1 j 1)α j ; b i = Δt ε 1 α j, i = 2,...,. j=2 j=i 2.1 Stability Considerations 84 The solution of the linear syste under the initial guess w = 0) is part of the Krylov 85 space K and can be represented by a atrix polynoial 86 w = p 1 I Δt Jyn ) ) f y n ). 87 The atrix-free LIE ethod applied to the Dahlquist test proble y = λ y, y0)=1, 88 gives the following solution: 89 y n+1 = y n + Δtw=1 + zp 1 1 z)) y n = Rz)y n, 90 with z = Δt λ. The stability function of the equivalent ERK ethod is the degree 91 polynoial Rz)=1 + zp 1 1 z). 92 Page 429
4 Adrian Sandu and Aik St-Cyr Theore 2. The stability region of the LIE ethod, with a Krylov atrix-free linear 93 solver, is necessarily finite. The unconditional stability of the original LIE ethod is 94 lost. 95 Siilar considerations hold for Krylov space ethods that use an orthogonal basis 96 of the Krylov space, built by Arnoldi iterations [7] Accuracy Considerations 98 The ethod accuracy is difficult to assess, as the coefficients depend on the tie 99 step. The relation between the finite difference scaling factor ε and the tie step Δt 100 is iportant in deterining accuracy. 101 Assue that the finite difference scaling factor is a constant fraction of the tie 102 step, ε/δt = const. This is a reasonable assuption: in order to increase accuracy 103 one decreases both Δt, to reduce the truncation error, and ε, to reduce the finite 104 difference error. Of course, for very sall ε the finite difference error becoes again 105 large due to roundoff.) Also assue that the coefficients α 1,...,α do not depend 106 on ε or Δt. 107 In this case the accuracy can be assessed using the classical approach. The order 108 conditions depend on the Krylov space coefficients α as follows: 109 Order 1: Order 2: b i = j=1 b i c i = α j = 1, i=2 7a) i 1)α i = b) Neither condition 7a) nor 7b) are autoatically satisfied by the Krylov iterative 110 ethods. In particular, 111 Lea 1. The first order accuracy of the atrix-free LIE is not autoatic when 112 ε/δt = const. Additional constraints need to be iposed on the Krylov solution 113 coefficients. 114 Consider now the case where ε is constant does not depend on Δt ). Assue 115 that the coefficients α 1,...,α do not depend on ε or Δt. A necessary condition for 116 the ethod to be accurate of order p is that its stability function approxiates the 117 exponential, Rz)=e z + O z p+1). The stability function does not depend on either 118 ε or Δt. The conditions 7a)and7b) on the Krylov solution coefficients α 1,...,α, 119 which are sufficient when ε = const Δt, see to be necessary in the case ε = const. 120 In the general case the Krylov solution coefficients α 1,...,α do depend on Δt. 121 For Δt 0wehavethatw f y n ) and therefore α 1 1, α 2,α 3, Asyp- 122 totically the condition 7a) holds. Moreover, the nuberof iterations also depends 123 on Δt through the convergence speed. Consequently, it is difficult to extend the clas- 124 sical accuracy analysis to atrix-free linearly iplicit ethods. It sees reasonable, 125 however, to odify the Krylov ethod and ipose at least condition 7a) onthe 126 Krylov coefficients. 127 Page 430
5 Stability Analysis of the Matrix-Free Linearly Iplicit Euler Method 3 Preconditioned Iterations 128 Consider the case where a preconditioner atrix M is used to speed up the iterations. 129 The linear syste 2) becoes 130 M 1 I ΔtJy n ) ) k = M 1 f y n ). 131 The Krylov space constructed in this case is 132 { K = span f y n )..., M 1 I ΔtJy n )) ) } 1 M 1 f y n ). In the atrix-free approach the following basis is constructed recursively 133 l 1 = M 1 f y n ), l i = M 1 l i 1 Δt ε 1 M 1 f y n + ε l i 1 )+Δt ε 1 l 1, i = 2,...,. Denote k 1 = Δt l 1 and k i = Δt l 1 ε l i for i = 2,,. Wehave 134 Mk 1 = Δt fy n ) 8) Mk i = Δt fy n + k 1 k i 1 )+k i 1 k 1, i = 2,...,. Consider, for coparison, a Rosenbrock-W ROW) ethod in the ipleentation- 135 friendly forulation [2, Sect. IV.7] 136 ) ] i 1 i 1 [I Δt γ Ĵ n k i = Δt γ f y n + a ij k j + γ c ij k j, y n+1 = y n + s j=1 j=1 i k i. 9) Here Ĵ n J y n ) is an approxiation of the exact Jacobian at the current step. We 137 identify the ethod coefficients γ = 1and 138 c i,1 = 1; c i,i 1 = 1; a i,1 = 1; a i,i 1 = 1, i = 2,,. Fro the solution w = α i l i = b i k i K we identify the weights 139 b 1 = α 1 Δt 1 + ε 1 α j ; j=2 b i = ε 1 α i, i = 2,...,. The preconditioner defines the Jacobian approxiation in the ROW ethod, 140 M = I Δt γ Ĵ n Ĵ n = Δt 1 I M). Theore 3. The preconditioned atrix-free LIE is equivalent to a linearly-iplicit 141 ROW ethod. The choice of the preconditioner, besides accelerating convergence, 142 iproves the stability of the atrix-free LIE ethod. The preconditioner defines the 143 Jacobian approxiation in the ROW ethod. 144 Note that the general approach can be applied to ROW ethods [2, Sect. IV.7] 145 by solving the linear syste of each stage with an iterative atrix free algorith. 146 The resulting schee is an explicit Runge Kutta ethod or a ROW ethod) with 147 s i stages. 148 Page 431
6 Adrian Sandu and Aik St-Cyr 4 Nuerical Results 149 Consider the one diensional scalar advection-diffusion equation 150 this figure will be printed in b/w u t +au) x = Du xx, ux,t = 0)=u 0 x). 10) A spectral discontinuous Galerkin spatial discretization is used with 20 eleents and 151 polynoials of order 8. The diffusive ter discretization is stabilized using the inter- 152 nal penalty ethod [8]. The LIE tie stepping is used with the atrix-free GMRES 153 solver [7]. 154 Iaginary Σ b i =0.60; M=4; CFL=2.05; Diff=1.0e Real a) D=0.001, =4 Iaginary Σ b i =0.94; M=10; CFL=2.05; Diff=1.0e Real b) D=0.001, =10 M=9; LIE Orders: Gauss=0.98 GMRES=0.98 MFree= CFL Error c) D=0:001; =9 Fig. 1. a) andb) The ERK stability regions for different nubers of GMRES iterations. c) The accuracy of the LIE schee using various approaches to invert the Jacobian atrix. The GMRES weights are restricted by 7b) such as to obtain a second order ethod. Advectiondiffusion equation 10), Δt = CFL, ε = 10 6 Δt In Fig.1a, b, the stability regions generated by the GMRES iterations are plotted 155 for a varying nuber of Krylov vectors. The regions grow quickly and encopass 156 the eigenvalues of the discrete advection-diffusion operator. Subsequent iterations 157 iprove solution accuracy but do not iprove linear stability. Additional experi- 158 ents not reported here due to space constraints) reveal that the stability region of 159 the resulting ERK ethod adapts to the eigenvalues of different discrete operators. 160 To verify the analysis in 7), we consider three different ways of coputing the 161 inverse of the linear Jacobian. The first is by Gauss eliination LU), the second 162 uses GMRES with the full Jacobian, and the third eploys atrix-free GMRES it- 163 erations. In the last approach the GMRES coefficients are restricted by 7b) such 164 as to obtain a second order tie discretization ethod. Figure 1c shows the work- 165 precision diagra for these approaches. The Gaussian eliination and traditional 166 GMRES solutions display first order converge, while the constrained GMRES solu- 167 tion displays second order convergence. 168 Gauss GMRES Mat.Free 5 Conclusions 169 Iplicit tie integration ethods are becoing widely used in the the siulation of 170 tie dependent PDEs, as they do not suffer fro CFL stability restrictions. While 171 Page 432
7 Stability Analysis of the Matrix-Free Linearly Iplicit Euler Method iplicit ethods can use uch larger tie steps than explicit ethods, their co- 172 putational cost per step is also higher. The coputational tie is doinated by the 173 solutions of non)linear systes of equations that define each stage of a linearly) 174 iplicit ethod. The iplicit code is ore effective only when the gains in step size 175 offset the extra cost. 176 To reduce the coputational overhead of LU decoposition, to alleviate eory 177 requireents, and to aid parallelization, iterative Krylov space ethods are used 178 to solve the large linear systes. A atrix-free ipleentation approxiates the 179 required Jacobian vector products by finite differences. 180 This paper studies the effect of the atrix-free iterative solutions on the proper- 181 ties of the nuerical integration ethod. The analysis reveals that atrix-free lin- 182 early iplicit ethods can be viewed as explicit Runge Kutta ethods. Their stabil- 183 ity region is finite, and the unconditional stability property of the original iplicit 184 ethod is lost. The equivalent Runge Kutta ethod is nonlinear, in the sense that 185 its weights depend on the tie step and on the stage vectors. This akes the ac- 186 curacy analysis difficult. Order conditions of the equivalent explicit Runge Kutta 187 ethod can be fulfilled by iposing additional conditions on the Krylov solution 188 coefficients. For preconditioned atrix-free iterations the overall tie stepping pro- 189 cess is equivalent to a Rosenbrock-W ethod, where the preconditioner deterines 190 the Jacobian approxiation. Future work will address the effect of a finite nuber 191 of Krylov iterations on the stability and accuracy of the overall schee, in the case 192 where an analytical Jacobian is used. 193 Acknowledgents This work was supported by the NSF projects DMS , CCF , OCI , PetaApps A. Sandu was also supported by NCAR s ASP 195 Faculty Fellowship Progra. 196 Bibliography 197 [1] A. Crivellini and F. Bassi. An iplicit atrix-free Discontinuous Galerkin solver 198 for viscous and turbulent aerodynaic siulations. Coputers and Fluids, ):81 93, [2] E. Hairer and G. Wanner. Solving Ordinary Differential Equations II. Stiff and 201 Differential-Algebraic Probles. 2Ed. Springer-Verlag, [3] E. Hairer, S.P. Norsett, and G. Wanner. Solving Ordinary Differential Equations 203 I. Nonstiff Probles. Springer-Verlag, Berlin, [4] Ti Kelley. Iterative Methods for Optiization. SIAM, [5] D.A. Knoll and D.E. Keyes. Jacobian-free Newton-Krylov ethods: a survey 206 of approaches and applications. Journal of Coputational Physics, 1932): , [6] H. Luo, J.D. Bau, and R. Lohner. A fast, atrix-free iplicit ethod for co- 209 pressible flows on unstructured grids. In Sixteenth international conference on 210 nuerical ethods in fluid dynaics, volue 515 of Lecture Notes in Physics, 211 pages 73 78, Page 433
8 Adrian Sandu and Aik St-Cyr [7] Y. Saad. Iterative Methods for Sparse Linear Systes. SIAM, 2 edition, [8] K. Shahbazi. An explicit expression for the penalty paraeter of the interior 214 penalty ethod. Journal of Coputational Physics, 205: , Page 434
Comparison 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 informationlecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 38: Linear Multistep Methods: Absolute Stability, Part II
lecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 3: Linear Multistep Methods: Absolute Stability, Part II 5.7 Linear ultistep ethods: absolute stability At this point, it ay well
More informationNumerical issues in the implementation of high order polynomial multidomain penalty spectral Galerkin methods for hyperbolic conservation laws
Nuerical issues in the ipleentation of high order polynoial ultidoain penalty spectral Galerkin ethods for hyperbolic conservation laws Sigal Gottlieb 1 and Jae-Hun Jung 1, 1 Departent of Matheatics, University
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 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 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 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 informationEMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS
EMPIRICAL COMPLEXITY ANALYSIS OF A MILP-APPROACH FOR OPTIMIZATION OF HYBRID SYSTEMS Jochen Till, Sebastian Engell, Sebastian Panek, and Olaf Stursberg Process Control Lab (CT-AST), University of Dortund,
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 informationlecture 36: Linear Multistep Mehods: Zero Stability
95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,
More informationDonald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles
More informationPOD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION EQUATION IN NEURO-MUSCULAR SYSTEM
6th European Conference on Coputational Mechanics (ECCM 6) 7th European Conference on Coputational Fluid Dynaics (ECFD 7) 1115 June 2018, Glasgow, UK POD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION
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 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 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 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 informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
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 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 informationP016 Toward Gauss-Newton and Exact Newton Optimization for Full Waveform Inversion
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
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 informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
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 informationFinding Rightmost Eigenvalues of Large Sparse. Non-symmetric Parameterized Eigenvalue Problems. Abstract. Introduction
Finding Rightost Eigenvalues of Large Sparse Non-syetric Paraeterized Eigenvalue Probles Applied Matheatics and Scientific Coputation Progra Departent of Matheatics University of Maryland, College Par,
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationRECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE
Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS
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 informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationNumerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term
Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2
More informationPHY307F/407F - Computational Physics Background Material for Expt. 3 - Heat Equation David Harrison
INTRODUCTION PHY37F/47F - Coputational Physics Background Material for Expt 3 - Heat Equation David Harrison In the Pendulu Experient, we studied the Runge-Kutta algorith for solving ordinary differential
More informationStability Ordinates of Adams Predictor-Corrector Methods
BIT anuscript No. will be inserted by the editor Stability Ordinates of Adas Predictor-Corrector Methods Michelle L. Ghrist Jonah A. Reeger Bengt Fornberg Received: date / Accepted: date Abstract How far
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
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 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 information. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe
PROPERTIES OF MULTIVARIATE HOMOGENEOUS ORTHOGONAL POLYNOMIALS Brahi Benouahane y Annie Cuyt? Keywords Abstract It is well-known that the denoinators of Pade approxiants can be considered as orthogonal
More informationOn the Communication Complexity of Lipschitzian Optimization for the Coordinated Model of Computation
journal of coplexity 6, 459473 (2000) doi:0.006jco.2000.0544, available online at http:www.idealibrary.co on On the Counication Coplexity of Lipschitzian Optiization for the Coordinated Model of Coputation
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 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 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 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 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 informationSD Numerical Simulation Technique for Hydrodynamic Flow Gas- Solids Mixing Mantilla Núñez, Irla* 1 S. De Vicente C. 2
SD Nuerical Siulation Technique for Hydrodynaic Flow Gas- Solids Mixing Mantilla Núñez, Irla* S. De Vicente C. National University of Engineering, Lia, Perú, Polytechnic University of Madrid. *Lia 5, irlan@uni.edu.pe.
More informationAN APPLICATION OF CUBIC B-SPLINE FINITE ELEMENT METHOD FOR THE BURGERS EQUATION
Aksan, E..: An Applıcatıon of Cubıc B-Splıne Fınıte Eleent Method for... THERMAL SCIECE: Year 8, Vol., Suppl., pp. S95-S S95 A APPLICATIO OF CBIC B-SPLIE FIITE ELEMET METHOD FOR THE BRGERS EQATIO by Eine
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 informationUfuk Demirci* and Feza Kerestecioglu**
1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,
More informationDESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *
IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS
More informationA BLOCK MONOTONE DOMAIN DECOMPOSITION ALGORITHM FOR A NONLINEAR SINGULARLY PERTURBED PARABOLIC PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volue 3, Nuber 2, Pages 211 231 c 2006 Institute for Scientific Coputing and Inforation A BLOCK MONOTONE DOMAIN DECOMPOSITION ALGORITHM FOR A NONLINEAR
More informationComparison of Charged Particle Tracking Methods for Non-Uniform Magnetic Fields. Hann-Shin Mao and Richard E. Wirz
42nd AIAA Plasadynaics and Lasers Conferencein conjunction with the8th Internati 27-30 June 20, Honolulu, Hawaii AIAA 20-3739 Coparison of Charged Particle Tracking Methods for Non-Unifor Magnetic
More informationNumerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda
Journal of Math-for-Industry, Vol. 3 (C-), pp. Nuerically repeated support splitting and erging phenoena in a porous edia equation with strong absorption To the eory of y friend Professor Nakaki. Kenji
More informationON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN
More informationROBUST REPRESENTATION OF DRY CELLS IN MODFLOW
Final Technical Report for ROBUST REPRESENTATION OF DRY CELLS IN MODFLOW SwRI Project 20-13003 Prepared for Edwards Aquifer Authority 1615 N. St. Mary's Street San Antonio, TX 78215 May 2007 SOUTHWEST
More informationQuantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search
Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths
More informationProjectile Motion with Air Resistance (Numerical Modeling, Euler s Method)
Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,
More informationHybrid System Identification: An SDP Approach
49th IEEE Conference on Decision and Control Deceber 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Hybrid Syste Identification: An SDP Approach C Feng, C M Lagoa, N Ozay and M Sznaier Abstract The
More informationarxiv: v1 [physics.comp-ph] 20 Oct 2017 Abstract
Convergence Properties & Extension to Non-Conservative Multidiensional ystes Haran Jackson Cavendish Laboratory, JJ Thoson Ave, Cabridge, UK, CB3 0HE arxiv:70.07603v [physics.cop-ph] 20 Oct 207 Abstract
More informationWAVE-BODY INTERACTIONS FOR A SURFACE-PIERCING BODY IN WATER OF FINITE DEPTH *
745 010,(6):745-75 DOI: 10.1016/S1001-6058(09)6011-8 WAVE-BODY INTERACTIONS FOR A SURFACE-PIERCING BODY IN WATER OF FINITE DEPTH * LI Yong, LIN Mian Institute of Mechanics, Chinese Acadey of Sciences,
More informationModel Predictive Control Approach to Design Practical Adaptive Cruise Control for Traffic Jam
Taku Takahaa et al./international Journal of Autootive Engineering Vol.9, No.3 (218) pp.99-14 Research Paper 218495 Model Predictive Control Approach to Design Practical Adaptive Cruise Control for Traffic
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 informationModeling Chemical Reactions with Single Reactant Specie
Modeling Cheical Reactions with Single Reactant Specie Abhyudai Singh and João edro Hespanha Abstract A procedure for constructing approxiate stochastic odels for cheical reactions involving a single reactant
More informationANALYSIS OF A NUMERICAL SOLVER FOR RADIATIVE TRANSPORT EQUATION
ANALYSIS OF A NUMERICAL SOLVER FOR RADIATIVE TRANSPORT EQUATION HAO GAO AND HONGKAI ZHAO Abstract. We analyze a nuerical algorith for solving radiative transport equation with vacuu or reflection boundary
More informationShort Papers. Test Data Compression and Decompression Based on Internal Scan Chains and Golomb Coding
IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 1, NO. 6, JUNE 00 715 Short Papers Test Data Copression and Decopression Based on Internal Scan Chains and Golob Coding
More informationAn Improved Particle Filter with Applications in Ballistic Target Tracking
Sensors & ransducers Vol. 72 Issue 6 June 204 pp. 96-20 Sensors & ransducers 204 by IFSA Publishing S. L. http://www.sensorsportal.co An Iproved Particle Filter with Applications in Ballistic arget racing
More informationANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE
DRAFT Proceedings of the ASME 014 International Mechanical Engineering Congress & Exposition IMECE014 Noveber 14-0, 014, Montreal, Quebec, Canada IMECE014-36371 ANALYTICAL INVESTIGATION AND PARAMETRIC
More informationNonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy
Storage Capacity and Dynaics of Nononotonic Networks Bruno Crespi a and Ignazio Lazzizzera b a. IRST, I-38050 Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I-38050 Povo (Trento) Italy INFN Gruppo
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationScalable Symbolic Model Order Reduction
Scalable Sybolic Model Order Reduction Yiyu Shi Lei He -J Richard Shi Electrical Engineering Dept, ULA Electrical Engineering Dept, UW Los Angeles, alifornia, 924 Seattle, WA, 985 {yshi, lhe}eeuclaedu
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 informationSpine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area
Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular
More informationAnalysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning
Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics
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 informationTHE POLYNOMIAL REPRESENTATION OF THE TYPE A n 1 RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n
THE POLYNOMIAL REPRESENTATION OF THE TYPE A n RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n SHEELA DEVADAS AND YI SUN Abstract. We study the polynoial representation of the rational Cherednik algebra
More informationSupplementary Material for Fast and Provable Algorithms for Spectrally Sparse Signal Reconstruction via Low-Rank Hankel Matrix Completion
Suppleentary Material for Fast and Provable Algoriths for Spectrally Sparse Signal Reconstruction via Low-Ran Hanel Matrix Copletion Jian-Feng Cai Tianing Wang Ke Wei March 1, 017 Abstract We establish
More informationEffective joint probabilistic data association using maximum a posteriori estimates of target states
Effective joint probabilistic data association using axiu a posteriori estiates of target states 1 Viji Paul Panakkal, 2 Rajbabu Velurugan 1 Central Research Laboratory, Bharat Electronics Ltd., Bangalore,
More informationConvex Programming for Scheduling Unrelated Parallel Machines
Convex Prograing for Scheduling Unrelated Parallel Machines Yossi Azar Air Epstein Abstract We consider the classical proble of scheduling parallel unrelated achines. Each job is to be processed by exactly
More informationOptimization of ripple filter for pencil beam scanning
Nuclear Science and Techniques 24 (2013) 060404 Optiization of ripple filter for pencil bea scanning YANG Zhaoxia ZHANG Manzhou LI Deing * Shanghai Institute of Applied Physics, Chinese Acadey of Sciences,
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 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 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 informationEfficient Numerical Solution of Diffusion Convection Problem of Chemical Engineering
Cheical and Process Engineering Research ISSN 4-7467 (Paper) ISSN 5-9 (Online) Vol, 5 Efficient Nuerical Solution of Diffusion Convection Proble of Cheical Engineering Bharti Gupta VK Kukreja * Departent
More informationProbability Distributions
Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples
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 informationHyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005
Hyperbolic Horn Helical Mass Spectroeter (3HMS) Jaes G Hageran Hageran Technology LLC & Pacific Environental Technologies April 5 ABSTRACT This paper describes a new type of ass filter based on the REFIMS
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 informationFurther applications of scheme for reducing numerical viscosity: 3D hypersonic flow
Journal of Physics: Conference Series Further applications of schee for reducing nuerical viscosity: 3D hypersonic flow To cite this article: S A Elaskar et al 2009 J. Phys.: Conf. Ser. 166 012018 View
More informationPaul M. Goggans Department of Electrical Engineering, University of Mississippi, Anderson Hall, University, Mississippi 38677
Evaluation of decay ties in coupled spaces: Bayesian decay odel selection a),b) Ning Xiang c) National Center for Physical Acoustics and Departent of Electrical Engineering, University of Mississippi,
More informationESTIMATE OF THE TRUNCATION ERROR OF FINITE VOLUME DISCRETISATION OF THE NAVIER-STOKES EQUATIONS ON COLOCATED GRIDS.
ESTIMATE OF THE TRUNCATION ERROR OF FINITE VOLUME DISCRETISATION OF THE NAVIER-STOKES EQUATIONS ON COLOCATED GRIDS. Alexandros Syrakos, Apostolos Goulas Departent of Mechanical Engineering, Aristotle University
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 informationThe Simplex Method is Strongly Polynomial for the Markov Decision Problem with a Fixed Discount Rate
The Siplex Method is Strongly Polynoial for the Markov Decision Proble with a Fixed Discount Rate Yinyu Ye April 20, 2010 Abstract In this note we prove that the classic siplex ethod with the ost-negativereduced-cost
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 informationChaotic Coupled Map Lattices
Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each
More informationKinematics and dynamics, a computational approach
Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for
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 informationCHAPTER 19: Single-Loop IMC Control
When I coplete this chapter, I want to be able to do the following. Recognize that other feedback algoriths are possible Understand the IMC structure and how it provides the essential control features
More informationAn Adaptive UKF Algorithm for the State and Parameter Estimations of a Mobile Robot
Vol. 34, No. 1 ACTA AUTOMATICA SINICA January, 2008 An Adaptive UKF Algorith for the State and Paraeter Estiations of a Mobile Robot SONG Qi 1, 2 HAN Jian-Da 1 Abstract For iproving the estiation accuracy
More informationThe Solution of One-Phase Inverse Stefan Problem. by Homotopy Analysis Method
Applied Matheatical Sciences, Vol. 8, 214, no. 53, 2635-2644 HIKARI Ltd, www.-hikari.co http://dx.doi.org/1.12988/as.214.43152 The Solution of One-Phase Inverse Stefan Proble by Hootopy Analysis Method
More informationFriction Induced Hunting Limit Cycles: An Event Mapping Approach
Friction Induced Hunting Liit Cycles: An Event Mapping Approach Ron H.A. Hensen, Marinus (René) J.G. van de Molengraft Control Systes Technology Group, Departent of Mechanical Engineering, Eindhoven University
More informationMatrix Inversion-Less Signal Detection Using SOR Method for Uplink Large-Scale MIMO Systems
Matrix Inversion-Less Signal Detection Using SOR Method for Uplin Large-Scale MIMO Systes Xinyu Gao, Linglong Dai, Yuting Hu, Zhongxu Wang, and Zhaocheng Wang Tsinghua National Laboratory for Inforation
More informationHaarOrthogonalFunctionsBasedParameterIdentification ofwastewatertreatmentprocesswithdistributedparameters
M. TODOROVA and T. PENCHEVA, Haar Orthogonal Functions Based, Che. Bioche. Eng. Q. 19 (3) 275 282 (2005) 275 HaarOrthogonalFunctionsBasedParaeterIdentification ofwastewatertreatentprocesswithdistributedparaeters
More informationarxiv: v1 [math.na] 10 Oct 2016
GREEDY GAUSS-NEWTON ALGORITHM FOR FINDING SPARSE SOLUTIONS TO NONLINEAR UNDERDETERMINED SYSTEMS OF EQUATIONS MÅRTEN GULLIKSSON AND ANNA OLEYNIK arxiv:6.395v [ath.na] Oct 26 Abstract. We consider the proble
More informationMolecular dynamics algorithm for multiple time scales: Systems with long range forces
Molecular dynaics algorith for ultiple tie scales: Systes with long range forces Mark E. Tuckeran@ Bruce J. Berne Departent of Cheistry, Colubia University, New York, New York 10027 Glenn J. Martyna Departent
More information