Abstract:Purpose: The finite Hilbert transform (FHT) or inverse finite Hilbert transform (IFHT) is recently found to have some important applications in computerized tomography (CT) arena [1-6], where they are used to filter the derivatives of back-projected data in the chord-line based CT reconstruction algorithms. In this paper, we implemented, improved and validated a fast numerical solution to the FHT via a double exponential (DE) integration scheme. A same strategy can be used to compute IFHT.Methods: To overcome the Underflow of floating-point numbers, we first determined the range of variable transformation from the minimum positive value of single or double precision floating point number, the integration step can be further determined by the range of variable transformation and the integration level. Two functions with their known analytical FHTs are used to validate the implementation of the FHT via DE scheme. The surface map and 2D contour of the FHT transformation error with respect to integration level and the range of the variable transformation are used to numerically determine the optimal numbers for a fast FHT.Results: Given a specific precision, the lowest integration level and the optimal range of variable transformation, which are used to transform a signal with a certain degree Of fluctuation, can be numerically determined by the surface map and 2D contour of the standard deviation of transformation error. These two numbers can then be taken to efficiently compute the FHT for other signals with the same or less degree of fluctuation.Conclusions: The FHT via DE scheme and the numerical method to determine the integration level and the range of transformation can be used for fast FHT in certain applications, such as data filtering in chord-line based CT reconstruction algorithms.

A Fast Precise Algorithm for Duhamel Integrals

Efficient Evaluation of Double Surface Integrals in Time-Domain Integral Equation Formulations

Semi-analytic Fresnel Diffraction Calculation with Polynomial Decomposition

Approximate Floating-Point FFT Design with Wide Precision-Range and High Energy Efficiency.

A Fast Numerical Solution Method for Two Dimensional Fredholm Integral Equations of the Second Kind

The fast Euler-Maruyama method for solving multiterm Caputo fractional stochastic delay integro-differential equations

Fast Finite Hilbert Transform Via DE Quadrature Scheme

Fast IMEX Time Integration of Nonlinear Stiff Fractional Differential Equations

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

A Fast Fourier-Galerkin Method Solving Boundary Integral Equations for the Helmholtz Equation with Exponential Convergence

Fast Q1 finite element for two-dimensional integral fractional Laplacian

A Multi-Domain Spectral Method for Time-Fractional Differential Equations

Computation of highly oscillatory integrals using a Fourier extension approximation

Fast numerical schemes for nonlinear space-fractional multidelay reaction-diffusion equations by implicit integration factor methods

A fast time integral finite difference method for a space-time fractional FitzHugh-Nagumo monodomain model in irregular domains

A Systematic and Efficient Method to Compute Multi-Loop Master Integrals

Efficient Exponential Integrator Finite Element Method for Semilinear Parabolic Equations

A Robust and Efficient Composite Time Integration Algorithm for Nonlinear Structural Dynamic Analysis

Spectral Collocation Method for Fractional Differential/Integral Equations with Generalized Fractional Operator

A Parareal exponential integrator finite element method for linear parabolic equations