Related papers: Increased Performance of Matsubara space calculati…
We present an algorithm which uses Fujiwara's inequality to bound algebraic functions over ellipses of a certain type, allowing us to concretely implement a rigorous Gauss-Legendre integration method for algebraic functions over a line…
In many astrophysical plasmas, the Coulomb collision is insufficient to maintain an isotropic temperature, and the system is driven to the anisotropic regime. In this case, magnetohydrodynamic (MHD) models with anisotropic pressure are…
Simple finite differencing of the anisotropic diffusion equation, where diffusion is only along a given direction, does not ensure that the numerically calculated heat fluxes are in the correct direction. This can lead to negative…
When used to accelerate the convergence of fixed-point iterative methods, such as the Picard method, which is a kind of nonlinear fixed-point iteration, polynomial extrapolation techniques can be very effective. The numerical solution of…
We obtain a comparison result for solutions to nonlinear fully anisotropic elliptic problems by means of anisotropic symmetrization. As consequence we deduce a priori estimates for norms of the relevant solutions.
In this paper, we present a high order finite difference solver for anisotropic diffusion problems based on the first-order hyperbolic system method. In particular, we demonstrate that the construction of a uniformly accurate fifth-order…
In this paper we propose fast solution methods for the Cauchy problem for the multidimensional Schr\"odinger equation. Our approach is based on the approximation of the data by the basis functions introduced in the theory of approximate…
A simple method for numerical analytic continuation is developed. It is designed to analytically continue the imaginary time (Matsubara frequency) quantum Monte Carlo simulation results to the real time (real frequency) domain. Such a…
We give a unified proof for the well-posedness of a class of linear half-space equations with general incoming data and construct a Galerkin method to numerically resolve this type of equations in a systematic way. Our main strategy in both…
There are three main types of numerical computations for the Bessel function of the second kind: series expansion, continued fraction, and asymptotic expansion. In addition, they are combined in the appropriate domain for each. However,…
Finite difference based micromagnetic simulations are a powerful tool for the computational investigation of magnetic structures. In this paper, we demonstrate how the discretization of continuous micromagnetic equations introduces a…
In this work, an efficient approximation scheme has been proposed for getting accurate approximate solution of nonlinear partial differential equations with constant or variable coefficients satisfying initial conditions in a series of…
A general synthetic iterative scheme is proposed to solve the Enskog equation within a Monte Carlo framework. The method demonstrates rapid convergence by reducing intermediate Monte Carlo evolution and preserves the asymptotic-preserving…
This is an expository version of our paper [arXiv:1902.07384]. Our aim is to present recent Macaulay2 algorithms for computation of mixed multiplicities of ideals in a Noetherian ring which is either local or a standard graded algebra over…
Efficient ab initio calculations of correlated materials at finite temperature require compact representations of the Green's functions both in imaginary time and Matsubara frequency. In this paper, we introduce a general procedure which…
High-temperature series are computed for a generalized $3d$ Ising model with arbitrary potential. Two specific ``improved'' potentials (suppressing leading scaling corrections) are selected by Monte Carlo computation. Critical exponents are…
The convergence analysis of a third-order scheme for the highly nonlinear Landau-Lifshitz-Gilbert equation with a non-convex constraint is considered. In this paper, we first present a fully discrete semi-implicit method for solving the…
The quasi-isotropic inhomogeneous solution of the Einstein equations near a cosmological singularity in the form of a series expansion in the synchronous system of reference, first found by Lifshitz and Khalatnikov in 1960, is generalized…
A second order accurate numerical scheme is proposed and implemented for the Landau-Lifshitz-Gilbert equation, which models magnetization dynamics in ferromagnetic materials, with large damping parameters. The main advantages of this method…
There are many numerical methods for solving partial different equations (PDEs) on manifolds such as classical implicit, finite difference, finite element, and isogeometric analysis methods which aim at improving the interoperability…