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).     

Semilocal convergence for a class of improved multi-step Chebyshev–Halley-like methods under extended conditions

Semilocal convergence for a class of improved multi-step Chebyshev–Halley-like methods is considered in this paper. Compared with the results for the Chebyshev method in Hernández (J Comput Appl Math 126:131–143, 2000), the R-order of convergence is heightened and the Hölder continuity of second derivative is also relaxed. Moreover, an existence-uniqueness theorem is proved under the extended conditions.     

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.     

Semilocal convergence of multi-point improved super-Halley-type methods without the second derivative under generalized weak condition

In this paper, we consider the semilocal convergence of multi-point improved super-Halley-type methods in Banach space. Different from the results of super-Halley method studied in reference Gutiérrez, J.M. and Hernández, M.A. (Comput. Math. Appl. 36,1–8, 1998) these methods do not require second derivative of an operator, the R-order is improved and the convergence condition is also relaxed. We prove a convergence theorem to show existence and uniqueness of the solution.     

Preconditioning of spectral methods via Birkhoff interpolation

High-order differentiation matrices as calculated in spectral collocation methods usually include a large round-off error and have a large condition number (Baltensperger and Berrut Computers and Mathematics with Applications 37(1), 41–48 1999; Baltensperger and Trummer SIAM J. Sci. Comput. 24(5), 1465–1487 2003; Costa and Don Appl. Numer. Math. 33(1), 151–159 2000). Wang et al. (Wang et al. SIAM J. Sci. Comput. 36(3), A907–A929 2014) present a method to precondition these matrices using Birkhoff interpolation. We generalize this method for all orders and boundary conditions and allowing arbitrary rows of the system matrix to be replaced by the boundary conditions. The preconditioner is an exact inverse of the highest-order differentiation matrix in the equation; thus, its product with that matrix can be replaced by the identity matrix. We show the benefits of the method for high-order differential equations. These include improved condition number and, more importantly, higher accuracy of solutions compared to other methods.     

A Double-Parameter Scaling Broyden–Fletcher–Goldfarb–Shanno Method Based on Minimizing the Measure Function of Byrd and Nocedal for Unconstrained Optimization

In this paper, the first two terms on the right-hand side of the Broyden–Fletcher–Goldfarb–Shanno update are scaled with a positive parameter, while the third one is also scaled with another positive parameter. These scaling parameters are determined by minimizing the measure function introduced by Byrd and Nocedal (SIAM J Numer Anal 26:727–739, 1989). The obtained algorithm is close to the algorithm based on clustering the eigenvalues of the Broyden–Fletcher–Goldfarb–Shanno approximation of the Hessian and on shifting its large eigenvalues to the left, but it is not superior to it. Under classical assumptions, the convergence is proved by using the trace and the determinant of the iteration matrix. By using a set of 80 unconstrained optimization test problems, it is proved that the algorithm minimizing the measure function of Byrd and Nocedal is more efficient and more robust than some other scaling Broyden–Fletcher–Goldfarb–Shanno algorithms, including the variants of Biggs (J Inst Math Appl 12:337–338, 1973), Yuan (IMA J Numer Anal 11:325–332, 1991), Oren and Luenberger (Manag Sci 20:845–862, 1974) and of Nocedal and Yuan (Math Program 61:19–37, 1993). However, it is less efficient than the algorithms based on clustering the eigenvalues of the iteration matrix and on shifting its large eigenvalues to the left, as shown by Andrei (J Comput Appl Math 332:26–44, 2018, Numer Algorithms 77:413–432, 2018).     

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.     

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.     

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).     

Absolute Cesàro Series Spaces and Matrix Operators

In this paper we derive a series space \(\vert C_{\lambda,\mu} \vert _{k}\) using the well known absolute Cesàro summability \(\vert C_{\lambda,\mu} \vert _{k}\) of Das (Proc. Camb. Philol. Soc. 67:321–326, 1970), compute its \(\beta\)-dual, give some algebraic and topological properties, and characterize some matrix operators defined on that space. So we generalize some results of Bosanquet (J. Lond. Math. Soc. 20:39–48, 1945), Flett (Proc. Lond. Math. Soc. 7:113–141, 1957), Mehdi (Proc. Lond. Math. Soc. (3)10:180–199, 1960), Mazhar (Tohoku Math. J. 23:433–451, 1971), Orhan and Sar?göl (Rocky Mt. J. Math. 23(3):1091–1097, 1993) and Sar?göl (Commun. Math. Appl. 7(1):11–22, 2016; Math. Comput. Model. 55:1763–1769, 2012).     

Numerical Studies on Asymptotics of European Option Under Multiscale Stochastic Volatility

Multiscale stochastic volatilities models relax the constant volatility assumption from Black-Scholes option pricing model. Such models can capture the smile and skew of volatilities and therefore describe more accurately the movements of the trading prices. Christoffersen et al. Manag Sci 55(2):1914–1932 (2009) presented a model where the underlying price is governed by two volatility components, one changing fast and another changing slowly. Chiarella and Ziveyi Appl Math Comput 224:283–310 (2013) transformed Christoffersen’s model and computed an approximate formula for pricing American options. They used Duhamel’s principle to derive an integral form solution of the boundary value problem associated to the option price. Using method of characteristics, Fourier and Laplace transforms, they obtained with good accuracy the American option prices. In a previous research of the authors (Canhanga et al. 2014), a particular case of Chiarella and Ziveyi Appl Math Comput 224:283–310 (2013) model is used for pricing of European options. The novelty of this earlier work is to present an asymptotic expansion for the option price. The present paper provides experimental and numerical studies on investigating the accuracy of the approximation formulae given by this asymptotic expansion. We present also a procedure for calibrating the parameters produced by our first-order asymptotic approximation formulae. Our approximated option prices will be compared to the approximation obtained by Chiarella and Ziveyi Appl Math Comput 224:283–310 (2013).     

Semilocal convergence on a family of root-finding multi-point methods in Banach spaces under relaxed continuity condition

In this paper, we consider the semilocal convergence on a family of root-finding multi-point methods. Compared with the results in reference (Hernández, M.A., Salanova, M.A., J. Comput. Appl. Math. 126, 131–143 3), these multi-point methods do not require the second derivative, Hölder continuity condition is relaxed, and the R-order is also enhanced. We prove an existence-uniqueness theorem of the solution. The R-order for these multi-point methods is at least 6 + q with relaxed continuous second derivative, where q∈[0,1].     

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.     

Multivariate Higher-Degree Stochastic Increasing Convexity

Building on the seminal work by Shaked and Shanthikumar (Adv Appl Probab 20:427–446, 1988a; Stoch Process Appl 27:1–20, 1988b), Denuit et al. (Eng Inf Sci 13:275–291, 1999; Methodol Comput Appl Probab 2:231–254, 2000; 2001) studied the stochastic s-increasing convexity properties of standard parametric families of distributions. However, the analysis is restricted there to a single parameter. As many standard families of distributions involve several parameters, multivariate higher-order stochastic convexity properties also deserve consideration for applications. This is precisely the topic of the present paper, devoted to stochastic \((s_1,s_2,\ldots ,s_d)\)-increasing convexity of distribution families indexed by a vector \((\theta _1,\theta _2,\ldots ,\theta _d)\) of parameters. This approach accounts for possible correlation in multivariate mixture models.     

Remarks on defective Fano manifolds

This note continues our previous work on special secant defective (specifically, conic connected and local quadratic entry locus) and dual defective manifolds. These are now well understood, except for the prime Fano ones. Here we add a few remarks on this case, completing the results in our papers (Russo in Math Ann 344:597–617, 2009; Ionescu and Russo in Compos Math 144:949–962, 2008; Ionescu and Russo in J Reine Angew Math 644:145–157, 2010; Ionescu and Russo in Am J Math 135:349–360, 2013; Ionescu and Russo in Math Res Lett 21:1137–1154, 2014); see also the recent book (Russo, On the Geometry of Some Special Projective Varieties, Lecture Notes of the Unione Matematica Italiana, Springer, 2016).     

Performance of convex underestimators in a branch-and-bound framework

The efficient determination of tight lower bounds in a branch-and-bound algorithm is crucial for the global optimization of models spanning numerous applications and fields. The global optimization method \(\alpha \)-branch-and-bound (\(\alpha \)BB, Adjiman et al. in Comput Chem Eng 22(9):1159–1179, 1998b, Comput Chem Eng 22(9):1137–1158, 1998a; Adjiman and Floudas in J Global Optim 9(1):23–40, 1996; Androulakis et al. J Global Optim 7(4):337–363, 1995; Floudas in Deterministic Global Optimization: Theory, Methods and Applications, vol. 37. Springer, Berlin, 2000; Maranas and Floudas in J Chem Phys 97(10):7667–7678, 1992, J Chem Phys 100(2):1247–1261, 1994a, J Global Optim 4(2):135–170, 1994), guarantees a global optimum with \(\epsilon \)-convergence for any \(\mathcal {C}^2\)-continuous function within a finite number of iterations via fathoming nodes of a branch-and-bound tree. We explored the performance of the \(\alpha \)BB method and a number of competing methods designed to provide tight, convex underestimators, including the piecewise (Meyer and Floudas in J Global Optim 32(2):221–258, 2005), generalized (Akrotirianakis and Floudas in J Global Optim 30(4):367–390, 2004a, J Global Optim 29(3):249–264, 2004b), and nondiagonal (Skjäl et al. in J Optim Theory Appl 154(2):462–490, 2012) \(\alpha \)BB methods, the Brauer and Rohn+E (Skjäl et al. in J Global Optim 58(3):411–427, 2014) \(\alpha \)BB methods, and the moment method (Lasserre and Thanh in J Global Optim 56(1):1–25, 2013). Using a test suite of 40 multivariate, box-constrained, nonconvex functions, the methods were compared based on the tightness of generated underestimators and the efficiency of convergence of a branch-and-bound global optimization algorithm.     

A general pinching principle for mean curvature flow and applications

We prove a sharp pinching estimate for immersed mean convex solutions of mean curvature flow which unifies and improves all previously known pinching estimates, including the umbilic estimate of Huisken (J Differ Geom 20(1):237–266, 1984), the convexity estimates of Huisken–Sinestrari (Acta Math 183(1):45–70, 1999) and the cylindrical estimate of Huisken–Sinestrari (Invent Math 175(1):137–221, 2009; see also Andrews and Langford in Anal PDE 7(5):1091–1107, 2014; Huisken and Sinestrari in J Differ Geom 101(2):267–287, 2015). Namely, we show that the curvature of the solution pinches onto the convex cone generated by the curvatures of any shrinking cylinder solutions admitted by the initial data. For example, if the initial data is \((m+1)\)-convex, then the curvature of the solution pinches onto the convex hull of the curvatures of the shrinking cylinders \(\mathbb {R}^m\times S^{n-m}_{\sqrt{2(n-m)(1-t)}}\), \(t<1\). In particular, this yields a sharp estimate for the largest principal curvature, which we use to obtain a new proof of a sharp estimate for the inscribed curvature for embedded solutions (Brendle in Invent Math 202(1):217–237, 2015; Haslhofer and Kleiner in Int Math Res Not 15:6558–6561, 2015; Langford in Proc Am Math Soc 143(12):5395–5398, 2015). Making use of a recent idea of Huisken–Sinestrari (2015), we then obtain a series of sharp estimates for ancient solutions. In particular, we obtain a convexity estimate for ancient solutions which allows us to strengthen recent characterizations of the shrinking sphere due to Huisken–Sinestrari (2015) and Haslhofer–Hershkovits (Commun Anal Geom 24(3):593–604, 2016).     

Asymptotic Behavior of Underlying NT Paths in Interior Point Methods for Monotone Semidefinite Linear Complementarity Problems

An interior point method (IPM) defines a search direction at each interior point of the feasible region. These search directions form a direction field, which in turn gives rise to a system of ordinary differential equations (ODEs). Thus, it is natural to define the underlying paths of the IPM as solutions of the system of ODEs. In Sim and Zhao (Math. Program. Ser. A 110:475–499, 2007), these off-central paths are shown to be well-defined analytic curves and any of their accumulation points is a solution to the given monotone semidefinite linear complementarity problem (SDLCP). In Sim and Zhao (Math. Program. Ser. A 110:475–499, 2007; J. Optim. Theory Appl. 137:11–25, 2008) and Sim (J. Optim. Theory Appl. 141:193–215, 2009), the asymptotic behavior of off-central paths corresponding to the HKM direction is studied. In particular, in Sim and Zhao (Math. Program. Ser. A 110:475–499, 2007), the authors study the asymptotic behavior of these paths for a simple example, while, in Sim and Zhao (J. Optim. Theory Appl. 137:11–25, 2008) and Sim (J. Optim. Theory Appl. 141:193–215, 2009), the asymptotic behavior of these paths for a general SDLCP is studied. In this paper, we study off-central paths corresponding to another well-known direction, the Nesterov-Todd (NT) direction. Again, we give necessary and sufficient conditions for these off-central paths to be analytic w.r.t. \(\sqrt{\mu}\) and then w.r.t. μ, at solutions of a general SDLCP. Also, as in Sim and Zhao (Math. Program. Ser. A 110:475–499, 2007), we present off-central path examples using the same SDP, whose first derivatives are likely to be unbounded as they approach the solution of the SDP. We work under the assumption that the given SDLCP satisfies a strict complementarity condition.     

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.     

O $$\varvec{}$$-orm preservers ad the Aleksadrov coservative $$\varvec{}$$-distace problem

The goal of this paper is to point out that the results obtained in the recent papers (Chen and Song in Nonlinear Anal 72:1895–1901, 2010; Chu in J Math Anal Appl 327:1041–1045, 2007; Chu et al. in Nonlinear Anal 59:1001–1011, 2004a, J. Math Anal Appl 289:666–672, 2004b) can be seriously strengthened in the sense that we can significantly relax the assumptions of the main results so that we still get the same conclusions. In order to do this first, we prove that for \(n \ge 3\) any transformation which preserves the n-norm of any n vectors is automatically plus-minus linear. This will give a re-proof of the well-known Mazur–Ulam-type result that every n-isometry is automatically affine (\(n \ge 2\)) which was proven in several papers, e.g. in Chu et al. (Nonlinear Anal 70:1068–1074, 2009). Second, following the work of Rassias and ?emrl (Proc Am Math Soc 118:919–925, 1993), we provide the solution of a natural Aleksandrov-type problem in n-normed spaces, namely, we show that every surjective transformation which preserves the unit n-distance in both directions (\(n\ge 2\)) is automatically an n-isometry.