Related papers: A Fast Potential and Self-Gravity Solver for Non-A…
Accelerating the solution of linear systems of equations is critical due to their central role in numerous applications, such as numerical simulations, data analytics, and machine learning. This paper presents an analog solver circuit…
The plane wave method is most widely used for solving the Kohn-Sham equations in first-principles materials science computations. In this procedure, the three-dimensional (3-dim) trial wave functions' fast Fourier transform (FFT) is a…
In the present work, we study how to develop an efficient solver for the fast resolution of large and sparse linear systems that occur while discretizing elliptic partial differential equations using isogeometric analysis. Our new approach…
We present accurate models of the gravitational potential produced by a radially exponential disk mass distribution. The models are produced by combining three separate Miyamoto-Nagai disks. Such models have been used previously to model…
Self-gravity plays an important role in the evolution of rotationally supported systems such as protoplanetary disks, accretion disks around black holes, or galactic disks, as it can both feed turbulence or lead to gravitational…
The density functional theory (DFT) in electronic structure calculations can be formulated as either a nonlinear eigenvalue or direct minimization problem. The most widely used approach for solving the former is the so-called…
We apply our method of indirect integration, described in Part I, at fourth order, to the radial fall affected by the self-force. The Mode-Sum regularisation is performed in the Regge-Wheeler gauge using the equivalence with the harmonic…
Exchange symmetry in acceleration partitions the configuration space of an N particle, one-dimensional, gravitational system into N! equivalent cells. We take advantage of the resulting small angular extent of each cell to construct a…
We have developed a computational package for the calculation of numerically exact internal and external gravitational potential, its functional derivatives and sensitivity kernels, in an aspherical, heterogeneous planet. We detail our…
We study and compare different numerical differential equation solvers on the basis of numerical complexity, energy conservation, and stable solution in phase-space for the Simple Harmonic Oscillation (SHM) problem. We conclude and show…
Axially symmetric solutions for f (R)-gravity can be derived starting from exact spherically sym- metric solutions achieved by Noether symmetries. The method takes advantage of a complex coordi- nate transformation previously developed by…
We present an efficient solver for massively-parallel direct numerical simulations of incompressible turbulent flows. The method uses a second-order, finite-volume pressure-correction scheme, where the pressure Poisson equation is solved…
We study the computation of equilibrium points of electrostatic potentials: locations in space where the electrostatic force arising from a collection of charged particles vanishes. This is a novel scenario of optimization in which…
We describe three approaches for computing a gravity signal from a density anomaly. The first approach consists of the classical "summation" technique, whilst the remaining two methods solve the Poisson problem for the gravitational…
I formulate a general finite element method (FEM) for self-gravitating stellar systems. I split the configuration space to finite elements, and express the potential and density functions over each element in terms of their nodal values and…
We experimentally study the mechanical pressure exerted by a set of respectively passive isotropic and self-propelled polar disks onto two different flexible unidimensional membranes. In the case of the isotropic disks, the mechanical…
In this paper, we present a systematic study of the force-free field equation for simple axisymmetric configurations in spherical geometry and apply it to the solar active regions. The condition of separability of solutions in the radial…
The construction of robust solvers for linear systems obtained from the discretization of partial differential equations using Isogeometric Analysis is challenging since the condition number of the system matrix not only grows with the…
This work focuses on a quasi-linear-in-complexity strategy for a hybrid surface-wire integral equation solver for the electroencephalography forward problem. The scheme exploits a block diagonally dominant structure of the wire self block…
A direct integration algorithm is described to compute the magnetostatic field and energy for given magnetization distributions on not necessarily uniform tensor grids. We use an analytically-based tensor approximation approach for…