Related papers: On Broyden's method for the solution of the multil…
In this paper, we propose novel, fully Bayesian non-parametric tests for one-sample and two-sample multivariate location problems. We model the underlying distribution using a Dirichlet process prior, and develop a testing procedure based…
The inversion of spectropolarimetric observations of the solar upper atmosphere is one of the most challenging goals in solar physics. If we account for all relevant ingredients of the spectral line formation process, such as the…
A finite element method for solving the resonance line transfer problem in moving media is presented. The algorithm works in three spatial dimensions on unstructured grids which are adaptively refined by means of an a posteriori error…
We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core…
We derive a nonlinear moment model for radiative transfer equation in 3D space, using the method to derive the nonlinear moment model for the radiative transfer equation in slab geometry. The resulted 3D HMPN model enjoys a list of…
In a transformation method the numerical solution of a given boundary value problem is obtained by solving one or more related initial value problems. This paper is concerned with the application of the iterative transformation method to…
We present a parallel version of the well-known Split-Step Fourier method (SSF) for solving the Nonlinear Schr\"odinger equation, a mathematical model describing wave packet propagation in fiber optic lines. The algorithm is implemented…
Here we present a new non-parametric approach to density estimation and classification derived from theory in Radon transforms and image reconstruction. We start by constructing a "forward problem" in which the unknown density is mapped to…
Stochastic differential equation (SDE in short) solvers find numerous applications across various fields. However, in practical simulations, we usually resort to using Ito-Taylor series-based methods like the Euler-Maruyama method. These…
A numerical scheme is proposed for the solution of the three-dimensional radiative transfer equation with variable optical depth. We show that time-dependent ray tracing is an attractive choice for simulations of astrophysical ionization…
The mutual information (MI) of Gaussian multi-input multi-output (MIMO) channels has been evaluated by utilizing random matrix theory (RMT) and shown to asymptotically follow Gaussian distribution, where the ergodic mutual information (EMI)…
This paper presents an overview of the radiative transfer problem of calculating the spectral line intensity and polarization that emerges from a (generally magnetized) astrophysical plasma composed of atoms and molecules whose excitation…
This paper presents a fast wavefield evaluation method for two-dimensional wave scattering problems. The proposed method is based on a modified version of proxy-surface-accelerated interpolative decomposition, making it effective even if…
3-D astrophysical atmospheres will have random velocity fields. We seek to combine the methods we have developed for solving the 1-D problem with arbitrary flows to those that we have developed for solving the fully 3-D relativistic…
A meshless method is presented to solve the radiative transfer equation in the even parity formulation of the discrete ordinates method in complex 2D and 3D geometries. Prediction results of radiative heat transfer problems obtained by the…
X-ray scatter has been a serious concern in computed tomography (CT), leading to image artifacts and distortion of CT values. The linear Boltzmann transport equation (LBTE) is recognized as a fast and accurate approach for scatter…
In this paper we examine iterative methods for solving the forward ($A{\bf x}={\bf b}$) and adjoint ($A^{T}{\bf y}={\bf g}$) systems of linear equations used to approximate the scattering amplitude, defined by ${\bf g}^{T}{\bf x}={\bf…
In this work we apply the Adomian decomposition method combined with the Laplace transform (LADM) in order to solve the 1-dimensional nonlinear Schrodinger equation whose nonlinear term presents a nonlinear defocusing strength that varies…
In this paper, we consider numerical simulations of the nonlocal optical response of metallic nanostructure arrays inside a dielectric host, which is of particular interest to the nanoplasmonics community due to many unusual properties and…
Numerical transfer matrices have been widely used in the study of wave propagation and scattering. These may be viewed as descretizations of a recently introduced fundamental notion of transfer matrix which admits a representation in terms…