Related papers: Iterative solution and preconditioning for the tan…
In this paper, fast numerical methods are established for solving a class of time distributed-order and Riesz space fractional diffusion-wave equations. We derive new difference schemes by the weighted and shifted Gr$\ddot{\rm{u}}$nwald…
Gradient schemes is a framework that enables the unified convergence analysis of many numerical methods for elliptic and parabolic partial differential equations: conforming and non-conforming Finite Element, Mixed Finite Element and Finite…
The current-induced magnetisation dynamics in a ferromagnet at elevated temperatures can be described by the Landau--Lifshitz--Bloch (LLB) equation with spin-torque terms. In this paper, we focus on the regime above the Curie temperature.…
Here we describe a development of computer algorithm to simulate the Time Dependent Ginzburg-Landau equation (TDGL) and its application to understand superconducting vortex dynamics in confined geometries. Our initial motivation to get…
Finite element methods for symmetric linear hyperbolic systems using unstructured advancing fronts (satisfying a causality condition) are considered in this work. Convergence results and error bounds are obtained for mapped tent pitching…
The importance of simulating pinning arrays in superconducting samples for the increase of critical currents has been increasing in the last few years. Since the Time Dependent Ginzburg Landau (TDGL) can be more accurate than alternative…
We consider linear iterative schemes for the time-discrete equations stemming from a class of nonlinear, doubly-degenerate parabolic equations. More precisely, the diffusion is nonlinear and may vanish or become multivalued for certain…
Combinatorial Optimization problems are widespread in domains such as logistics, manufacturing, and drug discovery, yet their NP-hard nature makes them computationally challenging. Recent Neural Combinatorial Optimization methods leverage…
The dynamics of the magnetic distribution in a ferromagnetic material is governed by the Landau-Lifshitz equation, which is a nonlinear geometric dispersive equation with a nonconvex constraint that requires the magnetization to remain of…
In this paper, we further investigate and refine the subspace-constrained preconditioning technique to enhance the theoretical and numerical convergence properties of randomized iterative methods for solving linear systems. In particular,…
The Lagrange-mesh method is a powerful method to solve eigenequations written in configuration space. It is very easy to implement and very accurate. Using a Gauss quadrature rule, the method requires only the evaluation of the potential at…
We present a scalable block preconditioning strategy for the trace system coming from the high-order hybridized discontinuous Galerkin (HDG) discretization of incompressible resistive magnetohydrodynamics (MHD). We construct the block…
We present a stationary iteration based upon a block splitting for a class of indefinite least squares problem. Convergence of the proposed method is investigated and optimal value of the involving parameter is used. The induced…
We present our open-source Python module Commics for the study of the magnetization dynamics in ferromagnetic materials via micromagnetic simulations. It implements state-of-the-art unconditionally convergent finite element methods for the…
We consider 4-dimensional spacetime manifolds that are piecewise Lorentzian, where the Lorentzian components of the manifold are separated by codimension-one planes (spacelike or timelike) on which the metric is degenerate. Such manifolds…
We consider linear systems arising from the use of the finite element method for solving scalar linear elliptic problems. Our main result is that these linear systems, which are symmetric and positive semidefinite, are well approximated by…
In this paper, a linearized Crank-Nicolson Galerkin finite element method (FEM) for generalized Ginzburg-Landau equation (GLE) is considered, in which, the difference method in time and the standard Galerkin FEM are employed. Based on the…
The convergence of GMRES for solving linear systems can be influenced heavily by the structure of the right hand side. Within the solution of eigenvalue problems via inverse iteration or subspace iteration, the right hand side is generally…
Sparse linear iterative solvers are essential for many large-scale simulations. Much of the runtime of these solvers is often spent in the implicit evaluation of matrix polynomials via a sequence of sparse matrix-vector products. A variety…
The dynamical equation of the magnetization has been reconsidered with enlarging the phase space of the ferromagnetic degrees of freedom to the angular momentum. The generalized Landau-Lifshitz-Gilbert equation that includes inertial terms,…