Abel maps and presentation schemes

We sharpen the two main tools used to treat the compactified Jacobian of a singular curve: Abel maps and presentation schemes. First we prove a smoothness theorem for bigraded Abel maps. Second we study the two complementary filtrations provided by the images of certain Abel maps and certain presentation schemes. Third we study a lifting of the Abel map of bidegree (m, 1) to the corresponding presentation scheme. Fourth we prove that, if a curve is blown up at a double point, then the corresponding presentation scheme is a ¹-bundle. Finally, using Abel maps of bidegree (m, 1), we characterize the curves having double points at worst.     

Homogeneous Stäckel-type systems

A class of Hamiltonian dynamic systems integrated by the variable separation method is considered. The integration for this class is the inversion of an Abel mapping on hyperelliptic curves. We prove that the derivative of the Abel mapping is the Stäckel matrix, which determines a diagonal Riemannian metric and curvilinear orthogonal coordinate systems in a flat space. Lax representations with the spectral parameter are constructed. The corresponding classicalr-matrices are dynamic. It is shown how the class of pointwise canonical transformations can be naturally generalized using the Abel integral reduction theory.     

Some analysis of a stochastic logistic growth model

We derive several new results on a well-known stochastic logistic equation. For the martingale case, we compute the distribution of the solution, mean passage times, and the distribution of hitting times, all in closed form. For the case of constant coefficients, we also find mean passage times and for the general equation we give the weak solution expressed in terms of stochastic quadratures. We also show how these quadratures may be considerably simplified using the results for the martingale case. As it turns out, the martingale case has a particularly elegant weak solution, and to a large degree its structure carries over to the general case.     

On the solution of the Abel equation of the second kind by the shifted Chebyshev polynomials

This paper presents a new approximate method of Abel differential equation. By using the shifted Chebyshev expansion of the unknown function, Abel differential equation is approximately transformed to a system of nonlinear equations for the unknown coefficients. A desired solution can be determined by solving the resulting nonlinear system. This method gives a simple and closed form of approximate solution of Abel differential equation. The solution is calculated in the form of a series with easily computable components. The numerical results show the effectiveness of the method for this type of equation. Comparing the methodology with some known techniques shows that the present approach is relatively easy and highly accurate.     

On Generalized Gaussian Quadratures for Exponentials and Their Applications

We introduce new families of Gaussian-type quadratures for weighted integrals of exponential functions and consider their applications to integration and interpolation of bandlimited functions.We use a generalization of a representation theorem due to Carathéodory to derive these quadratures. For each positive measure, the quadratures are parameterized by eigenvalues of the Toeplitz matrix constructed from the trigonometric moments of the measure. For a given accuracy ε, selecting an eigenvalue close to ε yields an approximate quadrature with that accuracy. To compute its weights and nodes, we present a new fast algorithm.These new quadratures can be used to approximate and integrate bandlimited functions, such as prolate spheroidal wave functions, and essentially bandlimited functions, such as Bessel functions. We also develop, for a given precision, an interpolating basis for bandlimited functions on an interval.     

Closed-form solutions of certain Abel equations of the first kind

By using appropriate transformations in combination with specific Abel equations solvable in closed form containing arbitrary functions, an implicit solution as well as the associated sufficient condition are derived for certain differential equations of the Abel class of the first kind.     

Generalized Abel type integral equations with localized fractional integrals and derivatives

Generalized Abel type integral equations with Gauss, Kummer's and Humbert's confluent hypergeometric functions in the kernel and generalized Abel type integral equations with localized fractional integrals are considered. The left-hand sides of these equations are inversed by using generalized fractional derivatives. Explicit solutions of the equations in the class of locally summable functions are obtained. They are represented in terms of hypergeometric functions. Asymptotic power exponential type expansions of the generalized and localized fractional integrals are obtained. The base solutions of the generalized Abel type integral equation are given in the form of asymptotic series.     

On computing quadrature-based bounds for the A-norm of the error in conjugate gradients

In their original paper, Golub and Meurant (BIT 37:687–705, 1997) suggest to compute bounds for the A-norm of the error in the conjugate gradient (CG) method using Gauss, Gauss-Radau and Gauss-Lobatto quadratures. The quadratures are computed using the (1,1)-entry of the inverse of the corresponding Jacobi matrix (or its rank-one or rank-two modifications). The resulting algorithm called CGQL computes explicitly the entries of the Jacobi matrix and its modifications from the CG coefficients. In this paper, we use the fact that CG computes the Cholesky decomposition of the Jacobi matrix which is given implicitly. For Gauss-Radau and Gauss-Lobatto quadratures, instead of computing the entries of the modified Jacobi matrices, we directly compute the entries of the Cholesky decompositions of the (modified) Jacobi matrices. This leads to simpler formulas in comparison to those used in CGQL.     

Higher-order quadratures for circulant preconditioned Wiener-Hopf equations

In this paper, we consider solving matrix systems arising from the discretization of Wiener-Hopf equations by preconditioned conjugate gradient (PCG) methods. Circulant integral operators as preconditioners have been proposed and studied. However, the discretization of these preconditioned equations by employing higher-order quadratures leads to matrix systems that cannot be solved efficiently by using fast Fourier transforms (FFTs). The aim of this paper is to propose new preconditioners for Wiener-Hopf equations. The discretization of these preconditioned operator equations by higher-order quadratures leads to matrix systems that involve only Toeplitz, circulant and diagonal matrix-vector multiplications and hence can be computed efficiently by FFTs in each iteration. We show that with the proper choice of kernel functions of Wiener-Hopf equations, the resulting preconditioned operators will have clustered spectra and therefore the PCG method converges very fast. Numerical examples are given to illustrate the fast convergence of the method and the improvement of the accuracy of the computed solutions with using higher-order quadratures.Research supported by the Cooperative Research Centre for Advanced Computational Systems.Research supported in part by Lee Ka Shing scholarship.     

级数的常规可和，Cesàro可和与Abel可和的几点讨论

讨论级数常规可和、Cesaro可和与Abel可和的关系．利用数学分析级数理论,证明Abel可和适用范围最广,Cesaro可和其次,级数常规可和适用范围最小．这个结论丰富了经典级数理论,为实际应用中选用合适可和提供依据．     

Computing Discrepancies of Smolyak Quadrature Rules

In recent years, Smolyak quadrature rules (also called quadratures on hyperbolic cross points or sparse grids) have gained interest as possible competition to number theoretic quadratures for high dimensional problems. A standard way of comparing the quality of multivariate quadrature formulas consists in computing theirL₂-discrepancy. Especially for larger dimensions, such computations are a highly complex task. In this paper we develop a fast recursive algorithm for computing theL₂-discrepancy (and related quality measures) of general Smolyak quadratures. We carry out numerical comparisons between the discrepancies of certain Smolyak rules and Hammersley and Monte Carlo sequences.     

An error analysis of Runge–Kutta convolution quadrature

An error analysis is given for convolution quadratures based on strongly A-stable Runge–Kutta methods, for the non-sectorial case of a convolution kernel with a Laplace transform that is polynomially bounded in a half-plane. The order of approximation depends on the classical order and stage order of the Runge–Kutta method and on the growth exponent of the Laplace transform. Numerical experiments with convolution quadratures based on the Radau IIA methods are given on an example of a time-domain boundary integral operator.     

High-order quadratures for integral operators with singular kernels

A numerical integration method that has rapid convergence for integrands with known singularities is presented. Based on endpoint corrections to the trapezoidal rule, the quadratures are suited for the discretization of a variety of integral equations encountered in mathematical physics. The quadratures are based on a technique introduced by Rokhlin (1990). The present modification controls the growth of the quadrature weights and permits higher-order rules in practice. Several numerical examples are included.     

Abel’s theorem and Bäcklund transformations for the Hamilton-Jacobi equations

We consider an algorithm for constructing auto-Bäcklund transformations for finitedimensional Hamiltonian systems whose integration reduces to the inversion of the Abel map. In this case, using equations of motion, one can construct Abel differential equations and identify the sought Bäcklund transformation with the well-known equivalence relation between the roots of the Abel polynomial. As examples, we construct Bäcklund transformations for the Lagrange top, Kowalevski top, and Goryachev–Chaplygin top, which are related to hyperelliptic curves of genera 1 and 2, as well as for the Goryachev and Dullin–Matveev systems, which are related to trigonal curves in the plane.     

Some relations between Nörlund and Abel summability

The definition of generalised Abel summability is extended to positive orders and a definition of strong generalised Abel summability is introduced. A result of Jurkat and Peyerimhoff concerning the implication between Nörlund and generalised Abel summabilities is extended to positive integral order. Analogous results for the corresponding strong methods are also given.

     

三维洛仑兹空间形式中的类光曲线

讨论三维洛仑兹空间形式中的类光曲线粒子模型,研究依赖于粒子轨道的Cartan曲率的作用,找出了沿着极值曲线的Killing场,通过Killing场构造合适的坐标系,用积分求出极值曲线的参数表达式.     

Quadrature formulas on the unit circle based on rational functions

Quadrature formulas on the unit circle were introduced by Jones in 1989. On the other hand, Bultheel also considered such quadratures by giving results concerning error and convergence. In other recent papers, a more general situation was studied by the authors involving orthogonal rational functions on the unit circle which generalize the well-known Szeg polynomials. In this paper, these quadratures are again analyzed and results about convergence given. Furthermore, an application to the Poisson integral is also made.     

The second Hamiltonian structure for a special case of the Lotka-Volterra equations

A special case of the Lotka-Volterra equations is considered for which it is possible to find the second Hamiltonian structure that is complementary to the known one. The form of the new Hamiltonian makes it possible to solve the equations by quadratures, which is the main feature of the case under examination. As a consequence, the period can also be represented by quadratures. In terms of the new variables, the equations of motion admit a mechanical analogy with the oscillations of a mass on a nonlinear spring.     

Convergence of Newton-Cotes quadratures for analytic functions

We determine (Theorem 3) the smallest closed region, containing the interva of integration, such that the analyticity of the integrand in this closed region implies the convergence of the Newton-Cotes quadratures. By considering, in particular, certain ellipses as regions of analyticity, we obtain (Theorem 4) an improvement of Davis' result on the convergence of Newton-Cotes quadratures for analytic functions.     

On quadratures related to Padé approximants for the inversion of the Laplace transform

Properties of the quadratures for the numerical inversion of the Laplace transform generated by Padé approximants of the exponential function are examined. In particular, quadratures of the highest possible degree of accuracy are considered.