Related papers: Preconditioned Bi-Conjugate Gradient Method for Ra…
In this paper, we present a structured solver based on the preconditioned conjugate gradient method (PCGM) for solving the linear quadratic (LQ) optimal control problem for $K \times N$ sub-systems connected in a two-dimensional (2D) grid…
Convolution-type integral equations arise from various fields, \textit{e.g.}, finite impulse response filters in signal processing and deblurring problems in image processing. When solving these equations, conventional numerical methods,…
Observations and magnetohydrodynamic simulations of solar and stellar atmospheres reveal an intermittent behavior or steep gradients in physical parameters, such as magnetic field, temperature, and bulk velocities. The numerical solution of…
Properties of Superiorized Preconditioned Conjugate Gradient (SupPCG) algorithms in image reconstruction from projections are examined. Least squares (LS) is usually chosen for measuring data-inconsistency in these inverse problems.…
A novel method for solving the linear radiative transport equation (RTE) in a three-dimensional homogeneous medium is proposed and illustrated with numerical examples. The method can be used with an arbitrary phase function A(s,s') with the…
We study the solution of large symmetric positive-definite linear systems in a matrix-free setting with a limited iteration budget. We focus on the preconditioned conjugate gradient (PCG) method with spectral preconditioning. Spectral…
This paper presents a spatial two-grid (STG) compact difference scheme for a two-dimensional (2D) nonlinear diffusion-wave equation with variable exponent, which describes, e.g., the propagation of mechanical diffusive waves in viscoelastic…
Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide essential insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area…
A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods - the polytropic cooling approximation as devised by Stamatellos et al (2007), and…
Statistical preconditioning enables fast methods for distributed large-scale empirical risk minimization problems. In this approach, multiple worker nodes compute gradients in parallel, which are then used by the central node to update the…
A structured preconditioned conjugate gradient (PCG) solver is developed for the Newton steps in second-order methods for a class of constrained network optimal control problems. Of specific interest are problems with discrete-time dynamics…
A method for the fast and accurate solution of the radiative transfer equation in plane-parallel media with coherent isotropic scattering is presented. This largely analytical approach uses the formalism of meromorphic functions in order to…
Accurate photometric and kinematic modelling of disc galaxies requires the inclusion of radiative transfer models. Due to the complexity of the radiative transfer equation (RTE), sophisticated techniques are required. Various techniques…
An implicit method for radiative transfer in SPH is described. The diffusion approximation is used, and the hydrodynamic calculations are performed by a fully three--dimensional SPH code. Instead of the energy equation of state for an ideal…
The atmospheres of planets (including Earth) and the outer layers of stars have often been treated in radiative transfer as plane-parallel media, instead of spherical shells, which can lead to inaccuracy, e.g. limb darkening. We give an…
Efficient numerical solvers for partial differential equations empower science and engineering. One of the commonly employed numerical solvers is the preconditioned conjugate gradient (PCG) algorithm which can solve large systems to a given…
We compare four different methods to calculate radiative transfer through a magnetized stellar atmosphere, and apply them to the case of magnetic white dwarfs. All methods are numerically stable enough to allow determination of the magnetic…
A new polynomial preconditioner for symmetric complex linear systems based on Hermitian and skew-Hermitian splitting (HSS) for complex symmetric linear systems is herein presented. It applies to Conjugate Orthogonal Conjugate Gradient…
We develop a new numerical scheme for solving the radiative transfer equation in a spherically symmetric system. This scheme does not rely on any kind of diffusion approximation and it is accurate for optically thin, thick, and intermediate…
This paper introduces a method to preform optical tomography, using 3D radiative transfer as the forward model. We use an iterative approach predicated on the Spherical Harmonics Discrete Ordinates Method (SHDOM) to solve the optimization…