Computational properties of three-term recurrence relations for Kummer functions

Several three-term recurrence relations for confluent hypergeometric functions are analyzed from a numerical point of view. Minimal and dominant solutions for complex values of the variable z are given, derived from asymptotic estimates of the Whittaker functions with large parameters. The Laguerre polynomials and the regular Coulomb wave functions are studied as particular cases, with numerical examples of their computation.     

<Emphasis Type="Italic">N</Emphasis>=1 supersymmetric conformal block recursion relations

We present explicit recursion relations for the four-point superconformal block functions that are essentially particular contributions of the given conformal class to the four-point correlation function. The approach is based on the analytic properties of the superconformal blocks as functions of the conformal dimensions and the central charge of the superconformal algebra. We compare the results with the explicit analytic expressions obtained for special parameter values corresponding to the truncated operator product expansion. These recursion relations are an efficient tool for numerically studying the four-point correlation function in superconformal field theory in the framework of the bootstrap approach, similar to that in the case of the purely conformal symmetry. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 152, No. 3, pp. 476–487, September, 2007.     

The Ackermann function. a theoretical,computational, and formula manipulative study

Ackermann's function is of highly recursive nature and of two arguments. It is here treated as a class of functions of one argument, where the other argument defines the member of the class. The first members are expressed with elementary functions, the higher members with a hierarchy of primitive recursive functions. The number of calls of the function needed in a straightforward recursive computation is given for the first members. The maximum depth in the recursion during the evaluation is investigated.Results from tests with the Ackermann function of recursive procedure implementations in ALGOL-60, ALGOL W, PL/I and SIMULA-67 on IBM 360/75 and CD 6600 are given.A SYMBAL formula manipulating program, that automatically solves recurrence relations for the first members of the function class and for the number of calls needed in their straightforward computation, is given.The Ackermann rating of programming languages is discussed.     

A stress recovery procedure for solving geometrically non-linear problems in the mechanics of a deformable solid by the finite element method

A stress recovery procedure is presented for non-linear and linearized problems, based on the determination of the forces at the mesh points using a stiffness matrix obtained by the finite element method for the Lagrange variational equation written in the initial configuration using an asymmetric Piola–Kirchhoff stress tensor. Vectors of the forces reduced to the mesh points are constructed using the displacements at the mesh points found by solving this equation and for the known stiffness matrices of the elements. On the other hand, these forces at the mesh points are defined in terms of unknown forces distributed over the surface of an element and given shape functions. As a result, a system of Fredholm integral equations of the first kind is obtained, the solution of which gives these distributed forces. The values of the Piola–Kirchhoff stress tensor of the first kind at the mesh points are determined using the values found for the distributed forces on the surfaces of the finite element mesh (including at the mesh points) using the Cauchy relations for the initial configuration. The linearized representation of this tensor enables all the derivatives of the increment in the strain vector with respect to the coordinates to be found without invoking the operation of differentiation. The particular features of the use of the stress recovery procedure are demonstrated for a plane problem in the non-linear theory of elasticity.     

Functional relations in lattice statistical mechanics, enumerative combinatorics, and discrete dynamical systems

We recall some non-trivial, non-linear functional relations appearing in various domains of mathematics and physics, such as lattice statistical mechanics, quantum mechanics, or enumerative combinatorics. We focus, more particularly, on the analyticity properties of the solutions of these functional relations. We then consider discrete dynamical systems corresponding to birational transformations. The rational expressions for dynamical zeta functions obtained for a particular two-dimensional birational mapping, depending on two parameters, are recalled, as well as some non-trivial functional relations satisfied by these dynamical zeta functions. We finally give some functional equations corresponding to some singled out orbits of this two-dimensional birational mapping for particular values of the two parameters. This example shows that functional equations associated with curves, for real values of the variables, are actually compatible with a chaotic dynamical system.     

A local search method for optimization problem with d.c. inequality constraints

This paper addresses a nonconvex optimization problem with the cost function and inequality constraints given by d.c. functions. The original problem is reduced to a problem without inequality constraints by the exact penalization procedure. A special local search method for the penalized problem is developed, which is based, first, on the linearization procedure with respect to the basic nonconvexity and, second, on the consecutive solutions of linearized convex problems. Convergence properties of the method are investigated. In particular, it is shown that a limit point of the sequence produced by the method is considerably stronger than the usual KKT-vector.In addition, the relations between an approximate solution of linearized convex problem and the KKT-vector of the original problem are established, and the various stopping criteria are substantiated. Besides, we established the relations among the Lagrange multipliers of the original problem, those ones of the linearized problem, and the value of the penalty parameter. Finally, a preliminary computational testing of the LSM developed has been carried out on several test problems taken from literature.     

Riemann type integrals for functions taking values in a locally convex space

The McShane and Kurzweil-Henstock integrals for functions taking values in a locally convex space are defined and the relations with other integrals are studied. A characterization of locally convex spaces in which Henstock Lemma holds is given.     

Two new multivariable generating relations

A multidimensional extension of Bailey's transform is utilised to deduce two new generating relations of quite a general character. These expressions are then specialised to give more practical formulae in terms of Karlsson's generalised Kampé de Fériet functions which embody very many generating relations. A number of interesting special cases are given in an appendix which includes results involving Lauricella polynomials, generalised hypergeometric polynomials and the polynomials of Meixner, Charlier and Laguerre.     

Recurrence relations for discrete hypergeometric functions

We present a general procedure for finding linear recurrence relations for the solutions of the second order difference equation of hypergeometric type. Applications to wave functions of certain discrete system are also given.     

Quasi-arithmetic means of covariance functions with potential applications to space–time data

The theory of quasi-arithmetic means represents a powerful tool in the study of covariance functions across space–time. In the present study we use quasi-arithmetic functionals to make inferences about the permissibility of averages of functions that are not, in general, permissible covariance functions. This is the case, e.g., of the geometric and harmonic averages, for which we obtain permissibility criteria. Also, some important inequalities involving covariance functions and preference relations as well as algebraic properties can be derived by means of the proposed approach. In particular, quasi-arithmetic covariances allow for ordering and preference relations, for a Jensen-type inequality and for a minimal and maximal element of their class. The general results shown in this paper are then applied to the study of spatial and spatio-temporal random fields. In particular, we discuss the representation and smoothness properties of a weakly stationary random field with a quasi-arithmetic covariance function. Also, we show that the generator of the quasi-arithmetic means can be used as a link function in order to build a space–time nonseparable structure starting from the spatial and temporal margins, a procedure that is technically sound for those working with copulas. Several examples of new families of stationary covariances obtainable with this procedure are shown. Finally, we use quasi-arithmetic functionals to generalise existing results concerning the construction of nonstationary spatial covariances, and discuss the applicability and limits of this generalisation.     

Two new multivariable generating relations

A multidimensional extension of Bailey's transform is utilised to deduce two new generating relations of quite a general character. These expressions are then specialised to give more practical formulae in terms of Karlsson's generalised Kampé de Fériet functions which embody very many generating relations. A number of interesting special cases are given in an appendix which includes results involving Lauricella polynomials, generalised hypergeometric polynomials and the polynomials of Meixner, Charlier and Laguerre.     

On the numerical solution of ill‐conditioned linear systems by regularization and iteration

We propose to reduce the (spectral) condition number of a given linear system by adding a suitable diagonal matrix to the system matrix, in particular by shifting its spectrum. Iterative procedures are then adopted to recover the solution of the original system. The case of real symmetric positive definite matrices is considered in particular, and several numerical examples are given. This approach has some close relations with Riley's method and with Tikhonov regularization. Moreover, we identify approximately the aforementioned procedure with a true action of preconditioning.     

A superharmonic approach to solving infinite horizon partially observable Markov decision problems

In this paper we use an approach which uses a superharmonic property of a sequence of functions generated by an algorithm to show that these functions converge in a non-increasing manner to the optimal value function for our problem, and bounds are given for the loss of optimality if the computational process is terminated at any iteration. The basic procedure is to add an additional linear term at each iteration, selected by solving a particular optimisation problem, for which primal and dual linear programming formulations are given.     

An Extremal Problem Arising in the Dynamics of Two-Phase Materials That Directly Reveals Information about the Internal Geometry

In two-phase materials, each phase having a non-local response in time, it has been found that for some driving fields the response somehow untangles at specific times, and allows one to directly infer useful information about the geometry of the material, such as the volume fractions of the phases. Motivated by this, and to obtain an algorithm for designing appropriate driving fields, we find approximate, measure independent, linear relations between the values that Markov functions take at a given set of possibly complex points, not belonging to the interval [-1,1] where the measure is supported. The problem is reduced to simply one of polynomial approximation of a given function on the interval [-1,1] and, to simplify the analysis, Chebyshev approximation is used. This allows one to obtain explicit estimates of the error of the approximation, in terms of the number of points and the minimum distance of the points to the interval [-1,1]. Assuming this minimum distance is bounded below by a number greater than 1/2, the error converges exponentially to zero as the number of points is increased. Approximate linear relations are also obtained that incorporate a set of moments of the measure. In the context of the motivating problem, the analysis also yields bounds on the response at any particular time for any driving field, and allows one to estimate the response at a given frequency using an appropriately designed driving field that effectively is turned on only for a fixed interval of time. The approximation extends directly to Markov-type functions with a positive semidefinite operator valued measure, and this has applications to determining the shape of an inclusion in a body from boundary flux measurements at a specific time, when the time-dependent boundary potentials are suitably tailored. © 2022 Wiley Periodicals, Inc.     

Indépendance algébrique sur les T-modules

We study the transcendence degree (of the values at algebraic points) of the coordinate functions of a given one parameter subgroup on a simple g dimensional T-module. In such a situation, we obtain a lower bound for this transcendence degree depending on g and on the growth order of the parameter subgroup. As a particular case, one gets the finite characteristic analogue of a result of Siegel on the algebraic independence for the values of the classical Bessel function. Résumé. Nous obtenons une minoration du degré de transcendance d'un ensemble de valeurs en des points algébriques des fonctions coordonnées d'un sous-groupe à un paramètre d'un T-module simple, en fonction de l'ordre de croissance de ces dernières et de la dimension du T-module. Un cas particulier de ce résultat permet d'obtenir en caractéristique finie un analogue d'un résultat de Siegel sur l'indépendance algébrique des valeurs de la fonction de Bessel usuelle.     

Stress recovery procedure for solving boundary value problems in the mechanics of a deformable solid by the finite element method

A stress recovery procedure, based on the determination of the forces at the mesh points using a stiffness matrix obtained by the finite element method for the variational Lagrange equation, is described. The vectors of the forces reduced to the mesh points are constructed for the known stiffness matrices of the elements using the displacements at the mesh points found from the solution of the problem. On the other hand, these mesh point forces are determined in terms of the unknown forces distributed over the surface of an element and given shape functions. As a result, a system of Fredholm integral equations of the first kind is obtained, the solution of which gives these distributed forces. The stresses at the mesh points are determined for the values of these forces found on the surfaces of the finite element mesh (including at the mesh points) using the Cauchy relations, which relate the forces, stresses and the normal to the surface. The special features of the use of the stress recovery procedure are demonstrated for a plane problem in the linear theory of elasticity.     

Nonlocal hydrodynamics in the quantum field modelϕ 4

The procedure of passing from quantum statistical mechanics to the hydrodynamics previously developed by the author is now applied to the quantum field model ⁴. For a certain class of external forces, the equations of many-body systems in quantum theory appear to be equivalent to the equations of nonlocal hydrodynamics. The hydrodynamic nonlocalities arising in constituent relations are expressed through the Green's functions for currents. Some properties of the nonlocal kernels, in particular, the conditions related to dissipation and T-invariance of the model ⁴ (an analogue of Onsager's relations), are deduced from the general symmetry properties. In hydrodynamics, nonlocality allows causality and dissipativity to be consistently combined. The connection between the classical transport coefficients and the hydrodynamic kernels is established. An algorithm for calculating constituent relations by perturbation theory, using the technique of temperature Green's functions, is described.Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 108, No. 1, pp. 50–68, July, 1996.     

Pontryagin-de Branges-Rovnyak spaces of slice hyperholomorphic functions

We study reproducing kernel Hilbert and Pontryagin spaces of slice hyperholomorphic functions. These are analogs of the Hilbert spaces of analytic functions introduced by de Branges and Rovnyak. In the first part of the paper, we focus on the case of Hilbert spaces and introduce, in particular, a version of the Hardy space. Then we define Blaschke factors and Blaschke products and consider an interpolation problem. In the second part of the paper, we turn to the case of Pontryagin spaces. We first prove some results from the theory of Pontryagin spaces in the quaternionic setting and, in particular, a theorem of Shmulyan on densely defined contractive linear relations. We then study realizations of generalized Schur functions and of generalized Carathéodory functions.     

Gaussian Quadrature Rules with Simple Node-Weight Relations

In this paper, we consider the symmetric Gaussian and L-Gaussian quadrature rules associated with twin periodic recurrence relations with possible variations in the initial coefficient. We show that the weights of the associated Gaussian quadrature rules can be given as rational functions in terms of the corresponding nodes where the numerators and denominators are polynomials of degree at most 4. We also show that the weights of the associated L-Gaussian quadrature rules can be given as rational functions in terms of the corresponding nodes where the numerators and denominators are polynomials of degree at most 5. Special cases of these quadrature rules are given. Finally, an easy to implement procedure for the evaluation of the nodes is described.