7,494 research outputs found
The fractional orthogonal derivative
This paper builds on the notion of the so-called orthogonal derivative, where
an n-th order derivative is approximated by an integral involving an orthogonal
polynomial of degree n. This notion was reviewed in great detail in a paper in
J. Approx. Theory (2012) by the author and Koornwinder. Here an approximation
of the Weyl or Riemann-Liouville fractional derivative is considered by
replacing the n-th derivative by its approximation in the formula for the
fractional derivative. In the case of, for instance, Jacobi polynomials an
explicit formula for the kernel of this approximate fractional derivative can
be given. Next we consider the fractional derivative as a filter and compute
the transfer function in the continuous case for the Jacobi polynomials and in
the discrete case for the Hahn polynomials. The transfer function in the Jacobi
case is a confluent hypergeometric function. A different approach is discussed
which starts with this explicit transfer function and then obtains the
approximate fractional derivative by taking the inverse Fourier transform. The
theory is finally illustrated with an application of a fractional
differentiating filter. In particular, graphs are presented of the absolute
value of the modulus of the transfer function. These make clear that for a good
insight in the behavior of a fractional differentiating filter one has to look
for the modulus of its transfer function in a log-log plot, rather than for
plots in the time domain.Comment: 32 pages, 7 figures. The section between formula (4.15) and (4.20) is
correcte
Mathematical computer programs: A compilation
Computer programs, routines, and subroutines for aiding engineers, scientists, and mathematicians in direct problem solving are presented. Also included is a group of items that affords the same users greater flexibility in the use of software
Differentiation by integration using orthogonal polynomials, a survey
This survey paper discusses the history of approximation formulas for n-th
order derivatives by integrals involving orthogonal polynomials. There is a
large but rather disconnected corpus of literature on such formulas. We give
some results in greater generality than in the literature. Notably we unify the
continuous and discrete case. We make many side remarks, for instance on
wavelets, Mantica's Fourier-Bessel functions and Greville's minimum R_alpha
formulas in connection with discrete smoothing.Comment: 35 pages, 3 figures; minor corrections, subsection 3.11 added;
accepted by J. Approx. Theor
Rotation method for accelerating multiple-spherical Bessel function integrals against a numerical source function
A common problem in cosmology is to integrate the product of two or more
spherical Bessel functions (sBFs) with different configuration-space arguments
against the power spectrum or its square, weighted by powers of wavenumber.
Naively computing them scales as with the number of
configuration space arguments and the grid size, and they cannot be
done with Fast Fourier Transforms (FFTs). Here we show that by rewriting the
sBFs as sums of products of sine and cosine and then using the product to sum
identities, these integrals can then be performed using 1-D FFTs with scaling. This "rotation" method has the potential to
accelerate significantly a number of calculations in cosmology, such as
perturbation theory predictions of loop integrals, higher order correlation
functions, and analytic templates for correlation function covariance matrices.
We implement this approach numerically both in a free-standing,
publicly-available \textsc{Python} code and within the larger,
publicly-available package \texttt{mcfit}. The rotation method evaluated with
direct integrations already offers a factor of 6-10 speed-up over the
naive approach in our test cases. Using FFTs, which the rotation method
enables, then further improves this to a speed-up of
over the naive approach. The rotation method should be useful in light of
upcoming large datasets such as DESI or LSST. In analysing these datasets
recomputation of these integrals a substantial number of times, for instance to
update perturbation theory predictions or covariance matrices as the input
linear power spectrum is changed, will be one piece in a Monte Carlo Markov
Chain cosmological parameter search: thus the overall savings from our method
should be significant
Formalization of Transform Methods using HOL Light
Transform methods, like Laplace and Fourier, are frequently used for
analyzing the dynamical behaviour of engineering and physical systems, based on
their transfer function, and frequency response or the solutions of their
corresponding differential equations. In this paper, we present an ongoing
project, which focuses on the higher-order logic formalization of transform
methods using HOL Light theorem prover. In particular, we present the
motivation of the formalization, which is followed by the related work. Next,
we present the task completed so far while highlighting some of the challenges
faced during the formalization. Finally, we present a roadmap to achieve our
objectives, the current status and the future goals for this project.Comment: 15 Pages, CICM 201
Noise in a Calorimeter Readout System Using Periodic Sampling
Fourier transform analysis of the calorimeter noise problem gives
quantitative results on a) the time-height correlation, b) the effect of
background on optimal shaping and on the ENC, c) sampling frequency
requirements, and d) the relation between sampling frequency and the required
quantization error
Quantum Computation of Scattering in Scalar Quantum Field Theories
Quantum field theory provides the framework for the most fundamental physical theories to be confirmed experimentally, and has enabled predictions of unprecedented precision. However, calculations of physical observables often require great computational complexity and can generally be performed only when the interaction strength is weak. A full understanding of the foundations and rich consequences of quantum field theory remains an outstanding challenge. We develop a quantum algorithm to compute relativistic scattering amplitudes in massive phi-fourth theory in spacetime of four and fewer dimensions. The algorithm runs in a time that is polynomial in the number of particles, their energy, and the desired precision, and applies at both weak and strong coupling. Thus, it offers exponential speedup over existing classical methods at high precision or strong coupling
Introduction to Quantum Algorithms for Physics and Chemistry
In this introductory review, we focus on applications of quantum computation
to problems of interest in physics and chemistry. We describe quantum
simulation algorithms that have been developed for electronic-structure
problems, thermal-state preparation, simulation of time dynamics, adiabatic
quantum simulation, and density functional theory.Comment: 44 pages, 5 figures; comments or suggestions for improvement are
welcom
- …