EFFICIENT DOMAIN DECOMPOSITION METHOD FOR ACOUSTIC SCATTERING IN MULTI-LAYERED MEDIA
|
|
- Jack Todd
- 6 years ago
- Views:
Transcription
1 European Conference on Computatonal Flud Dynamcs ECCOMAS CFD 2006 P. Wesselng, E. Oñate and J. Péraux (Eds) c TU Delft, The Netherlands, 2006 EFFICIENT DOMAIN DECOMPOSITION METHOD FOR ACOUSTIC SCATTERING IN MULTI-LAYERED MEDIA Kazufum Ito and Jar Tovanen Center for Research n Scentfc Computng, Box 8205 North Carolna State Unversty Ralegh, NC , USA e-mal: kto@ncsu.edu, jatovan@ncsu.edu Key words: Acoustc Scatterng, Geologcal survey, Helmholtz equaton, Doman Decomposton Method, Precondtoned Iteratve Method Abstract. An effcent doman decomposton method s developed for dscretzed Helmholtz equatons resultng from acoustcal geologcal survey problems. A Schwarz-type multplcatve precondtoner s constructed usng a fast drect method n subdomans by embeddng them to rectangular domans. The GMRES teratons are reduced to a subspace correspondng to nterfaces. Numercal experments wth several two-dmensonal problems demonstrate the effectveness and scalablty of the proposed method. 1 INTRODUCTION Acoustcal geologcal survey problems are often modeled usng the Helmholtz equaton wth varyng wave number. The dscretzaton of ths equatons requres suffcently many nodes per wavelength 7 and due to ths the resultng system of lnear equatons s very large. It s usual that grds have from 500 to a few thousands grd ponts n each drecton. Thus, the resultng problem can have several mllon unknowns for a two-dmensonal problem and several bllon unknowns for a three-dmensonal problem. The system matrx s ndefnte and complex-valued. Thus, many effcent soluton procedures avalable for postve defnte problems are not applcable. Iteratve methods wthout a good precondtoner converge slowly for Helmholtz problems whch leads to too hgh computatonal cost. Furthermore, whle drect methods lke LU decomposton mght be applcable for two-dmensonal problems they are computatonally too expensve for three-dmensonal problems. Thus, the research on effcent solvers for the Helmholtz equaton wth varyng wave number has concentrated on developng good precondtoners. The man technques employed to construct these precondtoners are doman decomposton methods, multgrd methods, and fast drect solvers. Plessx and Mulder 13 studed tensor product form precondtoners based on fast drect solvers. Whle for low frequency problems they were very effectve ther performance 1
2 deterorated severely when the frequency approached practcally nterestng level. The authors 6,8 used a smlar approach wth good success for almost perfectly vertcally layered meda. Larson and Holmgren 12 developed a doman decomposton precondtoner employng a fast drect solver and tested t on a parallel computer. Erlangga, Oosterlee, and Vuk 3 developed a multgrd precondtoner whch s farly effectve for geologcal survey problems. In ths paper, we develop a Schwarz-type multplcatve doman decomposton precondtoner. The subdoman precondtoners employ a fast drect solver va embeddng the subdomans nto rectangles. By fast drect methods, we mean drect methods based on FFT or cyclc reducton, for example, whch requre order of m log m or m(log m) 2 operatons to solve a system wth m unknowns. Wth a pecewse constant wave number, the precondtoner can be made to match the system matrx n the nteror of the subdomans. Ths has two mportant consequences: The resdual vanshed n the nteror of the subdomans and due to ths the teratons can be reduced on the nterfaces and ther near-by grd ponts. Thus, the soluton procedure can be consdered as a precondtoned teratve method on a sparse subspace 9,10. Due to ths the memory requrement of the GMRES method 15 s greatly reduced. Secondly, the condtonng of the precondtoned system s mproved due to the matchng n the nteror of the subdomans. The numercal examples demonstrate the effectveness of the developed approach. 2 MODEL PROBLEMS In the rectangular computatonal doman Π, the pressure feld p satsfes the Helmholtz partal dfferental equaton p k 2 p = g, (1) where k s the wave number gven by k = ω, ω s the angular frequency, c s the speed of c sound, and g corresponds to a sound source. We assume the speed of sound and, thus, also the wave number to be pecewse constant functons. On the boundares Π, we mpose a second-order absorbng boundary condton 1,8 p n = kp p on Π 2 s k s (2) p n = 3k 4 p at C, where n and s denote the unt outward normal and tangent vectors, respectvely, and C denotes the set of the corner ponts of Π. 3 NUMERICAL METHOD 3.1 Dscretzaton We dscretze the Helmholtz equaton on a unform grd wth the grd step sze h. The dscretzaton stencls for nteror grd ponts, boundary grd ponts, and corner grd ponts 2
3 1 4 (k(x + h/2, y + h/2) h) (3) Kazufum Ito and Jar Tovanen are gven by Fgure 1, where the dagonal weghts are defned by d = 4 1 (k(x ± h/2, y ± h/2) h) 2, 4 b = ( ) k(x, y ± h/2) h 1 (k(x + h/2, y ± h/2) h) 2, 2 k(x, y ± h/2) h 4 a = ( ) k(x, y + h/2) h + k(x + h/2, y) h 2 k(x, y + h/2) h k(x + h/2, y) h (1 + k(x,y+h/2) h ) 1 2 (1 + k(x,y+h/2) h ) 1 d 1 b 1 a 1 2 (1 + k(x+h/2,y) h ) (1 + k(x,y h/2) h ) Fgure 1: The dscretzaton stencl at (x, y) for an nteror grd pont (left), a grd pont on the left boundary (mddle), and the left lower corner pont (rght). The weghts d, b, and a are gven by (3). The dscretzaton leads to a system of lnear equatons 3.2 Iteratve soluton procedure Ap = b. (4) Instead of solvng the system of lnear equatons (4), we solve teratvely usng the GMRES method 15 a left precondtoned system AB 1 q = b, (5) where B s the precondtoner. Once we have obtaned q the soluton of the orgnal problem s gven by p = B 1 q. In order to descrbe our doman decomposton precondtoner B, we defne each subdoman as the unon of cells n whch the speed of sound s a gven constant. More precsely, let c j, j = 1,..., n, denote the dfferent speeds of sound. The closure of the jth subdoman Ω j s defned by Ω j = [x h/2, x + h/2] [y h/2, y + h/2], (6) c(x, y) = c j (x, y) M 3
4 where M s the set of cell mdponts. Let the rectangular matrx R j correspond to the restrcton operator nto the subdoman Ω j. Multplyng a vector by t results a vector contanng the components assocated to the subdoman Ω j and ts boundary. We denote the precondtoner n Ω j by B j and t wll be defned n Secton 3.3. We employ a Schwarz-type multplcatve doman decomposton precondtoner B = P n whch s defned recursvely as P 1 j = P 1 j 1 + RT j B 1 j R j (I AP 1 j 1 ), j = 2,..., n, (7) where P 1 1 = R T 1 B 1 1 R 1. Thus, we obtan y = y (n) = B 1 x by performng the followng sequence of operatons: 3.3 Subdoman precondtoners y (1) = R T 1 B 1 1 R 1 x, y (2) = y (1) + R T 2 B 1 2 R 2 (x Ay (1) ),. y (n) = y (n 1) + R T n B 1 n R n (x Ay (n 1) ). Systems of lnear equatons correspondng to general shaped subdomans cannot be solved usng a fast drect solver. For ths reason, we embed each subdoman nto a larger rectangular doman and we defned the subdoman precondtoner as a Schur complement matrx. In order to descrbe ths more precsely, we defne the mnmum and maxmum values of the coordnates for each subdoman Ω j as x j mn = argmn (x,y) Ω j x, x j max = argmax (x,y) Ω j x, y j mn = argmn (x,y) Ω j y, and y j max = argmax (x,y) Ω j y. Smlarly, we defne x mn, x max, y mn, and y max for the computatonal doman Π. Then the extended rectangular subdomans are defned by Ω j = [ max{x j mn lh, x mn}, mn{x j max + lh, x max } ] [ max{y j mn lh, y mn}, mn{y j max + lh, y max } ], where l s a nonnegatve nteger. The jth subdoman precondtoner s based on the Helmholtz equaton (8) (9) (10) p k 2 j p = g n Ω j, (11) where k j = ω/c j and the absorbng boundary condtons (2) are posed on the boundares Ω j. By usng the same grd and dscretzaton for the problem (11), we obtan a matrx C j whch has a block form ( ) Cj,dd C C j = j,de (12) C j,ed 4 C j,ee
5 once the unknowns correspondng to the subdoman Ω j are numbered frst (subscrpt d) and the unknowns correspondng to the extenson Ω j \ Ω j are numbered second (subscrpt e). The precondtoner B j for the jth subdoman s gven by a Schur complement matrx A system of lnear equatons B j y j = x j can solved as B j = C j,dd C j,de C 1 j,ee C j,ed. (13) y j = ( 1 0 ) C 1 j ( ) 1 x 0 j (14) usng a fast drect solver lke the cyclc reducton method 5,14. Ths soluton requres order of m j log m j floatng pont operatons, where m j s the sze of C j. 3.4 Sparse subspace teratons The GMRES method 15 forms and stores bass vectors for a Krylov subspace whch after l teratons s K l (v) = { v, AB 1 v, (AB 1 ) 2 v,..., (AB 1 ) l v }, (15) where v s the ntal resdual. For example, f the ntal guess for soluton n the GMRES method s zero then v = b. In general case, storng the bass vector requres a vast amount of memory for large problems unless the teratons converge already after a few teratons. The vectors v j, j = 1,..., l, defnng K l (v) are gven by v j = AB 1 v j 1 = [(A B) + B] B 1 v j 1 = (A B)B 1 v j 1 + v j 1, (16) where we have denoted v 0 = v. It s easy to see from (16) usng nducton that v j+1 X {v}, where X = range(a B). (17) The subspace X corresponds to the rows of the matrces A and B whch dffer. Due to the constructon of the precondtoner, these rows are assocated to the nterfaces between subdomans and ther near-by grd ponts 4,9,10. Thus, the subspace X s very sparse. For a two-dmensonal problem, the number of nonzero components n vectors on X s of order m, where m s the sze of the system of lnear equatons. If the ntal guess s chosen to be b then the ntal resdual v belongs to X and, thus, K l (v) X, that s, the Krylov subspaces are sparse. Alternatvely, X {v} s sparse when b s sparse and the ntal guess s zero. These observatons reduce the memory usage by a large factor and allows us to use the GMRES method wthout restarts. Furthermore, the so-called partal soluton technque 2,11 can be used to solve the systems of lnear equatons wth B on the sparse subspace 4,9. 5
6 4 NUMERICAL RESULTS In all experments, the stoppng crteron for the GMRES method was that the norm of the resdual s reduced by the factor In all embeddngs, we extended the subdomans by one mesh step, that s, l = 1 n (10). The sound source s always a pont source at the mdpont of the upper boundary. The ntal guess for the soluton n the GMRES method s zero for all problems. All experments are performed wthout GMRES restarts. 4.1 Perfectly vertcally layered medum These test problems are defned n the unt square and we use a grd to solve these problems. In the frst test problem, the speed of sound s constant c = 1 n the whole doman. However, we dvde the unt square nto two and four equally szed subdomans whch are extended accordng to (10) for precondtonng. The number of GMRES teratons for the resultng precondtoned systems are gven by the thrd and fourth column of Table 1. The next two problems have perfectly vertcally layered meda shown n Fgure 2. The number of GMRES teratons for these problems are gven by the two last columns of Table 1. The real part of the soluton at the frequency f = 25 s shown by Fgure 2. c = 1 c = 1 c = 2 c = 2 c = 3 c = 4 Fgure 2: Layered test problems n the unt square and the soluton for the four layer problem at f = 25. homog. homog. layered layered subdomans f grd ter. ter. ter. ter Table 1: The number of teratons for dfferent test problems and frequences f. 6
7 4.2 Wedge problem The speed of sound for the wedge test problem 3,13 s shown n Fgure 3. The doman decomposton s gven by Fgure 4 together wth the one grd step extensons for the subdomans. Table 2 gves the number of GMRES teratons and CPU tmes n seconds at dfferent frequences on a PC wth 3.8 GHz Xeon processor. The real parts of the solutons at f = 30 Hz and 50 Hz are plotted n Fgure m c = 2000 m/s 500 m c = 1500 m/s 400 m 100 m m m c = 3000 m/s 600 m Fgure 3: The wedge test problem and the solutons at f = 30 Hz and 50 Hz. subdoman Ω 1 subdoman Ω 2 subdoman Ω 3 Fgure 4: The subdomans Ω j for the wedge test problem (red and black), ther extensons Ω j \ Ω j (green), and the assocated parts of the sparse subspace X {v} (black). 4.3 Salt depost problem Ths test problem mmcs a smple salt depost. The length of the sde of the square computatonal doman s chosen to be 600 meters. The profle for the speed of sound s shown n Fgure 5. It s perfectly vertcally layered except for the salt dsk wth the radus 7
8 f [Hz] grd ter. tme Table 2: The number of teratons and CPU tmes n seconds for dfferent frequences f for the wedge problem. of 150 meters and the center at (300 m, 250 m) when the lower left corner s at orgn. We consder three frequences f = 30 Hz, 60 Hz, and 120 Hz and three grds , , and Table 3 gves the number of GMRES teratons for all possble combnatons of frequences and grds. Fgure 6 shows the real parts of the solutons at the consdered frequences. The doman decomposton s shown n Fgure m/s 2500 m/s 3500 m/s 4500 m/s 5500 m/s Fgure 5: The speed of sound for the salt depost problem. Fgure 6: The solutons for the salt depost problem at f = 30 Hz, 60 Hz, and 120 Hz. 8
9 subdoman Ω 1 subdoman Ω 2 subdoman Ω 3 subdoman Ω 4 subdoman Ω 5 Fgure 7: The subdomans Ω j for the salt depost problem (red and black), ther extensons Ω j \ Ω j (green), and the assocated parts of the sparse subspace X {v} (black). grd f [Hz] Table 3: The number of teratons for dfferent frequences f and grds for the salt depost problem. 5 CONCLUSIONS We proposed a Schwarz-type multplcatve doman decomposton precondtoner for the Helmholtz equaton wth a pecewse constant wave number. The computatonal effcency results from the use of fast drect solvers for subdoman precondtonng and the good condtonng of the precondtoned system. Furthermore, teratons can be reduced on the subdoman boundares and ther neghborng grd ponts. The number of GMRES teratons requred to reduce the resdual by a gven factor depends only mldly on the grd step sze and frequency. The numercal experments demonstrated the method to be effcent for two-dmensonal problems wth realstc frequences and jumps of the speed of sounds encountered n acoustcal geologcal survey problems. The proposed method can be generalzed n a straghtforward manner for three-dmensonal problems. ACKNOWLEDGEMENTS The research was supported by the Academy of Fnland grant # and the Offce of Naval Research grant N REFERENCES [1] A. Bamberger, P. Joly, and J E. Roberts, Second-order absorbng boundary condtons for the wave equaton: a soluton for the corner problem. SIAM J. Numer. Anal., 27, , (1990). [2] A. Banegas, Fast Posson solvers for problems wth sparsty, Math. Comp., 32, , (1978). 9
10 [3] Y. A. Erlangga, C. W. Oosterlee, and C. Vuk, A novel multgrd based precondtoner for heterogeneous Helmholtz problems. SIAM J. Sc. Comput., 27, , (2006). [4] E. Hekkola, Y. A. Kuznetsov, P. Nettaanmäk, and J. Tovanen, Fcttous doman methods for the numercal soluton of two-dmensonal scatterng problems. J. Comput. Phys., 145, , (1998). [5] E. Hekkola, T. Ross, and J. Tovanen, Fast drect soluton of the Helmholtz equaton wth a perfectly matched layer/an absorbng boundary condton. Internat. J. Numer. Methods Engrg., 57, , (2003). [6] Q. Huynh, K. Ito, and J. Tovanen, A fast Helmholtz solver for scatterng by a soundsoft target n sedment. Proceedngs of the 16th Internatonal Conference on Doman Decomposton Methods, To appear. [7] F. Ihlenburg, Fnte element analyss of acoustc scatterng. Sprnger-Verlag, New York, (1998). [8] K. Ito and J. Tovanen, A fast teratve solver for scatterng by elastc objects n layered meda, Appl. Numer. Math. To appear. [9] K. Ito and J. Tovanen, Precondtoned teratve methods on sparse subspaces, Appl. Math. Letters. To appear. [10] Y. A. Kuznetsov, Numercal methods n subspaces. Vychsltel nye Processy Sstemy II, Nauka, Moscow, , (1985). In Russan. [11] Y. A. Kuznetsov and A. M. Matsokn, On partal soluton of systems of lnear algebrac equatons, Sov. J. Numer. Anal. Math. Modellng, 4, , (1989). [12] E. Larsson and S. Holmgren, Parallel soluton of the Helmholtz equaton n a multlayer doman. BIT, 43, , (2003). [13] R. E. Plessx and W. A. Mulder, Separaton-of-varables as a precondtoner for an teratve Helmholtz solver. Appl. Numer. Math., 44, , (2003). [14] T. Ross and J. Tovanen, A parallel fast drect solver for block trdagonal systems wth separable matrces of arbtrary dmenson. SIAM J. Sc. Comput., 20, , (1999). [15] Y. Saad and M. H. Schultz, GMRES: a generalzed mnmal resdual algorthm for solvng nonsymmetrc lnear systems. SIAM J. Sc. Statst. Comput., 7, , (1986). 10
Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationVector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.
Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationPreconditioning techniques in Chebyshev collocation method for elliptic equations
Precondtonng technques n Chebyshev collocaton method for ellptc equatons Zh-We Fang Je Shen Ha-We Sun (n memory of late Professor Benyu Guo Abstract When one approxmates ellptc equatons by the spectral
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationDeriving the X-Z Identity from Auxiliary Space Method
Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationNonlinear Overlapping Domain Decomposition Methods
Nonlnear Overlappng Doman Decomposton Methods Xao-Chuan Ca 1 Department of Computer Scence, Unversty of Colorado at Boulder, Boulder, CO 80309, ca@cs.colorado.edu Summary. We dscuss some overlappng doman
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationPRECONDITIONING TECHNIQUES IN CHEBYSHEV COLLOCATION METHOD FOR ELLIPTIC EQUATIONS
ITERATIOAL JOURAL OF UMERICAL AALYSIS AD MODELIG Volume 15 umber 1-2 Pages 277 287 c 2018 Insttute for Scentfc Computng and Informaton PRECODITIOIG TECHIQUES I CHEBYSHEV COLLOCATIO METHOD FOR ELLIPTIC
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationA New Refinement of Jacobi Method for Solution of Linear System Equations AX=b
Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,
More informationAnnexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances
ec Annexes Ths Annex frst llustrates a cycle-based move n the dynamc-block generaton tabu search. It then dsplays the characterstcs of the nstance sets, followed by detaled results of the parametercalbraton
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More information[7] R.S. Varga, Matrix Iterative Analysis, Prentice-Hall, Englewood Clis, New Jersey, (1962).
[7] R.S. Varga, Matrx Iteratve Analyss, Prentce-Hall, Englewood ls, New Jersey, (962). [8] J. Zhang, Multgrd soluton of the convecton-duson equaton wth large Reynolds number, n Prelmnary Proceedngs of
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationAdditive Schwarz Method for DG Discretization of Anisotropic Elliptic Problems
Addtve Schwarz Method for DG Dscretzaton of Ansotropc Ellptc Problems Maksymlan Dryja 1, Potr Krzyżanowsk 1, and Marcus Sarks 2 1 Introducton In the paper we consder a second order ellptc problem wth dscontnuous
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More information5 The Rational Canonical Form
5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces
More informationLecture 3. Ax x i a i. i i
18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest
More informationSolutions to exam in SF1811 Optimization, Jan 14, 2015
Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
More informationThe Minimum Universal Cost Flow in an Infeasible Flow Network
Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationDeveloping an Improved Shift-and-Invert Arnoldi Method
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Vol. 5, Issue (June 00) pp. 67-80 (Prevously, Vol. 5, No. ) Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) Developng an
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationAcceleration of Gas Bubble-Free Surface Interaction Computation Using Basis Preconditioners
Acceleraton of Gas Bubble-Free Surface Interacton Computaton Usng Bass Precondtoners K. L. Tan, B. C. Khoo and J. K. Whte ABSTRACT The computaton of gas bubble-free surface nteracton entals a tme-steppng
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationYong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )
Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationNP-Completeness : Proofs
NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationPrimer on High-Order Moment Estimators
Prmer on Hgh-Order Moment Estmators Ton M. Whted July 2007 The Errors-n-Varables Model We wll start wth the classcal EIV for one msmeasured regressor. The general case s n Erckson and Whted Econometrc
More informationMin Cut, Fast Cut, Polynomial Identities
Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationA Fast Iterative Solver for Scattering by Elastic Objects in Layered Media
A Fast Iterative Solver for Scattering by Elastic Objects in Layered Media K. Ito and J. Toivanen Center for Research in Scientific Computation, North Carolina State University, Raleigh, North Carolina
More informationThe interface control domain decomposition (ICDD) method for the Stokes problem. (Received: 15 July Accepted: 13 September 2013)
Journal of Coupled Systems Multscale Dynamcs Copyrght 2013 by Amercan Scentfc Publshers All rghts reserved. Prnted n the Unted States of Amerca do:10.1166/jcsmd.2013.1026 J. Coupled Syst. Multscale Dyn.
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationTR A BDDC ALGORITHM FOR PROBLEMS WITH MORTAR DISCRETIZATION. September 4, 2005
TR2005-873 A BDDC ALGORITHM FOR PROBLEMS WITH MORTAR DISCRETIZATION HYEA HYUN KIM, MAKSYMILIAN DRYJA, AND OLOF B WIDLUND September 4, 2005 Abstract A BDDC balancng doman decomposton by constrants algorthm
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationON COMPUTING MAXIMUM/MINIMUM SINGULAR VALUES OF A GENERALIZED TENSOR SUM
Electronc Transactons on Numercal Analyss. Volume 43, pp. 244 254, 215. Copyrght c 215, Kent State Unversty. ISSN 168 9613. ETNA Kent State Unversty ON COMPUTING MAXIMUM/MINIMUM SINGULAR VALUES OF A GENERALIZED
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to
THE INVERSE POWER METHOD (or INVERSE ITERATION) -- applcaton of the Power method to A some fxed constant ρ (whch s called a shft), x λ ρ If the egenpars of A are { ( λ, x ) } ( ), or (more usually) to,
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationCalculation of time complexity (3%)
Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationSecond Order Analysis
Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to
More informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationTwo-level Parallel Computation for Approximate Inverse with AISM Method
Vol 49 No 6 2164 2168 (June 2008) AISM 2 1 2 3 2 AISM SGI Orgn 2400 MPI PE16 PE 13672 Two-level Parallel Computaton for Approxmate Inverse wth AISM Method Lnje Zhang, 1 Kentaro Morya 2 and Takash Nodera
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationOn a class of preconditioners for solving the Helmholtz equation
Appled Numercal athematcs 50 (2004) 409 425 www.elsever.com/locate/apnum On a class of precondtoners for solvng the Helmholtz equaton Y.A. Erlangga,C.Vuk,C.W.Oosterlee Department of Appled athematcal Analyss,
More informationThe FEAST Algorithm for Sparse Symmetric Eigenvalue Problems
The FEAST Algorthm for Sparse Symmetrc Egenvalue Problems Susanne Bradley Department of Computer Scence, Unversty of Brtsh Columba Abstract The FEAST algorthm s a method for computng extreme or nteror
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationMultigrid Methods and Applications in CFD
Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac
More informationDIVIDE AND CONQUER LOW-RANK PRECONDITIONING TECHNIQUES
DIVIDE AND CONQUER LOW-RANK PRECONDITIONING TECHNIQUES RUIPENG LI AND YOUSEF SAAD Abstract. Ths paper presents a precondtonng method based on a recursve multlevel lowrank approxmaton approach. The basc
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationRESTRICTED ADDITIVE SCHWARZ METHOD WITH HARMONIC OVERLAP
RESTRICTED ADDITIVE SCHWARZ METHOD WITH HARMONIC OVERLAP XIAO-CHUAN CAI, MAKSYMILIAN DRYJA, AND MARCUS SARKIS Abstract. In ths paper, we ntroduce a new Schwarz precondtoner and wth a new coarse space.
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationThe lower and upper bounds on Perron root of nonnegative irreducible matrices
Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationNumerical Properties of the LLL Algorithm
Numercal Propertes of the LLL Algorthm Frankln T. Luk a and Sanzheng Qao b a Department of Mathematcs, Hong Kong Baptst Unversty, Kowloon Tong, Hong Kong b Dept. of Computng and Software, McMaster Unv.,
More informationDFT with Planewaves pseudopotential accuracy (LDA, PBE) Fast time to solution 1 step in minutes (not hours!!!) to be useful for MD
LLNL-PRES-673679 Ths work was performed under the auspces of the U.S. Department of Energy by under contract DE-AC52-07NA27344. Lawrence Lvermore Natonal Securty, LLC Sequoa, IBM BGQ, 1,572,864 cores O(N)
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationSupplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso
Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed
More informationAsymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation
Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton
More information