Related papers: On Broyden's method for the solution of the multil…
On an example of the open nonlinear electrodynamic system - transverse non-homogeneous, isotropic, nonlinear (a Kerr-like dielectric nonlinearity) dielectric layer, the algorithms of solution of the diffraction problem of a plane wave on…
The Transfer Matrix Method (TMM) is a widely used technique for modeling linear propagation of electromagnetic waves through stratified layered media. However, since its extension to inhomogeneous and nonlinear systems is not…
The branching methods developed are effective methods to solve some semi linear PDEs and are shown numerically to be able to solve some full non linear PDEs. These methods are however restricted to some small coefficients in the PDE and…
Several possibilities exist to implement the propagation step of the lattice Boltzmann method. This paper describes common implementations which are compared according to the number of memory transfer operations they require per lattice…
Modern AI relies on huge matrix multiplications (MatMuls), whose computation poses a scalability problem for inference and training. We propose an alternative, GPU native bilinear operator to MatMuls in neural networks, which offers a…
Nowadays computational complexity of fast walsh hadamard transform and nonlinearity for Boolean functions and large substitution boxes is a major challenge of modern cryptography research on strengthening encryption schemes against linear…
An efficient linear solver plays an important role while solving partial differential equations (PDEs) and partial integro-differential equations (PIDEs) type mathematical models. In most cases, the efficiency depends on the stability and…
Solving stiff ordinary differential equations (StODEs) requires sophisticated numerical solvers, which are often computationally expensive. In general, traditional explicit time integration schemes with restricted time step sizes are not…
This paper studies system identification for nonlinear state-space models, a problem that arises across many fields yet remains challenging in practice. Focusing on maximum likelihood estimation, we employ Bayesian optimization (BayesOpt)…
Bayesian model comparison (BMC) offers a principled approach for assessing the relative merits of competing computational models and propagating uncertainty into model selection decisions. However, BMC is often intractable for the popular…
When integrating the radiative transfer equation for polarized light, the necessity of high-order numerical methods is well known. In fact, well-performing high-order formal solvers enable higher accuracy and the use of coarser spatial…
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…
We introduce a new methodology for a fast and reliable discrimination between ordered and chaotic orbits in multidimensional Hamiltonian systems which we call the Linear Dependence Index (LDI). The new method is based on the recently…
The fast multipole method (FMM) performs fast approximate kernel summation to a specified tolerance $\epsilon$ by using a hierarchical division of the domain, which groups source and receiver points into regions that satisfy local…
In this paper we explore the possibility for solving the 3D Maxwell's equations in the presence of nonlinear and/or inhomogeneous material response. We propose using a hybrid approach which combines a bound- ary integral representation with…
The purpose of the current work is to find numerical solutions of the steady state inhomogeneous Vlasov equation. This problem has a wide range of applications in the kinetic simulation of non-thermal plasmas. However, the direct…
We use variational methods to study the existence of nontrivial and radially symmetric solutions to the H\`enon-Lane-Emden system with weights, when the exponents involved lie on the "critical hyperbola". We also discuss qualitative…
We derive and analyze a symmetric interior penalty discontinuous Galerkin scheme for the approximation of the second-order form of the radiative transfer equation in slab geometry. Using appropriate trace lemmas, the analysis can be carried…
Recently, researchers have utilized neural networks to accurately solve partial differential equations (PDEs), enabling the mesh-free method for scientific computation. Unfortunately, the network performance drops when encountering a high…
We formulate the problem of a two-level system in a linearly polarized laser field in terms of a nonlinear Riccati-type differential equation and solve the equation analytically in time intervals much shorter than half the optical period.…