Search results for " Numerical analysis"
showing 10 items of 106 documents
Iterative Reconstruction of Signals on Graph
2020
We propose an iterative algorithm to interpolate graph signals from only a partial set of samples. Our method is derived from the well known Papoulis-Gerchberg algorithm by considering the optimal value of a constant involved in the iteration step. Compared with existing graph signal reconstruction algorithms, the proposed method achieves similar or better performance both in terms of convergence rate and computational efficiency.
Cell-average WENO with progressive order of accuracy close to discontinuities with applications to signal processing
2020
In this paper we translate to the cell-average setting the algorithm for the point-value discretization presented in S. Amat, J. Ruiz, C.-W. Shu, D. F. Y\'a\~nez, A new WENO-2r algorithm with progressive order of accuracy close to discontinuities, submitted to SIAM J. Numer. Anal.. This new strategy tries to improve the results of WENO-($2r-1$) algorithm close to the singularities, resulting in an optimal order of accuracy at these zones. The main idea is to modify the optimal weights so that they have a nonlinear expression that depends on the position of the discontinuities. In this paper we study the application of the new algorithm to signal processing using Harten's multiresolution. Se…
Multidomain spectral method for the Gauss hypergeometric function
2018
International audience; We present a multidomain spectral approach for Fuchsian ordinary differential equations in the particular case of the hypergeometric equation. Our hybrid approach uses Frobenius’ method and Moebius transformations in the vicinity of each of the singular points of the hypergeometric equation, which leads to a natural decomposition of the real axis into domains. In each domain, solutions to the hypergeometric equation are constructed via the well-conditioned ultraspherical spectral method. The solutions are matched at the domain boundaries to lead to a solution which is analytic on the whole compactified real line R∪∞, except for the singular points and cuts of the Rie…
Numerical study of soliton stability, resolution and interactions in the 3D Zakharov–Kuznetsov equation
2021
International audience; We present a detailed numerical study of solutions to the Zakharov-Kuznetsov equation in three spatial dimensions. The equation is a three-dimensional generalization of the Korteweg-de Vries equation, though, not completely integrable. This equation is L-2-subcritical, and thus, solutions exist globally, for example, in the H-1 energy space.We first study stability of solitons with various perturbations in sizes and symmetry, and show asymptotic stability and formation of radiation, confirming the asymptotic stability result in Farah et al. (0000) for a larger class of initial data. We then investigate the solution behavior for different localizations and rates of de…
On 1-Laplacian Elliptic Equations Modeling Magnetic Resonance Image Rician Denoising
2015
Modeling magnitude Magnetic Resonance Images (MRI) rician denoising in a Bayesian or generalized Tikhonov framework using Total Variation (TV) leads naturally to the consideration of nonlinear elliptic equations. These involve the so called $1$-Laplacian operator and special care is needed to properly formulate the problem. The rician statistics of the data are introduced through a singular equation with a reaction term defined in terms of modified first order Bessel functions. An existence theory is provided here together with other qualitative properties of the solutions. Remarkably, each positive global minimum of the associated functional is one of such solutions. Moreover, we directly …
WENO schemes applied to the quasi-relativistic Vlasov-Maxwell model for laser-plasma interaction
2014
Abstract In this paper we focus on WENO-based methods for the simulation of the 1D Quasi-Relativistic Vlasov–Maxwell (QRVM) model used to describe how a laser wave interacts with and heats a plasma by penetrating into it. We propose several non-oscillatory methods based on either Runge–Kutta (explicit) or Time-Splitting (implicit) time discretizations. We then show preliminary numerical experiments.
Annihilation Operators for Exponential Spaces in Subdivision
2022
We investigate properties of differential and difference operators annihilating certain finite-dimensional subspaces of exponential functions in two variables that are connected to the representation of real-valued trigonometric and hyperbolic functions. Although exponential functions appear in a variety of contexts, the motivation behind this work comes from considering subdivision schemes with the capability of preserving those exponential functions required for an exact description of surfaces parametrized in terms of trigonometric and hyperbolic functions.
On the analysis of Catalan thin vaults
2010
This paper is concerned with the identification, modelling and analysis of thin layered vaults, typical in the Catalan constructions of the XIX century. These special structures, also known as bovedas tabicadas, are characterized by very low thickness with respect to the medium surface dimensions and by the presence of different superimposed layers of bricks tied with mortar, are studied in order to individuate a coherent mechanical model for describing the material behaviour, to recognize the structural response utilizing as comparison adequate experimental results, and to extend the obtained results for the analysis of new vaults or for the restoration design of existing vaults. Firstly, …
Numerical study of blow-up in solutions to generalized Kadomtsev-Petviashvili equations
2013
We present a numerical study of solutions to the generalized Kadomtsev-Petviashvili equations with critical and supercritical nonlinearity for localized initial data with a single minimum and single maximum. In the cases with blow-up, we use a dynamic rescaling to identify the type of the singularity. We present a discussion of the observed blow-up scenarios.
Towards Stable Radial Basis Function Methods for Linear Advection Problems
2021
In this work, we investigate (energy) stability of global radial basis function (RBF) methods for linear advection problems. Classically, boundary conditions (BC) are enforced strongly in RBF methods. By now it is well-known that this can lead to stability problems, however. Here, we follow a different path and propose two novel RBF approaches which are based on a weak enforcement of BCs. By using the concept of flux reconstruction and simultaneous approximation terms (SATs), respectively, we are able to prove that both new RBF schemes are strongly (energy) stable. Numerical results in one and two spatial dimensions for both scalar equations and systems are presented, supporting our theoret…