Related papers: Increased Performance of Matsubara space calculati…
Surface integral equation (SIE) methods are of great interest for the numerical solution of Maxwell's equations in the presence of homogeneous objects. However, existing SIE algorithms have limitations, either in terms of scalability,…
An efficient method is described to handle mesh indexes in multidimensional problems like numerical integration of partial differential equations, lattice model simulations, and determination of atomic neighbor lists. By creating an…
This paper deals with the numerical study of a strongly anisotropic heat equation. The use of standard schemes in this situation leads to poor results, due to the high anisotropy. Furthermore, the recently proposed Asymptotic-Preserving…
In some cases, computational benefit can be gained by exploring the hyper parameter space using a deterministic set of grid points instead of a Markov chain. We view this as a numerical integration problem and make three unique…
In this paper we present a treatment of hypersingular integral equations, which have relevant applications in many problems of wave dynamics, elasticity and fluid mechanics with mixed boundary conditions. The main goal of the present work…
We present an extension to high-order of a first-order Lagrange-projection like method for the approximation of the Euler equations introduced in Coquel {\it et al.} (Math. Comput., 79 (2010), pp.~1493--1533). The method is based on a…
A new method is proposed to improve the numeri- cal simulation of time dependent problems when the initial and boundary data are not compatible. Unlike earlier methods limited to space dimension one, this method can be used for any space…
The aim of this work is to obtain new analitical solutions for Einstein equations in the anisotropical domain. This will be done via the minimal geometric deformation (MGD) approach, which is a simple and systematical method that allow us…
We study the Einstein-Maxwell system of equations in spherically symmetric gravitational fields for static interior spacetimes. The condition for pressure isotropy is reduced to a recurrence equation with variable, rational coefficients. We…
In order to perform accurate and stable long-term numerical integration of the Einstein equations, several hyperbolic systems have been proposed. We here report our numerical comparisons between weakly hyperbolic, strongly hyperbolic, and…
We introduce in this document a direct method allowing to solve numerically inverse type problems for linear parabolic equations. We consider the reconstruction of the full solution of the parabolic equation posed in $\Omega\times (0,T)$ -…
The numerical simulation of the 3D incompressible Euler equation is analyzed with respect to different integration methods. The numerical schemes we considered include spectral methods with different strategies for dealiasing and two…
We show how symmetry properties can be used to greatly increase the accuracy and efficiency in auxiliary-field quantum Monte Carlo (AFQMC) calculations of electronic systems. With the Hubbard model as an example, we study symmetry…
To improve the computational efficiencies of the real-space orbital-free density functional theory, this work develops a new single-grid solver by directly providing the closed-form solution to the inner iteration and using an improved…
We present an auxiliary space theory that provides a unified framework for analyzing various iterative methods for solving linear systems that may be semidefinite. By interpreting a given iterative method for the original system as an…
The analytic continuation of imaginary-time quantum Monte Carlo data to extract real-frequency spectra remains a key problem in connecting theory with experiment. Here we present a fast and efficient stochastic optimization method (FESOM)…
An approach is treated for numerical integration of ordinary differential equations systems of the first order with choice of a computation scheme, ensuring the required local precision. The treatment is made on the basis of schemes of…
We present sharp estimates for the extremal eigenvalues of the Schur complements arising in saddle point problems. These estimates are derived using the auxiliary space theory, in which a given iterative method is interpreted as an…
When solving partial differential equations using classical schemes such as finite difference or finite volume methods, sufficiently fine meshes and carefully designed schemes are required to achieve high-order accuracy of numerical…
In this paper the finite-difference time-domain general vector auxiliary differential equation method [Greene, J. H. and A. Taflove, Opt. Express 14, 8305 (2006)], nonlinear polarization vector, the nonlinear electric dipole moment per unit…