Substructuring Preconditioners for the Bidomain Extracellular Potential Problem
|
|
- Derrick Floyd
- 5 years ago
- Views:
Transcription
1 Substructuring Preconditioners for the Bidomain Extraceuar Potentia Probem Mico Pennacchio 1 and Vaeria Simoncini 2,1 1 IMATI - CNR, via Ferrata, 1, Pavia, Itay mico@imaticnrit 2 Dipartimento di Matematica, Università di Boogna, Piazza di Porta S Donato, 5, Boogna, and CIRSA Ravenna, Itay vaeria@dmuniboit Summary We study the efficient soution of the inear system arising from the discretization by the mortar method of mathematica modes in eectrocardioogy We focus on the bidomain extraceuar potentia probem and on the cass of substructuring preconditioners We verify that the condition number of the preconditioned matrix ony grows poyogarithmicay with the number of degrees of freedom as predicted by the theory and vaidated by numerica tests Moreover, we discuss the roe of the conductivity tensors in buiding the preconditioner 1 Introduction A macroscopic mode accounting for the excitation process in the myocardium is the bidomain mode that yieds the foowing Reaction-Diffusion (R-D) system of equations for the intra-, extraceuar and transmembrane potentia u i, u and v = u i u: find (v(x, t), u(x, t)), x Ω, t [0, T] such that c m t v div M i v + I(v) = div M i u + I app in Ω ]0, T[ div M u = div M i v in Ω ]0, T[ (1) with M i, M e, M = M i +M e conductivity tensors modeing the cardiac fibers, I app an appied current used to initiate the process, c m the surface capacitance of the membrane The function I(v) is the transmembrane ionic current which is assumed for simpicity to depend ony on v and to be a cubic poynomia (see [7]) The genera R-D system can be more compex, incuding additiona ordinary differentia equations that govern the evoution of v These modes are computationay chaenging because of the different space and time scaes invoved; reaistic three-dimensiona simuations with uniform grids yied discrete probems with more than O(10 7 ) unknowns at every time step To improve computationa efficiency, we consider a non conforming non overapping domain decomposition, within the mortar finite eement method This aows us to concentrate the computationa work ony in regions of high
2 2 Mico Pennacchio and Vaeria Simoncini eectrica activity; in addition, the matching of different discretizations on adjacent subdomains are weaky enforced In [8, 9] we compared this technique to the cassica conforming FEM verifying its better performance In this paper, we focus on the probem of the efficient soution of the inear system arising from this discretization and here for simpicity, we concentrate on the probem with the eiptic equation of (1): for each time instant t find u(x, t), soution of: { div M u = div Mi v in Ω n T M u = n T (2) M i v on Γ Such probem is of interest in its own right, as it represents a separate mode for the bidomain extraceuar potentia [8] We consider substructuring preconditioners and we report our numerica experience on soving probem (2) Our experiments confirm the theory depicting poyogarithmic bound for the condition number of the preconditioned matrix Moreover, attention is devoted to tuning the preconditioner so as to take into account the conductivity tensor M in (2) 2 Mortar Method The computationa domain Ω is decomposed as the union of L subdomains Ω 1,, Ω L (see, eg, [10]) We set Γ n = Ω n Ω, S = Γ n and we denote by γ (i), i = 1,, 4 the i-th side of the -th domain so that Ω = 4 i=1 γ(i) Here we consider ony geometricay conforming decomposition, ie each edge γ (i) coincides with Γ n for some n The Mortar Method is appied by choosing a spitting of the skeeton S as the disjoint union of a certain number of subdomain sides γ (i), caed mortar or save sides: we fix an index set I {1,,L} {1,,4} such that S = (,i) I γ(i) The index set corresponding to trace or master sides wi be denoted by I : I {1,, L} {1,,4}, I I = and S = (,i) I γ(i) Let the spaces X and T be X = H1 (Ω ), T = H1/2 ( Ω ) with the broken norms: u 2 X = u 2 1,Ω and η 2 T = η 2 1/2, Ω For each, et aso Vh be a famiy of finite dimensiona subspaces of H1 (Ω ) C 0 ( Ω ), depending on a parameter h = h > 0, X h = L =1 V h X, T h = V h Ω and T h = L =1 T h T Then we define two composite biinear forms a X, a i X : X X R as: a X (u, φ) = φ T M u dx, a i X(u, φ) = φ T M i u dx (3) Ω Ω Since these biinear forms are not coercive on X, we consider proper subspaces of X consisting of functions satisfying a suitabe weak continuity constraint, eading to the foowing constrained approximation and trace spaces
3 Substructuring Preconditioners in Eectrocardioogy 3 X h = {v h X h, [v h ]λds = 0, λ M h } (4) S T h = {η T h, [η]λds = 0, λ M h }, (5) S with M h a suitaby chosen finite dimensiona mutipier space We can write the discrete probem, whose soution existence was proved in [8, Theorem 31]: Probem 1 Find u h X h such that for a φ h X h a X (u h, φ h ) = a i X (v h, φ h ) (6) We remark that Probem 1 admits a soution unique up to an additive constant reated to the reference potentia chosen In this paper we consider as reference potentia the one given by the potentia at a reference point x 0 Ω 3 Substructuring Preconditioners A key aspect of substructuring preconditioners is that they distinguish among three types of degrees of freedom: interior (corresponding to basis functions vanishing on the skeeton and supported on one sub-domain), edge and vertex degrees of freedom [4, 1] Thus, each function u X h can be written as the sum of three suitaby defined components: u = u 0 +u E +u V More specificay, et w = (w ) =1,,L X h be any discrete function, then w = w 0 + R h (w), w 0 X 0 h, (7) with w 0 Xh 0 interior function and R h(w) a discrete ifting, ie R h (w) = (Rh (w )) =1,,K, where Rh (w ) is the unique eement in Vh satisfying R h (w ) = w on Γ and Ω i,j M x i Consequenty, the spaces X h and X h can be spit as: and it can be verified that x j R h (w )v h dx = 0, v h V h (8) X h = X 0 h R h(t h ) X h = X 0 h R h(t h ) a X(w, v) = a X(w 0, v 0 ) + a X(R h (w), R h (v)) = a X(w 0, v 0 ) + s(η(w), η(v)), (9) where the discrete Stekov-Poincaré operator s : T h T h R is defined by s(ξ, η) := (M(x) Rh (ξ)) R h (η) (10) Ω
4 4 Mico Pennacchio and Vaeria Simoncini Furthermore the space of constrained skeeton functions T h can be spit as the sum of vertex and edge functions More specificay, denoting by L L =1 H1/2 ( Ω ) the space L = {(η ) =1,,L, η is inear on each edge of Ω }, then we can define the space of constrained vertex functions as T V h = P h L (11) with P h the correction operator imposing the constraint We make the (not restrictive) assumption L T h, which yieds Th V T h, and we introduce the space of constrained edge functions Th E T h defined by T E h = {η = (η ) =1,,L T h, η (A) = 0, vertex A of Ω } (12) We can easiy verify that T h = Th V T h E Then we wi consider a bock Jacobi type preconditioner s : T h T h R defined as s(η, ξ) = b V (η V, ξ V ) + b E (η E, ξ E ) (13) with bocks reated to the foowing edge and vertex goba biinear forms b E : Th E h E R such that be (η E, η E ) s(η E, η E ) b V : Th V h V such that bv (η V, η V ) s(η V, η V ) (14) 31 Matrix form In this section we derive the matrix form of the discrete Stekov-Poincaré operator s in (10) Equation (6) yieds the foowing inear system of equations: Au = b with b = A i v, (15) where A, A i are the stiffness matrices associated to the discretization of a X, a i X defined in (3) It can be shown that the matrix A is positive semidefinite and the system is consistent We reorder the vector of unknowns as: u = (u 0,u E,u V,u S ) T, with u 0,u E,u V,u S interior, edge, vertex and save nodes, respectivey From the mortar condition, it foows that the interior nodes of the mutipier sides are not associated with genuine degrees of freedom in the FEM space Indeed, the vaue of the coefficients u S corresponding to basis functions iving on save sides is uniquey determined by the remaining coefficients through the jump (mortar) condition and can be eiminated from the goba vector u, ie C S u S = C E u E C V u V u S =: Q E u E + Q V u V (16) where Q E = C 1 S C E, Q V = C 1 S C V The entries of C S, C E, C V are given by c ij = γ m [φ j ]λ i ds, λ i M h with φ j corresponding to the different noda basis functions on the save and master side and associated with the vertices Since biorthogona basis functions are empoyed, the square matrix C S is
5 Substructuring Preconditioners in Eectrocardioogy 5 diagona and easiy invertibe (cf [11]) The reduction in (16) may be written in matrix form as u = Q u I u E with Q = 0 I E I u V (17) V 0 Q E Q V where Q is a goba switching matrix The resuting reduced system is thus given by Ãu M = b (18) with à = QT AQ and b = Q T b We note that the (1,1) bock in à is cheapy invertibe therefore, the Schur compement of the system reative to the (1,1) bock can readiy be obtained, yieding the further reduced system ( ) ( ) ue be S = u V b V The Schur compement S represents the matrix form of the Stekov-Poincaré operator s(, ) To obtain the matrix form of s(, ) we consider the space L of inear functions, used in the spitting of the trace space (11) Then, we introduce an interpoation map denoted by RH T (say piecewise interpoation) from the noda vaue on V (vertices) onto a nodes of S The matrix R H can be viewed as the weighted (( ) restriction ) map from S onto V By defining the IE square matrix J = R O H, with I E the n E n E identity matrix, we can derive the new Schur compement matrix S, after the vertex correction, ) ( SE S = J T SEV SJ = S EV T S V (19) 32 The preconditioner We describe a generaization of a known optima preconditioner for S, and some more computationay effective variants The matrix discretization of the form s yieds the foowing (bock-jacobi type) diagona preconditioner ( ) PE 0 P =, 0 P V where P E, P V are the matrix counterparts of the biinear forms b E and b V in (14), respectivey It can be verified that the preconditioned matrix P 1 S satisfies the theory deveoped in [1, 3] so that ( )) 2 H cond(p 1 S) (1 + og (20) h
6 6 Mico Pennacchio and Vaeria Simoncini with H size of the subdomains and h finest meshsize of the finite eement spaces used Moreover, if an auxiiary coarse mesh is chosen for the vertex bock with mesh size δ > h as studied in [3], then a simiar estimate can be obtained but with a factor ( 1 + og ( )) H 3 h The next three variants make the preconditioner above computationay more appeaing with no essentia oss of optimaity This goa is achieved by repacing either or both the edge and vertex bocks P E and P V with more convenient approximations In their construction, we were inspired by a simiar approach first proposed in [4, 1, 3] for eiptic probems For ater considerations, we reca here an important bound for the condition number of the preconditioned matrix, expressed in terms of the preconditioning quaity of the two diagona bocks More precisey, et P = diag(p 1, P 2 ) be a Jacobitype preconditioner, and et µ M = max{λ max (P1 1 min{λ min (P1 1 S E ), λ min (P2 1 S V )} Then cond(p 1 S) 1 + γ 1 γ S E ), λ max (P 1 2 S V )}, µ m = µ M µ m γ 1, (21) where 1 + γ is the argest eigenvaue of the preconditioned matrix obtained by using the bock diagona of S as preconditioner P (see, eg, [6]) Foowing [1, 4] a simpe approach consists in dropping a coupings between different edges and between edges and vertex points: P E is repaced by its bock diagona part with one bock for each mortar This simpification provides our first variant, ( ) P diag P 1 = E 0 0 P V Assembing the edge and vertex bock preconditioner with such a choice coud be quite expensive A more efficient preconditioner may be obtained by approximating the edge bock P E of P as P (R) E = αr where R is the square root of the stiffness matrix associated on each edge to the discretization of the operator d 2 /dx 2 with homogeneous Dirichet conditions at the extrema [4, 5, 10] The choice of the parameter α is discussed beow Thus our second variant is ( ) (R) P P 2 = E 0 0 P V It can be easiy verified that (cf, eg, [4, 5]) c 1 v T SE v v T Rv c 2 (1 + og ( )) 2 H v T SE v (22) h
7 Substructuring Preconditioners in Eectrocardioogy 7 where c 1, c 2 are independent of H, h but may depend on the coefficients of M(x) Since P V = S V, in (21) we obtain µ M = max{λ max ((αr) 1 SE ), 1} ( ( max{α 1 c og H )) 2 h, 1} Anaogousy, µm min{α 1 c 1, 1} Therefore, the determination of µ M, µ m is infuenced by the magnitude of c 1, c 2 and of α The anisotropic conductivity tensor M = M i +M e is given as M s = M s (x) = σt si + (σs σs t )aat, s = i, e, where a = a(x) is the unit vector tangent to the cardiac fiber at a point x Ω, I is the identity matrix and σ s, σs t for s = i, e are the conductivity coefficients aong and across fiber, in the (i) and (e) media, assumed constant with σ s > σt s > 0 As aready mentioned, the magnitude of c 1, c 2 depends on the conductivity coefficients Therefore, to minimize the bound on the condition number in (21), it is standard practice to seect α of the same order of magnitude as the conductivity coefficients [5] In our case, by choosing α as M 2σ t +σ =: α we optimize the upper bound µ M with respect to the conductivity coefficients Since α 1, this vaue of α usuay aso eads to the ower estimate µ m = 1 Numerica experiments vaidated this choice Our third variant copes with the aready mentioned fact that buiding P V becomes expensive when grid refinements are required Various choices have been discussed in the iterature [10]; for instance, in [3] the vertex preconditioner was chosen as the vertex bock of the Schur compement matrix on a fixed auxiiary coarse mesh, independent of the space discretization We thus approximate P V with the matrix P V c obtained with a fixed coarse mesh, yieding ( ) (R) P P 3 = E 0 0 P V c This variant eads to very moderate (cose to unit) vaues of λ min (P 1 S V c V ) and λ max (P 1 S V c V ) to achieve an estimate in (21) Therefore, we maintained the seection of α as discussed above In Tabe 1 we report numerica experiments with the preconditioners P 1, P 2 and P 3 for the Schur compement system associated with the matrix in (19) The resuts are in cose agreement with the theory: the condition number of the preconditioned matrix grows at most poyogarithmicay with the number of degrees of freedom per subdomain, as indicated by (20) The coumns with α = 1 refer to such parameter seection in P (R) E This corresponds to discarding information on the conductivity tensor M in (2) The worse convergence vaidates our choice and shows the importance of an appropriate choice of the parameter Note that the seected vaue of α was more effective than other choices of simiar magnitude Indeed choosing α = σ and α = σ t for N 2 = 256, n = 5, 10, 20, 40 (fourth row in the tabe) and P 3, we obtained 29, 31, 33, 34 and 42, 45, 48, 49, respectivey In summary, our experiments demonstrate that the proposed variants aow us to imit the computationa cost (the cost of forming P (R) E and P V c is much ower than that for their origina counterparts), with basicay no oss in convergence rate, for the appropriate scaing factor
8 8 Mico Pennacchio and Vaeria Simoncini Tabe 1 Number of conjugate gradient iterations needed to reduce the residua of a factor 10 5 with the preconditioners P 1, P 2, P 3 and P 2, P 3 with α = 1 K = N 2 : # of subdomains n 2 : # of eements per subdomain Symbo * : the preconditioner P 1 coud not be buit due to memory constraints P 1 P 2 P 2 α = 1 P 3 P 3 α = 1 N 2 \n * * * References 1 Achdou, Y, Maday, Y, Widund, O: Substructuring preconditioners for the mortar method in dimension two SIAM J Numer Ana, 36, (1999) 2 Bernardi, C, Maday, Y, Patera, AT: Domain decomposition by the mortar eement method In: Kaper, H et a (eds) Asymptotic and Numerica Method for Partia Differentia Equations with Critica Parameters Dordrecht: Reide, (1993) 3 Bertouzza, S, Pennacchio, M: Preconditioning the mortar method by substructuring: the high order case App Num Ana Comp Math, 1, (2004) 4 Brambe, J H, Pasciak, JE, Schatz, A H: The construction of preconditioners for eiptic probems by substructuring, I Math Comp, 47, (1986) 5 Chan, T F, and Mathew, T P, Domain Decomposition Agorithms, In: Acta Numerica 1994, Cambridge University Press, (1994) 6 Y Notay, Agebraic mutigrid and agebraic mutieve methods: a theoretica comparison Numer Lin Ag App, 12, (2005) 7 Pennacchio, M and Simoncini, V: Efficient agebraic soution of reaction diffusion systems for the cardiac excitation process J Comput App Math, 145 (1): (2002) 8 Pennacchio, M : The mortar finite eement method for the cardiac bidomain mode of extraceuar potentia, J Sci Comput, 20, n2, pag (2004) 9 Pennacchio, M: The Mortar Finite Eement Method for Cardiac Reaction- Diffusion Modes, Computers in Cardioogy 2004; IEEE Proc, 31: (2004) 10 Tosei, A, Widund, O Domain decomposition methods agorithms and theory, Springer Series in Computationa Mathematics, voume 34 Springer-Verag, Berin (2005) 11 Wohmuth, B: Discretization Methods and Iterative Sovers Based on Domain Decomposition, Lecture Notes in Computationa Science and Engineering, voume 17, Springer (2001)
Multiscale Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Problems UNCORRECTED PROOF
AQ Mutiscae Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Probems Yachin Efendiev, Juan Gavis, Raytcho Lazarov, Svetozar Margenov 2, 4 and Jun Ren 5 Department of Mathematics, TAMU,
More informationBCCS TECHNICAL REPORT SERIES
BCCS TCHNICAL RPORT SRIS The Crouzeix-Raviart F on Nonmatching Grids with an Approximate Mortar Condition Taa Rahman, Petter Bjørstad and Xuejun Xu RPORT No. 17 March 006 UNIFOB the University of Bergen
More informationA two-level Schwarz preconditioner for heterogeneous problems
A two-eve Schwarz preconditioner for heterogeneous probems V. Doean 1, F. Nataf 2, R. Scheich 3 and N. Spiane 2 1 Introduction Coarse space correction is essentia to achieve agorithmic scaabiity in domain
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationMultigrid Method for Elliptic Control Problems
J OHANNES KEPLER UNIVERSITÄT LINZ Netzwerk f ür Forschung, L ehre und Praxis Mutigrid Method for Eiptic Contro Probems MASTERARBEIT zur Erangung des akademischen Grades MASTER OF SCIENCE in der Studienrichtung
More informationA FETI-DP method for Crouzeix-Raviart finite element discretizations 1
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vo. 1 (2001), No. 1, pp. 1 19 c 2001 Editoria Board of the Journa Computationa Methods in Appied Mathematics Joint Co. Ltd. A FETI-DP method for Crouzeix-Raviart
More informationPartial permutation decoding for MacDonald codes
Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationCryptanalysis of PKP: A New Approach
Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vo. 36, No. 2, pp. 58 606 c 999 Society for Industria and Appied Mathematics OVERLAPPING NONMATCHING GRID MORTAR ELEMENT METHODS FOR ELLIPTIC PROBLEMS XIAO-CHUAN CAI, MAKSYMILIAN DRYJA,
More informationSUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS
ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia
More informationNumerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet
Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationConvergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems
Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,
More informationNumerische Mathematik
Numer. Math. (1999) 84: 97 119 Digita Object Identifier (DOI) 10.1007/s002119900096 Numerische Mathematik c Springer-Verag 1999 Mutigrid methods for a parameter dependent probem in prima variabes Joachim
More informationParallel algebraic multilevel Schwarz preconditioners for a class of elliptic PDE systems
Parae agebraic mutieve Schwarz preconditioners for a cass of eiptic PDE systems A. Borzì Vaentina De Simone Daniea di Serafino October 3, 213 Abstract Agebraic mutieve preconditioners for agebraic probems
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationHILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT
Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is
More informationFEMAG: A High Performance Parallel Finite Element Toolbox for Electromagnetic Computations
Internationa Journa of Energy and Power Engineering 2016; 5(1-1): 57-64 Pubished onine November 28, 2015 (http://www.sciencepubishinggroup.com/j/ijepe) doi: 10.11648/j.ijepe.s.2016050101.19 ISSN: 2326-957X
More informationBALANCING REGULAR MATRIX PENCILS
BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationModel Solutions (week 4)
CIV-E16 (17) Engineering Computation and Simuation 1 Home Exercise 6.3 Mode Soutions (week 4) Construct the inear Lagrange basis functions (noda vaues as degrees of freedom) of the ine segment reference
More informationIntroduction to Riemann Solvers
CO 5 BOLD WORKSHOP 2, 2 4 June Introduction to Riemann Sovers Oskar Steiner . Three major advancements in the numerica treatment of the hydrodynamic equations Three major progresses in computationa fuid
More informationOn Complex Preconditioning for the solution of the Kohn-Sham Equation for Molecular Transport
On Compex Preconditioning for the soution of the Kohn-Sham Equation for Moecuar Transport Danie Osei-Kuffuor Lingzhu Kong Yousef Saad James R. Cheikowsky December 23, 2007 Abstract This paper anayzes the
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationSmoothers for ecient multigrid methods in IGA
Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third
More informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique Vol. 35, N o 4, 2001, pp THE MORTAR METHOD IN THE WAVELET CONTEXT
Mathematica Modeing and Numerica Anaysis ESAIM: MAN Modéisation Mathématique et Anayse Numérique Vo. 35, N o 4, 001, pp. 647 673 THE MORTAR METHOD IN THE WAVELET CONTEXT Sivia Bertouzza 1 and Vaérie Perrier
More informationUnconditional security of differential phase shift quantum key distribution
Unconditiona security of differentia phase shift quantum key distribution Kai Wen, Yoshihisa Yamamoto Ginzton Lab and Dept of Eectrica Engineering Stanford University Basic idea of DPS-QKD Protoco. Aice
More informationCopyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU
Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationAlgorithms to solve massively under-defined systems of multivariate quadratic equations
Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Heilbronner Str. 10, D Berlin - Wilmersdorf
Konrad-Zuse-Zentrum für Informationstechnik Berin Heibronner Str. 10, D-10711 Berin - Wimersdorf Raf Kornhuber Monotone Mutigrid Methods for Eiptic Variationa Inequaities II Preprint SC 93{19 (Marz 1994)
More informationDistributed average consensus: Beyond the realm of linearity
Distributed average consensus: Beyond the ream of inearity Usman A. Khan, Soummya Kar, and José M. F. Moura Department of Eectrica and Computer Engineering Carnegie Meon University 5 Forbes Ave, Pittsburgh,
More informationPrimal and dual active-set methods for convex quadratic programming
Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:
More informationVTU-NPTEL-NMEICT Project
MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid
More informationWavelet Galerkin Solution for Boundary Value Problems
Internationa Journa of Engineering Research and Deveopment e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Voume, Issue 5 (May 24), PP.2-3 Waveet Gaerkin Soution for Boundary Vaue Probems D. Pate, M.K.
More informationAbsolute Value Preconditioning for Symmetric Indefinite Linear Systems
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Absoute Vaue Preconditioning for Symmetric Indefinite Linear Systems Vecharynski, E.; Knyazev, A.V. TR2013-016 March 2013 Abstract We introduce
More informationREVISITING AND EXTENDING INTERFACE PENALTIES FOR MULTI-DOMAIN SUMMATION-BY-PARTS OPERATORS
REVISITING AND EXTENDING INTERFACE PENALTIES FOR MULTI-DOMAIN SUMMATION-BY-PARTS OPERATORS Mark H. Carpenter, Jan Nordström David Gottieb Abstract A genera interface procedure is presented for muti-domain
More informationPreconditioned Locally Harmonic Residual Method for Computing Interior Eigenpairs of Certain Classes of Hermitian Matrices
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Preconditioned Locay Harmonic Residua Method for Computing Interior Eigenpairs of Certain Casses of Hermitian Matrices Vecharynski, E.; Knyazev,
More information6 Wave Equation on an Interval: Separation of Variables
6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.
More informationHaar Decomposition and Reconstruction Algorithms
Jim Lambers MAT 773 Fa Semester 018-19 Lecture 15 and 16 Notes These notes correspond to Sections 4.3 and 4.4 in the text. Haar Decomposition and Reconstruction Agorithms Decomposition Suppose we approximate
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationAppendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model
Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationResearch Article Numerical Range of Two Operators in Semi-Inner Product Spaces
Abstract and Appied Anaysis Voume 01, Artice ID 846396, 13 pages doi:10.1155/01/846396 Research Artice Numerica Range of Two Operators in Semi-Inner Product Spaces N. K. Sahu, 1 C. Nahak, 1 and S. Nanda
More informationA Robust Multigrid Method for Isogeometric Analysis using Boundary Correction. C. Hofreither, S. Takacs, W. Zulehner. G+S Report No.
A Robust Mutigrid Method for Isogeometric Anaysis using Boundary Correction C. Hofreither, S. Takacs, W. Zuehner G+S Report No. 33 Juy 2015 A Robust Mutigrid Method for Isogeometric Anaysis using Boundary
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationPreprint BUW-SC 11/1. Bergische Universität Wuppertal. Fachbereich C Mathematik und Naturwissenschaften. Mathematik
Preprint BUW-SC 11/1 Bergische Universität Wupperta Fachbereich C Mathematik und Naturwissenschaften Mathematik A. Brandt, J. Brannick, K. Kah, and I. Livshits Bootstrap AMG Januar 2011 http://www-ai.math.uni-wupperta.de/scicomp/
More informationSymbolic models for nonlinear control systems using approximate bisimulation
Symboic modes for noninear contro systems using approximate bisimuation Giordano Poa, Antoine Girard and Pauo Tabuada Abstract Contro systems are usuay modeed by differentia equations describing how physica
More information(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].
PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationincompete LU factorizations and the agebraic mutigrid method. These methods possess features of mutieve methods as we as some features of ILU factoriz
Adapting Agebraic Recursive Mutieve Sovers (ARMS) for soving CFD probems Yousef Saad Λ Azzeddine Souaimani y Ridha Touihri y Abstract This paper presents resuts using preconditioners that are based on
More informationSASIMI: Sparsity-Aware Simulation of Interconnect-Dominated Circuits with Non-Linear Devices
SASIMI: Sparsity-Aware Simuation of Interconnect-Dominated Circuits with Non-Linear Devices Jitesh Jain, Stephen Cauey, Cheng-Kok Koh, and Venkataramanan Baakrishnan Schoo of Eectrica and Computer Engineering
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationA HIERARCHICAL LOW-RANK SCHUR COMPLEMENT PRECONDITIONER FOR INDEFINITE LINEAR SYSTEMS
A HERARCHCAL LOW-RANK SCHUR COMPLEMENT PRECONDTONER FOR NDEFNTE LNEAR SYSTEMS GEOFFREY DLLON, YUANZHE X, AND YOUSEF SAAD Abstract. Nonsymmetric and highy indefinite inear systems can be quite difficut
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More information18-660: Numerical Methods for Engineering Design and Optimization
8-660: Numerica Methods for Engineering esign and Optimization in i epartment of ECE Carnegie Meon University Pittsburgh, PA 523 Side Overview Conjugate Gradient Method (Part 4) Pre-conditioning Noninear
More informationM. Aurada 1,M.Feischl 1, J. Kemetmüller 1,M.Page 1 and D. Praetorius 1
ESAIM: M2AN 47 (2013) 1207 1235 DOI: 10.1051/m2an/2013069 ESAIM: Mathematica Modeing and Numerica Anaysis www.esaim-m2an.org EACH H 1/2 STABLE PROJECTION YIELDS CONVERGENCE AND QUASI OPTIMALITY OF ADAPTIVE
More informationGeneral Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping
Journa of Mathematica Research with Appications Jan.,, Vo. 3, No., pp. 53 6 DOI:.377/j.issn:95-65...7 Http://jmre.dut.edu.cn Genera Decay of Soutions in a Viscoeastic Equation with Noninear Locaized Damping
More informationA New Modeling Approach for Planar Beams: Finite-Element Solutions based on Mixed Variational Derivations
A New Modeing Approach for Panar Beams: Finite-Eement Soutions based on Mixed Variationa Derivations Ferdinando Auricchio,a,b,c,d, Giuseppe Baduzzi a,e, Caro Lovadina c,d,e May 31, 2 a Dipartimento di
More informationc 2011 Society for Industrial and Applied Mathematics BOOTSTRAP AMG
SIAM J. SCI. COMPUT. Vo. 0, No. 0, pp. 000 000 c 2011 Society for Industria and Appied Mathematics BOOTSTRAP AMG A. BRANDT, J. BRANNICK, K. KAHL, AND I. LIVSHITS Abstract. We deveop an agebraic mutigrid
More informationMultilevel Monte Carlo Methods and Applications to Elliptic PDEs with Random Coefficients
CVS manuscript No. wi be inserted by the editor) Mutieve Monte Caro Methods and Appications to Eiptic PDEs with Random Coefficients K.A. Ciffe M.B. Gies R. Scheich A.L. Teckentrup Received: date / Accepted:
More informationASummaryofGaussianProcesses Coryn A.L. Bailer-Jones
ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe
More informationComponentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems
Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,
More informationSource and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems
Source and Reay Matrices Optimization for Mutiuser Muti-Hop MIMO Reay Systems Yue Rong Department of Eectrica and Computer Engineering, Curtin University, Bentey, WA 6102, Austraia Abstract In this paper,
More informationSequential Decoding of Polar Codes with Arbitrary Binary Kernel
Sequentia Decoding of Poar Codes with Arbitrary Binary Kerne Vera Miosavskaya, Peter Trifonov Saint-Petersburg State Poytechnic University Emai: veram,petert}@dcn.icc.spbstu.ru Abstract The probem of efficient
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna
More informationc 2000 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vo. 2, No. 5, pp. 909 926 c 2000 Society for Industria and Appied Mathematics A DEFLATED VERSION OF THE CONJUGATE GRADIENT ALGORITHM Y. SAAD, M. YEUNG, J. ERHEL, AND F. GUYOMARC H
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationA nonlinear multigrid for imaging electrical conductivity and permittivity at low frequency
INSTITUTE OF PHYSICS PUBLISHING INVERSE PROBLEMS Inverse Probems 17 (001) 39 359 www.iop.org/journas/ip PII: S066-5611(01)14169-9 A noninear mutigrid for imaging eectrica conductivity and permittivity
More informationStatistics for Applications. Chapter 7: Regression 1/43
Statistics for Appications Chapter 7: Regression 1/43 Heuristics of the inear regression (1) Consider a coud of i.i.d. random points (X i,y i ),i =1,...,n : 2/43 Heuristics of the inear regression (2)
More informationA Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem
Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,
More informationProbabilistic Graphical Models
Schoo of Computer Science Probabiistic Graphica Modes Gaussian graphica modes and Ising modes: modeing networks Eric Xing Lecture 0, February 0, 07 Reading: See cass website Eric Xing @ CMU, 005-07 Network
More informationCoupling of LWR and phase transition models at boundary
Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of
More informationEXPONENTIAL DECAY OF SOLUTIONS TO A VISCOELASTIC EQUATION WITH NONLINEAR LOCALIZED DAMPING
Eectronic Journa of Differentia Equations, Vo. 24(24), No. 88, pp. 1 1. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (ogin: ftp) EXPONENTIAL DECAY
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.
More informationMA 201: Partial Differential Equations Lecture - 11
MA 201: Partia Differentia Equations Lecture - 11 Heat Equation Heat conduction in a thin rod The IBVP under consideration consists of: The governing equation: u t = αu xx, (1) where α is the therma diffusivity.
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationP-MULTIGRID PRECONDITIONERS APPLIED TO HIGH-ORDER DG AND HDG DISCRETIZATIONS
6th European Conference on Computationa Mechanics (ECCM 6) 7th European Conference on Computationa Fuid Dynamics (ECFD 7) 11-15 June 2018, Gasgow, UK P-MULTIGRID PRECONDITIONERS APPLIED TO HIGH-ORDER DG
More informationMultilayer Kerceptron
Mutiayer Kerceptron Zotán Szabó, András Lőrincz Department of Information Systems, Facuty of Informatics Eötvös Loránd University Pázmány Péter sétány 1/C H-1117, Budapest, Hungary e-mai: szzoi@csetehu,
More informationMechanism Deduction from Noisy Chemical Reaction Networks
Mechanism Deduction from Noisy Chemica Reaction Networks arxiv:1803.09346v1 [physics.chem-ph] 25 Mar 2018 Jonny Proppe and Markus Reiher March 28, 2018 ETH Zurich, Laboratory of Physica Chemistry, Vadimir-Preog-Weg
More informationQuantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal
Adv. Studies Theor. Phys., Vo. 6, 01, no. 3, 19-133 Quantum Eectrodynamica Basis for Wave Propagation through Photonic Crysta 1 N. Chandrasekar and Har Narayan Upadhyay Schoo of Eectrica and Eectronics
More informationMARKOV CHAINS AND MARKOV DECISION THEORY. Contents
MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationKing Fahd University of Petroleum & Minerals
King Fahd University of Petroeum & Mineras DEPARTMENT OF MATHEMATICAL SCIENCES Technica Report Series TR 369 December 6 Genera decay of soutions of a viscoeastic equation Saim A. Messaoudi DHAHRAN 3161
More informationLecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling
Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest
More informationNumerical methods for elliptic partial differential equations Arnold Reusken
Numerica methods for eiptic partia differentia equations Arnod Reusken Preface This is a book on the numerica approximation of partia differentia equations. On the next page we give an overview of the
More informationEfficient numerical solution of Neumann problems on complicated domains. Received: July 2005 / Revised version: April 2006 Springer-Verlag 2006
CACOO 43, 95 120 (2006) DOI: 10.1007/s10092-006-0118-4 CACOO Springer-Verag 2006 Efficient numerica soution of Neumann probems on compicated domains Serge Nicaise 1, Stefan A. Sauter 2 1 Macs, Université
More informationData Mining Technology for Failure Prognostic of Avionics
IEEE Transactions on Aerospace and Eectronic Systems. Voume 38, #, pp.388-403, 00. Data Mining Technoogy for Faiure Prognostic of Avionics V.A. Skormin, Binghamton University, Binghamton, NY, 1390, USA
More informationA SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS
A SIPLIFIED DESIGN OF ULTIDIENSIONAL TRANSFER FUNCTION ODELS Stefan Petrausch, Rudof Rabenstein utimedia Communications and Signa Procesg, University of Erangen-Nuremberg, Cauerstr. 7, 958 Erangen, GERANY
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More informationEfficiently Generating Random Bits from Finite State Markov Chains
1 Efficienty Generating Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown
More information