Related papers: H-Matrix Accelerated Direct Matrix Solver for Maxw…
Large, sparse linear systems are pervasive in modern science and engineering, and Krylov subspace solvers are an established means of solving them. Yet convergence can be slow for ill-conditioned matrices, so practical deployments usually…
Computation of the trace of a matrix function plays an important role in many scientific computing applications, including applications in machine learning, computational physics (e.g., lattice quantum chromodynamics), network analysis and…
In this paper, we propose an interior-point method for linearly constrained optimization problems (possibly nonconvex). The method - which we call the Hessian barrier algorithm (HBA) - combines a forward Euler discretization of Hessian…
We present BiqBin, an exact solver for linearly constrained binary quadratic problems. Our approach is based on an exact penalty method to first efficiently transform the original problem into an instance of Max-Cut, and then to solve the…
This paper introduces a novel boundary integral equation (BIE) method for the numerical solution of problems of planewave scattering by periodic line arrays of two-dimensional penetrable obstacles. Our approach is built upon a direct BIE…
A novel tensor-based formula for solving the linear systems involving Kronecker sum is proposed. Such systems are directly related to the matrix and tensor forms of Sylvester equation. The new tensor-based formula demonstrates the…
This text proposes a fast, rapidly convergent Nystr\"{o}m method for the solution of the Lippmann-Schwinger integral equation that mathematically models the scattering of time-harmonic acoustic waves by inhomogeneous obstacles, while…
Considered in this paper the Maxwell-Bloch (MB) equations became known after Lamb [1-4]. In [5] Ablowitz, Kaup and Newell proposed the inverse scattering transform (IST) to the Maxwell-Bloch equations for studying a physical phenomenon…
This manuscript details the use of the rational Chebyshev transform for describing the transverse dynamics of high-power laser diodes, either broad area lasers, index guided lasers or monolithic master oscillator power amplifier devices.…
An efficient coarse-mesh nodal integral method (NIM), based on cell-centered variables and termed the cell-centered NIM (CCNIM), is developed and applied to solve multi-dimensional, time-dependent, nonlinear Burgers equations, extending the…
Emerging tensor network techniques for solutions of Partial Differential Equations (PDEs), known for their ability to break the curse of dimensionality, deliver new mathematical methods for ultrafast numerical solutions of high-dimensional…
The Half-Space Matching (HSM) method has recently been developed as a new method for the solution of 2D scattering problems with complex backgrounds, providing an alternative to Perfectly Matched Layers (PML) or other artificial boundary…
We apply boundary integral equations for the first time to the two-dimensional scattering of time-harmonic waves from a smooth obstacle embedded in a continuously-graded unbounded medium. In the case we solve the square of the wavenumber…
The paper is concerned with the three-dimensional electromagnetic scattering from a large open rectangular cavity that is embedded in a perfectly electrically conducting infinite ground plane. By introducing a transparent boundary…
A preconditioned, multipole-accelerated, Krylov-subspace iterative algorithm for the electromagnetic scattering analysis of three dimensional (3D), arbitrary shaped dielectric structures composed of single and multi-layered dielectric…
A semi-spectral Chebyshev method for solving numerically singular integral equations is presented and applied in the quarkonium bound-state problem in momentum space. The integrals containing both, logarithmic and Cauchy singular kernels,…
We present a fast direct solver for structured linear systems based on multilevel matrix compression. Using the recently developed interpolative decomposition of a low-rank matrix in a recursive manner, we embed an approximation of the…
We investigate high-order Convolution Quadratures methods for the solution of the wave equation in unbounded domains in two dimensions that rely on Nystr\"om discretizations for the solution of the ensemble of associated Laplace domain…
In this paper, based on an optimization problem, a sketch-and-project method for solving the linear matrix equation AXB = C is proposed. We provide a thorough convergence analysis for the new method and derive a lower bound on the…
We consider the problem of efficiently solving large-scale linear least squares problems that have one or more linear constraints that must be satisfied exactly. Whilst some classical approaches are theoretically well founded, they can face…