A NEW GEOMETRIC ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD
|
|
- Harvey Rolf Davis
- 6 years ago
- Views:
Transcription
1 Electronic Transactions on Numerical Analysis. Volume 45, , Coyright c 2016,. ISSN ETNA A NEW GEOMETRIC ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD WILLIAMS LÓPEZ Abstract. We develo a geometrical acceleration scheme for the von Neumann-Halerin alternating rojection method, when alied to the roblem of finding the rojection of a oint onto the intersection of a finite number of closed subsaces of a Hilbert sace. We study the convergence roerties of the new scheme. We also resent some encouraging reliminary numerical results to illustrate the erformance of the new scheme when comared with a well-nown geometrical acceleration scheme, and also with the original von Neumann-Halerin alternating rojection method. Key words. von Neumann-Haerin algorithm, alternating rojection methods, orthogonal rojections, acceleration schemes AMS subject classifications. 52A20, 46C07, 65H10, 47J25 1. Introduction. Let H be a Hilbert sace with inner roduct.,.. For a given 0 H and many closed subsaces M 1,..., M in H, we consider the best aroimation roblem: find the closest oint to 0 in M = i=1 M i, which can be stated as an otimization roblem as follows: 1.1) minimize 0 subject to M, where, for any z H, z 2 = z, z. The unique solution of roblem 1.1) is called the rojection of 0 onto M and is denoted as P M 0 ). In 1933, von Neumann solved roblem 1.1) for the articular case of two closed subsaces. THEOREM 1.1 von Neumann [28]). If M 1 and M 2 are closed subsaces in H, then for each 0 H, lim P M 2 P M1 ) 0 ) = P M1 M 2 0 ). Figure 1.1 shows the geometric interretation of Theorem 1.1. The etension to more than two subsaces was established in 1962 by Halerin. THEOREM 1.2 I. Halerin [23]). If M 1,..., M are closed subsaces in H, then for each 0 H, lim P M P M 1 P M1 ) 0 ) = P i=1 Mi 0). Theorem 1.2 suggests an algorithm, called the method of alternating rojections or MAP for short); see [11, 16], which can be described as follows: for any 0 H, set 1.2) 0 = 1 i = P Mi i 1) i = 1, 2,...,, for Z +, with initial value 0 = 0. The MAP is closely related to Kaczmarz alternating rojection method [25] for solving linear systems of equations. Received October 27, Acceted June 9, Published online on Setember 2, Recommended by Fiorella Sgallari. Deartamento de Física y Matemáticas, Núcleo Universiraio Rafael Rangel, Universidad de Los Andes, Trujillo, Venezuela wloez@ula.ve). 330
2 ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD 331 P M2 P M1 ) 0 ) M 2 0 P M1 M 2 0 ) M 1 P M1 0 ) FIG lim P M 2 P M1 ) 0 ) = P M1 M 2 0 ). THEOREM 1.3. For any i = 1, 2,...,, the sequence { i } generated by 1.2) converges to P M 0 ). Proof. Let i {1, 2,..., }. Then i = P Mi P Mi 1 P M1 ) ), for all 0, where = P M P M 1 P M1 ) 0 ). From Theorem 1.2 we have that P M 0 ) when and since P Mi, P Mi 1,..., P M1 are bounded linear oerators, then i = P Mi P Mi 1 P M1 ) ) P Mi P Mi 1 P M1 )P M 0 )) = P M 0 ) when. Alications of MAP in various areas of mathematics can be found in [16, 9, 10, 3, 7, 8, 22]. The MAP has an r-linear rate of convergence that can be very slow when the angles between the subsaces in the sense of Friedrichs [20]) are small; see, e.g., [16]. Franchetti and Light [19] showed that the convergence in Theorem 1.1 may be arbitrarily slow if M 1 + M 2 is not closed i.e., if the angle between M 1 and M 2 is zero). Consequently, several acceleration schemes have been roosed; see, e.g., [4, 21, 24, 12]. Motivated by a recent wor of Lóez and Raydan [26], in this wor we resent and analyze a geometrical acceleration scheme for MAP in 1.2), which solves roblem 1.1). The new scheme involves finding oints closest to a solution, ) in the roduct sace H H, where is the solution of 1.1). The rest of this aer is organized as follows. In Section 2 we rovide information about already eisting acceleration schemes for MAP. In Section 2.1 we develo the new acceleration scheme and discuss its theoretical roerties. In Section 3 we resent encouraging reliminary numerical results. 2. Accelerations for MAP. The sequences generated by the method of alternating rojections often converge very slowly. Such a slow convergence is observed, e.g., if M 1 and M 2 are two closed subsaces with angles close to 0 see, e.g., [19]). Since the method of alternating rojections has many ractical alications, any acceleration technique seems to be imortant. Several acceleration schemes with a geometrical flavor have been roosed to imrove the erformance of MAP; see, e.g., De Pierro and Iusem [13], Dos Santos [14], Gearhart and Koshy [21], Bausche et al. [4], Gubin, Polya and Rai [22], Martínez [27], Aleby and Smolarsi [1] and Censor [6]. In Section 2.1, we describe a geometrically aealing acceleration scheme following the resentation in [21]. Some other different acceleration ideas have also been develoed see, e.g., Echebest et al. [17, 18] and Scolni et al. [29, 30]) based on the so-called rojected aggregation method PAM). Other secialized acceleration scheme
3 332 W. LÓPEZ ideas have been develoed by Hernández-Ramos et. al. [24] based on the use of the conjugate gradient method for minimizing a related conve quadratic ma. Usually vector etraolation algorithms are used for accelerating the convergence of MAP; see, e.g., [3, 8, 7]. Dystra s algorithm [15, 5] solves roblem 1.1) when the involved sets are closed and conve not necessarily closed subsaces). Dystra s algorithm can be viewed as a natural etension to conve sets of von Neumann-Halerin s MAP for subsaces in a Hilbert sace. The sequences generated by Dystra s algorithm often converge very slowly. Such a slow convergence is observed, e.g., if M 1 is a hyerlane and M 2 is a ball which is tangent to M 1 ; see, e.g., [2, Eamle 5.3]. An acceleration scheme has been develoed for Dystra s algorithm [26]. In Section 2.1 we develo a new acceleration scheme for von Neumann-Halerin s MAP on closed subsaces, which is related to the scheme given in [26] A new acceleration scheme. We now discuss an acceleration scheme for MAP to solve roblem 1.1). For that we need to consider an auiliary sequence in the roduct sace H H that will be denoted as H 2. If H is a Hilbert sace with inner roduct.,., we will use the inner roduct in H 2 defined by, y), w, z) =, w + y, z, for all, y), w, z) H 2 and norm, y) 2 =, y),, y) = 2 + y 2, for all, y) H 2. Thus H 2 is a Hilbert sace. For each Z +, let us define 2.1) = 1, ) H 2, where 1 and are defined in 1.2). Let us also denote = P M 0 ), P M 0 )) H 2. It follows that { } is a sequence in H 2 that converges to. Our idea to accelerate the convergence of MAP consists in building another sequence in H 2 which converges faster to than { }. For that we need to define a suitable subsace Π of H 2 that contains, as follows: Π = {, ) : H}. Clearly, Π is a closed subsace in H 2 and Π. Let us denote by P Π ) the orthogonal rojection of H 2 onto Π. For each Z +, let us consider P Π ) Π, where { } is given by 2.1). THEOREM 2.1. For all 1, P Π ). Proof. For all 1, since P Π ) is orthogonal to Π, and since P Π ) and belong to Π, then P Π ) is orthogonal to P Π ). Hence 2 = P Π ) + P Π ) 2 = P Π ) 2 + P Π ) 2. Consequently, 2 P Π ) 2, thus P Π ). The sequence {P Π )} will lay a ey role in the develoment of the acceleration scheme for MAP. REMARK 2.2. If, y) H 2, then see, e.g., [26]) 2.2) P Π, y) = 1/2 + y, + y). Notice that comuting the rojection onto Π only requires us to comute the average of the two involved vectors, i.e., it requires a very ineensive calculation. Let { } be the sequence defined by 2.1). For any 1, with P Π ) P Π ), let ô Π be the rojection of onto the line that goes through P Π ) and P Π ) see
4 ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD 333 Figure 2.1). Then ô, P Π ) P Π ) = 0 and there eists a unique α R such that Therefore ô = P Π ) + α P Π ) P Π )). 2.3) P Π ) + α P Π ) P Π )), P Π ) P Π ) = 0. Solving for α, from 2.3), gives 2.4) α = P Π ), P Π ) P Π ) P Π ) P Π ) 2. In the following result we will give a formula for α in 2.4) that does not require nowledge of. Π ô P Π ) P Π ) FIG ô, P Π ) P Π ) = 0. THEOREM 2.3. Let { } be defined by 2.1). Then, P Π ) P Π ) = 0 for all 1. Proof. Let 1. Since Π and the rojection P Π is self-adjoint, we obtain, P Π ) P Π ) =, =, ), 1 1, ) =, 1 1 +,. On the other hand, since M = i=1 M i and the rojections P Mi are self-adjoint, we have, =,, =,, P M P... P M 1 M 1 ) ) =,, = 0. Similarly, we also obtain, 1 1 = 0. COROLLARY 2.4. Let α, for some, be given by 2.4). Then 2.5) α = P Π ), P Π ) P Π ) P Π ) P Π ) 2. Proof. This is a consequence of Theorem 2.3 and 2.4).
5 334 W. LÓPEZ Let us define the sequence {ô } in Π, for 1, as follow { P Π ) if P Π ) P Π ) = 0, 2.6) ô = P Π ) + α P Π ) P Π )) if P Π ) P Π ) 0, where α is given by 2.5). THEOREM 2.5. Let {ô } be given by 2.6). Then ô for all 1. Proof. Let 1. If P Π ) P Π ) = 0, then from Theorem 2.1 and 2.6) it follows that ô = P Π ). On the other hand, if P Π ) P Π ) 0, then ô is chosen to minimize the distance to along the line connecting P Π ) and P Π ). Therefore ô P Π ) and, from Theorem 2.1, it follows that ô. Theorem 2.5 guarantees than the sequence {ô } defined by 2.6) accelerates the convergence of the sequence { } defined by 2.1). REMARK 2.6. For the articular case of closed subsaces, the sequence {ô } defined by 2.6) accelerates the convergence of the sequence {ẑ } defined in [26]. For any 1, ẑ defined in [26] is the intersection of the subsace Π with the line that goes through and see Figure 2.2). Consequently ẑ = + α ) Π, for some α R. Since ẑ Π and P Π is a linear oerator, ẑ = P Π ẑ ) = P Π ) + α P Π ) P Π )). Therefore we have that ẑ belongs to the line connecting P Π ) and P Π ). Since ô is chosen to minimize the distance to along the line connecting P Π ) and P Π ), ô ẑ see Figure 2.2). Π P ẑ Π ) ô P Π ) FIG ô vs. ẑ. REMARK 2.7. Since the sequence {ô } defined by 2.6) belongs to Π, for any 1, there eists a unique o H such that 2.7) ô = o, o ) Π.
6 ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD 335 THEOREM 2.8. Let {o } the sequence defined by 2.7). Then, for each 1, o P M 0 ) ma{ 1 P M 0 ), P M 0 ) }. Proof. Let 1. From Theorem 2.5, we have that 2 o P M 0 ) = o P M 0 ) 2 + o P M 0 ) 2 = ô = 1 P M 0 ) 2 + P M 0 ) 2 2 ma{ 1 P M 0 ) 2, P M 0 ) 2 } = 2 ma{ 1 P M 0 ), P M 0 ) }, and the result is established. COROLLARY 2.9. The sequence {o } defined by 2.7) converges to P M 0 ). Proof. This is a consequence of Theorem 2.8 and Theorem 1.3. REMARK Observe that the sequence {ô } in H 2 accelerates the convergence of the sequence { } in H 2 see Theorem 2.5). However, Theorem 2.8 does not guarantee acceleration when the sequence {o } in H is comared to the original MAP sequence { } in H, since a maimal is involved. Nevertheless, based on numerical eerimentation, the sequence {o } almost always accelerates the convergence of the original MAP sequence { }. We will now describe elicitly the sequence {o } defined by 2.7). If P Π ) P Π ) = 0, then from 2.1) and 2.2) we have that ô = P Π ) = P Π 1, In this case o = 1/2 1 + then where ) = 1/2 ) = o, o ) Π. 1 +, 1 + ) H. If, on the other hand, PΠ ) P Π ) 0, ô = P Π ) + α P Π ) P Π )) = o, o ) Π, o = 1/2 + α From Corollary 2.9 it follows that o = 1/2 + α ) α 1 1)) H. ) α 1 1)) PM 0 ) when. Since P M is a bounded linear oerator where P M is the rojection oerator onto M ), [ P M 1/2 + α ) α 1 1))] PM P M 0 )) when. Consequently 1/2 [ + α ) + + α )] PM 0 )
7 336 W. LÓPEZ when. Therefore 2.8) + α ) PM 0 ) when. Equation 2.8) suggests how to define an accelerated sequence in H and a secialized algorithm for which convergence to the solution P M 0 ) will be later established. Nevertheless, to achieve this it is convenient to write α given by 2.5) as a function of the original MAP iterations in H. LEMMA Let α, for some, be given by 2.5). Then 2.9) α = 1 +, , Proof. Let 1. The numerator in 2.5) can be written as P Π ), P Π ) P Π ) = P Π ), = 1/2 ) 1 +, 1 +, 1 1, ) = 1/2 1 +, 1 ) ,. On the other hand P Π ) P Π ) = 1/2 1 +, 1 + ) 1/2 1 +, 1 + ) = 1/ , ). Therefore the denominator in 2.5) can be written as P Π ) P Π ) 2 = 1/ ) = 1/ , and the result holds. Let { i } be the MAP s iterates given by 1.2). Let us now define a new sequence {o } in H, for 1, as follows { 2.10) o = if = 0, + α ) if , where α is given by 2.9). THEOREM The sequence {o } defined by 2.10) converges to P M 0 ). Proof. This is a consequence of Corollary 2.9, equation 2.8), and Theorem 1.3. We are now ready to resent our acceleration scheme for MAP. Algorithm 1. Let M i, i = 1,...,, be closed subsaces in H. Given 0 H; set = 1. set 0 = 0 1 = P M 1... P M1 ) 1 ) = P M 1) for = 1, 2,... do 1 = P M 1... P M1 ) ) = P M 1 )
8 ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD 337 if then comute α using 2.9), and set o = + α ) else set o = end if end for In the following, we will comare our scheme, given in Algorithm 1, with the acceleration scheme given by Gearhart and Koshy [21]. Therefore, we will give a brief elanation of the Gearhart-Koshy scheme. Let us denote by 0 the given starting oint and by Q the comosition of the rojection oerators, i.e, Q = P M P M 1 P M1, where P Mi is the rojection oerator onto M i for all i. Let be the th iterate, and let Q be the net iterate after alying a swee of MAP. The idea is to search along the line through the oints and Q to obtain the oint closest to the solution = P i=1 Mi 0). Let us reresent any oint on this line as t) = tq + 1 t) = + tq ), for some real number t. Let t be the value of t for which this oint is the closest to. Then, 2.11) t ), Q = 0. Now, since i=1 M i and the rojections P Mi are self-adjoint, it follows, Q = P M1 P M2 P M, =,. Consequently,, Q = 0, and so can be eliminated from 2.11) to obtain t ), Q = 0. Solving for t gives To summarize, we have the following stes. t =, Q Q 2. Algorithm 2. Gearhart-Koshy) Let M i, i = 1,...,, be closed subsaces in H. Given 0 H. for = 1, 2,... do Q = Q t =, Q / Q 2 = t Q + 1 t ) end for 3. Numerical eeriments. We comare our acceleration scheme, given in Algorithm 1, with the acceleration scheme given by Gearhart and Koshy Algorithm 2), and with the original MAP with no acceleration given by 1.2). All comutations were erformed in MATLAB. For all eeriments we now the eact solution, and we sto each algorithm when the absolute error the distance from the -th iterate to ) is less than or equal to 10 6.
9 338 W. LÓPEZ For our first eeriment, we consider the following three subsaces of the sace of square real matrices R 3 3, with the Frobenius norm A 2 F = A, A = traceat A): M 1 = { A R 3 3 : A T = A }, M 2 = { A R 3 3 : a i,j+1 = 0, i = 1, 2, j = i + 1,..., 3 }, M 3 = { A R 3 3 : a 1,1 = a 1,3 = a 3,1 = a 3,3 }. If A = A ij ) R 3 3, then P M1 A) = A T + A)/2, A P M2 A) = A 21 A 22 0, and P M3 A) = A 31 A 32 A 33 where P = A 11 + A 13 + A 31 + A 33 )/4. We choose as the initial given oint. Then since A 0 = P M1 M 2 M 3 A 0 ) = P A 12 P A 21 A 22 A 23 P A 32 P, A : R = M 1 M 2 M 3, A, A 0 P M1 M 2 M 3 A 0 ) = trace The erformance of each method is shown in Figure 3.1. T , = 0. Consider the roblem of solving the linear system of equations A = 0, where A R m n and R n. This roblem can be generalized to any Hilbert sace H by considering the following formulation: find in the intersection of m closed subsaces given by H i = { H : a i, = 0}, for every i = 1,..., m. Here a i denotes the ith row of A or, in general, a fied given vector in H. If z H, then P Hi z) = z a i, z / a i, a i )a i. In the following two test roblems we will solve the linear equations A = 0, where A is a square nonsingular matri which is randomly generated in MATLAB we verified that deta) 0 after randomly generating the matri A in MATLAB). Therefore the redetermined solutions is always = 0. In both eeriments we chose 0 = 100, 100,..., 100) as the initial given oint. Figure 3.2 shows the erformance for each methods for A R 5 5 while Figure 3.3 shows the erformance for each method for A R Our reliminary numerical eeriments seem to indicate that our acceleration scheme, when comared to the other two cometitors, benefits from a reduction in the number of cycles, as well as in the comutational wor.
10 ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD 339 FIG Acceleration for three subsaces in R 3 3. FIG Acceleration for solving A = 0, where A R 5 5. Acnowledgements. I than the three anonymous referees for constructive remars and additional references. REFERENCES [1] G. APPLEBY AND D. C. SMOLARSKI, A linear acceleration row action method for rojecting onto subsaces, Electron. Trans. Numer. Anal, ), htt://etna.ricam.oeaw.ac.at/vol / dir/ df [2] H. H. BAUSCHKE AND J. M. BORWEIN, Dystra s alternating rojection algorithm for two sets, J. Aro.
11 340 W. LÓPEZ FIG Acceleration for solving A = 0, where A R Theory, ), [3] H. H. BAUSCHKE, P. L. COMBETTES, AND S. G. KRUK, Etraolation algorithm for affine-conve feasibility roblems, Numer. Algorithms, ), [4] H. H. BAUSCHKE, F. DEUTSCH, H. HUNDAL, AND S. H. PARK, Accelerating the convergence of the method of alternating rojections, Trans. Amer. Math. Soc., ), [5] J. P. BOYLE AND L. DYKSTRA, A method for finding rojections onto the intersection of conve sets in Hilbert saces, in Advances in Order Restricted Statistical Inference, R. Dystra, T. Robertson, and F. T. Wright, eds., Lecture Notes in Statistics, 37, Sringer, Berlin, 1986, [6] Y. CENSOR, Comutational acceleration of rojection algorithms for the linear best aroimation roblem, Linear Algebra Al., ), [7] P. L. COMBETTES, Conve set theoretic image recovery by etraolated iterations of arallel subgradient rojections, IEEE Trans. Image Process., ), [8], The conve feasibility roblem in image recovery, in Advances in Imaging and Electron Physics, 95 P. W. Hawes, ed., Academic Press, New Yor, 1996, [9] F. DEUTSCH, Alications of von Neumann s alternating rojections algorithm, in Mathematical Methods in Oerations Research, P. Kenderov, ed., Bulgarian Academy of Science, Sofia, 1983, [10], The method of alternating orthogonal rojections, in Aroimation Theory, Sline Functions and Alications, S. P. Singh, ed., Kluwer, Netherlands, 1992, [11], Best Aroimation in Inner Product Sace, Sringer, New Yor, [12], Accelerating the convergence of the method of alternating rojections via a line search: a brief survey, in Inherently Parallel Algorithms in Feasibility and Otimization and their Alications, D. Butnariu, Y. Censor, S. Reich, eds., Studies in Comutational Mathematics, 8, North-Holland Publishing, Amsterdam, 2001, [13] A. DE PIERRO AND A. N. IUSEM, A arallel rojection method for finding a common oint of a family of conve sets, Pesquisa Oeracional, ), [14] L. T. DOS SANTOS, A arallel subgradient rojections method for the conve feasibility roblem, J. Comut. Al. Math., ), [15] R. L. DYKSTRA, An algorithm for restricted least-squares regression, J. Amer. Statist. Assoc., ), [16] R. ESCALANTE AND M. RAYDAN, Alternating Projection Methods, SIAM, Philadelhia, [17] N. ECHEBEST, M. T. GUARDARUCCI, H. SCOLNIK, AND M. C. VACCHINO, An acceleration scheme for solving conve feasibility roblems using incomlete rojection algorithms, Numer. Algorithms, ), [18], An accelerated iterative method with diagonally scaled oblique rojections for solving linear feasibility roblems, Ann. Oer. Res., ), [19] C. FRANCHETTI AND W. LIGHT, On the von Neumann alternating algorithm in Hilbert sace, J. Math. Anal.
12 ACCELERATION OF THE VON NEUMANN-HALPERIN PROJECTION METHOD 341 Al., ), [20] K. FRIEDRICHS, On certain inequalities and characteristic value roblems for analytic functions and for functions of two variables, Trans. Amer. Math. Soc., ), [21] W. GEARHART AND M. KOSHY, Acceleration schemes for the method of alternating rojection, J. Comut. Al. Math., ), [22] L. G. GUBIN, B. T. POLYAK, AND E. V. RAIK, The method of rojections for finding the common oint of conve sets, USSR Comut. Math. Math. Phys., ), [23] I. HALPERIN, The roduct of rojection oerators, Acta Sci. Math. Szeged), ), [24] L. M. HERNÁNDEZ-RAMOS, R. ESCALANTE, AND M. RAYDAN, Unconstrained otimization techniques for the acceleration of alternating rojection methods, Numer. Funct. Anal. Otim., ), [25] S. KACZMARZ, Angenäherte Auflösung von Systemen linearer Gleichungen, Bull. Internat. Acad. Pol. Sci. Lett. Ser. A, ), [26] W. LÓPEZ AND M. RAYDAN, An acceleration scheme for Dystra s algorithm, Comut. Otim. Al., ), [27] J. M. MARTÍNEZ, An accelerated successive orthogonal rojection method for solving large-scale linear feasibility roblems, Comut. Math. Al., ), [28] J. VON NEUMANN, Functional oerators. Vol II. The Geometry of Orthogonal Saces, Annals of Math. Studies, vol. 22, Princeton University Press, Princeton, [29] H. SCOLNIK, N. ECHEBEST, M. T. GUARDARUCCI, AND M. C. VACCHINO, Acceleration scheme for arallel rojected aggregation methods for solving large linear systems, Ann. Oer. Res., ), [30], A class of otimized row rojection methods for solving large nonsymmetric linear systems, Al. Numer. Math., ),
TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES
Khayyam J. Math. DOI:10.22034/kjm.2019.84207 TRACES OF SCHUR AND KRONECKER PRODUCTS FOR BLOCK MATRICES ISMAEL GARCÍA-BAYONA Communicated by A.M. Peralta Abstract. In this aer, we define two new Schur and
More informationNONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS. The Goldstein-Levitin-Polyak algorithm
- (23) NLP - NONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS The Goldstein-Levitin-Polya algorithm We consider an algorithm for solving the otimization roblem under convex constraints. Although the convexity
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationRIEMANN-STIELTJES OPERATORS BETWEEN WEIGHTED BERGMAN SPACES
RIEMANN-STIELTJES OPERATORS BETWEEN WEIGHTED BERGMAN SPACES JIE XIAO This aer is dedicated to the memory of Nikolaos Danikas 1947-2004) Abstract. This note comletely describes the bounded or comact Riemann-
More informationConvex Analysis and Economic Theory Winter 2018
Division of the Humanities and Social Sciences Ec 181 KC Border Conve Analysis and Economic Theory Winter 2018 Toic 16: Fenchel conjugates 16.1 Conjugate functions Recall from Proosition 14.1.1 that is
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationThe Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment
he Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment William Glunt 1, homas L. Hayden 2 and Robert Reams 2 1 Department of Mathematics and Computer Science, Austin Peay State
More informationA Dykstra-like algorithm for two monotone operators
A Dykstra-like algorithm for two monotone operators Heinz H. Bauschke and Patrick L. Combettes Abstract Dykstra s algorithm employs the projectors onto two closed convex sets in a Hilbert space to construct
More informationL p Norms and the Sinc Function
L Norms and the Sinc Function D. Borwein, J. M. Borwein and I. E. Leonard March 6, 9 Introduction The sinc function is a real valued function defined on the real line R by the following eression: sin if
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. III Stability Theory - Peter C. Müller
STABILITY THEORY Peter C. Müller University of Wuertal, Germany Keywords: Asymtotic stability, Eonential stability, Linearization, Linear systems, Lyaunov equation, Lyaunov function, Lyaunov stability,
More informationThe Method of Alternating Projections
Matthew Tam Variational Analysis Session 56th Annual AustMS Meeting 24th 27th September 2012 My Year So Far... Closed Subspaces Honours student supervised by Jon Borwein. Thesis topic: alternating projections.
More informationGENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS
International Journal of Analysis Alications ISSN 9-8639 Volume 5, Number (04), -9 htt://www.etamaths.com GENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS ILYAS ALI, HU YANG, ABDUL SHAKOOR Abstract.
More informationSpecialized and hybrid Newton schemes for the matrix pth root
Secialized and hybrid Newton schemes for the matrix th root Braulio De Abreu Marlliny Monsalve Marcos Raydan June 15, 2006 Abstract We discuss different variants of Newton s method for comuting the th
More informationChapter 10. Supplemental Text Material
Chater 1. Sulemental Tet Material S1-1. The Covariance Matri of the Regression Coefficients In Section 1-3 of the tetbook, we show that the least squares estimator of β in the linear regression model y=
More informationEfficient algorithms for the smallest enclosing ball problem
Efficient algorithms for the smallest enclosing ball roblem Guanglu Zhou, Kim-Chuan Toh, Jie Sun November 27, 2002; Revised August 4, 2003 Abstract. Consider the roblem of comuting the smallest enclosing
More informationE-companion to A risk- and ambiguity-averse extension of the max-min newsvendor order formula
e-comanion to Han Du and Zuluaga: Etension of Scarf s ma-min order formula ec E-comanion to A risk- and ambiguity-averse etension of the ma-min newsvendor order formula Qiaoming Han School of Mathematics
More informationApplicable Analysis and Discrete Mathematics available online at HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS
Alicable Analysis and Discrete Mathematics available online at htt://efmath.etf.rs Al. Anal. Discrete Math. 4 (010), 3 44. doi:10.98/aadm1000009m HENSEL CODES OF SQUARE ROOTS OF P-ADIC NUMBERS Zerzaihi
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 9,. 29-36, 25. Coyright 25,. ISSN 68-963. ETNA ASYMPTOTICS FOR EXTREMAL POLYNOMIALS WITH VARYING MEASURES M. BELLO HERNÁNDEZ AND J. MíNGUEZ CENICEROS
More informationOn Erdős and Sárközy s sequences with Property P
Monatsh Math 017 18:565 575 DOI 10.1007/s00605-016-0995-9 On Erdős and Sárközy s sequences with Proerty P Christian Elsholtz 1 Stefan Planitzer 1 Received: 7 November 015 / Acceted: 7 October 016 / Published
More informationA viability result for second-order differential inclusions
Electronic Journal of Differential Equations Vol. 00(00) No. 76. 1 1. ISSN: 107-6691. URL: htt://ejde.math.swt.edu or htt://ejde.math.unt.edu ft ejde.math.swt.edu (login: ft) A viability result for second-order
More informationImproving AOR Method for a Class of Two-by-Two Linear Systems
Alied Mathematics 2 2 236-24 doi:4236/am22226 Published Online February 2 (htt://scirporg/journal/am) Imroving AOR Method for a Class of To-by-To Linear Systems Abstract Cuixia Li Shiliang Wu 2 College
More information#A45 INTEGERS 12 (2012) SUPERCONGRUENCES FOR A TRUNCATED HYPERGEOMETRIC SERIES
#A45 INTEGERS 2 (202) SUPERCONGRUENCES FOR A TRUNCATED HYPERGEOMETRIC SERIES Roberto Tauraso Diartimento di Matematica, Università di Roma Tor Vergata, Italy tauraso@mat.uniroma2.it Received: /7/, Acceted:
More informationGRACEFUL NUMBERS. KIRAN R. BHUTANI and ALEXANDER B. LEVIN. Received 14 May 2001
IJMMS 29:8 2002 495 499 PII S06720200765 htt://immshindawicom Hindawi Publishing Cor GRACEFUL NUMBERS KIRAN R BHUTANI and ALEXANDER B LEVIN Received 4 May 200 We construct a labeled grah Dn that reflects
More informationLecture 21: Quantum Communication
CS 880: Quantum Information Processing 0/6/00 Lecture : Quantum Communication Instructor: Dieter van Melkebeek Scribe: Mark Wellons Last lecture, we introduced the EPR airs which we will use in this lecture
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationON FREIMAN S 2.4-THEOREM
ON FREIMAN S 2.4-THEOREM ØYSTEIN J. RØDSETH Abstract. Gregory Freiman s celebrated 2.4-Theorem says that if A is a set of residue classes modulo a rime satisfying 2A 2.4 A 3 and A < /35, then A is contained
More informationMultiplicity of weak solutions for a class of nonuniformly elliptic equations of p-laplacian type
Nonlinear Analysis 7 29 536 546 www.elsevier.com/locate/na Multilicity of weak solutions for a class of nonuniformly ellitic equations of -Lalacian tye Hoang Quoc Toan, Quô c-anh Ngô Deartment of Mathematics,
More informationIMPROVED BOUNDS IN THE SCALED ENFLO TYPE INEQUALITY FOR BANACH SPACES
IMPROVED BOUNDS IN THE SCALED ENFLO TYPE INEQUALITY FOR BANACH SPACES OHAD GILADI AND ASSAF NAOR Abstract. It is shown that if (, ) is a Banach sace with Rademacher tye 1 then for every n N there exists
More informationarxiv:math/ v1 [math.fa] 5 Dec 2003
arxiv:math/0323v [math.fa] 5 Dec 2003 WEAK CLUSTER POINTS OF A SEQUENCE AND COVERINGS BY CYLINDERS VLADIMIR KADETS Abstract. Let H be a Hilbert sace. Using Ball s solution of the comlex lank roblem we
More informationProducts of Composition, Multiplication and Differentiation between Hardy Spaces and Weighted Growth Spaces of the Upper-Half Plane
Global Journal of Pure and Alied Mathematics. ISSN 0973-768 Volume 3, Number 9 (207),. 6303-636 Research India Publications htt://www.riublication.com Products of Comosition, Multilication and Differentiation
More informationA CLASS OF ALGEBRAIC-EXPONENTIAL CONGRUENCES MODULO p. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXXI, (2002),. 3 7 3 A CLASS OF ALGEBRAIC-EXPONENTIAL CONGRUENCES MODULO C. COBELI, M. VÂJÂITU and A. ZAHARESCU Abstract. Let be a rime number, J a set of consecutive integers,
More informationResearch Article An iterative Algorithm for Hemicontractive Mappings in Banach Spaces
Abstract and Alied Analysis Volume 2012, Article ID 264103, 11 ages doi:10.1155/2012/264103 Research Article An iterative Algorithm for Hemicontractive Maings in Banach Saces Youli Yu, 1 Zhitao Wu, 2 and
More information2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.
-D Analysis for Iterative Learning Controller for Discrete-ime Systems With Variable Initial Conditions Yong FANG, and ommy W. S. Chow Abstract In this aer, an iterative learning controller alying to linear
More informationR.A. Wibawa-Kusumah and H. Gunawan 1
TWO EQUIVALENT n-norms ON THE SPACE OF -SUMMABLE SEQUENCES RA Wibawa-Kusumah and H Gunawan 1 Abstract We rove the (strong) equivalence between two known n-norms on the sace l of -summable sequences (of
More informationFINDING BEST APPROXIMATION PAIRS RELATIVE TO TWO CLOSED CONVEX SETS IN HILBERT SPACES
FINDING BEST APPROXIMATION PAIRS RELATIVE TO TWO CLOSED CONVEX SETS IN HILBERT SPACES Heinz H. Bauschke, Patrick L. Combettes, and D. Russell Luke Abstract We consider the problem of finding a best approximation
More informationOn Total Convexity, Bregman Projections and Stability in Banach Spaces
Journal of Convex Analysis Volume 11 (2004), No. 1, 1 16 On Total Convexity, Bregman Projections and Stability in Banach Spaces Elena Resmerita Department of Mathematics, University of Haifa, 31905 Haifa,
More informationVariational inequalities for fixed point problems of quasi-nonexpansive operators 1. Rafał Zalas 2
University of Zielona Góra Faculty of Mathematics, Computer Science and Econometrics Summary of the Ph.D. thesis Variational inequalities for fixed point problems of quasi-nonexpansive operators 1 by Rafał
More informationAN EIGENVALUE STUDY ON THE SUFFICIENT DESCENT PROPERTY OF A MODIFIED POLAK-RIBIÈRE-POLYAK CONJUGATE GRADIENT METHOD S.
Bull. Iranian Math. Soc. Vol. 40 (2014), No. 1, pp. 235 242 Online ISSN: 1735-8515 AN EIGENVALUE STUDY ON THE SUFFICIENT DESCENT PROPERTY OF A MODIFIED POLAK-RIBIÈRE-POLYAK CONJUGATE GRADIENT METHOD S.
More informationWAVELETS, PROPERTIES OF THE SCALAR FUNCTIONS
U.P.B. Sci. Bull. Series A, Vol. 68, No. 4, 006 WAVELETS, PROPERTIES OF THE SCALAR FUNCTIONS C. PANĂ * Pentru a contrui o undină convenabilă Ψ este necesară şi suficientă o analiză multirezoluţie. Analiza
More informationWEAK CONVERGENCE OF RESOLVENTS OF MAXIMAL MONOTONE OPERATORS AND MOSCO CONVERGENCE
Fixed Point Theory, Volume 6, No. 1, 2005, 59-69 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.htm WEAK CONVERGENCE OF RESOLVENTS OF MAXIMAL MONOTONE OPERATORS AND MOSCO CONVERGENCE YASUNORI KIMURA Department
More informationGalois representations on torsion points of elliptic curves NATO ASI 2014 Arithmetic of Hyperelliptic Curves and Cryptography
Galois reresentations on torsion oints of ellitic curves NATO ASI 04 Arithmetic of Hyerellitic Curves and Crytograhy Francesco Paalardi Ohrid, August 5 - Setember 5, 04 Lecture - Introduction Let /Q be
More informationChapter 10. Classical Fourier Series
Math 344, Male ab Manual Chater : Classical Fourier Series Real and Comle Chater. Classical Fourier Series Fourier Series in PS K, Classical Fourier Series in PS K, are aroimations obtained using orthogonal
More informationReferences 1. Aleyner, A., & Reich, S. (2008). Block-iterative algorithms for solving convex feasibility problems in Hilbert and Banach spaces. Journa
References 1. Aleyner, A., & Reich, S. (2008). Block-iterative algorithms for solving convex feasibility problems in Hilbert and Banach spaces. Journal of Mathematical Analysis and Applications, 343, 427
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationSTRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2
STRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2 ANGELES ALFONSECA Abstract In this aer we rove an almost-orthogonality rincile for
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationInexact Spectral Projected Gradient Methods on Convex Sets
Inexact Spectral Projected Gradient Methods on Convex Sets Ernesto G. Birgin José Mario Martínez Marcos Raydan March 26, 2003 Abstract A new method is introduced for large scale convex constrained optimization.
More informationGradient Based Iterative Algorithms for Solving a Class of Matrix Equations
1216 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 50, NO. 8, AUGUST 2005 Exanding (A3)gives M 0 8 T +1 T +18 =( T +1 +1 ) F = I (A4) g = 0F 8 T +1 =( T +1 +1 ) f = T +1 +1 + T +18 F 8 T +1 =( T +1 +1 )
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationLINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL
LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL Mohammad Bozorg Deatment of Mechanical Engineering University of Yazd P. O. Box 89195-741 Yazd Iran Fax: +98-351-750110
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationA Projection Algorithm for the Quasiconvex Feasibility Problem 1
International Mathematical Forum, 3, 2008, no. 28, 1375-1381 A Projection Algorithm for the Quasiconvex Feasibility Problem 1 Li Li and Yan Gao School of Management, University of Shanghai for Science
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationarxiv: v1 [quant-ph] 20 Jun 2017
A Direct Couling Coherent Quantum Observer for an Oscillatory Quantum Plant Ian R Petersen arxiv:76648v quant-h Jun 7 Abstract A direct couling coherent observer is constructed for a linear quantum lant
More informationCONGRUENCES CONCERNING LUCAS SEQUENCES ZHI-HONG SUN
Int. J. Number Theory 004, no., 79-85. CONGRUENCES CONCERNING LUCAS SEQUENCES ZHI-HONG SUN School of Mathematical Sciences Huaiyin Normal University Huaian, Jiangsu 00, P.R. China zhihongsun@yahoo.com
More informationA Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression
Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi
More informationConvergence rate examples and theory
Convergence rate examples and theory Simen Gaure Abstract. If you use lfe for various tasks, you will notice that some estimations converge fast, whereas others converge slowly. Convergence rate of the
More informationAlaa Kamal and Taha Ibrahim Yassen
Korean J. Math. 26 2018), No. 1,. 87 101 htts://doi.org/10.11568/kjm.2018.26.1.87 ON HYPERHOLOMORPHIC Fω,G α, q, s) SPACES OF QUATERNION VALUED FUNCTIONS Alaa Kamal and Taha Ibrahim Yassen Abstract. The
More information1 1 c (a) 1 (b) 1 Figure 1: (a) First ath followed by salesman in the stris method. (b) Alternative ath. 4. D = distance travelled closing the loo. Th
18.415/6.854 Advanced Algorithms ovember 7, 1996 Euclidean TSP (art I) Lecturer: Michel X. Goemans MIT These notes are based on scribe notes by Marios Paaefthymiou and Mike Klugerman. 1 Euclidean TSP Consider
More informationUnsupervised Hyperspectral Image Analysis Using Independent Component Analysis (ICA)
Unsuervised Hyersectral Image Analysis Using Indeendent Comonent Analysis (ICA) Shao-Shan Chiang Chein-I Chang Irving W. Ginsberg Remote Sensing Signal and Image Processing Laboratory Deartment of Comuter
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationDEPARTMENT OF ECONOMICS ISSN DISCUSSION PAPER 20/07 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES
DEPARTMENT OF ECONOMICS ISSN 1441-549 DISCUSSION PAPER /7 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES ZuXiang Wang * & Russell Smyth ABSTRACT We resent two new Lorenz curve families by using the basic
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationConvolution Properties for Certain Meromorphically Multivalent Functions
Filomat 31:1 (2017), 113 123 DOI 10.2298/FIL1701113L Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: htt://www.mf.ni.ac.rs/filomat Convolution Proerties for Certain
More informationLEIBNIZ SEMINORMS IN PROBABILITY SPACES
LEIBNIZ SEMINORMS IN PROBABILITY SPACES ÁDÁM BESENYEI AND ZOLTÁN LÉKA Abstract. In this aer we study the (strong) Leibniz roerty of centered moments of bounded random variables. We shall answer a question
More informationA Certain Subclass of Multivalent Analytic Functions Defined by Fractional Calculus Operator
British Journal of Mathematics & Comuter Science 4(3): 43-45 4 SCIENCEDOMAIN international www.sciencedomain.org A Certain Subclass of Multivalent Analytic Functions Defined by Fractional Calculus Oerator
More informationAnalysis of the convergence rate for the cyclic projection algorithm applied to semi-algebraic convex sets
Analysis of the convergence rate for the cyclic projection algorithm applied to semi-algebraic convex sets Liangjin Yao University of Newcastle June 3rd, 2013 The University of Newcastle liangjin.yao@newcastle.edu.au
More informationL p -CONVERGENCE OF THE LAPLACE BELTRAMI EIGENFUNCTION EXPANSIONS
L -CONVERGENCE OF THE LAPLACE BELTRAI EIGENFUNCTION EXPANSIONS ATSUSHI KANAZAWA Abstract. We rovide a simle sufficient condition for the L - convergence of the Lalace Beltrami eigenfunction exansions of
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More informationOn Doob s Maximal Inequality for Brownian Motion
Stochastic Process. Al. Vol. 69, No., 997, (-5) Research Reort No. 337, 995, Det. Theoret. Statist. Aarhus On Doob s Maximal Inequality for Brownian Motion S. E. GRAVERSEN and G. PESKIR If B = (B t ) t
More informationOn the Properties for Iteration of a Compact Operator with Unstructured Perturbation
On the Proerties for Iteration of a Comact Oerator with Unstructured Perturbation Jarmo Malinen Helsinki University of Technology Institute of Mathematics Research Reorts A360 (1996) 1 2 Jarmo Malinen:
More informationASYMPTOTIC BEHAVIOR FOR THE BEST LOWER BOUND OF JENSEN S FUNCTIONAL
75 Kragujevac J. Math. 25 23) 75 79. ASYMPTOTIC BEHAVIOR FOR THE BEST LOWER BOUND OF JENSEN S FUNCTIONAL Stojan Radenović and Mirjana Pavlović 2 University of Belgrade, Faculty of Mechanical Engineering,
More informationExistence and Approximation of Fixed Points of. Bregman Nonexpansive Operators. Banach Spaces
Existence and Approximation of Fixed Points of in Reflexive Banach Spaces Department of Mathematics The Technion Israel Institute of Technology Haifa 22.07.2010 Joint work with Prof. Simeon Reich General
More informationA numerical implementation of a predictor-corrector algorithm for sufcient linear complementarity problem
A numerical imlementation of a redictor-corrector algorithm for sufcient linear comlementarity roblem BENTERKI DJAMEL University Ferhat Abbas of Setif-1 Faculty of science Laboratory of fundamental and
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationSOME INEQUALITIES FOR (α, β)-normal OPERATORS IN HILBERT SPACES. 1. Introduction
SOME INEQUALITIES FOR (α, β)-normal OPERATORS IN HILBERT SPACES SEVER S. DRAGOMIR 1 AND MOHAMMAD SAL MOSLEHIAN Abstract. An oerator T is called (α, β)-normal (0 α 1 β) if α T T T T β T T. In this aer,
More informationGeneralized Coiflets: A New Family of Orthonormal Wavelets
Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University
More informationConvex Feasibility Problems
Laureate Prof. Jonathan Borwein with Matthew Tam http://carma.newcastle.edu.au/drmethods/paseky.html Spring School on Variational Analysis VI Paseky nad Jizerou, April 19 25, 2015 Last Revised: May 6,
More informationFor q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i
Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:
More informationExamples from Elements of Theory of Computation. Abstract. Introduction
Examles from Elements of Theory of Comutation Mostafa Ghandehari Samee Ullah Khan Deartment of Comuter Science and Engineering University of Texas at Arlington, TX-7609, USA Tel: +(87)7-5688, Fax: +(87)7-784
More informationAn Algorithm for Total Variation Minimization and Applications
Journal of Mathematical Imaging and Vision 0: 89 97, 004 c 004 Kluwer Academic Publishers. Manufactured in The Netherlands. An Algorithm for Total Variation Minimization and Alications ANTONIN CHAMBOLLE
More informationSUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS *
Journal of Comutational Mathematics Vol.8, No.,, 48 48. htt://www.global-sci.org/jcm doi:.48/jcm.9.-m6 SUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationRotations in Curved Trajectories for Unconstrained Minimization
Rotations in Curved rajectories for Unconstrained Minimization Alberto J Jimenez Mathematics Deartment, California Polytechnic University, San Luis Obiso, CA, USA 9407 Abstract Curved rajectories Algorithm
More informationClass Numbers and Iwasawa Invariants of Certain Totally Real Number Fields
Journal of Number Theory 79, 249257 (1999) Article ID jnth.1999.2433, available online at htt:www.idealibrary.com on Class Numbers and Iwasawa Invariants of Certain Totally Real Number Fields Dongho Byeon
More informationON THE RESIDUE CLASSES OF (n) MODULO t
#A79 INTEGERS 3 (03) ON THE RESIDUE CLASSES OF (n) MODULO t Ping Ngai Chung Deartment of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts briancn@mit.edu Shiyu Li Det of Mathematics,
More information#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS
#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt
More informationAn interior-point gradient method for large-scale totally nonnegative least squares problems
An interior-point gradient method for large-scale totally nonnegative least squares problems Michael Merritt and Yin Zhang Technical Report TR04-08 Department of Computational and Applied Mathematics Rice
More informationCommutators on l. D. Dosev and W. B. Johnson
Submitted exclusively to the London Mathematical Society doi:10.1112/0000/000000 Commutators on l D. Dosev and W. B. Johnson Abstract The oerators on l which are commutators are those not of the form λi
More informationDIVISIBILITY CRITERIA FOR CLASS NUMBERS OF IMAGINARY QUADRATIC FIELDS
IVISIBILITY CRITERIA FOR CLASS NUMBERS OF IMAGINARY QUARATIC FIELS PAUL JENKINS AN KEN ONO Abstract. In a recent aer, Guerzhoy obtained formulas for certain class numbers as -adic limits of traces of singular
More informationTranspose of the Weighted Mean Matrix on Weighted Sequence Spaces
Transose of the Weighted Mean Matri on Weighted Sequence Saces Rahmatollah Lashkariour Deartment of Mathematics, Faculty of Sciences, Sistan and Baluchestan University, Zahedan, Iran Lashkari@hamoon.usb.ac.ir,
More informationAdditive results for the generalized Drazin inverse in a Banach algebra
Additive results for the generalized Drazin inverse in a Banach algebra Dragana S. Cvetković-Ilić Dragan S. Djordjević and Yimin Wei* Abstract In this aer we investigate additive roerties of the generalized
More informationNumerical Methods: Structured vs. unstructured grids. General Introduction: Why numerical methods? Numerical methods and their fields of application
Numerical Methods: Structured vs. unstructured grids The goals o this course General : Why numerical methods? Numerical methods and their ields o alication Review o inite dierences Goals: Understanding
More informationBest approximation by linear combinations of characteristic functions of half-spaces
Best aroximation by linear combinations of characteristic functions of half-saces Paul C. Kainen Deartment of Mathematics Georgetown University Washington, D.C. 20057-1233, USA Věra Kůrková Institute of
More informationOn the Distribution of Perfect Powers
2 3 47 6 23 Journal of Integer Sequences, Vol. 4 (20), rticle.8.5 On the Distribution of Perfect Powers Rafael Jakimczuk División Matemática Universidad Nacional de Luján Buenos ires rgentina jakimczu@mail.unlu.edu.ar
More informationIntroduction to MVC. least common denominator of all non-identical-zero minors of all order of G(s). Example: The minor of order 2: 1 2 ( s 1)
Introduction to MVC Definition---Proerness and strictly roerness A system G(s) is roer if all its elements { gij ( s)} are roer, and strictly roer if all its elements are strictly roer. Definition---Causal
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More information