18,701 research outputs found
Shifted Laplacian multigrid for the elastic Helmholtz equation
The shifted Laplacian multigrid method is a well known approach for
preconditioning the indefinite linear system arising from the discretization of
the acoustic Helmholtz equation. This equation is used to model wave
propagation in the frequency domain. However, in some cases the acoustic
equation is not sufficient for modeling the physics of the wave propagation,
and one has to consider the elastic Helmholtz equation. Such a case arises in
geophysical seismic imaging applications, where the earth's subsurface is the
elastic medium. The elastic Helmholtz equation is much harder to solve than its
acoustic counterpart, partially because it is three times larger, and partially
because it models more complicated physics. Despite this, there are very few
solvers available for the elastic equation compared to the array of solvers
that are available for the acoustic one. In this work we extend the shifted
Laplacian approach to the elastic Helmholtz equation, by combining the complex
shift idea with approaches for linear elasticity. We demonstrate the efficiency
and properties of our solver using numerical experiments for problems with
heterogeneous media in two and three dimensions
An Elimination Method for Solving Bivariate Polynomial Systems: Eliminating the Usual Drawbacks
We present an exact and complete algorithm to isolate the real solutions of a
zero-dimensional bivariate polynomial system. The proposed algorithm
constitutes an elimination method which improves upon existing approaches in a
number of points. First, the amount of purely symbolic operations is
significantly reduced, that is, only resultant computation and square-free
factorization is still needed. Second, our algorithm neither assumes generic
position of the input system nor demands for any change of the coordinate
system. The latter is due to a novel inclusion predicate to certify that a
certain region is isolating for a solution. Our implementation exploits
graphics hardware to expedite the resultant computation. Furthermore, we
integrate a number of filtering techniques to improve the overall performance.
Efficiency of the proposed method is proven by a comparison of our
implementation with two state-of-the-art implementations, that is, LPG and
Maple's isolate. For a series of challenging benchmark instances, experiments
show that our implementation outperforms both contestants.Comment: 16 pages with appendix, 1 figure, submitted to ALENEX 201
Implicitization of curves and (hyper)surfaces using predicted support
We reduce implicitization of rational planar parametric curves and (hyper)surfaces to linear algebra, by interpolating the coefficients of the implicit equation.
For predicting the implicit support, we focus on methods that exploit input and output structure in the sense of sparse (or toric) elimination theory, namely by computing the Newton polytope of the implicit polynomial, via sparse resultant theory.
Our algorithm works even in the presence of base points but, in this case, the implicit equation shall be obtained as a factor of the produced polynomial.
We implement our methods on Maple, and some on Matlab as well, and study their numerical stability and efficiency on several classes of curves and surfaces.
We apply our approach to approximate implicitization,
and quantify the accuracy of the approximate output,
which turns out to be satisfactory on all tested examples; we also relate our measures to Hausdorff distance.
In building a square or rectangular matrix, an important issue is (over)sampling the given curve or surface: we conclude that unitary complexes offer the best tradeoff between speed and accuracy when numerical methods are employed, namely SVD, whereas for exact kernel computation random integers is the method of choice.
We compare our prototype to existing software and find that it is rather competitive
Compositional Performance Modelling with the TIPPtool
Stochastic process algebras have been proposed as compositional specification formalisms for performance models. In this paper, we describe a tool which aims at realising all beneficial aspects of compositional performance modelling, the TIPPtool. It incorporates methods for compositional specification as well as solution, based on state-of-the-art techniques, and wrapped in a user-friendly graphical front end. Apart from highlighting the general benefits of the tool, we also discuss some lessons learned during development and application of the TIPPtool. A non-trivial model of a real life communication system serves as a case study to illustrate benefits and limitations
Track 3: Computations in theoretical physics -- techniques and methods
Here, we attempt to summarize the activities of Track 3 of the 17th
International Workshop on Advanced Computing and Analysis Techniques in Physics
Research (ACAT 2016).Comment: 10 pages, 3 figures, to appear in the proceedings of ACAT 201
- âŠ