Perturbation Analyses for the Cholesky Factorization with Backward Rounding Errors
|
|
- Joseph Payne
- 5 years ago
- Views:
Transcription
1 Perturbation Analyses for the holesky Fatorization with Bakward Rounding Errors Xiao-Wen hang Shool of omputer Siene, MGill University, Montreal, Quebe, anada, H3A A7 Abstrat. This paper gives perturbation analyses of the holesky fatorization with the form of perturbations we ould expet from the equivalent bakward error in A resulting from numerially stable omputations. The analyses more aurately reflet the sensitivity of the problem than previous suh results. Both numerial results and an analysis show the standard method of symmetri pivoting usually improves the ondition of the problem. It follows that the omputed R will usually have more auray when we use the standard symmetri pivoting strategy. Introdution Let A R n n be a symmetri positive definite matrix. Then A has a unique holesky fatorization of the form A = R T R, where R is upper triangular with positive diagonal elements and is alled the holesky fator of A. Suppose there is a perturbation A in A and A + A is also symmetri positive definite, then A + A has a unique holesky fatorization A + A = (R + R) T (R + R). The goal of the perturbation analysis is to give a bound on R (or R ) in terms of A (or A ). There have been several papers dealing with the perturbation analysis for the holesky fator for a general perturbation A, see [], [8], and [0] for normwise analyses, and [, ] for omponentwise analyses. Reently perturbation results of a different flavor were presented in [4]. For A orresponding to the bakward rounding error in A resulting from standard numerial stable omputations, a nie norm-based perturbation result was presented in [4]. Note [, ] also inluded omponent perturbation bounds for a different and somewhat ompliated form of the bakward rounding error in A. The purpose of this paper is to establish new first-order perturbation bounds for A orresponding to the bakward rounding error, whih are sharper than the orresponding first-order result in [4]. Now we introdue notation. To simplify the presentation, for any n n matrix X, we define the upper triangular matrix x x x n 0 up(x) x x n......, () 0 0 x nn
2 and the diagonal matrix diag(x) diag(x, x,..., x nn ) diag(x ii ). For any matrix ( ij ) [,..., n ] R n n, denote by (i) j the vetor of the leading i elements of j. Using u to denote upper, we define uve() [ () T, () T,, n (n) T ] T, () the vetor formed by staking the olumns of the upper triangular part of into one long vetor. Note that for any upper triangular X, uve(x) = X F. To help our norm analysis, for any matrix R n n we define duve() diag (,,,...,,,...,, )uve(). (3) }{{}}{{} n Thus for any symmetri matrix X we have duve(x) = X F. Finally let D n be the set of all n n real positive-definite diagonal matries. At the end of this setion we summarize the strongest norm-based results for a general perturbation A presented in [], as those results will be ompared with the new results for A having the form of bakward rounding error. Theorem. (see []) Let A R n n be symmetri positive definite with the holesky fatorization A = R T R. Let A be symmetri. If ǫ A F / A satisfies κ (A)ǫ <, where κ (A) A A, then A+ A has the holesky fatorization A + A = (R + R) T (R + R), suh that where R F R κ (A)ǫ + O(ǫ ), (4) κ/ (A) κ (A) κ (A) κ (A), (5) κ (A) W R A /, (6) κ (A) inf D D n κ (A, D), κ (A, D) κ (R)κ (D R), (7) with W R R n(n+) n(n+) defined by r r r r r W R = rn r. (8) rn rn r r r n r n r nn
3 If standard symmetri pivoting is used in omputing the holesky fatorization PAP T = R T R, we have κ/ (A) κ (PAP T ) κ (PAP T ) κ / (A) n(n + )(4 n + 6n )/6. In [] we showed the bound (4) is attainable to first order in ǫ. So κ (A) is the ondition number for the holesky fator R. It is expensive to estimate or evaluate κ (A) diretly by using the usual methods. Fortunately we an use κ (A) as an approximation of κ (A). Aording to the well-known result of van de Sluis [3], if D is hosen to equilibrate the rows of R, i.e., all rows of D R have unit length, then κ (A, D) with suh a D is at most a fator of n off the infimum κ (A). Notie that κ (A, D) an be estimated by a standard ondition estimator in O(n ) flops. Main results In this setion we first derive a tight perturbation bound, leading to the ondition number χ (A) for perturbation having the form of equivalent bakward rounding error. Also we derive a pratial estimate χ (A) of χ (A). Then we ompare χ (A) with κ (A), and χ (A) with κ (A). Finally we show how standard pivoting improves the ondition number χ (A). Before proeeding we introdue the following result presented in [3], see also [6, Theorems 0.5 and 0.7]. Lemma. (see [3]) Let A D HD R n n be a symmetri positive definite floating point matrix, where D diag(a) /. If nǫ H <, (9) where ǫ (n+)u/( (n+)u) with u being the unit round-off, then holesky fatorization applied to A sueeds (barring underflow and overflow) and produes a nonsingular R, whih satisfies where d i = a / ii. A + A = R T R, A ǫ dd T, (0) This lemma is appliable to any standard algorithms for the holesky fatorization (see, for example, [5] for the standard algorithms). Based on this, we establish the following result. Theorem 3. Suppose all of the assumptions in Lemma hold. Let A = R T R be the holesky fatorization of A. Set ˆR RD. Then for the perturbation A and result R in (0) we have with R R R that R F R χ (A)ǫ + O(ǫ ), () χ (A) χ (A) n H, () χ (A) H, (3)
4 where χ (A) n D W ˆR / R, χ (A) inf D D n χ (A, D), χ (A, D) n ˆR ˆR D D R / R, with D diag (a /, a/, a/ } {{ },..., a / nn, a / nn,..., a / nn ) R n(n+) }{{} n n(n+), (4) and W ˆR being just W R in (8) with eah entry r ij replaed by ˆr ij. The bound in () is approximately attainable to first-order in ǫ. Proof. Let G A/ǫ. By (0) and (9) it is easy to show A + tg is symmetri positive definite for all t [0, ǫ], and so it has the holesky fatorization A + tg = R T (t)r(t), t ǫ, (5) with R(0) = R and R(ǫ) = R R + R. Differentiating (5) and setting t = 0 in the result gives G = R T Ṙ(0) + Ṙ(0)T R, (6) whih, ombined with R = ˆRD, gives ˆR T Ṙ(0)D + D Ṙ T (0) ˆR = D GD, (7) The upper and lower triangular parts of (7) ontain idential information, and it is easy to show that by using () and (3) the upper triangular part of (7) an be rewritten as the following matrix-vetor equation form W ˆR uve(ṙ(0)d ) = duve(d GD ). (8) Notie from () that uve(ṙ(0)d ) = D uve(ṙ(0)) with D as in (4), then from (8) we have uve(ṙ(0)) = D W ˆR whih with G = A/ǫ and (0) gives duve(d GD ), (9) Ṙ(0) F D W ˆR D dd T D F = n D W ˆR. (0) The Taylor expansion of R(t) about t = 0 gives at t = ǫ R + R R = R(ǫ) = R(0) + ǫṙ(0) + O(ǫ ), () whih, ombined with (0), gives (). Obviously there exists a symmetri matrix F R n n suh that D W ˆR duve(d FD ) = D W ˆR D FD F.
5 Then by taking G = (min fij 0 d i d j / f ij )F, we have A ǫ dd T and from (9) that Ṙ(0) F = ( min d id j / f ij ) D W f ij 0 ˆR D FD F D W ˆR, whih shows the first-order bound in () is approximately attained for suh G. It remains to prove () and (3). From (6) with R = ˆRD, it follows that Ṙ(0)R + R T Ṙ T (0) = ˆR T D GD ˆR. From this we see with notation up (see ()) that Ṙ(0) = up( ˆR T D GD ˆR )R, () so with (9) and the fat that up(xd) = up(x)d and up(x) F X F for any X R n n and D D n we have D W ˆR duve(d GD ) = up( ˆR T D GD ˆR )R F ˆR D GD F ˆR D D R. Atually this holds for any symmetri G R n n sine it was essentially obtained from the matrix equation R T X + X T R = G with X triangular. Notie duve(d GD ) = D GD F, thus we must have D W ˆR ˆR ˆR D D R. Sine this is true for any D D n, we see that the first inequality in () holds. The seond inequality in () follows from χ (A, I) = n ˆR = n H. Similarly we an prove (3) by using the fat that up(x) F X F for any symmetri X R n n. Sine the first-order bound () is approximately attainable, χ (A) an be regarded as the ondition number for the holesky fatorization with the form of perturbation error satisfying (0). Reently Drma et al. [4] gave a perturbation analysis with the same form of perturbation error. Their asymptoti perturbation bound an be presented as follows: R R ( + log n )n H ǫ + O(ǫ ). (3) This suggests that H an be regarded as a onditioning measure. From () we see our bound () is sharper [ than (3). ] In fat [ it ] an be muh sharper. δ δ As an illustration suppose A = δ δ + δ 4, R = 0 δ. Then it is easy to show for small δ > 0, χ (A) χ (A) = O(/δ), H = O(/δ ). Thus (3) may overestimate the true relative error of the omputed holesky fator, and the approximation χ (A) to the ondition number χ (A) is a signifiant improvement to H. Furthermore it is easy to see H is invariant if pivoting is used in omputing the holesky fatorization of A, whereas χ (A)
6 and χ (A) depend on any pivoting. Thus χ (A) and χ (A) more losely reflets the true sensitivity of the holesky fatorization than H. As far as we an see, it is expensive to ompute or approximate χ (A) diretly by the usual approah. Fortunately χ (A) is quite easy to estimate. By the result of van der Sluis [3], ˆR D D R will be nearly minimum when the rows of D R are equilibrated. Then a proedure for obtaining an O(n ) ondition estimator for the holesky fatorization with bakward rounding errors is to hoose D = D r diag( R(i, :) ), and then use a standard norm estimator to estimate all fators in χ (A, D). Numerial experiments suggest usually χ (A) is a reasonable approximation of χ (A). But χ (A) an still be very larger than χ (A), even though it an be muh smaller than n H. For the example above we have χ (A) = O(). Numerial tests suggest that usually χ (A) is smaller or muh smaller than κ (A), defined by (6). This is not surprising sine (0) provides more information about the perturbation in the data. In [] we proved the following results. Theorem 4. n χ maxi aii (A) κ (A) min i a ii χ (A), (4) n χ (A) κ maxi aii (A) min i a ii χ (A). The first inequality in (4) is attainable, sine equality will hold by taking A = I with > 0. The seond inequality is at least nearly attainable. In fat taking A = R T R, where R = diag(δ n, δ n,...,δ, ) + e e T n with small δ > 0, we easily obtain max i a ii κ (A) = O( ), χ δn (A) = min i a ii δ n O() = O( ). δn This example also suggests that possibly κ (A) is muh larger than χ (A) if the maximum element is muh larger than the minimum one on the diagonal of A. The standard symmetri pivoting strategy an usually improve χ (A), just as it an usually improve κ (A). In [] we showed the following theorem. Theorem5. Let A R n n be symmetri positive definite with the holesky fatorization PAP T = R T R when the standard pivoting strategy is used. Then χ (PAP T ) χ (PAP T ) H / n n(n + )(4 n + 6n )/6, (5) where PAP T D HD with D diag(pap T ) /. One may not be impressed by the 4 n fator in the upper bound, and may wonder if it an signifiantly be improved. In fat we an prove for any n the upper bound an nearly be approximated by a parametrized family of matries (f. [, pp.5 6]). But suh examples are rare in pratie. Our numerial experiments onfirms that χ (PAP T ) is usually (muh) smaller than χ (A). Thus holesky fatorization with standard symmetri pivoting will usually give more aurate R. By following the approah of [7], it is straightforward, but detailed and lengthy, to extend the first-order results to provide rigorous perturbation bounds. For suh results, see [].
7 Table. Results for Pasal matries without pivoting. n χ (A) χ (A, D r) κ (A) κ (A, D r) 5.0e 0.0e e 0.0e+00 7.e 0.e e+00.e e e e e+0 9.7e e+0 3.6e e+0.3e+0 5.5e+0 4.8e+0 3.0e+0 5.7e+0.4e e+0 4.9e+03.5e e+0 9.e+03.5e+03 5.e+04.e e e+04.9e e+05.0e e e+05.5e e+06.7e e+04.3e+06.3e e+07.5e e+05.5e+07.e e+08.4e+08.7e+06 9.e e e+09.e e e e+08.0e+.e e e e+09.e+ 9.8e+0 4.0e+08.e+0 7.e+0.4e+3 8.8e e+08.3e+ 6.5e+.6e+4 7.9e+ Table. Results for Pasal matries with pivoting, Ã PAP T. n χ (Ã) χ(ã, Dr) κ(ã) κ(ã, Dr) 5.0e 0.0e e 0.0e+00 7.e 0.6e e+00.5e+00 4.e e e+00.8e+0 5.e+00.6e+0 3.6e+0 4.3e+0 6.e+0.e+0 8.0e+0 3.0e e+0.e+0 8.3e+0 3.3e+0.5e e+0 6.7e+0.5e+0.3e+03.e e+0.e e+0 5.e+03.0e e+0 7.8e e+03.4e+04.7e+06 9.e+03.7e+04.6e+04.0e+05.5e e e e e+05.4e e+03.7e+05.4e+05.8e+06.e+09.8e+04 9.e e+05 8.e+06.e e+04.9e+06 3.e+06 3.e e e e+06.3e+07.e e+ 5.e+05.8e e e e+ 3 Numerial experiments In setion we presented a new first-order perturbation bound for the holesky fator with the hange aused by bakward rounding errors, defined χ (A) n D W ˆR / A / as the ondition number, and suggested χ (A) ould be estimated in pratie by χ (A, D r ) n ˆR ˆR D r Dr R / R with D r = diag( R(i, :) ), whih an be estimated by standard norm estimators
8 in O(n ) flops. Our new first-order result is (muh) better than the previous orresponding result. Also we ompare χ (A) with κ (A), the ondition number for general perturbation A, and ompare the orresponding estimates χ (A) with κ (A) as well. Now we give a set of examples to show our findings. The matries are n n Pasal matries, (with elements a j = a i =, a ij = a i,j + a i,j ), n =,,..., n. The results are shown in Table without pivoting and in Table with pivoting. Note in Tables and how n H an be worse than χ (A). In Table pivoting is seen to give a signifiant improvement to χ (A). Also we observe from both the tables that χ (A) is a reasonable approximation of χ (A). We see χ (A) is smaller than κ (A) for n >. Aknowledgement. I would like to thank hris Paige for his valuable omments and suggestions. Referenes. hang, X.-W.: Perturbation analysis of some matrix fatorizations. Ph.D thesis, MGill University, February 997. hang, X.-W., Paige,.., Stewart, G. W.: New perturbation analyses for the holesky fatorization. IMA J. Numer. Anal., 6 (996) Demmel, J. W.: On floating point errors in holesky. Tehnial Report S-89-87, Department of omputer Siene, University of Tennessee, 989. LAPAK Working Note 4 4. Drma, Z., Omladi, M., Veselić, K.: On the perturbation of the holesky fatorization. SIAM J. Matrix Anal. Appl., 5 (994) Golub, G. H., Van Loan,. F.: Matrix omputations. Third Edition. The Johns Hopkins University Press, Baltimore, Maryland, Higham, N. J.: Auray and Stability of Numerial Algorithms. Soiety for Industrial and Applied Mathematis, Philadelphia, PA, Stewart, G. W.: Error and perturbation bounds for subspaes assoiated with ertain eigenvalue problems. SIAM Rev., 5 (973) Stewart, G. W.: Perturbation bounds for the QR fatorization of a matrix. SIAM J. Numer. Anal., 4 (977) Stewart, G. W.: On the perturbation of LU, holesky, and QR fatorizations. SIAM J. Matrix Anal. Appl., 4 (993) Sun, J.-G.: Perturbation bounds for the holesky and QR fatorization. BIT, 3 (99) Sun, J.-G.: Rounding-error and perturbation bounds for the holesky and LDL T fatorizations. Linear Algebra and Appl., 73 (99) Sun, J.-G.: omponentwise perturbation bounds for some matrix deompositions. BIT, 3 (99) van der Sluis, A.: ondition numbers and equilibration of matries. Numerishe Mathematik, 4 (966) 4 3 This artile was proessed using the L A TEX maro pakage with LLNS style
On the Perturbation of the Q-factor of the QR Factorization
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. ; :1 6 [Version: /9/18 v1.] On the Perturbation of the Q-factor of the QR Factorization X.-W. Chang McGill University, School of Comptuer
More informationMULTIPLICATIVE PERTURBATION ANALYSIS FOR QR FACTORIZATIONS. Xiao-Wen Chang. Ren-Cang Li. (Communicated by Wenyu Sun)
NUMERICAL ALGEBRA, doi:10.3934/naco.011.1.301 CONTROL AND OPTIMIZATION Volume 1, Number, June 011 pp. 301 316 MULTIPLICATIVE PERTURBATION ANALYSIS FOR QR FACTORIZATIONS Xiao-Wen Chang School of Computer
More informationMultiplicative Perturbation Analysis for QR Factorizations
Multiplicative Perturbation Analysis for QR Factorizations Xiao-Wen Chang Ren-Cang Li Technical Report 011-01 http://www.uta.edu/math/preprint/ Multiplicative Perturbation Analysis for QR Factorizations
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationc 2010 Society for Industrial and Applied Mathematics
SIAM J MATRIX ANAL APPL Vol 31, No 5, pp 841 859 c 010 Society for Industrial and Applied Mathematics RIGOROUS PERTURBATION BOUNDS OF SOME MATRIX FACTORIZATIONS X-W CHANG AND D STEHLÉ Abstract This article
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationREFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS. 1. Introduction
Version of 5/2/2003 To appear in Advanes in Applied Mathematis REFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS MATTHIAS BECK AND SHELEMYAHU ZACKS Abstrat We study the Frobenius problem:
More informationA NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM
NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum
More informationLecture 7: Sampling/Projections for Least-squares Approximation, Cont. 7 Sampling/Projections for Least-squares Approximation, Cont.
Stat60/CS94: Randomized Algorithms for Matries and Data Leture 7-09/5/013 Leture 7: Sampling/Projetions for Least-squares Approximation, Cont. Leturer: Mihael Mahoney Sribe: Mihael Mahoney Warning: these
More informationInstitute for Advanced Computer Studies. Department of Computer Science. On the Perturbation of. LU and Cholesky Factors. G. W.
University of Maryland Institute for Advanced Computer Studies Department of Computer Science College Park TR{95{93 TR{3535 On the Perturbation of LU and Cholesky Factors G. W. Stewart y October, 1995
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationWe first repeat some well known facts about condition numbers for normwise and componentwise perturbations. Consider the matrix
BIT 39(1), pp. 143 151, 1999 ILL-CONDITIONEDNESS NEEDS NOT BE COMPONENTWISE NEAR TO ILL-POSEDNESS FOR LEAST SQUARES PROBLEMS SIEGFRIED M. RUMP Abstract. The condition number of a problem measures the sensitivity
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationMOLECULAR ORBITAL THEORY- PART I
5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationA variant of Coppersmith s Algorithm with Improved Complexity and Efficient Exhaustive Search
A variant of Coppersmith s Algorithm with Improved Complexity and Effiient Exhaustive Searh Jean-Sébastien Coron 1, Jean-Charles Faugère 2, Guénaël Renault 2, and Rina Zeitoun 2,3 1 University of Luxembourg
More informationAssessing the Performance of a BCI: A Task-Oriented Approach
Assessing the Performane of a BCI: A Task-Oriented Approah B. Dal Seno, L. Mainardi 2, M. Matteui Department of Eletronis and Information, IIT-Unit, Politenio di Milano, Italy 2 Department of Bioengineering,
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationNonreversibility of Multiple Unicast Networks
Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast
More informationA Backward Stable Hyperbolic QR Factorization Method for Solving Indefinite Least Squares Problem
A Backward Stable Hyperbolic QR Factorization Method for Solving Indefinite Least Suares Problem Hongguo Xu Dedicated to Professor Erxiong Jiang on the occasion of his 7th birthday. Abstract We present
More informationAyan Kumar Bandyopadhyay
Charaterization of radiating apertures using Multiple Multipole Method And Modeling and Optimization of a Spiral Antenna for Ground Penetrating Radar Appliations Ayan Kumar Bandyopadhyay FET-IESK, Otto-von-Guerike-University,
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationQ2. [40 points] Bishop-Hill Model: Calculation of Taylor Factors for Multiple Slip
27-750, A.D. Rollett Due: 20 th Ot., 2011. Homework 5, Volume Frations, Single and Multiple Slip Crystal Plastiity Note the 2 extra redit questions (at the end). Q1. [40 points] Single Slip: Calulating
More informationSimplified Buckling Analysis of Skeletal Structures
Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation
More information7.1 Roots of a Polynomial
7.1 Roots of a Polynomial A. Purpose Given the oeffiients a i of a polynomial of degree n = NDEG > 0, a 1 z n + a 2 z n 1 +... + a n z + a n+1 with a 1 0, this subroutine omputes the NDEG roots of the
More informationStability of alternate dual frames
Stability of alternate dual frames Ali Akbar Arefijamaal Abstrat. The stability of frames under perturbations, whih is important in appliations, is studied by many authors. It is worthwhile to onsider
More informationarxiv: v2 [math.pr] 9 Dec 2016
Omnithermal Perfet Simulation for Multi-server Queues Stephen B. Connor 3th Deember 206 arxiv:60.0602v2 [math.pr] 9 De 206 Abstrat A number of perfet simulation algorithms for multi-server First Come First
More informationGeneral Equilibrium. What happens to cause a reaction to come to equilibrium?
General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember
More informationLAPACK-Style Codes for Pivoted Cholesky and QR Updating
LAPACK-Style Codes for Pivoted Cholesky and QR Updating Sven Hammarling 1, Nicholas J. Higham 2, and Craig Lucas 3 1 NAG Ltd.,Wilkinson House, Jordan Hill Road, Oxford, OX2 8DR, England, sven@nag.co.uk,
More informationCoding for Random Projections and Approximate Near Neighbor Search
Coding for Random Projetions and Approximate Near Neighbor Searh Ping Li Department of Statistis & Biostatistis Department of Computer Siene Rutgers University Pisataay, NJ 8854, USA pingli@stat.rutgers.edu
More informationCopyright 2018 Society of Photo-Optical Instrumentation Engineers (SPIE). One print or electronic copy may be made for personal use only.
Copyright 018 Soiety of Photo-Optial Instrumentation Engineers (SPIE) One print or eletroni opy may be made for personal use only Systemati reprodution and distribution, dupliation of any material in this
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationLAPACK-Style Codes for Pivoted Cholesky and QR Updating. Hammarling, Sven and Higham, Nicholas J. and Lucas, Craig. MIMS EPrint: 2006.
LAPACK-Style Codes for Pivoted Cholesky and QR Updating Hammarling, Sven and Higham, Nicholas J. and Lucas, Craig 2007 MIMS EPrint: 2006.385 Manchester Institute for Mathematical Sciences School of Mathematics
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationStudy of EM waves in Periodic Structures (mathematical details)
Study of EM waves in Periodi Strutures (mathematial details) Massahusetts Institute of Tehnology 6.635 partial leture notes 1 Introdution: periodi media nomenlature 1. The spae domain is defined by a basis,(a
More informationProceedings of Meetings on Acoustics
Proeedings of Meetings on Aoustis Volume 19, 2013 http://aoustialsoiety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Arhitetural Aoustis Session 4aAAa: Room Aoustis Computer Simulation I 4aAAa8.
More informationQCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines
DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany
More informationSPLINE ESTIMATION OF SINGLE-INDEX MODELS
SPLINE ESIMAION OF SINGLE-INDEX MODELS Li Wang and Lijian Yang University of Georgia and Mihigan State University Supplementary Material his note ontains proofs for the main results he following two propositions
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationBerry s phase for coherent states of Landau levels
Berry s phase for oherent states of Landau levels Wen-Long Yang 1 and Jing-Ling Chen 1, 1 Theoretial Physis Division, Chern Institute of Mathematis, Nankai University, Tianjin 300071, P.R.China Adiabati
More informationSINCE Zadeh s compositional rule of fuzzy inference
IEEE TRANSACTIONS ON FUZZY SYSTEMS, VOL. 14, NO. 6, DECEMBER 2006 709 Error Estimation of Perturbations Under CRI Guosheng Cheng Yuxi Fu Abstrat The analysis of stability robustness of fuzzy reasoning
More informationAn Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances
An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory
More informationEvaluation of effect of blade internal modes on sensitivity of Advanced LIGO
Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationQuasi-Monte Carlo Algorithms for unbounded, weighted integration problems
Quasi-Monte Carlo Algorithms for unbounded, weighted integration problems Jürgen Hartinger Reinhold F. Kainhofer Robert F. Tihy Department of Mathematis, Graz University of Tehnology, Steyrergasse 30,
More informationFrequency Domain Analysis of Concrete Gravity Dam-Reservoir Systems by Wavenumber Approach
Frequeny Domain Analysis of Conrete Gravity Dam-Reservoir Systems by Wavenumber Approah V. Lotfi & A. Samii Department of Civil and Environmental Engineering, Amirkabir University of Tehnology, Tehran,
More informationarxiv:math/ v1 [math.ca] 27 Nov 2003
arxiv:math/011510v1 [math.ca] 27 Nov 200 Counting Integral Lamé Equations by Means of Dessins d Enfants Sander Dahmen November 27, 200 Abstrat We obtain an expliit formula for the number of Lamé equations
More informationIMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION
09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE
More informationarxiv:gr-qc/ v2 6 Feb 2004
Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this
More informationAn iterative least-square method suitable for solving large sparse matrices
An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationDense LU factorization and its error analysis
Dense LU factorization and its error analysis Laura Grigori INRIA and LJLL, UPMC February 2016 Plan Basis of floating point arithmetic and stability analysis Notation, results, proofs taken from [N.J.Higham,
More informationQuantum Mechanics: Wheeler: Physics 6210
Quantum Mehanis: Wheeler: Physis 60 Problems some modified from Sakurai, hapter. W. S..: The Pauli matries, σ i, are a triple of matries, σ, σ i = σ, σ, σ 3 given by σ = σ = σ 3 = i i Let stand for the
More informationBäcklund Transformations: Some Old and New Perspectives
Bäklund Transformations: Some Old and New Perspetives C. J. Papahristou *, A. N. Magoulas ** * Department of Physial Sienes, Helleni Naval Aademy, Piraeus 18539, Greee E-mail: papahristou@snd.edu.gr **
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS
ON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS MARIO LEFEBVRE and JEAN-LUC GUILBAULT A ontinuous-time and ontinuous-state stohasti proess, denoted by {Xt), t }, is defined from a proess known as
More information(a) We desribe physics as a sequence of events labelled by their space time coordinates: x µ = (x 0, x 1, x 2 x 3 ) = (c t, x) (12.
2 Relativity Postulates (a) All inertial observers have the same equations of motion and the same physial laws. Relativity explains how to translate the measurements and events aording to one inertial
More informationRobust Recovery of Signals From a Structured Union of Subspaces
Robust Reovery of Signals From a Strutured Union of Subspaes 1 Yonina C. Eldar, Senior Member, IEEE and Moshe Mishali, Student Member, IEEE arxiv:87.4581v2 [nlin.cg] 3 Mar 29 Abstrat Traditional sampling
More informationHILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES
HILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES L ERBE, A PETERSON AND S H SAKER Abstrat In this paper, we onsider the pair of seond-order dynami equations rt)x ) ) + pt)x
More informationModel-based mixture discriminant analysis an experimental study
Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,
More informationAverage Rate Speed Scaling
Average Rate Speed Saling Nikhil Bansal David P. Bunde Ho-Leung Chan Kirk Pruhs May 2, 2008 Abstrat Speed saling is a power management tehnique that involves dynamially hanging the speed of a proessor.
More informationComputer Science 786S - Statistical Methods in Natural Language Processing and Data Analysis Page 1
Computer Siene 786S - Statistial Methods in Natural Language Proessing and Data Analysis Page 1 Hypothesis Testing A statistial hypothesis is a statement about the nature of the distribution of a random
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationGeneral probability weighted moments for the three-parameter Weibull Distribution and their application in S-N curves modelling
General probability weighted moments for the three-parameter Weibull Distribution and their appliation in S-N urves modelling Paul Dario Toasa Caiza a,, Thomas Ummenhofer a a KIT Stahl- und Leihtbau, Versuhsanstalt
More informationSensitivity Analysis in Markov Networks
Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores
More informationSensor management for PRF selection in the track-before-detect context
Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationPlanning with Uncertainty in Position: an Optimal Planner
Planning with Unertainty in Position: an Optimal Planner Juan Pablo Gonzalez Anthony (Tony) Stentz CMU-RI -TR-04-63 The Robotis Institute Carnegie Mellon University Pittsburgh, Pennsylvania 15213 Otober
More informationAnalysis of Block LDL T Factorizations for Symmetric Indefinite Matrices
Analysis of Block LDL T Factorizations for Symmetric Indefinite Matrices Haw-ren Fang August 24, 2007 Abstract We consider the block LDL T factorizations for symmetric indefinite matrices in the form LBL
More informationCMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017
CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,
More informationNumerische Mathematik
Numer. Math. (2001) 88: 319 345 Digital Object Identifier (DOI) 10.1007/s002110000236 Numerische Mathematik Componentwise perturbation analyses for the Q factorization Xiao-Wen Chang, Christopher C. Paige
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematis A NEW ARRANGEMENT INEQUALITY MOHAMMAD JAVAHERI University of Oregon Department of Mathematis Fenton Hall, Eugene, OR 97403. EMail: javaheri@uoregon.edu
More informationFor δa E, this motivates the definition of the Bauer-Skeel condition number ([2], [3], [14], [15])
LAA 278, pp.2-32, 998 STRUCTURED PERTURBATIONS AND SYMMETRIC MATRICES SIEGFRIED M. RUMP Abstract. For a given n by n matrix the ratio between the componentwise distance to the nearest singular matrix and
More informationEffective Resistances for Ladder Like Chains
Effetive Resistanes for Ladder Like Chains Carmona, A., Eninas, A.M., Mitjana, M. May 9, 04 Abstrat Here we onsider a lass of generalized linear hains; that is, the ladder like hains as a perturbation
More informationCONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION. Smail Mahdi
Serdia Math. J. 30 (2004), 55 70 CONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION Smail Mahdi Communiated by N. M. Yanev Abstrat. A two-sided onditional onfidene interval
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationBEHAVIOR OF SQUARE CONCRETE-FILLED TUBULAR COLUMNS UNDER ECCENTRIC COMPRESSION WITH DOUBLE CURVATURE DEFLECTION
Otober 2-7, 28, Beijing, China BEHAVIOR OF SQARE CONCRETE-FILLED TBLAR COLNS NDER ECCENTRIC COPRESSION WITH DOBLE CRVATRE DEFLECTION T. Fujinaga, H. Doi 2 and Y.P. Sun 3 Assoiate Professor, Researh Center
More informationV. Interacting Particles
V. Interating Partiles V.A The Cumulant Expansion The examples studied in the previous setion involve non-interating partiles. It is preisely the lak of interations that renders these problems exatly solvable.
More informationA total positivity property of the Marchenko-Pastur Law
Eletroni Journal of Linear Algebra Volume 30 Volume 30 (205) Artile 7 205 A total positivity property of the Marhenko-Pastur Law Ana Maro Universidad de Alala, Alala de Henares, ana.maro@uah.es Jose-Javier
More informationLecture 3 - Lorentz Transformations
Leture - Lorentz Transformations A Puzzle... Example A ruler is positioned perpendiular to a wall. A stik of length L flies by at speed v. It travels in front of the ruler, so that it obsures part of the
More informationPHYSICS 104B, WINTER 2010 ASSIGNMENT FIVE SOLUTIONS N { N 2 + 1, N 2 + 2, N
PHYSICS 4B, WINTER 2 ASSIGNMENT FIVE SOLUTIONS [.] For uniform m =. and springs g = 2. the eigenvalues have the analyti form, The ode jaobi test yields λ k = 2g ( os k) m k = k n = 2πn N = 2π N { N 2 +,
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume, No 4, 01 Copyright 010 All rights reserved Integrated Publishing servies Researh artile ISSN 0976 4399 Strutural Modelling of Stability
More informationLOGISTIC REGRESSION IN DEPRESSION CLASSIFICATION
LOGISIC REGRESSIO I DEPRESSIO CLASSIFICAIO J. Kual,. V. ran, M. Bareš KSE, FJFI, CVU v Praze PCP, CS, 3LF UK v Praze Abstrat Well nown logisti regression and the other binary response models an be used
More informationEstimating the probability law of the codelength as a function of the approximation error in image compression
Estimating the probability law of the odelength as a funtion of the approximation error in image ompression François Malgouyres Marh 7, 2007 Abstrat After some reolletions on ompression of images using
More informationThe Second Postulate of Euclid and the Hyperbolic Geometry
1 The Seond Postulate of Eulid and the Hyperboli Geometry Yuriy N. Zayko Department of Applied Informatis, Faulty of Publi Administration, Russian Presidential Aademy of National Eonomy and Publi Administration,
More informationModeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers
A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation
More informationVariation Based Online Travel Time Prediction Using Clustered Neural Networks
Variation Based Online Travel Time Predition Using lustered Neural Networks Jie Yu, Gang-Len hang, H.W. Ho and Yue Liu Abstrat-This paper proposes a variation-based online travel time predition approah
More informationRemark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.
Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need
More informationHamacher Sum and Hamacher Product of Fuzzy Matrices
Intern J Fuzzy Maematial Arhive Vol 13, No, 017, 191-198 ISSN: 30 34 (P), 30 350 (online) Published on Deember 017 wwwresearhmasiorg DOI: http://dxdoiorg/10457/fmav13na9 International Journal of amaher
More informationTHE computation of Fourier transforms (see Definition
IAENG International Journal of Applied Mathematis, 48:, IJAM_48 6 The Computation of Fourier Transform via Gauss-Regularized Fourier Series Qingyue Zhang Abstrat In this paper, we onstrut a new regularized
More informationPerforming Two-Way Analysis of Variance Under Variance Heterogeneity
Journal of Modern Applied Statistial Methods Volume Issue Artile 3 5--003 Performing Two-Way Analysis of Variane Under Variane Heterogeneity Sott J. Rihter University of North Carolina at Greensboro, sjriht@ung.edu
More informationSensitivity analysis for linear optimization problem with fuzzy data in the objective function
Sensitivity analysis for linear optimization problem with fuzzy data in the objetive funtion Stephan Dempe, Tatiana Starostina May 5, 2004 Abstrat Linear programming problems with fuzzy oeffiients in the
More informationThe transition between quasi-static and fully dynamic for interfaces
Physia D 198 (24) 136 147 The transition between quasi-stati and fully dynami for interfaes G. Caginalp, H. Merdan Department of Mathematis, University of Pittsburgh, Pittsburgh, PA 1526, USA Reeived 6
More informationAdvances in Radio Science
Advanes in adio Siene 2003) 1: 99 104 Copernius GmbH 2003 Advanes in adio Siene A hybrid method ombining the FDTD and a time domain boundary-integral equation marhing-on-in-time algorithm A Beker and V
More informationAn Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems
An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,
More informationSolving Constrained Lasso and Elastic Net Using
Solving Constrained Lasso and Elasti Net Using ν s Carlos M. Alaı z, Alberto Torres and Jose R. Dorronsoro Universidad Auto noma de Madrid - Departamento de Ingenierı a Informa tia Toma s y Valiente 11,
More information2 The Bayesian Perspective of Distributions Viewed as Information
A PRIMER ON BAYESIAN INFERENCE For the next few assignments, we are going to fous on the Bayesian way of thinking and learn how a Bayesian approahes the problem of statistial modeling and inferene. The
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More information