10,487 research outputs found
Fast hyperbolic Radon transform represented as convolutions in log-polar coordinates
The hyperbolic Radon transform is a commonly used tool in seismic processing,
for instance in seismic velocity analysis, data interpolation and for multiple
removal. A direct implementation by summation of traces with different moveouts
is computationally expensive for large data sets. In this paper we present a
new method for fast computation of the hyperbolic Radon transforms. It is based
on using a log-polar sampling with which the main computational parts reduce to
computing convolutions. This allows for fast implementations by means of FFT.
In addition to the FFT operations, interpolation procedures are required for
switching between coordinates in the time-offset; Radon; and log-polar domains.
Graphical Processor Units (GPUs) are suitable to use as a computational
platform for this purpose, due to the hardware supported interpolation routines
as well as optimized routines for FFT. Performance tests show large speed-ups
of the proposed algorithm. Hence, it is suitable to use in iterative methods,
and we provide examples for data interpolation and multiple removal using this
approach.Comment: 21 pages, 10 figures, 2 table
Probability Distributions and Hilbert Spaces: Quantum and Classical Systems
We use the fact that some linear Hamiltonian systems can be considered as
``finite level'' quantum systems, and the description of quantum mechanics in
terms of probabilities, to associate probability distributions with this
particular class of linear Hamiltonian systems.Comment: LATEX,13pages,accepted by Physica Scripta (1999
Review of Summation-by-parts schemes for initial-boundary-value problems
High-order finite difference methods are efficient, easy to program, scales
well in multiple dimensions and can be modified locally for various reasons
(such as shock treatment for example). The main drawback have been the
complicated and sometimes even mysterious stability treatment at boundaries and
interfaces required for a stable scheme. The research on summation-by-parts
operators and weak boundary conditions during the last 20 years have removed
this drawback and now reached a mature state. It is now possible to construct
stable and high order accurate multi-block finite difference schemes in a
systematic building-block-like manner. In this paper we will review this
development, point out the main contributions and speculate about the next
lines of research in this area
- …