Related papers: Preconditioned Bi-Conjugate Gradient Method for Ra…
Space-time modulation adds another powerful degree of freedom to the manipulation of classical wave systems. It opens the door for complex control of wave behavior beyond the reach of stationary systems, such as nonreciprocal wave transport…
Rectangular spectral collocation (RSC) methods have recently been proposed to solve linear and nonlinear differential equations with general boundary conditions and/or other constraints. The involved linear systems in RSC become extremely…
We present an extension of TRAPHIC, the method for radiative transfer of ionising radiation in smoothed particle hydrodynamics simulations that we introduced in Pawlik & Schaye (2008). The new version keeps all advantages of the original…
The conjugate gradient (CG) method is a classic Krylov subspace method for solving symmetric positive definite linear systems. We introduce an analogous semi-conjugate gradient (SCG) method for unsymmetric positive definite linear systems.…
There exist two major problems in application of the conventional block BiCGSTAB method to the O(a)-improved Wilson-Dirac equation with multiple right-hand-sides: One is the deviation between the true and the recursive residuals. The other…
Birefringent crystals are extensively used to manipulate polarized light. The generalized transfer matrix developed allows efficient calculation of the full polarization state of light transmitted through and reflected by a stack of…
Diffusion models have shown remarkable promise for image restoration by leveraging powerful priors. Prominent methods typically frame the restoration problem within a Bayesian inference framework, which iteratively combines a denoising step…
The use of the Preconditioned Conjugate Gradient (PCG) method for computing the Generalized Least Squares (GLS) estimator of the General Linear Model (GLM) is considered. The GLS estimator is expressed in terms of the solution of an…
In this paper, we propose two mixed precision algorithms for Block-Jacobi preconditioner(BJAC): a fixed low precision strategy and an adaptive precision strategy. We evaluate the performance improvement of the proposed mixed precision BJAC…
The conjugate gradient method is a widely used algorithm for the numerical solution of a system of linear equations. It is particularly attractive because it allows one to take advantage of sparse matrices and produces (in case of infinite…
This paper is devoted to deal with some mathematical and numerical aspects of the radiative integral transfer equations. First, the properties of the raidative integral operators are analyzed. Based on these results, the existence and…
A new method, 'Configurationally-Resolved-Super-Transition-Arrays', for calculation of the spectral absorption coefficient in hot plasmas is presented. In the new method, the spectrum of each Super-Transition-Array is evaluated as the…
The linearly-polarized solar limb spectrum that is produced by scattering processes contains a wealth of information on the physical conditions and magnetic fields of the solar outer atmosphere, but the modeling of many of its strongest…
Reconstruction of Attosecond Beating By Interference of Two-photon Transitions (RABBITT) is a technique that can be used to determine the phases of atomic transition elements in photoionization processes. In the traditional RABBITT scheme,…
In several natural and engineering systems, changes in control parameters can trigger bifurcations that lead to sustained or growing periodic oscillations, indicating the onset of oscillatory instabilities. Such emergent behaviour often…
We present a variant of the s-step Preconditioned Conjugate Gradient (PCG) method that combines a Chebyshev-stabilized Krylov basis with a Forward Gauss-Seidel (FGS) iteration for the solution of the reduced Gram systems. In s-step…
Kinematic equations for the motion of slowly propagating, weakly curved fronts in bistable media are derived. The equations generalize earlier derivations where algebraic relations between the normal front velocity and its curvature are…
We present a general formalism for computing self-consistent, numerical solutions to the time-dependent radiative transfer equation in low velocity, multi-level ions undergoing radiative interactions. Recent studies of time-dependent…
The formulation of the on-surface radiation condition (OSRC) is extended to handle wave scattering problems in the presence of multiple obstacles. The new multiple-OSRC simultaneously accounts for the outgoing behavior of the wave fields,…
We derive new equations using the mixed-frame approach for one- and two-dimensional (axisymmetric) time-dependent radiation transport and the associated couplings with matter. Our formulation is multi-group and multi-angle and includes…