A contribution to perturbation analysis for total least squares problems

In this note, we present perturbation analysis for the total least squares (Tls) problems under the genericity condition. We review the three condition numbers proposed respectively by Zhou et al. (Numer. Algorithm, 51 (2009), pp. 381–399), Baboulin and Gratton (SIAM J. Matrix Anal. Appl. 32 (2011), pp. 685–699), Li and Jia (Linear Algebra Appl. 435 (2011), pp. 674–686). We also derive new perturbation bounds.     

Finite-Codimensional Compressions of Symmetric and Self-Adjoint Linear Relations in Krein Spaces

Theorems due to Stenger (Bull Am Math Soc 74:369–372, 1968) and Nudelman (Int Equ Oper Theory 70:301–305, 2011) in Hilbert spaces and their generalizations to Krein spaces in Azizov and Dijksma (Int Equ Oper Theory 74(2):259–269, 2012) and Azizov et al. (Linear Algebra Appl 439:771–792, 2013) generate additional questions about properties a finite-codimensional compression \({T_0}\) of a symmetric or self-adjoint linear relation \({T}\) may or may not inherit from \({T}\). These questions concern existence of invariant maximal nonnegative subspaces, definitizability, singular critical points and defect indices.     

Local convergence and the dynamics of a two-point four parameter Jarratt-like method under weak conditions

We present a local convergence analysis of a two-point four parameter Jarratt-like method of high convergence order in order to approximate a locally unique solution of a nonlinear equation. In contrast to earlier studies such us (Amat et al. Aequat. Math. 69(3), 212–223 2015; Amat et al. J. Math. Anal. Appl. 366(3), 24–32 2010; Behl, R. 2013; Bruns and Bailey Chem. Eng. Sci. 32, 257–264 1977; Candela and Marquina. Computing 44, 169–184 1990; Candela and Marquina. Computing 45(4), 355–367 1990; Chun. Appl. Math. Comput. 190(2), 1432–1437 2007; Cordero and Torregrosa. Appl. Math. Comput. 190, 686–698 2007; Deghan. Comput. Appl Math. 29(1), 19–30 2010; Deghan. Comput. Math. Math. Phys. 51(4), 513–519 2011; Deghan and Masoud. Eng. Comput. 29(4), 356–365 15; Cordero and Torregrosa. Appl. Math. Comput. 190, 686–698 2012; Deghan and Masoud. Eng. Comput. 29(4), 356–365 2012; Ezquerro and Hernández. Appl. Math. Optim. 41(2), 227–236 2000; Ezquerro and Hernández. BIT Numer. Math. 49, 325–342 2009; Ezquerro and Hernández. J. Math. Anal. Appl. 303, 591–601 2005; Gutiérrez and Hernández. Comput. Math. Appl. 36(7), 1–8 1998; Ganesh and Joshi. IMA J. Numer. Anal. 11, 21–31 1991; González-Crespo et al. Expert Syst. Appl. 40(18), 7381–7390 2013; Hernández. Comput. Math. Appl. 41(3-4), 433–455 2001; Hernández and Salanova. Southwest J. Pure Appl. Math. 1, 29–40 1999; Jarratt. Math. Comput. 20(95), 434–437 1966; Kou and Li. Appl. Math. Comput. 189, 1816–1821 2007; Kou and Wang. Numer. Algor. 60, 369–390 2012; Lorenzo et al. Int. J. Interact. Multimed. Artif. Intell. 1(3), 60–66 2010; Magreñán. Appl. Math. Comput. 233, 29–38 2014; Magreñán. Appl. Math. Comput. 248, 215–224 2014; Parhi and Gupta. J. Comput. Appl. Math. 206(2), 873–887 2007; Rall 1979; Ren et al. Numer. Algor. 52(4), 585–603 2009; Rheinboldt Pol. Acad. Sci. Banach Ctr. Publ. 3, 129–142 1978; Sicilia et al. J. Comput. Appl. Math. 291, 468–477 2016; Traub 1964; Wang et al. Numer. Algor. 57, 441–456 2011) using hypotheses up to the fifth derivative, our sufficient convergence conditions involve only hypotheses on the first Fréchet-derivative of the operator involved. The dynamics of the family for choices of the parameters such that it is optimal is also shown. Numerical examples are also provided in this study     

Improved local convergence analysis of the Gauss–Newton method under a majorant condition

We present a local convergence analysis of Gauss-Newton method for solving nonlinear least square problems. Using more precise majorant conditions than in earlier studies such as Chen (Comput Optim Appl 40:97–118, 2008), Chen and Li (Appl Math Comput 170:686–705, 2005), Chen and Li (Appl Math Comput 324:1381–1394, 2006), Ferreira (J Comput Appl Math 235:1515–1522, 2011), Ferreira and Gonçalves (Comput Optim Appl 48:1–21, 2011), Ferreira and Gonçalves (J Complex 27(1):111–125, 2011), Li et al. (J Complex 26:268–295, 2010), Li et al. (Comput Optim Appl 47:1057–1067, 2004), Proinov (J Complex 25:38–62, 2009), Ewing, Gross, Martin (eds.) (The merging of disciplines: new directions in pure, applied and computational mathematics 185–196, 1986), Traup (Iterative methods for the solution of equations, 1964), Wang (J Numer Anal 20:123–134, 2000), we provide a larger radius of convergence; tighter error estimates on the distances involved and a clearer relationship between the majorant function and the associated least squares problem. Moreover, these advantages are obtained under the same computational cost.     

Reformulation of Hensel’s Lemma and extension of a theorem of Ore

In this paper, we extend the theorem of Ore regarding factorization of polynomials over p-adic numbers to henselian valued fields of arbitrary rank thereby generalizing the main results of Khanduja and Kumar (J Pure Appl Algebra 216:2648–2656, 2012) and Cohen et al. (Mathematika 47:173–196, 2000). As an application, we derive the analogue of Dedekind’s Theorem regarding splitting of rational primes in algebraic number fields as well as of its converse for general valued fields extending similar results proved for discrete valued fields in Khanduja and Kumar (Int J Number Theory 4:1019–1025, 2008). The generalized version of Ore’s Theorem leads to an extension of a result of Weintraub dealing with a generalization of Eisenstein Irreducibility Criterion (cf. Weintraub in Proc Am Math Soc 141:1159–1160, 2013). We also give a reformulation of Hensel’s Lemma for polynomials with coefficients in henselian valued fields which is used in the proof of the extended Ore’s Theorem and was proved in Khanduja and Kumar (J Algebra Appl 12:1250125, 2013) in the particular case of complete rank one valued fields.     

Central Limit Theorem for Lipschitz–Killing Curvatures of Excursion Sets of Gaussian Random Fields

Our interest in this paper is to explore limit theorems for various geometric functionals of excursion sets of isotropic Gaussian random fields. In the past, asymptotics of nonlinear functionals of Gaussian random fields have been studied [see Berman (Sojourns and extremes of stochastic processes, Wadsworth & Brooks, Monterey, 1991), Kratz and León (Extremes 3(1):57–86, 2000), Kratz and León (J Theor Probab 14(3):639–672, 2001), Meshenmoser and Shashkin (Stat Probab Lett 81(6):642–646, 2011), Pham (Stoch Proc Appl 123(6):2158–2174, 2013), Spodarev (Chapter in modern stochastics and applications, volume 90 of the series Springer optimization and its applications, pp 221–241, 2013) for a sample of works in such settings], the most recent addition being (Adler and Naitzat in Stoch Proc Appl 2016; Estrade and León in Ann Probab 2016) where a central limit theorem (CLT) for Euler integral and Euler–Poincaré characteristic, respectively, of the excursions set of a Gaussian random field is proven under some conditions. In this paper, we obtain a CLT for some global geometric functionals, called the Lipschitz–Killing curvatures of excursion sets of Gaussian random fields, in an appropriate setting.     

Spectral Properties of Unbounded Jacobi Matrices with Almost Monotonic Weights

We present a unified framework to identify spectra of Jacobi matrices. We give applications of the long-standing problem of Chihara (Mt J Math 21(1):121–137, 1991, J Comput Appl Math 153(1–2):535–536, 2003) concerning one-quarter class of orthogonal polynomials, to the conjecture posed by Roehner and Valent (SIAM J Appl Math 42(5):1020–1046, 1982) concerning continuous spectra of generators of birth and death processes, and to spectral properties of operators studied by Janas and Moszyńki (Integral Equ Oper Theory 43(4):397–416, 2002) and Pedersen (Proc Am Math Soc 130(8):2369–2376, 2002).     

Discrete derived categories I: homomorphisms,autoequivalences and t-structures

Discrete derived categories were studied initially by Vossieck (J Algebra 243:168–176, 2001) and later by Bobiński et al. (Cent Eur J Math 2:19–49, 2004). In this article, we describe the homomorphism hammocks and autoequivalences on these categories. We classify silting objects and bounded t-structures.     

On manifolds of tensors of fixed TT-rank

Recently, the format of TT tensors (Hackbusch and Kühn in J Fourier Anal Appl 15:706–722, 2009; Oseledets in SIAM J Sci Comput 2009, submitted; Oseledets and Tyrtyshnikov in SIAM J Sci Comput 31:5, 2009; Oseledets and Tyrtyshnikov in Linear Algebra Appl 2009, submitted) has turned out to be a promising new format for the approximation of solutions of high dimensional problems. In this paper, we prove some new results for the TT representation of a tensor \({U \in \mathbb{R}^{n_1\times \cdots\times n_d}}\) and for the manifold of tensors of TT-rank \({\underline{r}}\) . As a first result, we prove that the TT (or compression) ranks r _i of a tensor U are unique and equal to the respective separation ranks of U if the components of the TT decomposition are required to fulfil a certain maximal rank condition. We then show that the set \({\mathbb{T}}\) of TT tensors of fixed rank \({\underline{r}}\) locally forms an embedded manifold in \({\mathbb{R}^{n_1\times\cdots\times n_d}}\) , therefore preserving the essential theoretical properties of the Tucker format, but often showing an improved scaling behaviour. Extending a similar approach for matrices (Conte and Lubich in M2AN 44:759, 2010), we introduce certain gauge conditions to obtain a unique representation of the tangent space \({\mathcal{T}_U\mathbb{T}}\) of \({\mathbb{T}}\) and deduce a local parametrization of the TT manifold. The parametrisation of \({\mathcal{T}_{U}\mathbb{T}}\) is often crucial for an algorithmic treatment of high-dimensional time-dependent PDEs and minimisation problems (Lubich in From quantum to classical molecular dynamics: reduced methods and numerical analysis, 2008). We conclude with remarks on those applications and present some numerical examples.     

Extensions of the Noncommutative Integration

In this paper we will continue the analysis undertaken in Bagarello et al. (Rend Circ Mat Palermo (2) 55:21–28, 2006), Bongiorno et al. (Rocky Mt J Math 40(6):1745–1777, 2010), Triolo (Rend Circ Mat Palermo (2) 60(3):409–416, 2011) on the general problem of extending the noncommutative integration in a *-algebra of measurable operators. As in Aiena et al. (Filomat 28(2):263–273, 2014), Bagarello (Stud Math 172(3):289–305, 2006) and Bagarello et al. (Rend Circ Mat Palermo (2) 55:21–28, 2006), the main problem is to represent different types of partial *-algebras into a *-algebra of measurable operators in Segal’s sense, provided that these partial *-algebras posses a sufficient family of positive linear functionals (states) (Fragoulopoulou et al., J Math Anal Appl 388(2):1180–1193, 2012; Trapani and Triolo, Stud Math 184(2):133–148, 2008; Trapani and Triolo, Rend Circolo Mat Palermo 59:295–302, 2010; La Russa and Triolo, J Oper Theory, 69:2, 2013; Triolo, J Pure Appl Math, 43(6):601–617, 2012). In this paper, a new condition is given in an attempt to provide a extension of the non commutative integration.     

Superoscillating Sequences in Several Variables

In a series of papers (J Phys A 44:365304, 2011; Complex Anal Oper Theory 7:1299–1310, 2013; J Math Pures Appl 99:165–173, 2013; J Math Pures Appl 103:522–534, 2015), we have investigated some mathematical properties of superoscillating sequences in one variable, and their persistence in time. In this paper we study the notion of superoscillation in several variables and we show how to construct examples of sequences that exhibit this property.     

On <Emphasis Type="Italic">m</Emphasis>-ovoids of regular near polygons

We generalise the work of Segre (Ann Mat Pura Appl 4(70):1–201, 1965), Cameron et al. (J Algebra 55(2):257–280, 1978), and Vanhove (J Algebr Comb 34(3):357–373, 2011) by showing that nontrivial m-ovoids of the dual polar spaces \(\mathsf {DQ}(2d, q)\), \(\mathsf {DW}(2d-1,q)\) and \(\mathsf {DH}(2d-1,q^2)\) (\(d\geqslant 3\)) are hemisystems. We also provide a more general result that holds for regular near polygons.     

On the solution of systems of equations with constant rank derivatives

The famous for its simplicity and clarity Newton–Kantorovich hypothesis of Newton’s method has been used for a long time as the sufficient convergence condition for solving nonlinear equations. Recently, in the elegant study by Hu et al. (J Comput Appl Math 219:110–122, 2008), a Kantorovich-type convergence analysis for the Gauss–Newton method (GNM) was given improving earlier results by Häubler (Numer Math 48:119–125, 1986), and extending some results by Argyros (Adv Nonlinear Var Inequal 8:93–99, 2005, 2007) to hold for systems of equations with constant rank derivatives. In this study, we use our new idea of recurrent functions to extend the applicability of (GNM) by replacing existing conditions by weaker ones. Finally, we provide numerical examples to solve equations in cases not covered before (Häubler, Numer Math 48:119–125, 1986; Hu et al., J Comput Appl Math 219:110–122, 2008; Kontorovich and Akilov 2004).     

Data dependence for the amplitude equation of surface waves

We consider the amplitude equation for nonlinear surface wave solutions of hyperbolic conservation laws. This is an asymptotic nonlocal, Hamiltonian evolution equation with quadratic nonlinearity. For example, this equation describes the propagation of nonlinear Rayleigh waves (Hamilton et al. in J Acoust Soc Am 97:891–897, 1995), surface waves on current-vortex sheets in incompressible MHD (Alì and Hunter in Q Appl Math 61(3):451–474, 2003; Alì et al. in Stud Appl Math 108(3):305–321, 2002) and on the incompressible plasma–vacuum interface (Secchi in Q Appl Math 73(4):711–737, 2015). The local-in-time existence of smooth solutions to the Cauchy problem for the amplitude equation in noncanonical variables was shown in Hunter (J Hyperbolic Differ Equ 3(2):247–267, 2006), Secchi (Q Appl Math 73(4):711–737, 2015). In the present paper we prove the continuous dependence in strong norm of solutions on the initial data. This completes the proof of the well-posedness of the problem in the classical sense of Hadamard.     

Polynomial,rational and analytic first integrals for a family of 3-dimensional Lotka-Volterra systems

We extend the study of the integrability done by Leach and Miritzis (J Nonlinear Math Phys 13:535–548, 2006) on the classical model of competition between three species studied by May and Leonard (SIAM J Appl Math 29:243–256, 1975), to all real values of the parameters. Additionally, our results provide all polynomial, rational and analytic first integrals of this extended model. We also classify all the invariant algebraic surfaces of these models.     

Positivity for Generalized Cluster Variables of Affine Quivers

It has been proved in Lee and Schiffler, Ann. of Math. 182(1) 73–125 2015 that cluster variables of all skew-symmetric cluster algebras are positive. i.e., every cluster variable as a Laurent polynomial in the cluster variables of any fixed cluster has positive coefficients. We prove that every regular generalized cluster variable of an affine quiver is positive. As a corollary, we obtain that generalized cluster variables of affine quivers are positive and we also construct various positive bases. This generalizes the results in Dupont, J. Algebra Appl. 11(4) 19 2012 and Ding et al. Algebr. Represent. Theory 16(2) 491–525 2013.     

On the Semicircular Law of Large-Dimensional Random Quaternion Matrices

It is well known that the Gaussian symplectic ensemble is defined on the space of \(n\times n\) quaternion self-dual Hermitian matrices with Gaussian random elements. There is a huge body of literature regarding this kind of matrices based on the exact known form of the density function of the eigenvalues (see Erd?s in Russ Math Surv 66(3):507–626, 2011; Erd?s in Probab Theory Relat Fields 154(1–2):341–407, 2012; Erd?s et al. in Adv Math 229(3):1435–1515, 2012; Knowles and Yin in Probab Theory Relat Fields, 155(3–4):543–582, 2013; Tao and Vu in Acta Math 206(1):127–204, 2011; Tao and Vu in Electron J Probab 16(77):2104–2121, 2011). Due to the fact that multiplication of quaternions is not commutative, few works about large-dimensional quaternion self-dual Hermitian matrices are seen without normality assumptions. As in natural, we shall get more universal results by removing the Gaussian condition. For the first step, in this paper, we prove that the empirical spectral distribution of the common quaternion self-dual Hermitian matrices tends to the semicircular law. The main tool to establish the universal result is given as a lemma in this paper as well.     

Condition Length and Complexity for the Solution of Polynomial Systems

Smale’s 17th problem asks for an algorithm which finds an approximate zero of polynomial systems in average polynomial time (see Smale in Mathematical problems for the next century, American Mathematical Society, Providence, 2000). The main progress on Smale’s problem is Beltrán and Pardo (Found Comput Math 11(1):95–129, 2011) and Bürgisser and Cucker (Ann Math 174(3):1785–1836, 2011). In this paper, we will improve on both approaches and prove an interesting intermediate result on the average value of the condition number. Our main results are Theorem 1 on the complexity of a randomized algorithm which improves the result of Beltrán and Pardo (2011), Theorem 2 on the average of the condition number of polynomial systems which improves the estimate found in Bürgisser and Cucker (2011), and Theorem 3 on the complexity of finding a single zero of polynomial systems. This last theorem is similar to the main result of Bürgisser and Cucker (2011) but relies only on homotopy methods, thus removing the need for the elimination theory methods used in Bürgisser and Cucker (2011). We build on methods developed in Armentano et al. (2014).     

A new iterative method for solving a class of complex symmetric system of linear equations

We present a new stationary iterative method, called Scale-Splitting (SCSP) method, and investigate its convergence properties. The SCSP method naturally results in a simple matrix splitting preconditioner, called SCSP-preconditioner, for the original linear system. Some numerical comparisons are presented between the SCSP-preconditioner and several available block preconditioners, such as PGSOR (Hezari et al. Numer. Linear Algebra Appl. 22, 761–776, 2015) and rotate block triangular preconditioners (Bai Sci. China Math. 56, 2523–2538, 2013), when they are applied to expedite the convergence rate of Krylov subspace iteration methods for solving the original complex system and its block real formulation, respectively. Numerical experiments show that the SCSP-preconditioner can compete with PGSOR-preconditioner and even more effective than the rotate block triangular preconditioners.