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.     

On the Gauss–Newton method

We provide a new semilocal convergence analysis of the Gauss–Newton method (GNM) for solving nonlinear equation in the Euclidean space. Using a combination of center-Lipschitz, Lipschitz conditions, and our new idea of recurrent functions, we provide under the same or weaker hypotheses than before (Ben-Israel, J. Math. Anal. Appl. 15:243–252, 1966; Chen and Nashed, Numer. Math. 66:235–257, 1993; Deuflhard and Heindl, SIAM J. Numer. Anal. 16:1–10, 1979; Guo, J. Comput. Math. 25:231–242, 2007; Häußler, Numer. Math. 48:119–125, 1986; Hu et al., J. Comput. Appl. Math. 219:110–122, 2008; Kantorovich and Akilov, Functional Analysis in Normed Spaces, Pergamon, Oxford, 1982), a finer convergence analysis. The results can be extended in case outer or generalized inverses are used. Numerical examples are also provided to show that our results apply, where others fail (Ben-Israel, J. Math. Anal. Appl. 15:243–252, 1966; Chen and Nashed, Numer. Math. 66:235–257, 1993; Deuflhard and Heindl, SIAM J. Numer. Anal. 16:1–10, 1979; Guo, J. Comput. Math. 25:231–242, 2007; Häußler, Numer. Math. 48:119–125, 1986; Hu et al., J. Comput. Appl. Math. 219:110–122, 2008; Kantorovich and Akilov, Functional Analysis in Normed Spaces, Pergamon, Oxford, 1982).     

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

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.     

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

Fast convergence of an inexact interior point method for horizontal complementarity problems

In Andreani et al. (Numer. Algorithms 57:457–485, 2011), an interior point method for the horizontal nonlinear complementarity problem was introduced. This method was based on inexact Newton directions and safeguarding projected gradient iterations. Global convergence, in the sense that every cluster point is stationary, was proved in Andreani et al. (Numer. Algorithms 57:457–485, 2011). In Andreani et al. (Eur. J. Oper. Res. 249:41–54, 2016), local fast convergence was proved for the underdetermined problem in the case that the Newtonian directions are computed exactly. In the present paper, it will be proved that the method introduced in Andreani et al. (Numer. Algorithms 57:457–485, 2011) enjoys fast (linear, superlinear, or quadratic) convergence in the case of truly inexact Newton computations. Some numerical experiments will illustrate the accuracy of the convergence theory.     

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

New characterizations of the Takagi function via functional equations

We provide two new characterizations of the Takagi function as the unique bounded solution of some systems of two functional equations. The results are independent of those obtained by Kairies (Wy? Szko? Ped Krakow Rocznik Nauk Dydakt Prace Mat 196:73–82, 1998), Kairies (Aequ Math 53:207–241, 1997), Kairies (Aequ Math 58:183–191, 1999) and Kairies et al. (Rad Mat 4:361–374, 1989; Errata, Rad Mat 5:179–180, 1989).     

On spectral analysis in varieties containing the solutions of inhomogeneous linear functional equations

The aim of the paper is to investigate the solutions of special inhomogeneous linear functional equations using spectral analysis in a translation invariant closed linear subspace of additive/multiadditive functions containing the restrictions of the solutions to finitely generated fields. The application of spectral analysis in some related varieties is a new and important trend in the theory of functional equations; especially they have successful applications in the case of homogeneous linear functional equations. The foundations of the theory can be found in Kiss and Varga (Aequat Math 88(1):151–162, 2014) and Kiss and Laczkovich (Aequat Math 89(2):301–328, 2015). We are going to adopt the main theoretical tools to solve some inhomogeneous problems due to Kocl?ga-Kulpa and Szostok (Ann Math Sylesianae 22:27–40, 2008), see also Kocl?ga-Kulpa and Szostok (Georgian Math J 16:725–736, 2009; Acta Math Hung 130(4):340–348, 2011). They are motivated by quadrature rules of approximate integration.     

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

A Complex Analogue of Toda’s Theorem

Toda (SIAM J. Comput. 20(5):865–877, 1991) proved in 1989 that the (discrete) polynomial time hierarchy, PH, is contained in the class P ^#P , namely the class of languages that can be decided by a Turing machine in polynomial time given access to an oracle with the power to compute a function in the counting complexity class #P. This result, which illustrates the power of counting, is considered to be a seminal result in computational complexity theory. An analogous result (with a compactness hypothesis) in the complexity theory over the reals (in the sense of Blum–Shub–Smale real machines (Blum et al. in Bull. Am. Math. Soc. 21(1):1–46, 1989) was proved in Basu and Zell (Found. Comput. Math. 10(4):429–454, 2010). Unlike Toda’s proof in the discrete case, which relied on sophisticated combinatorial arguments, the proof in Basu and Zell (Found. Comput. Math. 10(4):429–454, 2010) is topological in nature; the properties of the topological join are used in a fundamental way. However, the constructions used in Basu and Zell (Found. Comput. Math. 10(4):429–454, 2010) were semi-algebraic—they used real inequalities in an essential way and as such do not extend to the complex case. In this paper, we extend the techniques developed in Basu and Zell (Found. Comput. Math. 10(4):429–454, 2010) to the complex projective case. A key role is played by the complex join of quasi-projective complex varieties. As a consequence, we obtain a complex analogue of Toda’s theorem. The results of this paper, combined with those in Basu and Zell (Found. Comput. Math. 10(4):429–454, 2010), illustrate the central role of the Poincaré polynomial in algorithmic algebraic geometry, as well as in computational complexity theory over the complex and real numbers: the ability to compute it efficiently enables one to decide in polynomial time all languages in the (compact) polynomial hierarchy over the appropriate field.     

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.     

Computationally efficient modular nonlinear filter stabilization for high Reynolds number flows

The nonlinear filter based stabilization proposed in Layton et al. (J. Math. Fluid Mech. 14(2), 325–354 2012) allows to incorporate an eddy viscosity model into an existing laminar flow codes in a modular way. However, the proposed nonlinear filtering step requires the assembly of the associated matrix at each time step and solving a linear system with an indefinte matrix. We propose computationally efficient version of the filtering step that only requires the assembly once, and the solution of two symmetric, positive definite systems at each time step. We also test a new indicator function based on the entropy viscosity model of Guermond (Int. J. Numer. Meth. Fluids. 57(9), 1153–1170 2008); Guermond et al. (J. Sci. Comput. 49(1), 35–50 2011).     

On a rigidity result of singularities

The aim of this note is to prove, in the spirit of a rigidity result for isolated singularities of Schlessinger see Schlessinger (Invent Math 14:17–26, 1971) or also Kleiman and Landolfi (Compositio Math 23:407–434, 1971), a variant of a rigidity criterion for arbitrary singularities (Theorem 2.1 below). The proof of this result does not use Schlessinger’s Deformation Theory [Schlessinger (Trans Am Math Soc 130:208–222, 1968) and Schlessinger (Invent Math 14:17–26, 1971)]. Instead it makes use of Local Grothendieck-Lefschetz Theory, see (Grothendieck 1968, Éxposé 9, Proposition 1.4, page 106) and a Lemma of Zariski, see (Zariski, Am J Math 87:507–536, 1965, Lemma 4, page 526). I hope that this proof, although works only in characteristic zero, might also have some interest in its own.     

Ax–Schanuel for linear differential equations

We generalise the exponential Ax–Schanuel theorem to arbitrary linear differential equations with constant coefficients. Using the analysis of the exponential differential equation by Kirby (The theory of exponential differential equations, 2006, Sel Math 15(3):445–486, 2009) and Crampin (Reducts of differentially closed fields to fields with a relation for exponentiation, 2006) we give a complete axiomatisation of the first order theories of linear differential equations and show that the generalised Ax–Schanuel inequalities are adequate for them.     

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

The Existence of Fixed Points via the Measure of Noncompactness and its Application to Functional-Integral Equations

In the present paper, using the concept of measure of non-compactness, we introduce the concept of a new contraction on a Banach space and obtain few generalizations of Darbo’s fixed-point theorem and extend some recent results of (Aghajani et al., J. Comput. Appl. Math. 260:68–77, 2014) and (Aghajani et al., Bull. Belg. Math. Soc. Simon Stevin 20:345–358, 2013). Also we show the applicability of obtained results to the theory of integral equations. A concrete example illustrating the mentioned applicability is also included.     

On an improved local convergence analysis for the Secant method

We provide a local convergence analysis for the Secant method in a Banach space setting under Hölder continuous conditions. Using more precise estimates, and under the same computational cost, we enlarge the radius of convergence obtained in Ren and Wu (J Comput Appl Math 194:284–293, 2006).