scholarly journals Fast algorithms for Jacobi expansions via nonoscillatory phase functions

2019 ◽  
Vol 40 (3) ◽  
pp. 2019-2051
Author(s):  
James Bremer ◽  
Haizhao Yang

Abstract We describe a suite of fast algorithms for evaluating Jacobi polynomials, applying the corresponding discrete Sturm–Liouville eigentransforms and calculating Gauss–Jacobi quadrature rules. Our approach, which applies in the case in which both of the parameters $\alpha $ and $\beta $ in Jacobi’s differential equation are of magnitude less than $1/2$, is based on the well-known fact that in this regime Jacobi’s differential equation admits a nonoscillatory phase function that can be loosely approximated via an affine function over much of its domain. We illustrate this with several numerical experiments, the source code for which is publicly available.

Mathematics ◽  
2021 ◽  
Vol 9 (13) ◽  
pp. 1573
Author(s):  
Waleed Mohamed Abd-Elhameed ◽  
Badah Mohamed Badah

This article deals with the general linearization problem of Jacobi polynomials. We provide two approaches for finding closed analytical forms of the linearization coefficients of these polynomials. The first approach is built on establishing a new formula in which the moments of the shifted Jacobi polynomials are expressed in terms of other shifted Jacobi polynomials. The derived moments formula involves a hypergeometric function of the type 4F3(1), which cannot be summed in general, but for special choices of the involved parameters, it can be summed. The reduced moments formulas lead to establishing new linearization formulas of certain parameters of Jacobi polynomials. Another approach for obtaining other linearization formulas of some Jacobi polynomials depends on making use of the connection formulas between two different Jacobi polynomials. In the two suggested approaches, we utilize some standard reduction formulas for certain hypergeometric functions of the unit argument such as Watson’s and Chu-Vandermonde identities. Furthermore, some symbolic algebraic computations such as the algorithms of Zeilberger, Petkovsek and van Hoeij may be utilized for the same purpose. As an application of some of the derived linearization formulas, we propose a numerical algorithm to solve the non-linear Riccati differential equation based on the application of the spectral tau method.


2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Zohreh Zeinalabedini Charandabi ◽  
Hakimeh Mohammadi ◽  
Shahram Rezapour ◽  
Hashem Parvaneh Masiha

AbstractThe Sturm–Liouville differential equation is one of interesting problems which has been studied by researchers during recent decades. We study the existence of a solution for partial fractional Sturm–Liouville equation by using the α-ψ-contractive mappings. Also, we give an illustrative example. By using the α-ψ-multifunctions, we prove the existence of solutions for inclusion version of the partial fractional Sturm–Liouville problem. Finally by providing another example and some figures, we try to illustrate the related inclusion result.


2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Guanglu Zhou ◽  
Boying Wu ◽  
Wen Ji ◽  
Seungmin Rho

This study presents numerical schemes for solving a parabolic partial differential equation with a time- or space-dependent coefficient subject to an extra measurement. Through the extra measurement, the inverse problem is transformed into an equivalent nonlinear equation which is much simpler to handle. By the variational iteration method, we obtain the exact solution and the unknown coefficients. The results of numerical experiments and stable experiments imply that the variational iteration method is very suitable to solve these inverse problems.


2019 ◽  
Vol 16 (09) ◽  
pp. 1950140 ◽  
Author(s):  
L. C. N. Santos ◽  
C. C. Barros

In this paper, we study the Dirac equation in the Rindler spacetime. The solution of the wave equation in an accelerated reference frame is obtained. The differential equation associated to this wave equation is mapped into a Sturm–Liouville problem of a Schrödinger-like equation. We derive a compact expression for the energy spectrum associated with the Dirac equation in an accelerated reference. It is shown that the noninertial effect of the accelerated reference frame mimics an external potential in the Dirac equation and, moreover, allows the formation of bound states.


2020 ◽  
Author(s):  
Eddie C. Davis

This research presents an intermediate compiler representation that is designed for optimization, and emphasizes the temporary storage requirements and execution schedule of a given computation to guide optimization decisions. The representation is expressed as a dataflow graph that describes computational statements and data mappings within the polyhedral compilation model. The targeted applications include both the regular and irregular scientific domains. The intermediate representation can be integrated into existing compiler infrastructures. A specification language implemented as a domain specific language in C++ describes the graph components and the transformations that can be applied. The visual representation allows users to reason about optimizations. Graph variants can be translated into source code or other representation. The language, intermediate representation, and associated transformations have been applied to improve the performance of differential equation solvers, or sparse matrix operations, tensor decomposition, and structured multigrid methods.


2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Zihan Li ◽  
Xiao-Bao Shu ◽  
Tengyuan Miao

AbstractIn this article, we consider the existence of solutions to the Sturm–Liouville differential equation with random impulses and boundary value problems. We first study the Green function of the Sturm–Liouville differential equation with random impulses. Then, we get the equivalent integral equation of the random impulsive differential equation. Based on this integral equation, we use Dhage’s fixed point theorem to prove the existence of solutions to the equation, and the theorem is extended to the general second order nonlinear random impulsive differential equations. Then we use the upper and lower solution method to give a monotonic iterative sequence of the generalized random impulsive Sturm–Liouville differential equations and prove that it is convergent. Finally, we give two concrete examples to verify the correctness of the results.


2008 ◽  
Vol 45 (2) ◽  
pp. 347-362 ◽  
Author(s):  
Saul C. Leite ◽  
Marcelo D. Fragoso

This paper is concerned with the characterization of weak-sense limits of state-dependent G-networks under heavy traffic. It is shown that, for a certain class of networks (which includes a two-layer feedforward network and two queues in tandem), it is possible to approximate the number of customers in the queue by a reflected stochastic differential equation. The benefits of such an approach are that it describes the transient evolution of these queues and allows the introduction of controls, inter alia. We illustrate the application of the results with numerical experiments.


Sign in / Sign up

Export Citation Format

Share Document