12 research outputs found
Numerical methods for the computation of the confluent and Gauss hypergeometric functions
The two most commonly used hypergeometric functions are the confluent hypergeometric function and the Gauss hypergeometric function. We review the available techniques for accurate, fast, and reliable computation of these two hypergeometric functions in different parameter and variable regimes. The methods that we investigate include Taylor and asymptotic series computations, Gauss-Jacobi quadrature, numerical solution of differential equations, recurrence relations, and others. We discuss the results of numerical experiments used to determine the best methods, in practice, for each parameter and variable regime considered. We provide 'roadmaps' with our recommendation for which methods should be used in each situation
On the computation of confluent hypergeometric functions for large imaginary part of parameters b and z
The final publication is available at http://link.springer.com/chapter/10.1007%2F978-3-319-42432-3_30We present an efficient algorithm for the confluent hypergeometric functions when the imaginary part of b and z is large. The algorithm is based on the steepest descent method, applied to a suitable representation of the confluent hypergeometric functions as a highly oscillatory integral, which is then integrated by using various quadrature methods. The performance of the algorithm is compared with open-source and commercial software solutions with arbitrary precision, and for many cases the algorithm achieves high accuracy in both the real and imaginary parts. Our motivation comes from the need for accurate computation of the characteristic function of the Arcsine distribution or the Beta distribution; the latter being required in several financial applications, for example, modeling the loss given default in the context of portfolio credit risk.Peer ReviewedPostprint (author's final draft
GPU-SUPPORTED SIMULATION FOR ABEP AND QOS ANALYSIS OF A COMBINED MACRO DIVERSITY SYSTEM IN A GAMMA-SHADOWED K-µ FADING CHANNEL
In this paper we have analyzed macro-diversity (MD) system with one macro SC diversity (MD SC) receiver and two micro MRC (mD MRC) receivers over correlated Gamma-shadowed k-µ fading channel. The average bit error probability (ABEP) is calculated using the moment generating function (MGF) approach for BDPSK and BPSK modulations. Graphical representation of the results illustrates the effects of different parameters of the system on its performance as well as the improvements due to the benefits of a combined micro and macro diversity. The obtained analytical expressions are used for the GPU-enabled mobile network modeling, planning and simulation environment to determine the value of Quality of Service (QoS) parameter. Finally, linear optimization is proposed as an approach to improve the QoS parameter of the fading-affected system observed in this paper
Computing hypergeometric functions rigorously
We present an efficient implementation of hypergeometric functions in
arbitrary-precision interval arithmetic. The functions , ,
and (or the Kummer -function) are supported for
unrestricted complex parameters and argument, and by extension, we cover
exponential and trigonometric integrals, error functions, Fresnel integrals,
incomplete gamma and beta functions, Bessel functions, Airy functions, Legendre
functions, Jacobi polynomials, complete elliptic integrals, and other special
functions. The output can be used directly for interval computations or to
generate provably correct floating-point approximations in any format.
Performance is competitive with earlier arbitrary-precision software, and
sometimes orders of magnitude faster. We also partially cover the generalized
hypergeometric function and computation of high-order parameter
derivatives.Comment: v2: corrected example in section 3.1; corrected timing data for case
E-G in section 8.5 (table 6, figure 2); adjusted paper siz
Non-Gaussian Geostatistical Modeling using (skew) t Processes
We propose a new model for regression and dependence analysis when addressing
spatial data with possibly heavy tails and an asymmetric marginal distribution.
We first propose a stationary process with marginals obtained through scale
mixing of a Gaussian process with an inverse square root process with Gamma
marginals. We then generalize this construction by considering a skew-Gaussian
process, thus obtaining a process with skew-t marginal distributions. For the
proposed (skew) process we study the second-order and geometrical
properties and in the case, we provide analytic expressions for the
bivariate distribution. In an extensive simulation study, we investigate the
use of the weighted pairwise likelihood as a method of estimation for the
process. Moreover we compare the performance of the optimal linear predictor of
the process versus the optimal Gaussian predictor. Finally, the
effectiveness of our methodology is illustrated by analyzing a georeferenced
dataset on maximum temperatures in Australi
Improving on adjusted R-squared
Multivariate analysis of psychological dat