Related papers: A parallel code for multiprecision computations of…
We present a novel approach to the parallelization of the parabolic fast multipole method for a space-time boundary element method for the heat equation. We exploit the special temporal structure of the involved operators to provide an…
While not generally a conservation law, any symmetry of the equations of motion implies a useful reduction of any second-order equationto a first-order equation between invariants, whose solutions (first integrals) can then be integrated by…
We present a new parallel numerical method for solving the non-stationary Schr\"odinger equation with linear nonlocal condition and time-dependent potential which does not commute with the stationary part of the Hamiltonian. The given…
We consider a linear elliptic partial differential equation (PDE) with a generic uniformly bounded parametric coefficient. The solution to this PDE problem is approximated in the framework of stochastic Galerkin finite element methods. We…
The so-called 'direct' approach to separation of variables in linear PDEs is applied to the hydrodynamic stability problem. Calculations are made for the complete linear stability equations in cylindrical coordinates. Several classes of the…
This paper presents a highly-parallelizable parallel-in-time algorithm for efficient solution of nonlinear time-periodic problems. It is based on the time-periodic extension of the Parareal method, known to accelerate sequential…
In this work numerical methods for solving Einstein's equations are developed and applied to the study of inhomogeneous cosmological models. A two-dimensional computer code is described which implements two advanced numerical methods:…
We develop a "hybrid approximative scheme" in the framework of the post-Newtonian approximation for computing general-relativistic polytropic models simulating neutron stars in critical rigid rotation. We treat the differential equations…
It is quite common that a nonlinear partial differential equation (PDE) admits multiple distinct solutions and each solution may carry a unique physical meaning. One typical approach for finding multiple solutions is to use the Newton…
The approximation of invariant measures for nonlinear ergodic stochastic differential equations (SDEs) is a central problem in scientific computing, with important applications in stochastic sampling, physics, and ecology. We first propose…
Linear differential equations of arbitrary order with polynomial coefficients are considered. Specifically, necessary and sufficient conditions for the existence of polynomial solutions of a given degree are obtained for these equations. An…
The existence and multiplicity of solutions to a quasilinear, elliptic partial differential equation (PDE) with singular non-linearity is analyzed. The PDE is a recently derived variant of a canonical model used in the modeling of…
We propose in this paper the Wavelet-based Edge Multiscale Parareal (WEMP) Algorithm to efficiently solve parabolic equations with heterogeneous coefficients. This algorithm combines the advantages of multiscale methods that can deal with…
We derive computable error estimates for finite element approximations of linear elliptic partial differential equations (PDE) with rough stochastic coefficients. In this setting, the exact solutions contain high frequency content that…
The purpose of the research is to find the numerical solutions to the system of time dependent nonlinear parabolic partial differential equations (PDEs) utilizing the Modified Galerkin Weighted Residual Method (MGWRM) with the help of…
Deviations from kinetic equilibrium of massive particles caused by the universe expansion are calculated analytically in the Boltzmann approximation. For the case of an energy independent amplitude of elastic scattering, an exact partial…
Developing algorithms for solving high-dimensional partial differential equations (PDEs) has been an exceedingly difficult task for a long time, due to the notoriously difficult problem known as the "curse of dimensionality". This paper…
A numerical method to solve linear integro-differential equations is presented. This method has been used to solve the QCD Altarelli-Parisi evolution equations within the H1 Collaboration at DESY-Hamburg. Mathematical aspects and numerical…
Accurate numerical simulations of interaction between fluid and solid play an important role in applications. The task is challenging in practical scenarios as the media are usually highly heterogeneous with very large contrast. To overcome…
The TREE method has been widely used for long-range interaction {\it N}-body problems. We have developed a parallel TREE code for two-component classical plasmas with open boundary conditions and highly non-uniform charge distributions. The…